Changes in epigenetic profiles throughout early childhood and their relationship to the response to pneumococcal vaccination
Clinical Epigenetics volume 13, Article number: 29 (2021)
Pneumococcal infections are a major cause of morbidity and mortality in young children and immaturity of the immune system partly underlies poor vaccine responses seen in the young. Emerging evidence suggests a key role for epigenetics in the maturation and regulation of the immune system in health and disease. The study aimed to investigate epigenetic changes in early life and to understand the relationship between the epigenome and antigen-specific antibody responses to pneumococcal vaccination.
The epigenetic profiles from 24 healthy children were analyzed at 12 months prior to a booster dose of the 13-valent pneumococcal conjugate vaccine (PCV-13), and at 24 months of age, using the Illumina Methylation 450 K assay and assessed for differences over time and between high and low vaccine responders.
Our analysis revealed 721 significantly differentially methylated positions between 12 and 24 months (FDR < 0.01), with significant enrichment in pathways involved in the regulation of cell–cell adhesion and T cell activation. Comparing high and low vaccine responders, we identified differentially methylated CpG sites (P value < 0.01) associated with HLA-DPB1 and IL6.
These data imply that epigenetic changes that occur during early childhood may be associated with antigen-specific antibody responses to pneumococcal vaccines.
Young children are at higher risk of developing life-threatening diseases compared with adults, which is thought to be due to the immaturity of their immune system. They also show reduced immune responses to vaccines . Pneumococcal infections are a major public health problem worldwide causing high morbidity and mortality in young children who suffer pneumonia, meningitis and septicemia, leading up to 1 million deaths in children under the age of five years [2, 3]. Despite many years of research, the specific characteristics of the immune system in early childhood that underpins increased susceptibility to infection has still not been fully elucidated . Furthermore, there is great variability in the immune response to vaccination in children, with some individuals achieving more than 100-fold higher antibody levels post-immunization than others [5, 6]. The reason for this huge variation is not understood, and recent studies have begun to reveal the genetic determinants of immune responses to vaccines [7, 8].
The term epigenetics or ‘the epigenome’ describes the summary of heritable changes in expression (i.e., function) of genes that are not encoded in the DNA sequence and form a mechanism whereby past environmental events may be encoded and modify future cellular responses. Epigenetic mechanisms have repressive and permissive effects on gene expression and thus play a key role in cellular differentiation. DNA methylation forms the archetypal, heritable, epigenetic mark and involves the chemical modification of cytosine and cytosine:guanine dinucleotides, resulting in profound, typically repressive, effects . DNA methylation patterns can be measured across the whole genome covering 99% of currently annotated genes in a high-throughput manner . They differ between individuals, even between identical twins, and their change over time is influenced by intrauterine and postnatal environmental factors such as nutrition, toxins or drugs and illnesses . Recent studies have demonstrated that the levels of DNA methylation of certain CpG sites are also highly correlated with human chronological age . Moreover, these DNA methylation features have been proposed as a surrogate of the “biological age,” which may actually be more relevant than chronological age in predicting age-related health risks . Different conditions have been reported to accelerate age-associated DNA methylation profiles such as HIV infection and Down syndrome [14, 15].
There is increasing evidence demonstrating a key role of epigenetics in regulation of human immunity in health and disease [16, 17]. Consistent with this, DNA methylation changes during early childhood have been observed in genes implicated in inflammatory processes, encoded histone modifiers and chromatin remodeling factors . However, the effects of these changes on the immune system are unknown. Importantly, there is limited information on the relationship between epigenetic changes in infancy and the response to vaccination. Here, we sought to investigate epigenetic changes in early life in a well-studied cohort of healthy children, focusing on the characterization of the relationship between epigenetics and age in paired samples collected at 12 and 24 months of age, and the role of epigenetics on the strength of antigen-specific antibody responses to pneumococcal vaccination.
Twenty-four children (Table 1) who had received the 13-valent pneumococcal conjugate vaccine (PCV-13) at 12 months of age, as part of a previously published clinical trial [19,20,21], were selected based on their aggregated IgG antibody response and categorized into high and low responders (Fig. 1 and Additional file 1: Figure 1). Of note, this categorization was also valid for post-booster serotype-specific opsonophagocytic assay (OPA) titers and memory B cells responses (Additional file 1: Figures 2 and 3), therefore truly grouping individuals into high and low vaccine responders. Paired blood samples taken at 12 and 24 months of age were assessed by DNA epigenetic microarray analysis.
Widespread changes in the epigenome occur between 12 and 24 months of life
The first step of our analysis focused on the methylation changes between the two timepoints (12 vs. 24 months). After correcting our linear model for cell-type composition, we found 721 significant DMPs, associated with 421 unique genes. Considerable inflation could be observed in the QQ plot (Additional file 1: Figure 4), reflecting substantial age-related changes to the blood epigenome early in life. Of the significant DMPs, 314 (43.6%) were hypomethylated with the remaining 407 (56.4%) being hypermethylated, in samples collected at 12 months of age (Additional file 7: Table 1). The results without adjusting for cell compositions are shown in Additional file 8: Table 2. There was a clear separation between samples from 12 and 24 months of age using PCA and hierarchical clustering of the 721 DMPs (Fig. 2).
Further analyzing the 721 DMPs, we observed several DMRs where multiple adjacent CpGs probes were found together showing differential methylation between the two different timepoints. For instance, 10 probes were associated with nuclear factor I X (NFIX), and 12 CpGs were observed within the same island region, chr6:32118101–32118544, of the proline-rich transmembrane protein 1 (PRRT1). Hypomethylated CpGs were mainly found in CpG islands, whereas open sea, shore and shelf regions were mostly hypermethylated (Additional file 5: Figure 5a). Most DMPs were found in the gene body, TSS1500 and 5′ UTRs and dominated by hypomethylation patterns (Additional file 1: Figure 5b).
We further analyzed DMPs, which showed the highest age-dependent differences in methylation by applying a delta beta (difference in the mean β value of each group) threshold of ≥ 0.10. In this way, we detected 24 DMPs, considered as the DMPs most separating the two age groups (Table 2); of those, 18 CpGs exhibited a higher methylation level at 24 months, while the remaining 6 had a lower methylation pattern at 24 months.
Age-associated changes in the epigenome are enriched for pathways linked to T cell regulation and activation
To investigate which molecular functions or biological processes were associated with the DMPs between 12 and 24 months of age, gene set enrichment analysis was performed. We detected significant enrichment (FDR < 0.05) in 151 GO pathways, with the most significant pathways being related to regulation of T cell activation, cell–cell adhesion and regulation of cytokine production (Fig. 3; Additional file 9: Table 3).
Among DMPs, we observed 86 CpGs associated with 69 unique genes within immune system processes according to GO terms. Of these, 53 CpGs showed a decrease in methylation with age, while the other 32 exhibited an increase in methylation from 12 to 24 months (Fig. 4). Some of the genes were represented by multiple DMPs, such as HOXB6, NOD2 and LAG3, in the latter we found both hypo- and hypermethylated CpGs.
In order to understand if the different responsiveness to vaccination was associated with the “biological age” of the subjects, we used the epigenetic age estimator developed by Horvath in 2013  to calculate biological age based on methylation patterns. First, we applied age predictors to the results obtained from samples at the two time points and found that chronological and predicted biological ages were correlated (Fig. 5a). We did not observe significant differences in predicted biological age between high and low vaccine responders (Fig. 5b).
Hypomethylation of HLA-DPB1 and hypermethylation of IL-6 is associated with more robust responses to pneumococcal infant vaccination
We next assessed whether differences in antibody responses against vaccine antigens were associated with methylation patterns by comparing pre-vaccination blood epigenome profiles of high responders to those of low responders. After applying a FDR < 0.01 correction, no statistically significant DMPs were detected. However, using a less stringent statistical threshold (uncorrected P value < 0.01), we found 4067 differentially methylated CpG sites (Additional file 10: Table 4) associated with 2797 unique genes that distinguished high from low responders in a PCA (Fig. 6a) and an unsupervised clustering method (Fig. 6b). Similar results were obtained for a group comparison at 24 months (Additional file 6: Figure 6).
We further investigated among the 4067 DMPs those which showed the highest difference in methylation between high and low vaccine responders (mean β value) and assessed regions containing multiple DMPs in genes involved in immune system processes. Among the top DMPs identified, there were 4 hypomethylated positions in high vaccine responders within the chr6:33048416–33048814 island region of the HLA class II histocompatibility antigen DP(W2) beta chain (HLA-DPB1). There were also two CpGs associated with interleukin 6 (IL-6) that showed hypermethylation in high compared with low vaccine responders (Fig. 6c).
In the present study, we demonstrated age-associated DNA methylation changes between the first and second year of life that were largely related to T cell regulation and activation. In addition, using detailed antibody measurements following infant pneumococcal conjugate booster vaccination in a clinical vaccine trial, we found hypomethylation of HLA-DPB1 and hypermethylation of IL6 to be correlated with a stronger response to pneumococcal vaccination.
A total of 721 CpGs were found to be significantly differentially methylated between 12 and 24 months of age. Multiple DMPs within several DMRs showed similar methylation patterns suggesting a general epigenetic remodeling of these genetic loci. For example, within the chr6:32118101–32118544 locus, there were twelve DMPs that exhibited a significantly higher methylation pattern at 24 compared with 12 months of age. These DMPs are localized within the body of PRRT1 and could potentially alter the expression of PRRT1, as previously demonstrated for other hypermethylated positions in genes outside of the promoter area . Little is known about the function of PRRT1 although aberrant methylation of this gene has been related to neurodevelopmental disorders  and to hepatic tumorigenesis . Furthermore, we observed a total of 118 CpGs that were previously reported to be altered between 3 and 60 months in blood leukocytes . In agreement with this previous work, we observed 57 CpGs exhibiting an increase in methylation between 12 and 24 months of age, and 60 CpGs showing a decreased methylation with age .
Through an enrichment pathway analysis of the age dependent DMPs, we observed 86 CpGs that were associated with 69 genes involved in several immune system processes. Genes that were hypomethylated at 24 compared with 12 months included IL12A, previously shown to be repressed at chromatin level in neonatal mononuclear compared with adult cells [26,27,28], as well as BCL10 and NOD2, also previously demonstrated to decrease methylation with age . In contrast, CpGs with an increase in methylation between 12 and 24 months included genes involved in the regulation of innate immune responses such as tripartite motif containing 15 (TRIM15) , and the PML Nuclear Body Scaffold (PML) . A particular pattern of methylation was observed for lymphocyte activating 3 (LAG3), also required for T cell regulation . LAG3 was hypomethylated at two positions within the chr12:6882855–6883184 region and hypermethylated in an open sea CpG at 24 months, suggesting that diverging methylation patterns may regulate expression of this gene.
There was a clear separation of predicted biological age between samples taken at the two study time points. Predicted biological age was moderately higher than chronological age, most likely due to the small range of errors for whole blood samples that the age predictor exhibits . However, predicted biological age was not significantly different between high and low vaccine responders at any time point suggesting no differences attributable to “epigenetic age” between the groups.
The second part of our research study focused on the relationship between the epigenome at 12 months of age before immunization and the strength of the antigen-specific IgG response to vaccination. We identified a high number of CpGs that separated high from low vaccine responders using a less stringent statistical threshold (uncorrected P value < 0.01). These CpGs included multiple probes mapping to genes involved in immune system processes and were associated with large differences in methylation between the two groups. In this way, we observed two particularly interesting immune genes (HLA-DPB1 and IL6) that contained multiple DMPs differing between high and low vaccine responders. Hypomethylation of HLA-DPB1 as found in high compared with low vaccine responders may indeed strengthen the overall immune response as previously demonstrated by the contribution of HLA to cellular host immunity in response to anthrax vaccine adsorbed ; rubella vaccine ; hepatitis B vaccine; and the inactivated Japanese encephalitis vaccine . Additional evidence connecting HLA genotype and expression to vaccine immune responses was found in measles seropositive and seronegative individuals . Two CpGs in the body of IL6 showed higher methylation in high compared with low vaccine responders. These CpGs in open sea positions (i.e., isolated in the genome ) are considered to be “predicted enhancer elements” , therefore possibly impacting IL6 expression. Il-6 is a proinflammatory cytokine involved in a variety of immune processes and also plays an important role in mediating innate and adaptive immune responses ; therefore, differential expression of IL6 may impact vaccine immune responses. However, in contrast to our finding of hypermethylated (and possibly reduced expression of) IL6 in high vaccine responders, previous work suggested that IL-6 promotes stimulation of humoral and cellular immune responses [40,41,42,43].
There are several limitations to our study. First, several samples had to be excluded due to unavailability of genetic material, consent for genetic analysis or missing serum, which may have affected the results. In addition, high and low responders were categorized using the aggregated serotype-specific IgG response 1-month post-booster vaccination, which may not accurately reflect the strength of the overall immune response to immunization. Reassuringly, however, high vaccine responders showed similar patterns in OPA and memory B-cell responses, indicating that this categorization also seems to hold up in different measurements of the immune response.
In the present study, we have studied changes in methylation patterns of healthy children between 12 and 24 months of age and found that these involved regions dominated by immune pathways such as T cell regulation and activation. In addition, we found a number of CpGs associated with the strength of the antigen-specific IgG antibody response to a pneumococcal conjugate booster vaccine including differences in methylation patterns in HLA-DPB1 and IL6. These latter findings suggest that epigenetic patterns may influence antibody responses in young children and may potentially add to their susceptibility to infection.
Study subjects and sample selection
In a previous clinical trial, 178 children were randomized to receive a booster dose of either the 10- or the 13-valent pneumococcal conjugate vaccine (PCV-13) at 12 months of age following receipt of PCV-13 at 2 and 4 months of age [19,20,21]. Ethical approval for this clinical vaccine trial was obtained from the Oxfordshire Research Ethics Committee (Reference number 11/SC/0473). Ethical permission to store and study genetic material was obtained prior to enrolment into the study, in accordance with the Oxford Vaccine Centre Biobank (10/H0504/25). Detailed immune responses were measured before and at 1 and 12 months after the booster. At the three study time points, anti-polysaccharide serum Immunoglobulin G (IgG) and opsonophagocytic assay (OPA) titers to all PCV-13 serotypes were determined, and antigen-specific memory B cell responses were measured for selected pneumococcal polysaccharides. For the purpose of this study, 24 individuals who had received a PCV-13 booster dose at 12 months of age (n = 90) (baseline group) were selected based on their aggregated IgG antibody response and categorized into high and low responders. For this categorization, the strength of the polysaccharide-specific IgG response was ranked for each PCV-13 recipient and vaccine serotype at 1 month following the booster. The overall IgG antibody response was then calculated by determining the median rank for all serotypes. After excluding participants with missing serum or DNA samples, 38 individuals remained, of whom the 12 high and 12 low responders were selected (Fig. 1). This categorization was also valid for post-booster serotype-specific OPA and memory B cells responses. We used Infinium Methylation 450 K to profile (1) the epigenetic differences between two timepoints: 12 months (baseline group) and 24 months and (2) to elucidate the methylation differences between high and low PCV-13 booster vaccine responders, in order to identify pathways that show significant changes over time and those that are associated with a more robust immune response.
Polysaccharide-specific antibody measurements were performed as previously described [19, 44]. DNA was extracted from blood clots collected from the same participants at 12 and 24 months of age. Whole blood DNA was bisulfite converted using the EZ-96 DNA Bisulfite Zymo Research conversion protocol according to the manufacturer’s instructions. The treated DNA was then hybridized to the Illumina Infinium Methylation450K BeadChip, and the array was imaged using the Illumina iScan system in which the percent methylation state of each CpG site was quantified for the entire study group.
The epigenetic profile from study participants was measured at both time points using the Illumina Methylation 450 K assay covering > 485,000 methylation sites per sample at single-nucleotide resolution. DNA methylation analyses were conducted using R software (v.3.6.1) and several Bioconductor packages. Raw intensity files were preprocessed and transformed into β and M values using the R package minfi. The intensities of β values were estimated from the intensity ratio of the methylated signals over the total (methylated and unmethylated) signals for each site that represents the percentage of methylation at a given cytosine for an individual across the blood cells. The β value ranges between 0 (position not methylated in any of the cells) and 1 (position methylated in all cells in the sample). The M value instead is defined as the log2 ratio of the intensities of methylated vs. unmethylated probes and indicates hypermethylation (if positive) and hypomethylation (if negative). M values close to 0 represent a similar intensity between methylated and unmethylated probes.
All raw data passed quality control and normalized using the function preprocessQuantile of the minfi package. Subsequently, to avoid artifactual data, probes with a detection P value > 0.01, probes located in sex chromosomes, probes in sites containing SNPs or with a minor allele frequency < 5% and probes known to have cross-reaction were removed. This filtering removed 28,999 leaving 456,513 probes for subsequent analysis. As the DNA samples were obtained from peripheral blood and because of the extreme cell-type specificity of DNA methylation , cell composition heterogeneity was evaluated  using the package FlowSorted.Blood.450 K and the estimateCellCounts function. Batch effect was removed using Combat from SVA package.
Identification of differentially methylated probes/positions (DMPs) between the study timepoint and high- and low-responders to the vaccine was performed using limma package  fitting a linear regression model to each CpGs site. M-values are more homoscedastic than β-values and were thus used for all statistical analyses because of their better detection sensitivity . However, we report β-values in text and figures to facilitate interpretation. Through limma package, it was possible to include covariates in the specification of the model (M values ~ Group + Individuals + Responder + Sex + Cell Composition). Output of the linear model was further analyzed using an empirical Bayes method to moderate standard errors. P values were corrected for multiple testing using the Benjamini–Hochberg false discovery rate (FDR) method , and DMPs were defined as all positions with an FDR P value < 0.01. CpGs were annotated using Bioconductor package lluminaHumanMethylation450kanno.ilmn12.hg19 (v1.2) to genomic regions classified in regions 200 bp upstream of the transcription start site (TSS; TSS200), 1500 bp upstream of the TSS (TSS1500), the 5′ and 3′ untranslated regions (UTRs), 1st exon and gene body. Regions are further subdivided according to their relation to CpG-rich DNA sequences, the so-called CpG islands (CGI), in CpG N or S shelves (within ~ 4 kb north or south from CGI), CpG N or S shores (within ~ 2 kb north or south from CGI) and open sea (region > 4 kb from CGI). Differentially methylated regions (DMRs) refer to genomic regions with multiple differentially methylated positions close together and with different methylation patterns across samples.
Principal component analysis (PCA) was performed to visualize global epigenome patterns comparing 12- and 24-month samples as well as samples from high and low responders. Pathway analysis was performed employing the methylglm function of the methylGSA package . This function adjusts for the number of CpGs by implementing a logistic regression model using a functional class scoring method; all CpGs ranked by differential methylation test P values were considered in gene set testing since each CpGs could be of interest. Pathways enrichment was considered significant if false discovery rate (FDR) < 0.05.
Availability of data and materials
Methylation data have been deposited on the Gene Expression Omnibus (GSE165170; https://www.ncbi.nlm.nih.gov/geo/).
Blanchard-Rohner G, Pollard AJ. Long-term protection after immunization with protein–polysaccharide conjugate vaccines in infancy. Expert Rev Vaccines. 2011;10(5):673–84.
O’Brien KL, Wolfson LJ, Watt JP, Henkle E, Deloria-Knoll M, McCall N, et al. Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. The Lancet. 2009;374(9693):893–902.
World Health Organization. Pneumococcal conjugate vaccine for childhood immunization—WHO position paper. Wkly Epidemiol Rec. 2007;82(12):93–104.
Kollmann TR, Kampmann B, Mazmanian SK, Marchant A, Levy O. Protecting the newborn and young infant from infectious diseases: lessons from immune ontogeny. Immunity. 2017;46(3):350–63.
Westerink MJ, Schroeder HW Jr, Nahm MH. Immune responses to pneumococcal vaccines in children and adults: rationale for age-specific vaccination. Aging Dis. 2012;3(1):51.
Lipsitch M, Li LM, Patterson S, Trammel J, Juergens C, Gruber WC, et al. Serotype-specific immune responses to pneumococcal conjugate vaccine among children are significantly correlated by individual: analysis of randomized controlled trial data. Vaccine. 2018;36(4):473–8.
Mentzer AJ, O’Connor D, Pollard AJ, Hill AV. Searching for the human genetic factors standing in the way of universally effective vaccines. Philos Trans R Soc Lond B Biol Sci. 2015;370(1671):20140341.
O’Connor D, Png E, Khor CC, Snape MD, Hill AVS, van der Klis F, et al. Common genetic variations associated with the persistence of immunity following childhood immunization. Cell Rep. 2019;27(11):3241–53.
Jones PA, Takai D. The role of DNA methylation in mammalian epigenetics. Science. 2001;293(5532):1068–70.
Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, et al. High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006;16(3):383–93.
Perera F, Herbstman J. Prenatal environmental exposures, epigenetics, and disease. Reprod Toxicol. 2011;31(3):363–73.
Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, Sommer IE, et al. The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS ONE. 2009;4(8):e6767.
Bell CG, Lowe R, Adams PD, Baccarelli AA, Beck S, Bell JT, et al. DNA methylation aging clocks: challenges and recommendations. Genome Biol. 2019;20(1):249.
Boulias K, Lieberman J, Greer EL. An epigenetic clock measures accelerated aging in treated HIV infection. Mol Cell. 2016;62(2):153–5.
Horvath S, Garagnani P, Bacalini MG, Pirazzini C, Salvioli S, Gentilini D, et al. Accelerated epigenetic aging in down syndrome. Aging Cell. 2015;14(3):491–5.
Fernandez-Morera JL, Calvanese V, Rodriguez-Rodero S, Menendez-Torre E, Fraga MF. Epigenetic regulation of the immune system in health and disease. Tissue Antigens. 2010;76(6):431–9.
Kondilis-Mangum HD, Wade PA. Epigenetics and the adaptive immune response. Mol Asp Med. 2013;34(4):813–25.
Acevedo N, Reinius LE, Vitezic M, Fortino V, Soderhall C, Honkanen H, et al. Age-associated DNA methylation changes in immune genes, histone modifiers and chromatin remodeling factors within 5 years after birth in human blood leukocytes. Clin Epigenet. 2015;7:34.
Truck J, Jawad S, Goldblatt D, Roalfe L, Snape MD, Voysey M, et al. The antibody response following a booster with either a 10- or 13-valent pneumococcal conjugate vaccine in toddlers primed with a 13-valent pneumococcal conjugate vaccine in early infancy. Pediatr Infect Dis J. 2016;35(7):787–93.
Truck J, Kelly S, Jawad S, Snape MD, Voysey M, Pollard AJ. Differences in immunization site pain in toddlers vaccinated with either the 10- or the 13-valent pneumococcal conjugate vaccine. Pediatr Infect Dis J. 2018;37(4):e103–6.
Truck J, Snape MD, Tatangeli F, Voysey M, Yu LM, Faust SN, et al. Pneumococcal serotype-specific antibodies persist through early childhood after infant immunization: follow-up from a randomized controlled trial. PLoS ONE. 2014;9(3):e91413.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):3156.
Arechederra M, Daian F, Yim A, Bazai SK, Richelme S, Dono R, et al. Hypermethylation of gene body CpG islands predicts high dosage of functional oncogenes in liver cancer. Nat Commun. 2018;9(1):3164.
Ladd-Acosta C, Hansen KD, Briem E, Fallin MD, Kaufmann WE, Feinberg AP. Common DNA methylation alterations in multiple brain regions in autism. Mol Psych. 2014;19(8):862–71.
Maschietto M, Rodrigues TC, Kashiwabara AY, de Araujo ÉSS, Aguiar TFM, da Costa CML, et al. DNA methylation landscape of hepatoblastomas reveals arrest at early stages of liver differentiation and cancer-related alterations. Oncotarget. 2017;8(58):97871.
Henry CJ, Ornelles DA, Mitchell LM, Brzoza-Lewis KL, Hiltbold EM. IL-12 produced by dendritic cells augments CD8+ T cell activation through the production of the chemokines CCL1 and CCL17. J Immunol. 2008;181(12):8576–84.
Martinez GJ, Zhang Z, Chung Y, Reynolds JM, Lin X, Jetten AM, et al. Smad3 differentially regulates the induction of regulatory and inflammatory T cell differentiation. J Biol Chem. 2009;284(51):35283–6.
Scharschmidt E, Wegener E, Heissmeyer V, Rao A, Krappmann D. Degradation of Bcl10 induced by T-cell activation negatively regulates NF-κB signaling. Mol Cell Biol. 2004;24(9):3860–73.
DeVries A, Wlasiuk G, Miller SJ, Bosco A, Stern DA, Lohman IC, et al. Epigenome-wide analysis links SMAD3 methylation at birth to asthma in children of asthmatic mothers. J Allergy Clin Immunol. 2017;140(2):534–42.
Kawai T, Akira S. Regulation of innate immune signalling pathways by the tripartite motif (TRIM) family proteins. EMBO Mol Med. 2011;3(9):513–27.
Lunardi A, Gaboli M, Giorgio M, Rivi R, Bygrave A, Antoniou M, et al. A role for PML in innate immunity. Genes Cancer. 2011;2(1):10–9.
Huang CT, Workman CJ, Flies D, Pan X, Marson AL, Zhou G, et al. Role of LAG-3 in regulatory T cells. Immunity. 2004;21(4):503–13.
Ovsyannikova IG, Pankratz VS, Vierkant RA, Pajewski NM, Quinn CP, Kaslow RA, et al. Human leukocyte antigens and cellular immune responses to anthrax vaccine adsorbed. Infect Immun. 2013;81(7):2584–91.
Lambert ND, Haralambieva IH, Kennedy RB, Ovsyannikova IG, Pankratz VS, Poland GA. Polymorphisms in HLA-DPB1 are associated with differences in rubella virus-specific humoral immunity after vaccination. J Infect Dis. 2015;211(6):898–905.
Yao Y, Yang H, Shi L, Liu S, Li C, Chen J, et al. HLA class II genes HLA-DRB1, HLA-DPB1, and HLA-DQB1 are associated with the antibody response to inactivated Japanese encephalitis vaccine. Front Immunol. 2019;10:428.
Dhiman N, Bonilla RG, Jacobson RM, O’Kane D, Poland GA. Differential HLA gene expression in measles vaccine seropositive and seronegative subjects: a pilot study. Scand J Infect Dis. 2003;35(5):332–6.
Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, et al. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011;6(6):692–702.
Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.
Jones SA. Directing transition from innate to acquired immunity: defining a role for IL-6. J Immunol. 2005;175(6):3463–8.
Su L-K, Yu F, Li Z-F, Zeng C, Xu Q-A, Fan M-W. Intranasal co-delivery of IL-6 gene enhances the immunogenicity of anti-caries DNA vaccine. Acta Pharmacol Sin. 2014;35(5):592–8.
Su B, Wang J, Wang X, Jin H, Zhao G, Ding Z, et al. The effects of IL-6 and TNF-α as molecular adjuvants on immune responses to FMDV and maturation of dendritic cells by DNA vaccination. Vaccine. 2008;26(40):5111–22.
Kurtz SL, Foreman O, Bosio CM, Anver MR, Elkins KL. Interleukin-6 is essential for primary resistance to Francisella tularensis live vaccine strain infection. Infect Immun. 2013;81(2):585–97.
Velazquez-Salinas L, Verdugo-Rodriguez A, Rodriguez LL, Borca MV. The role of interleukin 6 during viral infections. Front Microbiol. 2019;10:1057.
Truck J, Mitchell R, Jawad S, Clutterbuck EA, Snape MD, Kelly DF, et al. Divergent memory B cell responses in a mixed infant pneumococcal conjugate vaccine schedule. Pediatr Infect Dis J. 2017;36(5):e130–5.
Ziller MJ, Gu H, Muller F, Donaghey J, Tsai LT, Kohlbacher O, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013;500(7463):477–81.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13:86.
Smyth GK. Limma: linear models for microarray data. Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer; 2005. p. 397–420.
Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinform. 2010;11:587.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodol). 1995;57(1):289–300.
Ren X, Kuan PF. methylGSA: a Bioconductor package and Shiny app for DNA methylation data length bias adjustment in gene set testing. Bioinformatics. 2019;35(11):1958–9.
We thank the High-Throughput Genomics Group at the Wellcome Trust Centre for Human Genetics (Funded by Wellcome Trust Grant reference is 090532/Z/09/Z and MRC Hub Grant G0900747 91070) for the generation of data.
Funding to perform methylation analysis was obtained through a startup research support scheme from the medical research fund of the University of Oxford. During analysis and manuscript writing, JT was supported by an Ambizione-SCORE grant of the Swiss National Science Foundation (PZ00P3_161147; PZ00P3_183777). BPF is funded by a Wellcome Intermediate Clinical Fellowship (201488/Z/16/Z). The research was supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC). The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health.
Ethics approval and consent to participate
Ethical approval for the clinical vaccine trial was obtained from the Oxfordshire Research Ethics Committee (Reference number 11/SC/0473). Ethical permission to store and study genetic material was obtained prior to enrolment into the study, in accordance with the Oxford Vaccine Centre Biobank (10/H0504/25).
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Boxplots of IgG antibody concentrations against PCV-13 serotypes (serotype indicated above each box) in low and high responders (n = 12 each; based on aggregated IgG antibody concentration 1 month post-vaccination) as well as middle responders (n = 14) and children without available DNA (n = 36).
. Boxplots of OPA antibody titers against PCV-13 serotypes (serotype indicated above each box) in low and high responders (n = 12 each; based on aggregated IgG antibody concentration 1 month post-vaccination) as well as middle responders (n = 14) and children without available DNA (n = 36).
. Boxplots of IgG antibody-secreting cell frequencies against selected PCV-13 serotypes (serotype indicated above each box) in low and high responders (n = 12 each; based on aggregated IgG antibody concentration 1 month post-vaccination) as well as middle responders (n = 14) and children without available DNA (n = 36).
. QQ plot of the observed p-values against the expected p-values.
. Distributions of the 721 age-dependent DMPs (a) in the island context (island, shore, shelf, open sea) and (b) in the genomic regions (1st exon, 3’ UTR, 5’ UTR, body, TSS1500, TSS200).
. (a) Principal component analysis and (b) unsupervised clustering and heatmap of 5233 differentially methylated CpG sites detected using a threshold of an uncorrected p-value < 0.01 at 24 months of age and grouping by high and low vaccine responders.
. Significantly differentially methylated positions (FDR adjusted p-value <0.01) between 12 and 24 months of age, adjusted for individual, responder, gender and cell composition.
. Differentially methylated positions found to be significant (FDR adjusted p-value < 0.01) between 12- and 24-months groups, without adjusting the model for cell composition.
. List of 152 significant enriched pathways in GO (FDR < 0.05).
. Differentially methylated positions (p-value < 0.01) found to be differentially methylated between high and low responders at 12 months of age, adjusted for gender and cell composition.
About this article
Cite this article
Pischedda, S., O’Connor, D., Fairfax, B.P. et al. Changes in epigenetic profiles throughout early childhood and their relationship to the response to pneumococcal vaccination. Clin Epigenet 13, 29 (2021). https://doi.org/10.1186/s13148-021-01012-w