Inflammatory and myeloid-associated gene expression before and one day after infant vaccination with MVA85A correlates with induction of a T cell response
BMC Infectious Diseases volume 14, Article number: 314 (2014)
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.
ClinicalTrials.gov number NCT00953927
Tuberculosis (TB) is a major global health problem with an estimated 8.6 million new cases and 1.3 million deaths in 2012 . 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 . 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 . Twelve novel TB vaccines are currently in clinical trials, including two in Phase 2b efficacy studies . The results of the first efficacy trial of a novel vaccine, Modified Vaccinia virus Ankara expressing antigen 85A (MVA85A) were published in early 2013 .
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 [5–8]. The vaccine shows a quantitatively lower immunogenicity in African adults and in younger children and infants in Africa compared to UK adults . In all trials to date, MVA85A induces antigen-specific Th1 and Th17 cells, believed to be important in protection against tuberculosis [9–11]. If, as has been suggested , 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 [13–16]. 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 . 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 . 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 . 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 , though a recent study found no decrease in T cell responses in older adults following a novel MVA-vectored influenza vaccine . 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 , 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, ClinicalTrials.gov 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 . 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.
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 . 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 (http://genome.ucsc.edu/), 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.
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.
The R package arrayQualityMetrics  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 [25–27]).
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 [28–30]. 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. . 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.
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.
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.
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 .
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.
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 . 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 . 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 . 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.
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 [37–40] 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 , 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.
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 . 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.
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.
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 . 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.
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.
World Health Organisation: Global Tuberculosis Report 2013. 2013, Available at: http://www.who.int/tb/publications/global_report/en/,
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.
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.
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.
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.
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.
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-
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.
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.
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.
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.
Dye C, Fine PE: A major event for new tuberculosis vaccines. Lancet. 2013, 6736: 3-5.
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.
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.
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.
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-
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.
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.
Goronzy JJ, Weyand CM: Understanding immunosenescence to improve responses to vaccines. Nat Immunol. 2013, 14: 428-436. 10.1038/ni.2588.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Merekoulias G, Alexopoulos EC, Belezos T, Panagiotopoulou E: Lymphocyte to monocyte ratio as a screening tool for influenza. Plos Curr. 2010, 2: RRN1154-
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.
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.
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.
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.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/314/prepub
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.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
About this article
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). https://doi.org/10.1186/1471-2334-14-314