Respiratory infection- and asthma-prone, low vaccine responder children demonstrate distinct mononuclear cell DNA methylation pathways

Abstract

Background: Infants with frequent viral and bacterial respiratory infections exhibit compromised immunity to routine immunizations. They are also more likely to develop chronic respiratory diseases in later childhood. This study investigated the feasibility of epigenetic profiling to reveal endotype-specific molecular pathways with potential for early identification and immuno-modulation. Peripheral blood mononuclear cells from respiratory infection allergy/asthma-prone (IAP) infants and non-infection allergy/asthma prone (NIAP) were retrospectively selected for genome-wide DNA methylation and single nucleotide polymorphism analysis. The IAP infants were enriched for the low vaccine responsiveness (LVR) phenotype (Fisher's exact p-value = 0.02).

Results: An endotype signature of 813 differentially methylated regions (DMRs) comprising 238 lead CpG associations (FDR < 0.05) emerged, implicating pathways related to asthma, mucin production, antigen presentation and inflammasome activation. Allelic variation explained only a minor portion of this signature. Stimulation of mononuclear cells with monophosphoryl lipid A (MPL), a TLR agonist, partially reversed this signature at a subset of CpGs, suggesting the potential for epigenetic remodeling.

Conclusions: This proof-of-concept study establishes a foundation for precision endotyping of IAP children and highlights the potential for immune modulation strategies using adjuvants for future investigation.

Fig. 1

Epigenome-wide association study of IAP and NIAP infants demonstrates differentially methylated regions enriched with immune response genes. A Volcano plot of differentially methylated regions significantly associated with IAP status. X-axis represents the average methylation change from NIAP group (delta methylation ratio), and the y-axis shows the log region size in base pairs. Points are colored according to maximum observed methylation difference in region. N = 30 B Circular visualization of significant SNP-CpG association pairs showing chromosomal location of associated markers. C Dot chart of summary statistics from gene ontology enrichment analysis of the infection-prone associated methylation signature highlighting significantly enriched pathways. The y-axis shows gene set nomenclature, point size reflects total count of differentially methylated genes, points are ranked by the percentage of genes in the set covered, and colored according to false discovery rate adjusted P-value. Diff Differentially, DMRs differentially methylated regions, IAP infection allergy/asthma prone, max maximum, bp base pairs, P.DE p-value for over-representation of the GO or KEGG term, Coverage percentage of genes identified within specified pathway

Fig. 2

Stimulation with MPL adjuvant modifies a subset of endotype-associated CpGs. A Dot chart of summary statistics from gene ontology enrichment analysis for the MPL-associated differentially methylated regions highlighting significantly enriched pathways (n = 30). The y-axis shows gene set nomenclature, point size reflects total count of differentially methylated genes, points are ranked by the percentage of genes in the set covered, and colored according to false discovery rate adjusted P-value. B Venn diagram of overlapping regions showing counts of unique and overlapping regions. C Bland–Altman visualization of methylation difference after MPL stimulation for 183 MPL-sensitive CpGs. Y-values represent the difference in log2 fold-change (NIAP – IAP) between PBS and MPL condition. X-axis shows mean methylation value for each CpG. Upper and lower 95% confidence intervals are shown as red lines. D Boxplot of median and quartile range for the four most significantly ranked endotype-associated CpGs modified by MPL treatment. DMRs differentially methylated regions, MPL monophosphoryl lipid A, P.DE p-value for over-representation of the GO or KEGG term, Coverage percentage of genes identified within specified pathway. Hypermeth hypermethylated regions. Hypometh hypomethylated regions.*** = P < 0.001, ** = P < 0.05–0.05, ns = P > 0.05, two-sided t-test