Skip to main content

Inflammatory and myeloid-associated gene expression before and one day after infant vaccination with MVA85A correlates with induction of a T cell response



Tuberculosis (TB) remains a global health problem, with vaccination likely to be a necessary part of a successful control strategy. Results of the first Phase 2b efficacy trial of a candidate vaccine, MVA85A, evaluated in BCG-vaccinated infants were published last year. Although no improvement in efficacy above BCG alone was seen, cryopreserved samples from this trial provide an opportunity to study the immune response to vaccination in this population.


We investigated blood samples taken before vaccination (baseline) and one and 28 days post-vaccination with MVA85A or placebo (Candin). The IFN-γ ELISpot assay was performed at baseline and on day 28 to quantify the adaptive response to Ag85A peptides. Gene expression analysis was performed at all three timepoints to identify early gene signatures predictive of the magnitude of the subsequent adaptive T cell response using the significance analysis of microarrays (SAM) statistical package and gene set enrichment analysis.


One day post-MVA85A, there is an induction of inflammatory pathways compared to placebo samples. Modules associated with myeloid cells and inflammation pre- and one day post-MVA85A correlate with a higher IFN-γ ELISpot response post-vaccination. By contrast, previous work done in UK adults shows early inflammation in this population is not associated with a strong T cell response but that induction of regulatory pathways inversely correlates with the magnitude of the T cell response. This may be indicative of important mechanistic differences in how T cell responses develop in these two populations following vaccination with MVA85A.


The results suggest the capacity of MVA85A to induce a strong innate response is key to the initiation of an adaptive immune response in South African infants but induction of regulatory pathways may be more important in UK adults. Understanding differences in immune response to vaccination between populations is likely to be an important aspect of developing successful vaccines and vaccination strategies.

Trial registration number NCT00953927

Peer Review reports


Tuberculosis (TB) is a major global health problem with an estimated 8.6 million new cases and 1.3 million deaths in 2012 [1]. Effective vaccination is likely to be necessary for the long-term control of the TB epidemic however Bacille Calmette-Guerin (BCG), the only currently licensed vaccine, provides variable protection against pulmonary disease [2]. Despite high BCG coverage, the incidence of TB remains high in endemic countries. Research efforts into new TB vaccines have focused largely on two strategies; either to modify BCG or replace it with an attenuated strain of Mycobacterium tuberculosis (M.tb), or to improve the protection provided by BCG through prime-boost regimes, often using viral vectors expressing TB antigens to enhance the pool of circulating memory cells primed by vaccination with BCG [3]. Twelve novel TB vaccines are currently in clinical trials, including two in Phase 2b efficacy studies [1]. The results of the first efficacy trial of a novel vaccine, Modified Vaccinia virus Ankara expressing antigen 85A (MVA85A) were published in early 2013 [4].

Although the efficacy of boosting BCG with MVA85A was not superior to that of BCG alone, the study demonstrated that a trial of a novel TB vaccine is feasible in a high burden setting. Furthermore, the collection of blood samples from all infants during the trial will enable research into the mechanisms of disease risk and response to vaccination in this setting. Previous studies with MVA85A have shown it to be safe and immunogenic in several diverse populations including adults without or with latent TB infection in the UK; healthy, latently infected and HIV-infected adults in Africa and healthy adolescents, children and infants in Africa [58]. The vaccine shows a quantitatively lower immunogenicity in African adults and in younger children and infants in Africa compared to UK adults [4]. In all trials to date, MVA85A induces antigen-specific Th1 and Th17 cells, believed to be important in protection against tuberculosis [911]. If, as has been suggested [12], low T cell responses to vaccination in this trial contributed to the lack of vaccine efficacy, understanding the mechanisms determining the magnitude of the response to vaccination is important to the development of an improved vaccine.

Several studies published in the last five years have demonstrated the power of genomics approaches in understanding the molecular mechanisms of the immune response to vaccination [1316]. Work using the yellow fever vaccine, YF-17D, identified a gene expression signature in circulating leukocytes of vaccinated volunteers shortly after vaccination which could predict the magnitude of the subsequent CD8+ T cell response [13]. Follow-up studies have yielded further mechanistic insight, showing that activation of the nutrient sensor GCN2 in dendritic cells following vaccination leads to increased antigen presentation and the development of a stronger immune response [17]. Similar approaches have been used by other groups and are beginning to reveal some of the factors contributing to the variability of the human immune response. The data show the importance of innate pathways in determining the magnitude of subsequent adaptive immune responses with a role for the stress response and gut microbiota in particular [18]. The cellular environment and its modification by vaccines and adjuvants are determined by many factors, underscoring the large variation seen in the immune response to vaccines across individuals and populations.

Understanding the differences in immune responses in different groups is key to developing targeted approaches to vaccination. The immune system changes with age, with a decrease in response to vaccination often noted in the elderly [19], though a recent study found no decrease in T cell responses in older adults following a novel MVA-vectored influenza vaccine [20]. Inflammation, apoptosis and immune senescence have all been linked to lower responses to vaccination in this age group [16, 21]. As the population of the developed world continues to age, understanding these factors will be important in developing effective vaccination strategies. At the other end of the spectrum are the immature immune systems of infants, which also differ to those of adults [8, 22]. Although children and infants, particularly in the developing world, are the target population for many vaccines against infectious diseases, the factors underpinning the immune response to vaccination in this population remain poorly characterized. Understanding the immune response to vaccination in infants living in areas with a high burden of disease, and how this differs from the immune response of healthy, young adults living in areas of lower disease prevalence –but in whom early testing of vaccines is usually carried out- is therefore a crucial component in the development and early selection of many of the vaccines in development.

In this study, we have analysed gene expression signatures pre- and post-vaccination in infants from the MVA85A Phase 2b efficacy trial who did not develop TB disease during the trial (non-cases) and correlated these changes to the antigen-specific T cell response to vaccination, as measured by IFN-γ ELISpot to Ag85A peptides, in an effort to understand the variability in response to vaccination in this setting. Finally, we compare these findings with previous work performed in UK adults receiving the same vaccine [23], in an effort to characterise some of the differences between these populations.


Origin of samples

Samples used in these experiments were cryopreserved peripheral blood mononuclear cells (PBMC) or whole blood in RNA lysis buffer from a double-blind, randomised, placebo-controlled Phase 2b efficacy trial of a candidate TB vaccine, MVA85A, in BCG-vaccinated, HIV-negative South African infants (South African National Clinical Trials Register DOH-27-0109-2654, NCT00953927). Infants were randomized to receive either one dose of MVA85A (1 × 108 plaque forming units in 0.06 mL) or an equal volume of Candida skin test antigen (Candin, AllerMed, USA) as placebo at 4–6 months of age [4]. The trial was approved by the University of Cape Town Faculty of Health Sciences Human Research Ethics Committee, Oxford University Tropical Research Ethics Committee, and the Medicines Control Council of South Africa. Parents or legal guardians provided written, informed consent. Storage of samples for exploratory immunological analyses was fully ethically approved.

The samples used in this study were selected from a subset of 100 infants who had a small blood sample taken one day post-vaccination. The samples were selected to exclude cases (infants who went on to develop TB disease) and controls (selected to demographically match the cases), which will be used in a future study looking at correlates of risk of TB disease. Therefore none of the samples used in this study were taken from infants diagnosed with TB during the course of the trial. Cells were collected 0–7 days pre- and 28 days post-vaccination with MVA85A/placebo in cell preparation tubes with sodium heparin (CPT; Vacutainer; BD) and PBMC separated and cryopreserved. The cells were thawed and stimulated as detailed below. One day post-vaccination, 50-300 μL of whole blood was collected directly via heel prick into a tube filled with RLT buffer (RNeasy kit, Qiagen) containing 10 μL/mL β-mercaptoethanol using a BD Quikheel Lancet. The sample was immediately frozen and RNA extracted as detailed below.

Cell thawing

Cryopreserved PBMC were rapidly thawed in a 37°C waterbath and transferred to a 15 mL Falcon tube containing 10 mL R10 (RPMI, 10% FCS, 1% L-glutamine, 1% Pen-Strep and 1% sodium pyruvate). PBMC were pelleted, supernatants discarded and resuspended in 10 mL R10 with 20 μL Benzonase (Merck Chemicals Ltd.) and rested overnight at 37°C, 5% CO2. PBMC were counted on a Casy Counter (Roche) and split into appropriate volumes for each assay. Not all assays were performed on all samples.

Ex-vivoIFN-γ ELISpot assay

The ex-vivo IFN-γ ELISpot assay was performed on thawed PBMC samples collected pre-vaccination and 28 days post-vaccination as previously described [23]. 3 × 105 PBMC were stimulated in triplicate with a pool of Ag85A peptides, consisting of 66 15mers, overlapping by 10 amino acids (2 μg/ml) (Peptide Protein Research).

Gene expression assays

PBMC from pre-vaccination and 28 days post-vaccination were incubated for 12 hours with either R10 media alone (unstimulated) or pooled Ag85A peptides as described for the ELISpot (2 μg/mL). After 12 hours supernatant was removed and the PBMC resuspended in 350uL RLT buffer (Qiagen) containing 10 μL/mL β-mercaptoethanol and frozen at -20C.

Blood in RLT buffer was thawed and RNA extracted using the RNeasy kit (Qiagen) according to manufacturer’s instructions, including the optional protocol for DNA digest (RNase-free DNase kit, Qiagen). The protocol was modified in the following way for the heelprick samples due to the small volume of whole blood collected: 80% ethanol was added to precipitate the RNA (rather than the recommended 70%) and an extra wash with 350 μL RW1 buffer was performed prior to DNA digest.

Messenger RNA was amplified from the total RNA using the Illumina Totalprep kit (Ambion) according to manufacturer’s instructions. RNA quantity and quality was assessed using a Nanodrop ND-1000 Spectrophotometer and an Agilent Bioanalyser (Agilent RNA 6000 Nano Kit).

750 ng amplified cRNA was labeled and hybridized to Illumina Human HT-12 v4 beadchips as specified in the manufacturer’s instructions. Beadchips were scanned on an Illumina iScan machine and data extracted using the GenomeStudio software.

RNA-Seq was performed on pre- and post-vaccination PBMC from 6 infants (12 samples). Total RNA was sent to the Beijing Genomics Institute (BGI). Libraries were constructed using the TruSeq kit (Illumina) and sequenced on a HiSeq2000 sequencer, using paired-end reads of 90 bp and 30 M sequencing depth. Quality control was performed and the reads aligned to the genome hg19, downloaded from UCSC (, using SOAPaligner2.21 with the following constraints: maximum number of mismatches allowed on one read is 5 bp, no gaps allowed, only repeat hits reported.

Data analysis

Ex vivo IFN-γ ELISpot assay

Phytohaemagglutinin (PHA) (Sigma) was used as a positive control and unstimulated wells were used as a measure of background IFN-γ production. Results are reported as spot forming cells (SFC) per million PBMC, calculated by subtracting the mean of the unstimulated wells from the mean of triplicate antigen wells and correcting for the number of PBMC in the well. An ELISpot response was deemed positive if the average count in the positive control wells was at least twice that in the negative control wells and at least 5 spots more than the negative control wells.

Illumina microarray

The R package arrayQualityMetrics [24] was used to assess sample quality. 82 heelprick samples and all PBMC samples passed quality control. The R package limma was used to perform background correction and normalization and the gene list was filtered using the gene filter package to remove genes with an expression IQR < 0.3 (log2 transformed). Lists of differentially expressed genes were generated using limma (p-value cut-off of 0.05 after Benjamini-Hochberg correction [2527]).

The package Significance Analysis of Microarrays (SAM) was used to rank genes correlating with the IFN-γ ELISpot response according to the strength of the correlation [2830]. The ranked gene list was then inputted into Broad Institute gene set enrichment analysis programme (GSEA) [31, 32] as an externally supplied preranked list. The reference gene set used was the Blood Transcription Modules compiled by Li et al. [33]. The significance of module enrichment was assessed by permutation in the GSEA program.

All heatmaps were generated in R, using Euclidean distance and average linkage as methods to calculate the distance matrix and hierarchical clustering respectively. Where correlations are shown, these use the Pearson product–moment correlation coefficient.

RNA sequencing

Genes for which the reads per kilobase per million value (rpkm) was <1 in over 40% of samples were excluded. RPKM was calculated as (10^9*C)/(N*L) where C = number of reads uniquely mapped to transcript, N = total number of uniquely mapped reads in sample and L = maximum length of transcript. Raw counts were then analysed in the R package limma and converted to log2 transformed counts per million. These values were compared to the gene expression values for the equivalent samples obtained by microarray analysis.

Accession codes

Gene expression omnibus: GSE56559 (day 1 heelprick, South African infants), GSE56561 (PBMC, South African infants) GSE40719 (UK adults).

Results and discussion

Table 1 shows the number of samples used in each assay in this study.

Table 1 Samples processed for each assay as part of this study

MVA85A induces an inflammatory signature one day post-vaccination

Illumina microarray gene expression analysis of 82 whole blood samples taken one day post-vaccination (37 MVA85A, 45 Candin placebo) identified 32 differentially expressed genes. These genes were largely associated with the immune response (differentially expressed genes are shown in Table 2 with genes associated with the immune response highlighted in bold). Hierarchical clustering showed that infants fall into three clusters with a mixed MVA85A/Candin cluster exhibiting a moderate level of expression of inflammatory genes (Figure 1). This analysis shows MVA85A-vaccinated infants exhibit a more inflammatory gene expression profile than those in the placebo group however the range is large in both groups and there is an overlap between the two groups. This overlap could be due to the immunomodulatory properties of Candin, which induces inflammation and may lead to functional reprogramming of monocytes associated with protection from subsequent infection [34].

Table 2 Differentially expressed genes 1 day post-vaccination: MVA85A vs Candin
Figure 1

Heatmap of differential gene expression 1 day post vaccination. Heatmap showing changes in gene expression one day post-vaccination with either MVA85A or a Candin placebo. Colors at the top show vaccination: MVA85A (blue) or Candin (red). Genes were selected on the basis of differential expression between the two groups (fdr < 0.05). Clustering using euclidean distance and average clustering methods. Red indicates up-regulated mRNA, blue indicates downregulated.

The immune response to stimulation with Ag85A peptides

The antigen-specific immune response to Ag85A was assessed by IFN-γ ELISpot and Illumina microarray gene expression analysis. Infants vaccinated with MVA85A had a significantly higher post-vaccination Ag85A-specific ELISpot response than the Candin group (Figure 2a). Differentially expressed genes between the vaccine and placebo groups, and pre- and post-vaccination time points, are shown in Table 3. The genes induced following Ag85A peptide stimulation are all associated with the STAT1 pathway and exhibit a highly correlated pattern of expression. Furthermore, upregulation of this pathway occurred in infants who also had a detectable Ag85A-specific T cell response by IFN-γ ELISpot assay but not in unstimulated cells or infants who received the candin placebo (Figure 2b,c). This observation suggests that, in this population, gene expression analysis does not add substantial information to that measured by the IFN-γ ELISpot assay in capturing the response to Ag85A peptide stimulation following MVA85A vaccination.

Figure 2

Stimulation of PBMC with Ag85A peptide pools. (a) IFN-γ ELISpot of PBMC taken pre- or post-vaccination in MVA85A or candin-vaccinated infants. Responses are significantly higher in the MVA85A group, 4 weeks post-vaccination (Wilcoxon test, p < 0.05). b,c. STAT1 expression in unstimulated (b) and 85A peptide-stimulated (c) PBMC. Squares = pre-vaccination, circles = 28 days post-vaccination; red = candin, blue = MVA85A. Dashed lines indicate a positive ELISpot response (x = 15) and elevated expression of STAT1 (y = 9.5).

Table 3 Differentially expressed genes following Ag85A stimulation post-MVA85A vaccination

A proportion of infants did not respond to antigenic stimulation with Ag85A peptides following MVA85A vaccination. This lack of response was observed both by IFN-γ ELISpot and gene expression analysis. It has been suggested that low or absent responses to MVA85A may be one explanation for the lack of efficacy observed in the trial [12]. Further analysis of cases and controls is under way to address this question however, in this smaller study, we next investigated some of the mechanisms underlying the magnitude of the adaptive immune response which develops following MVA85A vaccination.

Myeloid cells and inflammation are associated with a higher ELISpot response

The following analyses were performed using only samples from infants vaccinated with MVA85A, to further investigate the mechanisms underlying the immune response to this vaccine. The R package Significance Analysis of Microarrays (SAMR) was used to identify genes whose expression correlated with the Ag85A-specific T cell frequencies measured 28 days post-vaccination by IFN-γ ELISpot assay. SAM identifies genes significantly correlating with a continuous response variable, in this case the IFN-γ ELISpot, and outputs a positive and negative set of genes based on the strength of the correlation of each gene with higher (positive) or lower (negative) values of the response phenotype [29, 35]. Since the IFN-γ ELISpot responses in this study were very low, we have subsequently defined infants as responders or non-responders. An ELISpot response was deemed positive if the average count in the positive control wells was at least twice that in the negative control wells and at least 5 spots more than the negative control wells [36]. Sets of correlating genes pre-vaccination and one day post-vaccination were generated and ranked according to their significance score (Figure 3a). This ranked list of genes was then analysed using the Broad Institute Gene Set Enrichment Analysis PreRanked function, using the Blood Transcription Modules defined by Li et al. as a reference gene set [33]. Pre-vaccination, responder infants have an over-representation of genes enriched in monocytes, activated dendritic cells and neutrophils as well chemokines and inflammatory pathways (Figure 3b). One day post-vaccination, there is a negative association between lymphoid cells and the subsequent development of an ELISpot response (Figure 3c). Moreover, there is positive enrichment of gene sets associated with myeloid cells, inflammation and an antiviral response. The expression pattern of these clusters and their relationship to the ELISpot response is shown in heatmaps in Figure 4. Higher expression of genes associated myeloid cells and inflammation pre- and 1 day post-vaccination are both associated with the development of an antigen-specific T cell response to vaccination with MVA85A, suggesting the ability of MVA to induce a strong innate response is key to its function as a vaccine vector in this population. Additionally, higher expression of genes associated with activation of lymphoid cells such as NK cells and cytotoxic T cells one day post-vaccination is associated with an absent ELISpot response 28 days later.

Figure 3

Modular analysis of genes associated with a higher response to vaccination. Genes whose expression correlate with the IFNg ELISpot response were identified using the R package SAM (a) and ranked in order of their score. The ranked list was then analysed in GSEA using the Blood Transcription Modules compiled by Li et al. as the reference gene set [33]. The results of the analysis for genes expressed in PBMC taken pre-vaccination or whole blood taken one day post-vaccination are shown (b). Length of the bar shows Normalised enrichment score for each module, number in the bars indicates genes in the test list present ineach reference set. Colour saturation indicates genes present as a percentage of total genes within the module (signal). Red modules are positively associated with ELISpot response, blue modules are negatively associated.

Figure 4

Heatmaps of BTMs associated with a higher response to vaccination. Heatmaps show expression of genes for MVA85A-vaccinated infants in unstimulated PBMC taken pre-vaccination (a) or whole blood taken one day post-vaccination (b). The colour coding along the top of the heatmap shows responder (blue) and non-responder (black) infants as measured by IFN-γ ELISpot. Genes are arranged by module (right).

The ratio of myeloid to lymphoid cells has previously been associated with differences in TB disease risk in HIV-infected South African adults and susceptibility to malaria and influenza in other cohorts [3740] and this may be another example of an outcome associated with this ratio. The gene expression profiles associated with responder infants are present pre-vaccination and show a strong overlap pre- and one day post vaccination. This suggests the baseline inflammatory profile of the infant, including the differing proportion of circulating leukocytes, is key to determining the response to vaccination. This may be influenced by genetic influences on innate immunity or environmental exposure preceding vaccination, including the response to BCG, which all infants received at birth. As MVA preferentially infects myeloid cells [41], a higher proportion of myeloid cells may lead to overall increased viral expression of Ag85A protein in infants with higher frequencies of myeloid cells. Conversely, killing of infected myeloid cells by cytotoxic T cells and NK cells may decrease antigen expression, inhibiting the development of a response to Ag85A.In a subset of samples, the RNA extracted from the PBMC collected pre- and 28 days post-MVA85A was also measured by RNA Sequencing, as part of a pilot project for future studies. Gene expression values as measured independently by these two methods are highly correlated (Figure 5), providing a technical validation for this result.

Figure 5

Comparison of gene expression values measured by beadchip microarray and RNA sequencing. Plot of median values from the RNA-Seq and microarray data of 12 unstimulated PBMC samples for 74 genes used in the modular analysis. Significance assessed using Pearson’s correlation.

Gene expression signatures correlating with immunogenicity differ in South African infants and UK adults

We have previously described changes in gene expression in unstimulated PBMC from UK adults vaccinated with the same regime: a BCG prime followed by an MVA85A boost [23, 42]. We therefore wanted to compare the observations made in these two different populations. In the adults, there was no placebo and PBMC were collected for gene expression analysis pre-, two and seven days post-MVA85A. Comparing gene expression to the magnitude of the induced T cell response showed a positive association with TLR1 expression at baseline and a negative correlation with regulatory genes including STAT5B and CTLA4 at day 2 [42]. Table 4 shows the genes differentially expressed in infants one day post-MVA85A (compared to a placebo group at the equivalent timepoint). Genes in bold were also differentially expressed in UK adults two days post-MVA85A compared to baseline and fold changes are shown for both comparisons. In both populations, CXCL10 had the highest fold change at the innate timepoint (one or two days post-MVA85A). One interesting difference is the gene encoding indoleamine 2,3-dioxygenase (probes INDO and IDO1) which is highly differentially expressed in the infants but not the adults.

Table 4 Comparison of differentially expressed genes in South African infants and UK adults post-MVA85A

Further to this, the relationships between these genes show interesting differences between the two populations (Figure 6). In the South African infants, expression of both CXCL10 and IDO1 one day post-vaccination correlates with the IFN-γ ELISpot response and with each other. In the UK adults however, expression of these genes correlates neither with each other nor with the IFN-γ ELISpot response. In the UK adults, we have previously described a negative association between expression of CTLA4 post-vaccination and the IFN-γ ELISpot response. However, this association was not found in South African infants. Furthermore, CTLA4 correlated positively with IDO1 in the UK adults but negatively in the South African infants.

Figure 6

Differences in inflammation and regulation: RSA infants and UK adults. a. Expression in unstimulated whole blood 1 day post-MVA85A of CXCL10 and IDO1 but not CTLA4 correlates with ELISpot response in RSA infants. CTLA4 expression in unstimulated PBMC 2 days post-MVA85A is inversely correlated with ELISpot response in UK adults but IDO1 and CXCL10 show no significant association. b. In RSA infants, IDO1 expression correlates positively with CXCL10 and negatively with CTLA4 expression. In UK adults, CTLA4 and IDO1 expression show a positive association. All tests are Pearson’s correlations.

Indoleamine 2,3-dioxygenase (IDO1) is an enzyme catalyzing the first and rate-limiting step in tryptophan catabolism. This enzyme has multiple physiological effects, including immunosuppression and regulation of T cells. In UK adults, expression of IDO1 in PBMC correlates with that of the inhibitory co-receptor CTLA4. A recent study showing IDO1 is a critical resistance mechanism in antitumor T cell immunotherapy targeting CTLA-4 in a mouse model of melanoma suggests the interaction between these two genes is an important component to consider when inducing an immune response for therapeutic purposes [43]. In the South African infants however, we did not observe a correlation between expression of IDO1 and CTLA4 but IDO1 correlated with CXCL10 instead, suggesting either a different role for this enzyme or that inflammatory and regulatory responses are closely coupled in this population, perhaps to protect against excessive immune responses.

These disparities suggest there may be important differences in the immune response to vaccination in these two populations and in the pathways that determine the adaptive immune response. However, it is impossible to determine whether the differences observed are attributable to age or to genetic or environmental differences between the populations. Furthermore, the post-vaccination timepoint differed in both cell type (whole blood in the infants, PBMC in adults) and time taken (one or two days post-MVA85A). PBMC were used in the pre-vaccination timepoint in both studies. Despite this, the similarities in the list of differentially expressed genes suggest these comparisons are meaningful and warrant further investigation.

UK adults have a stronger innate response to the vaccine, however the magnitude of this response does not predict the magnitude of the antigen-specific response. By contrast, in the South African infants the magnitude of inflammatory gene expression one day post-MVA85A correlated with the ELISpot response at week 4. Further investigation of these mechanisms and how they differ between the different populations in which vaccines are tested and deployed may prove important in future vaccine development strategies.


This study has shown an association between high levels of inflammation and myeloid signatures both pre- and post-vaccination and the development of the antigen-specific T cell response to MVA85A in BCG-vaccinated South African infants. In this population, the capacity of the vaccine to induce a strong innate response appears key to its ability to initiate an adaptive immune response. Furthermore, we describe differences in the response to vaccination with MVA85A in UK adults and South African infants, suggesting different immune pathways may determine immune responses in these two very different cohorts. This study has focused on investigating the mechanisms underlying vaccine immunogenicity, not vaccine-induced protection. Ultimately, an understanding of both these aspects of vaccination and how they differ across individuals and populations is likely to be necessary in achieving protection in the diverse groups and areas of the world still plagued by TB.

Author contributions

MM, SAH, IS, LS, RT and MO carried out the experiments; MT, HMcS, MH, TJS, WAH and HAF were involved in study design and sample collection; MM analysed the data and drafted the manuscript. All authors reviewed and approved the final manuscript.


  1. 1.

    World Health Organisation: Global Tuberculosis Report 2013. 2013, Available at:,

    Google Scholar 

  2. 2.

    Colditz GA, Brewer TF, Berkey CS, Burdick E, Fineberg HV, Mosteller F: Efficacy of BCG vaccine in the prevention of tuberculosis. JAMA. 1994, 271: 698-702. 10.1001/jama.1994.03510330076038.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Ottenhoff THM, Kaufmann SHE: Vaccines against tuberculosis: where are we and where do we need to go?. PLoS Pathog. 2012, 8: e1002607-10.1371/journal.ppat.1002607.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Tameris MD, Hatherill M, Landry BS, Scriba TJ, Snowden MA, Lockhart S, Shea JE, McClain JB, Hussey GD, Hanekom WA, Mahomed H, McShane H: Safety and efficacy of MVA85A, a new tuberculosis vaccine, in infants previously vaccinated with BCG: a randomised, placebo-controlled phase 2b trial. Lancet. 2013, 6736: 1-8.

    Google Scholar 

  5. 5.

    Hawkridge T, Scriba TJ, Gelderbloem S, Smit E, Tameris M, Moyo S, Lang T, Veldsman A, Hatherill M, Van Der Merwe L, Fletcher HA, Mahomed H, Hill AVS, Hanekom WA, Hussey GD, McShane H: Safety and immunogenicity of a new tuberculosis vaccine, MVA85A, in healthy adults in South Africa. J Infect Dis. 2008, 198: 544-552. 10.1086/590185.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Sander CR, Pathan AA, Beveridge NER, Poulton I, Minassian A, Alder N, Van Wijgerden J, Hill AVS, Gleeson FV, Davies RJO, Pasvol G, McShane H: Safety and immunogenicity of a new tuberculosis vaccine, MVA85A, in Mycobacterium tuberculosis-infected individuals. Am J Respir Crit Care Med. 2009, 179: 724-733. 10.1164/rccm.200809-1486OC.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Minassian AM, Rowland R, Beveridge NER, Poulton ID, Satti I, Harris S, Poyntz H, Hamill M, Griffiths K, Sander CR, Ambrozak DR, Price DA, Hill BJ, Casazza JP, Douek DC, Koup RA, Roederer M, Winston A, Ross J, Sherrard J, Rooney G, Williams N, Lawrie AM, Fletcher HA, Pathan AA, McShane H: A Phase I study evaluating the safety and immunogenicity of MVA85A, a candidate TB vaccine, in HIV-infected adults. BMJ Open. 2011, 1: e000223-

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Scriba TJ, Tameris M, Mansoor N, Smit E, van der Merwe L, Mauff K, Hughes EJ, Moyo S, Brittain N, Lawrie A, Mulenga H, de Kock M, Gelderbloem S, Veldsman A, Hatherill M, Geldenhuys H, Hill AVS, Hussey GD, Mahomed H, Hanekom WA, McShane H: Dose-finding study of the novel tuberculosis vaccine, MVA85A, in healthy BCG-vaccinated infants. J Infect Dis. 2011, 203: 1832-1843. 10.1093/infdis/jir195.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Scriba TJ, Tameris M, Mansoor N, Smit E, van der Merwe L, Isaacs F, Keyser A, Moyo S, Brittain N, Lawrie A, Gelderbloem S, Veldsman A, Hatherill M, Hawkridge A, Hill AVS, Hussey GD, Mahomed H, McShane H, Hanekom WA: Modified vaccinia Ankara-expressing Ag85A, a novel tuberculosis vaccine, is safe in adolescents and children, and induces polyfunctional CD4+ T cells. Eur J Immunol. 2010, 40: 279-290.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Griffiths KL, Pathan AA, Minassian AM, Sander CR, Beveridge NER, Hill AVS, Fletcher HA, McShane H: Th1/Th17 cell induction and corresponding reduction in ATP consumption following vaccination with the novel Mycobacterium tuberculosis vaccine MVA85A. PLoS One. 2011, 6: e23463-10.1371/journal.pone.0023463.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    De Cassan SC, Pathan AA, Sander CR, Minassian A, Rowland R, Hill AVS, McShane H, Fletcher HA: Investigating the induction of vaccine-induced Th17 and regulatory T cells in healthy, Mycobacterium bovis BCG-immunized adults vaccinated with a new tuberculosis vaccine, MVA85A. Clin Vaccine Immunol. 2010, 17: 1066-1073. 10.1128/CVI.00047-10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Dye C, Fine PE: A major event for new tuberculosis vaccines. Lancet. 2013, 6736: 3-5.

    Google Scholar 

  13. 13.

    Querec TD, Akondy RS, Lee EK, Cao W, Nakaya HI, Teuwen D, Pirani A, Gernert K, Deng J, Marzolf B, Kennedy K, Wu H, Bennouna S, Oluoch H, Miller J, Vencio RZ, Mulligan M, Aderem A, Ahmed R, Pulendran B: Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans. Nat Immunol. 2009, 10: 116-125.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Gaucher D, Therrien R, Kettaf N, Angermann BR, Boucher G, Filali-Mouhim A, Moser JM, Mehta RS, Drake DR, Castro E, Akondy R, Rinfret A, Yassine-Diab B, Said EA, Chouikh Y, Cameron MJ, Clum R, Kelvin D, Somogyi R, Greller LD, Balderas RS, Wilkinson P, Pantaleo G, Tartaglia J, Haddad EK, Sékaly R-P: Yellow fever vaccine induces integrated multilineage and polyfunctional immune responses. J Exp Med. 2008, 205: 3119-3131. 10.1084/jem.20082292.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Nakaya HI, Wrammert J, Lee EK, Racioppi L, Marie-Kunze S, Haining WN, Means AR, Kasturi SP, Khan N, Li G-M, McCausland M, Kanchan V, Kokko KE, Li S, Elbein R, Mehta AK, Aderem A, Subbarao K, Ahmed R, Pulendran B: Systems biology of vaccination for seasonal influenza in humans. Nat Immunol. 2011, 12: 786-795. 10.1038/ni.2067.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Furman D, Jojic V, Kidd B, Shen-Orr S, Price J, Jarrell J, Tse T, Huang H, Lund P, Maecker HT, Utz PJ, Dekker CL, Koller D, Davis MM: Apoptosis and other immune biomarkers predict influenza vaccine responsiveness. Mol Syst Biol. 2013, 9: 659-

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Ravindran R, Khan N, Nakaya HI, Li S, Loebbermann J, Maddur MS, Park Y, Jones DP, Chappert P, Davoust J, Weiss DS, Virgin HW, Ron D, Pulendran B: Vaccine activation of the nutrient sensor GCN2 in dendritic cells enhances antigen presentation. Science. 2013, 2 (December): 1-8.

    Google Scholar 

  18. 18.

    Pulendran B, Oh JZ, Nakaya HI, Ravindran R, Kazmin DA: Immunity to viruses: learning from successful human vaccines. Immunol Rev. 2013, 255: 243-255. 10.1111/imr.12099.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Goronzy JJ, Weyand CM: Understanding immunosenescence to improve responses to vaccines. Nat Immunol. 2013, 14: 428-436. 10.1038/ni.2588.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Antrobus RD, Lillie PJ, Berthoud TK, Spencer AJ, McLaren JE, Ladell K, Lambe T, Milicic A, Price DA, Hill AVS, Gilbert SC: A T Cell-Inducing Influenza Vaccine for the Elderly: Safety and Immunogenicity of MVA-NP + M1 in Adults Aged over 50 Years. PLoS One. 2012, 7: e48322-10.1371/journal.pone.0048322.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    McElhaney JE, Zhou X, Talbot HK, Soethout E, Bleackley RC, Granville DJ, Pawelec G: The unmet need in the elderly: how immunosenescence, CMV infection, co-morbidities and frailty are a challenge for the development of more effective influenza vaccines. Vaccine. 2012, 30: 2060-2067. 10.1016/j.vaccine.2012.01.015.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Bodewes R, Fraaij PLA, Osterhaus ADME, Rimmelzwaan GF: Pediatric influenza vaccination: understanding the T-cell response. Expert Rev Vaccines. 2012, 11: 963-971. 10.1586/erv.12.69.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Meyer J, Harris SA, Satti I, Poulton ID, Poyntz HC, Tanner R, Rowland R, Griffiths KL, Fletcher HA, McShane H: Comparing the safety and immunogenicity of a candidate TB vaccine MVA85A administered by intramuscular and intradermal delivery. Vaccine. 2013, 31: 1026-1033. 10.1016/j.vaccine.2012.12.042.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Kauffmann A, Huber W: Microarray data quality control improves the detection of differentially expressed genes. Genomics. 2010, 95: 138-142. 10.1016/j.ygeno.2010.01.003.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Smyth GK: Limma: linear models for microarray data. Bioinforma Comput Biol Solut using R Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.

    Google Scholar 

  26. 26.

    Shi W, Oshlack A, Smyth GK: Optimizing the noise versus bias trade-off for Illumina whole genome expression BeadChips. Nucleic Acids Res. 2010, 38: e204-10.1093/nar/gkq871.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Ritchie ME, Diyagama D, Neilson J, van Laar R, Dobrovic A, Holloway A, Smyth GK: Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006, 7: 261-10.1186/1471-2105-7-261.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Zhang S: A comprehensive evaluation of SAM, the SAM R-package and a simple modification to improve its performance. BMC Bioinformatics. 2007, 8: 230-10.1186/1471-2105-8-230.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Shen-Orr SS, Tibshirani R, Khatri P, Bodian DL, Staedtler F, Perry NM, Hastie T, Sarwal MM, Davis MM, Butte AJ: Cell type-specific gene expression differences in complex tissues. Nat Methods. 2010, 7: 287-289. 10.1038/nmeth.1439.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstråle M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34: 267-273. 10.1038/ng1180.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Li S, Rouphael N, Duraisingham S, Romero-Steiner S, Presnell S, Davis C, Schmidt DS, Johnson SE, Milton A, Rajam G, Kasturi S, Carlone GM, Quinn C, Chaussabel D, Palucka AK, Mulligan MJ, Ahmed R, Stephens DS, Nakaya HI, Pulendran B: Molecular signatures of antibody responses derived from a systems biology study of five human vaccines. Nat Immunol. 2013, 15: 1-12.

    Article  Google Scholar 

  34. 34.

    Quintin J, Saeed S, Martens JHA, Giamarellos-Bourboulis EJ, Ifrim DC, Logie C, Jacobs L, Jansen T, Kullberg B-J, Wijmenga C, Joosten LAB, Xavier RJ, van der Meer JWM, Stunnenberg HG, Netea MG: Candida albicans infection affords protection against reinfection via functional reprogramming of monocytes. Cell Host Microbe. 2012, 12: 223-232. 10.1016/j.chom.2012.06.006.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Zang S, Guo R, Zhang L, Lu Y: Integration of statistical inference methods and a novel control measure to improve sensitivity and specificity of data analysis in expression profiling studies. J Biomed Inform. 2007, 40: 552-560. 10.1016/j.jbi.2007.01.002.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    McShane H, Pathan AA, Sander CR, Keating SM, Gilbert SC, Huygen K, Fletcher HA, Hill AV: Recombinant modified vaccinia virus Ankara expressing antigen 85A boosts BCG-primed and naturally acquired antimycobacterial immunity in humans. Nat Med. 2004, 10: 1240-1244. 10.1038/nm1128.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Naranbhai V, Hill AVS, Abdool Karim SS, Naidoo K, Abdool Karim Q, Warimwe GM, McShane H, Fletcher H: Ratio of monocytes to lymphocytes in peripheral blood identifies adults at risk of incident tuberculosis among HIV-infected adults initiating antiretroviral therapy. J Infect Dis. 2013, 209: 1-10.

    Google Scholar 

  38. 38.

    Warimwe GM, Murungi LM, Kamuyu G, Nyangweso GM, Wambua J, Naranbhai V, Fletcher HA, Hill AVS, Bejon P, Osier FHA, Marsh K: The ratio of monocytes to lymphocytes in peripheral blood correlates with increased susceptibility to clinical malaria in Kenyan children. PLoS One. 2013, 8: e57320-10.1371/journal.pone.0057320.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Merekoulias G, Alexopoulos EC, Belezos T, Panagiotopoulou E: Lymphocyte to monocyte ratio as a screening tool for influenza. Plos Curr. 2010, 2: RRN1154-

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Warimwe GM, Fletcher HA, Olotu A, Agnandji ST, Hill AVS, Marsh K, Bejon P: Peripheral blood monocyte-to-lymphocyte ratio at study enrollment predicts efficacy of the RTS,S malaria vaccine: analysis of pooled phase II clinical trial data. BMC Med. 2013, 11: 184-10.1186/1741-7015-11-184.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Berthoud TK, Fletcher H, Porter D, Thompson F, Hill AVS, Todryk SM: Comparing human T cell and NK cell responses in viral-based malaria vaccine trials. Vaccine. 2009, 28: 21-27. 10.1016/j.vaccine.2009.09.132.

    Article  PubMed  Google Scholar 

  42. 42.

    Matsumiya M, Stylianou E, Griffiths K, Lang Z, Meyer J, Harris SA, Rowland R, Minassian AM, Pathan AA, Fletcher H, Mcshane H: Roles for treg expansion and HMGB1 signaling through the TLR1-2-6 axis in determining the magnitude of the antigen-specific immune response to MVA85A. PLoS One. 2013, 8: 1-14.

    Article  Google Scholar 

  43. 43.

    Holmgaard RB, Zamarin D, Munn DH, Wolchok JD, Allison JP: Indoleamine 2,3-dioxygenase is a critical resistance mechanism in antitumor T cell immunotherapy targeting CTLA-4. J Exp Med. 2013, 210: 1389-1402. 10.1084/jem.20130066.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


The authors would like to thank the trial participants and their families and the MVA85A 020 Trial Study team.


This work was funded by Aeras, the Wellcome Trust and the Oxford Emergent Tuberculosis Consortium and NEWTBVAC (EC FP7); HMcS is a Wellcome Trust Senior Clinical Reseach Fellow.

Author information



Corresponding authors

Correspondence to Magali Matsumiya or Helen A Fletcher.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Matsumiya, M., Harris, S.A., Satti, I. et al. Inflammatory and myeloid-associated gene expression before and one day after infant vaccination with MVA85A correlates with induction of a T cell response. BMC Infect Dis 14, 314 (2014).

Download citation


  • Tuberculosis
  • Vaccine
  • Innate immunity
  • Transcriptomics
  • MVA85A