Skip to main content

Maternal asthma and newborn DNA methylation

Abstract

Background

Prenatal exposure to maternal asthma may influence DNA methylation patterns in offspring, potentially affecting their susceptibility to later diseases including asthma.

Objective

To investigate the relationship between parental asthma and newborn blood DNA methylation.

Methods

Epigenome-wide association analyses were conducted in 13 cohorts on 7433 newborns with blood methylation data from the Illumina450K or EPIC array. We used fixed effects meta-analyses to identify differentially methylated CpGs (DMCs) and comb-p to identify differentially methylated regions (DMRs) associated with maternal asthma during pregnancy and maternal asthma ever. Paternal asthma was analyzed for comparison. Models were adjusted for covariates and cell-type composition. We examined whether implicated sites related to gene expression analyses in publicly available data for childhood blood and adult lung.

Results

We identified 27 CpGs associated with maternal asthma during pregnancy at False Discovery Rate < 0.05 but none for maternal asthma ever. Two distinct CpGs were associated with paternal asthma. We observed 5 DMRs associated with maternal asthma during pregnancy 3 associated with maternal asthma ever and 13 DMRs associated with paternal asthma. Gene expression analysis using data in blood from 832 children and lung from 424 adults showed associations between identified DMCs using maternal asthma and expression of several genes, including HLA genes and HOXA5, previously implicated in asthma or lung function.

Conclusion

Parental asthma, especially maternal asthma during pregnancy, may be associated with alterations in newborn DNA methylation. These findings might shed light on underlying mechanisms for asthma susceptibility.

Introduction

Asthma is the most common chronic disease in childhood, leading to decreased quality of life for affected families and large costs to society due reduced productivity and missed school days [1]. Childhood asthma is a highly heritable disease with parental asthma being the strongest known risk factor for asthma in the offspring and heritability estimates between 35 and 95% [2]. Nevertheless, a significant portion of the asthma risk in offspring remains unexplained by genetic factors alone [3] and early life risk factors related to the perinatal environment [4] and comorbidity factors [5] have been reported. Information on epigenetic markers such as DNA methylation collected at birth, might provide insight into the mechanisms of prenatal programming of childhood asthma. Furthermore, clinical studies have suggested that there are sex differences in vulnerability to asthma. Such differences may be related to sex-specific DNA methylation levels [6, 7]. Therefore, we speculated that sex-specific effects in the association of parental asthma with offspring methylation might be observed.

Large-scale meta-analyses of DNA methylation have demonstrated that epigenetic alterations at birth and later in childhood are associated with childhood asthma, suggesting a potential role of epigenetic mechanisms in asthma development [8,9,10,11]. Some studies have found that maternal asthma is a stronger risk factor for childhood asthma than paternal asthma [12], implying a potential prenatal programming effect on the oocyte or fetus from a maternal asthma-associated milieu, possibly through epigenetic mechanisms. However, the potential role of parental asthma and the difference between maternal and paternal effects on offspring methylation has not been well explored.

Here, our main objective was to perform a large-scale meta-analysis of maternal asthma active during pregnancy, hypothesizing that epigenetic effects could be stronger from intrauterine “exposure” to active maternal asthma. We additionally examined associations between DNA methylation and maternal asthma history in relation to DNA methylation from various blood sources at birth in more than 7000 newborns from 13 cohorts. Given the interest in paternal exposures on offspring health, we similarly analyzed associations with paternal asthma, hypothesizing that a higher number of differentially methylated CpGs in the maternal compared to the paternal analysis would further support a direct exposure effect on offspring DNA methylation. Further, given the suggestions in the asthma literature about sex-specific effects of impacts of parental asthma in offspring, we investigated the role of offspring sex in each of these analyses by performing sex-stratified analyses and interaction testing. We conducted our analyses by testing associations of DNA methylation from various blood sources at both individual CpG sites and at differently methylated regions (DMR). We evaluated the potential functional impact of findings by integrating gene expression data from blood and lung tissue.

Methods

Study design

The study used DNA methylation array data from cohorts within the pregnancy and childhood epigenetics (PACE) Consortium. PACE is an international consortium of cohorts with DNA methylation data available at birth, in childhood and/or in adolescence using either the Illumina450K array or Illumina EPIC arrays [13, 14]. We evaluated a maternal diagnosis of asthma (ever and during pregnancy) in relation to DNA methylation data from various blood sources from newborns in a total of 13 cohorts (Avon Longitudinal Study of Parents and Children [ALSPAC], Children’s Health Study [CHS], Drakenstein Child Health Study [DCHS], Etude des Déterminants pré et post natals du développement et de la santé de l’Enfant [EDEN], Generation R [GENR], INfancia y Medio Ambiente [INMA], Isle of Wight 3rd Generation Birth Cohort [IOWF2], Lifestyle and environmental factors and their influence on the newborn allergy risk [LiNA], Father and Child Cohort Study the Multigenerational Familial and Environmental Risk for Autism [MINERvA] sample within The Integrative Psychiatric Research (iPSYCH) cohort, Norwegian Mother and Child cohorts [MoBa1] and [MoBa2], The NorthPop Birth Cohort Study [NorthPop] and the Upstate KIDS study [UpstateKIDS]). We also conducted analyses examining paternal asthma versus paternal asthma never in 11 cohorts with this information (ALSPAC, CHS, EDEN, GENR, INMA, IOW, LiNA, MoBa1, MoBa2, NorthPop and UpstateKIDS). A full list of all cohorts including case and control numbers across models is given in Table S1, and cohort-specific study descriptions and inclusion criteria are given in supplementary material.

Asthma definitions

Asthma during pregnancy was defined by maternally self-reported asthma and/or use of asthma medication and/or a doctor’s diagnosis during gestation. For all but one cohort, the asthma diagnosis was identified using questionnaires, whereas for the MINERvA cohort an asthma diagnosis was based on registries. Asthma ever was defined as self-reported asthma ever and/or use of asthma medication ever. By our definition, individuals classified as having asthma in pregnancy are also included in analyses of ever asthma. More detailed phenotype definitions for each cohort are given in Supplementary material. Asthma during pregnancy was defined in seven cohorts (EDEN, IOWF2, MINERvA, MoBA1, MoBA2, NorthPop and UpstateKIDS) with available data for 3899 individuals.

Methylation data measurement, quality control and annotation

Methylation was assessed using either the Illumina 450 K BeadChip platform or the Illumina EPIC 850 K chip. For all cohorts, the minimum recommended DNA amount of 500 ng was provided to the laboratories running the 450 K or EPIC arrays. Cohorts individually performed quality control, normalization and analyses of untransformed β values. Cross-reactive probes, probes located on X and Y chromosomes as well as probes that overlapped with known SNPs were excluded after meta-analysis [15]. Methylation beta values were trimmed using the 3*IQR trimming method as has been done previously [16], where beta values three times the interquartile range below the 25th percentile or above the 75th percentile for each CpG were removed [17].

Annotation of DNA methylation sites

We used the gene annotation provided in the Illumina annotation files for both DNA methylation chips. All annotations use the human GRCh37/hg19 assembly.

Cohort-specific statistical analyses

Each cohort ran the association between asthma types and DNA methylation using robust linear regression. Covariates included infant sex, gestational age as a continuous measure, mode of delivery with two categories: vaginal delivery and cesarian section delivery, maternal age as a continuous measure and socioeconomic status (cohort-specific definition, but in general maternal education and income). In this analysis, we adjusted for maternal smoking during pregnancy in three categories: none, quit early in pregnancy and those who smoke across pregnancy. Prior work in PACE cohorts has shown that the greatest impact of maternal smoking is seen for smoking that is sustained across the pregnancy not in the approximately half of smokers who quit early in pregnancy [18, 19]. The MoBa study found no associations for smoking by the mother that ended before pregnancy [18]. Cohorts were adjusted for batch effects by using ComBat [20] or by including a batch covariate in their models. The MINERvA cohort adjusted for DNA methylation smoking score at birth as a surrogate for maternal smoking [21]. If a selection factor was employed, cohorts additionally adjusted for this, for instance if a cohort contained cases and controls selected based on a condition or characteristic (see cohort-specific description given in Supplementary Material). Maternal BMI was not available for all included cohorts and was thus not included as a covariate; to accommodate this, we performed a lookup in a large PACE meta-analysis of maternal BMI and found no overlap with our findings [22].

Cell-type composition was adjusted by including all 7 estimated proportions of cells using the cord blood reference panel [23] calculated by the Houseman method [24] using the FlowSorted.Blood package available for minfi [25].

Meta-analyses

We meta-analyzed study-specific results with inverse variance weighting in METAL [26]. The meta-analysis was redone by an independent group using the same method and the results were compared to minimize the likelihood of human error. For the sex-stratified analyses, we restricted the meta-analyses to studies where there were at least 15 newborns of each sex exposed to the parental asthma condition under study (see Table S1) and used the resulting number of studies in the meta-analysis where sex and parental asthma were investigated using an interaction term. For maternal asthma ever, we used maternal asthma during pregnancy as a surrogate for the MINERvA and EDEN cohorts. We performed analyses restricted to either the 450 K chip including CpGs, after QC filtering, that overlapped between 450 K and the EPIC chip (424,403 CpGs) or were exclusive to the EPIC chip (321,034 CpGs). In total 10,922 CpGs were removed from the 450 K analysis because they were not captured on the EPIC chip. We restricted the analysis to CpGs with available data from at least three studies for probes on the 450 K chip and two studies for probes available on the EPIC chip data (as we only included two studies with such data) and accounted for multiple testing by controlling the false discovery rate (FDR) using a threshold of 5% for each chip-specific analysis along with a more strict FDR threshold of 0.025 as some would argue for using this threshold when separating CpGs into chip-specific analyses. We calculated if the observed effect sizes were homogeneous (I2 value) across cohorts using METAL [23]. We show forest plots for significantly differentially methylated CpGs (DMCs) including effect estimates and 95% confidence intervals for each cohort.

Analyses of differentially methylated regions

We identified differentially methylated regions (DMRs) using comb-p [27] as this method tends to be more conservative than DMRcate [28]. Comb-p corrects multiple comparisons through a one-step Šidák correction [27]. We identified significant DMRs using an adjusted FDR p value below 0.05, required at least three probes with a maximum distance of 500 bp. DMRs were annotated to the nearest gene, regulatory regions and proximity to CpG islands if present on the Illumina annotation file for the hg19 reference genome.

Identification of drug targets and mQTL associations using ChEMBL and GoDMC

We looked in the ChEMBL database (version 31, https://www.ebi.ac.uk/chembl/) to identify genes implicated in our analyses of DMCs or DMRs that previously have been targets of approved drugs or drugs in development. We did lookup of DMCs in the GoDMC database (http://mqtldb.godmc.org.uk/search) to look for potential methylation quantitative loci (mQTL) associated with asthma.

Correlation of DNA methylation and gene expression

We examined if DNA methylation at significant individual or DMR CpGs was related to gene expression using lookup in 39,749 significant expression quantitative trait methylation (eQTM) pairs (FDR < 0.05) of blood DNA methylation probes from the 450 K chip array and blood gene expression data from 832 children available in the HELIX consortium [29]. eQTMs were identified using linear regression of methylation levels in relation to expression at nearby genes (using a 1mb window centered on the TSS) [29]. In addition, we also included 8,646 significant eQTM pairs (and FDR < 0.05) of DNA methylation probes from the EPIC chip array and gene expression data from adult blood and lung tissue from 424 individuals in GTEx (https://www.gtexportal.org/home/). We did this by lookup in the summary statistic data from the cell-type-adjusted HELIX data and GTEx data separately and considered significance based on FDR p values below 0.05.

Data availability

Genome-wide meta-analysis results will be available in the following link upon publication: https://doiorg.publicaciones.saludcastillayleon.es/10.5281/zenodo.13219057.

Results

Demographic description

There were 3899 individuals across 7 cohorts available for maternal asthma during pregnancy and 7433 individuals across 13 cohorts for maternal asthma ever. Supplementary material contains cohort-specific asthma definitions and distributions of parental asthma phenotypes. The prevalence of self-reported maternal asthma ranged from 2 to 35%, sustained smoking ranged between 0.4 and 30%. Most cohorts used the 450 K chip. Two cohorts (NorthPop and UpstateKIDS) had DNA methylation measured using the EPIC chip. Participants were primarily of European descent (Table 1). We did not see any overlap between our findings and a recently published EWAS on maternal BMI [22].

Table 1 Characteristics of participating study cohorts

An overview of included analyses, primary quality control filters and analyses are shown in Fig. 1.

Fig. 1
figure 1

Flowchart describing the number of samples, the primary quality control filters as well as the analyses included

Maternal asthma during pregnancy and newborn DNA methylation

The meta-analysis of newborn DNA methylation in relation to maternal asthma during pregnancy included 278 exposed and 3621 non-exposed participants from 7 cohorts: IOW, EDEN, MoBA1, MoBA2, NorthPop, UpstateKIDS and MINERvA. We identified 1 significant DMC for the 450 K chip (λ = 1.09), cg26963854 within the south shelf of a CpG island, which was not annotated to a specific gene (FDR 5%). We also identified 26 DMCs exclusive to the EPIC chip (FDR < 5%) using 76 exposed and 675 non-exposed participants from NorthPop and UpstateKIDS (λ = 1.18) (Fig. 2, Table 2). Using a stricter P value threshold of (FDR < 2.5%), we observe 6 DMCs pertaining only to the EPIC chip (Table S2). We identified 5 DMRs (encompassing 28 CpGs) in relation to maternal asthma during pregnancy for the 450 K chip and none for probes exclusive to the EPIC chip Table S3). Forest plots with cohort-specific beta values and 95% confidence intervals for the identified CpGs are shown in Fig. S1. For the CpG site (cg26963854) identified using the 450 K chip, the effect estimate was lower in the MINERvA cohort compared to that of the other included cohorts with the 450 K data, but without evidence of heterogeneity (Pheterogeneity = 0.14). For the EPIC chip, we saw no evidence of heterogeneity (Fig. S1).

Fig. 2
figure 2

Manhattan plot of epigenome-wide association of maternal asthma during pregnancy and offspring methylation. Models include the following covariates: child sex, maternal smoking during pregnancy, gestational age, mode of delivery, maternal age at childbirth, maternal socioeconomic status, estimated cell type, batch covariates and ancestry. Panel A) shows analysis restricted to 435,329 probes on the 450 K chip for 7 cohorts using 202 exposed and 2946 non-exposed individuals (lambda = 1.09), and panel B) shows 321,034 probes exclusive to the EPIC chip for 2 cohorts using 76 exposed and 675 non-exposed individuals (lambda = 1.18). The red line indicates 5% FDR significance

Table 2 Twenty-seven differentially methylated CpGs (FDR < 0.05) from the meta-analysis of maternal asthma during pregnancy in relation to newborn methylation

Maternal asthma ever and newborn DNA methylation

The meta-analysis of newborn methylation and maternal asthma ever included 726 exposed and 6707 non-exposed participants from 13 cohorts: ALSPAC, CHS, DCHS, EDEN, GENR, INMA, IOW, LiNA, MINERvA, MoBA1, MoBA2, NorthPop and UpstateKIDS, and identified no associated probes (FDR < 0.05) for either the 450 K or EPIC chip (Fig. S2. We identified 3 DMRs (comprising 28 CpGs) for the 450 K chip but none for the EPIC chip (Table S4). We did not observe any significant probes (FDR < 0.05) if we restrict the maternal asthma ever analysis to the 7 cohorts in the maternal asthma during pregnancy analysis.

Paternal asthma and newborn DNA methylation

The meta-analysis of newborn methylation and paternal asthma ever included 440 exposed and 4479 non-exposed participants from 9 cohorts: ALSPAC, GENR, INMA, IOW, LiNA, MoBA1, MoBA2, NorthPop and UpstateKids, and identified no probes associated with paternal asthma (FDR < 0.05) for the 450 K chip and 2 associated probes for the EPIC chip, namely cg08311378 in the gene body of RPS6KA2 and cg07462855 in the gene body of FAM160B1 (Fig. S3, Table S5).These DMCs were not among the 26 DMCs identified using the EPIC chip for the analysis of maternal asthma during pregnancy nor were they at least nominally significant in the maternal asthma (active or ever) analysis. Forest plots for the identified DMCs are in Fig. S4.

We identified 11 DMRs encompassing 99 CpGs in relation to paternal asthma ever diagnosis for the 450 K and 2 for the EPIC chip, encompassing 9 CpGs (Tables S6). Some overlap was detected with DMRs associated with maternal pregnancy asthma status (in gene PPT2; PRRT1) and maternal ever asthma (HOXA genes).

Correlation analysis of newborn methylation across maternal and paternal asthma

We correlated methylation effect sizes across 450 K CpGs from the meta-analyses results from maternal asthma during pregnancy (MAP), maternal asthma ever (MAE), maternal asthma ever without using MAP data (MAE exclusive) as well as paternal asthma ever (PAE) (Fig. 3). We observed a moderately strong positive correlation between MAP and MAE effect sizes (rho = 0.55, P value < 0.001), as well as positive correlation between MAP and MAE exclusive effect sizes (rho = 0.47, P value < 0.001). We did not observe a correlation between MAP and PAE (rho = 0.002, P value = 0.15), but observed a weak positive correlation between MAE and PAE (rho = 0.05, P value < 0.001) and p between MAE exclusive and PAE (rho = 0.06, P value < 0.001).

Fig. 3
figure 3

Spearman correlations of effects from meta-analyses of all included models for A) CpGs pertaining to the overlap between the 450 K and EPIC chip and B) the CpGs unique to the EPIC chip. MAP, maternal asthma during pregnancy; MAE, maternal asthma ever while; PAE, paternal asthma ever. Sex-stratified models are defined by suffixes. ‘Int’ represent interaction models testing for differences between the sexes

For the EPIC chip, MAP and MAE effect sizes we again found moderate positive correlated (rho = 0.47, P value < 0.001). MAP and PAE were weakly negatively correlated (rho = − 0.04, P value < 0.001) as were MAE and PAE (rho = − 0.03, P value < 0.001) (Fig. 3).

The single DMC found on the 450 K chip, cg26963854 located on chromosome 14, that passed the FDR threshold of 5% in maternal asthma during pregnancy had a similar direction of effect across the maternal asthma ever and paternal asthma ever models but was only nominally significantly associated (P < 0.05) in the maternal asthma ever meta-analysis.

Of the 26 DMCs identified in maternal asthma during pregnancy meta-analysis on the EPIC chip, we focused on the 21 that were available across the maternal and paternal asthma ever models. We observed similar direction of effect in the maternal asthma during pregnancy and the maternal asthma ever meta-analysis for 19 of these 21 CpGs, and among these 19, we observed 13 nominally significant p values (Table S7). In contrast, none of the CpGs identified in the maternal asthma during pregnancy analysis had significant p values in the paternal asthma ever model nor consistent directions of effect (Table S7). Of the two DMCs, cg08311378 and cg07462855, identified in paternal asthma ever meta-analysis on the EPIC chip, we did not observe a similar direction of effect when compared to effect sizes in MAP (Table S8). We also did not observe nominally significant p values in the MAP meta-analysis results. Compared to the MAE analysis effect estimates, we saw that cg07462855 had a similar direction of effect but neither of the DMCs were nominally significant (Table S8).

CpGs associated in the literature with childhood asthma or pulmonary function

We uploaded the top CpGs on the 450 K chip identified for MAP (Table S9), MAE (Table S10) and PAE (Table S11) to the EWAS toolkit platform [30] to investigate enrichment in previous DNAm results. For the trait enrichment analyses using MAP-associated CpGs, we observed enrichment in the following traits (which were also among the top10 most associated traits): asthma, smoking, atopy and maternal smoking (all had enrichment P values < 1.98 × 10–13) (Table S12).

Similarly, we observed that asthma and smoking were the top two traits among all traits when using MAE-associated CpGs (Table S13). Using CpGs associated with paternal asthma ever, we also observed significant association with asthma and smoking (Table S14).

In addition, we investigated enrichment for DMCs identified in the literature specifically for childhood asthma and lung function [8, 11, 31, 32]. We included DMCs identified in blood in an investigation of neonates developing asthma and among children with a clinical diagnosis of asthma in Reese et al. 2019 [8], DMCs in whole blood from childhood asthma in Xu et al. 2018 [31], DMCs identified in nasal epithelial cells from Qi et al. 2020 [32] and unique DMCs identified across DMRs for FEV1, FEV1/FVC and FEF75 in cord blood [11]. In total, 766 previously identified DMCs were used as the enrichment target, and we considered significant enrichment using a P value cutoff of 5% from Fisher’s exact test and used as input the CpGs in our analyses with P values below 0.005 (Table S9, S10 and S11). We did not observe any enrichment across maternal asthma during pregnancy (P value = 0.57), maternal asthma ever (P value = 0.84) or paternal asthma ever (P value = 0.85).

Sex-specific analyses

Sex-stratified analyses were performed for all 3 main phenotypes. For the sex-stratified meta-analysis of newborn methylation and maternal asthma during pregnancy, we included 102/943 exposed/non-exposed for boys and 83/781 exposed/non-exposed for girls from 3 cohorts: MoBA1, MoBA2 and NorthPop. We observed 6 DMCs for boys and 325 DMCs for girls related to maternal asthma during pregnancy at FDR < 0.05. None of the 6 DMCs identified in boys gave evidence of significant interaction with sex (FDR < 0.05) (Table S15. Investigating the 325 DMCs identified for girls, 154 had nominal significance for the sex-specific interaction; however, none was statistically significant (FDR < 0.05). Of these 325, 19 were also nominally significant (14 had same direction of effect) in the boys and not among the 6 DMCs identified in boys alone (Table S16).

For the sex-stratified meta-analysis of newborn methylation and maternal asthma ever, we included 309/2347 exposed/non-exposed for boys and 278/2206 exposed/non-exposed for girls from 7 cohorts: ALSPAC, GENR, LiNA, MoBA1, MoBA2, NorthPop and UpstateKIDS; we identified 1 DMC in boys and 25 in girls; none gave evidence of interaction with sex (FDR < 0.05) but 15 of the DMCs identified for girls met nominal significance (Table S17). For paternal asthma exposure, we included 183/1920 exposed/non-exposed for boys and 182/1834 exposed/non-exposed for girls from 7 cohorts: ALSPAC, GENR, LiNA, MoBA1, MoBA2, NorthPop and UpstateKIDS, and we observed 100 DMCs in boys and 95 DMCs in girls (Table S18, S19). Among the 100 DMCs identified in boys only, we saw nominally significant evidence of interaction for 5 sites (Table S18). Among the 95 DMCs found in girls, we observed ten DMCs, three with known annotation near genes HCCA2, C1orf198 and PNMT, with statistically significant interaction (FDR < 0.05) (Table S19). For meta-analysis of sex interaction for paternal asthma exposure, we identified 12 DMCs (Table S20), where we observed in general lower methylation and a stronger effect in girls compared to boys.

Differential DNA methylation and gene expression in blood and lung

To investigate whether differently methylated sites may be associated with gene expression, we analyzed eQTM pairs for 832 blood samples available from the HELIX consortium for child blood [29] and eQTM pairs for 424 lung samples from the GTEx consortium [33].

Among the 402 DMCs identified across all models using the 450 K chip, we observed 15 unique DMCs with significant associations with gene expression in blood (Table S21). Among these associations, we observed a decreased expression of FAM43A with increased methylation of cg02072170, and this gene is associated with eosinophil counts and thus associated with asthma etiology [34]. We also observed increased expression of LTBP1 with increased cg15772133 methylation, and this gene has been associated with FEV1/FVC in adults [35, 36].

Among the 441 unique DMCs that are encompassed by identified DMRs across all models using the 450 K chip, we found 188 significant unique methylation and gene expression pairs using the HELIX data (Table S22). We observed CpGs that annotated to genes that previously have been associated with asthma in adults including HLA genes [37].

Examining 56 DMCs identified from DMRs in models with maternal asthma during pregnancy and ever as the primary exposures, we found 35 unique DMCs were associated with gene expression of 4 genes. One gene of particular interest is HOXA5, which has previously been associated with organogenesis [38], lung function in adults [39, 40] and mental disorder phenotypes [41] (Table S22). The two other identified genes (KDM2B and KCTD11) were found to be involved with neurodevelopmental disorders [42] and cancer [43], respectively.

Among the 402 unique DMCs identified across all models using the 450 K chip, we observed one association between increased methylation of cg20810675 and decreased expression of C4orf27 in lung tissue in GTEx (Table S23). Among the 411 DMCs that are encompassed in DMRs across all models, we identified substantially more associations with gene expression compared to the 402 unique single DMCs identified (Table S24). Specifically, these associations included many HLA gene variants and also the NOTCH4 gene which has been associated with schizophrenia [44], psoriasis [45] and asthma [46].

Druggable targets

We identified differential methylation in regions related to parental asthma involving the HOXA5 and HLA genes. HOXA5 is also the target of CHEMBL4224852, a lysine demethylase, which may implicate epigenetic regulation in asthma development. Several HLA genes were also identified as targets for drugs, including HLA-C (target of CHEMBL4680046) and HLA-DRB1 (target of CHEMBL2109447).

Discussion

We investigated parental asthma and newborn DNA methylation using data from 13 cohorts in the PACE consortium and found evidence that parental asthma is associated with differential DNA methylation in newborns. We identified 27 differentially methylated CpGs and 5 differently methylation regions associated with maternal asthma during pregnancy and none for the maternal asthma ever. These results suggest the relative importance of active maternal asthma on offspring methylation patterns. These CpGs were enriched for published associations with asthma and related phenotypes. While we found many more differentially methylated CpGs in girls, no sex interactions were significant at the epigenome-wide level and few reached nominal significance.

Interpretation

A stronger effect of maternal asthma during pregnancy compared to maternal asthma with respect to numbers of significant DMCs suggests that the timing of the disease exposure is important. In addition, exposure to maternal asthma during pregnancy may suffer from less misclassification that maternal asthma ever, because exposure was more recent and during pregnancy therefore more likely to be physician diagnosed. Also, we cannot rule out that the timing of asthma in the asthma ever exposure may be a relevant factor and thus “dilute” this exposure compared to maternal asthma during pregnancy. Secondly, we cannot rule out that the observed differential methylation is caused by the mother taking asthma medication during pregnancy. However, we also observed a high correlation of effects between the two maternal exposures, suggesting that at least part of the effect can also be seen for an ever diagnosis.

For paternal asthma, findings were enriched for those identified for maternal asthma during pregnancy, but enrichment appeared to be less strong than for the maternal asthma ever exposure- and chip-specific. The two CpGs identified for paternal asthma were also identified as significant in the maternal asthma analyses, indicating that these may represent more general processes not unique to maternal asthma exposure. We note that the power in terms of included cohorts and individuals for the paternal asthma ever was like the analysis of maternal asthma ever, suggesting that the difference in findings is not due to sample size. Taken together, differential DNA methylation was more pronounced if the exposure is active maternal asthma during pregnancy as compared to paternal asthma, substantiating the specificity of the prenatal window of exposure.

Our findings for parental asthma in relation to methylation in newborns were enriched for overlap with those for asthma-relevant traits in prior EWAS. Specifically, we observed enrichment in asthma, smoking and atopy as well as maternal smoking for CpGs identified using maternal asthma during pregnancy, while maternal asthma ever showed enrichment for asthma and smoking. Taken in concert, we also believe that the enrichment in previous asthma-related CpGs, but not CpGs identified specifically for childhood asthma, underlines that the risk induced by maternal asthma exposure pertains more to a T2 type asthma characterized by more inflammation and a later debut [47]. These findings need to be confirmed in future research on maternal health and its impact on child development and disease predisposition.

We also identified differential methylation in regions related to parental asthma involving the HOXA5 and HLA genes. Both genes have been targets for the drug. Altered expression of HLA genes has previously been linked to asthma and allergy [37] while perturbed expression of HOXA5 has been associated with impaired lung function in children [11] and in adults [39, 40] highlighting the importance of future studies investigating their role in asthma pathogenesis. Notably, HLA genes and HOXA5 have also been associated with mental disorders and the latter may be important in development [41].

There is considerable interest in differential impacts of parental factors on child DNA methylation depending on sex where the effect of parental allergy on childhood allergic diseases has been demonstrated to depend on the sex of the child [48]. More DMCs were identified in girls compared to boys, suggesting that girls might be more susceptible to methylation changes from exposure to maternal asthma. However, we did not observe any statistically significant sex interactions in our analyses (P < 0.05). Thus, the findings presented here should be interpreted with caution. Other studies have suggested a differential effect of maternal asthma on female children in terms of microbiome composition [49] on fetal growth [50] and resulting lower birth weight [51]. Future studies would be needed to determine lack of sex-specific effects.

Strengths and limitations

The major strength of this study is the large sample size, inclusion of multiple cohorts from different populations and the investigation of epigenetic marks in samples collected at the same developmental stage, which enhances the generalizability of our findings. Furthermore, the study adjusted for potential confounders and included cord blood-specific cell-type adjustment [24], which is crucial given that DNA methylation patterns can vary significantly between different cell types. An additional strength of our study is that it used both data from the Illumina MethylationEPIC BeadChip and the Illumina HumanMethylation450K BeadChip, thus increasing the number of CpGs investigated in relation to parental asthma.

The study also has limitations. First, most of the single DMCs were identified in EPIC unique meta-analyses that included only two cohorts. And secondly, none of the 27 identified CpGs has previously been reported to be association with childhood asthma, allergy or lung function in previous large-scale EWAS meta-analyses [8, 11, 31, 32], but we did observe enrichment datasets previously associated with asthma, smoking and atopy. Thirdly, because methylation was measured in various blood sources, the findings may not directly translate to lung tissue, which is a primary organ affected in asthma [52]. However, a previous study showed a high level of agreement between DNAm in blood and bronchial epithelial cells in functional relevant regions [53]. Also, blood for DNA methylation analyses came primarily from cord blood but in one study, newborn blood spots and in another newborn peripheral blood was used. Reference panels are available for cell deconvolution only for cord blood. This could be a source of between study heterogeneity, but we see little evidence for this. Another limitation is that asthma was based on self-reported or, for parental asthma reports in most cohorts and no report on the timing of disease for the parental ever exposure. This may have led to misclassification or noise in the exposure. Finally, to maximize the number of studies and overall sample size for analyses, we focused on probes overlapping between the 450 K and EPIC arrays. However, this approach inherently excludes the few probes present on the 450 K array that are absent on the EPIC array, potentially omitting relevant findings from these non-overlapping probes.

Our findings have important research implications. Identifying differentially methylated CpGs and regions linked to maternal and paternal asthma may improve understanding of asthma mechanisms and reveal potential therapeutic targets if DNAm changes are causal. While no SNPs were associated with the identified DMCs, asthma-related SNPs affecting DNAm cannot be excluded. Stronger findings for maternal asthma during pregnancy highlight the importance of exposure timing in DNAm changes.

Future studies should study the timing of the asthma diagnosis for both maternal and paternal diagnosis in predicting offspring DNAm changes including if both mother and father had asthma at birth and/or conception. Future research should also explore the functional consequences of any DNAm changes and their use as predictive markers for asthma in children. Longitudinal research on the persistence of these patterns and their impact on disease outcomes could provide further insight into DNAm’s role in asthma.

Conclusion

Our findings suggest that parental asthma is associated with DNA methylation patterns of newborns, with several DMCs and DMRs identified, specifically in relation to maternal asthma during pregnancy.

Availability of data and materials

Genome-wide meta-analysis results is available in the following link: https://doiorg.publicaciones.saludcastillayleon.es/10.5281/zenodo.13219057

References

  1. Eder W, Ege MJ, von Mutius E. The asthma epidemic. N Engl J Med. 2006;355:2226–35.

    Article  PubMed  CAS  Google Scholar 

  2. Ober C, Yao T-C. The genetics of asthma and allergic disease: a 21st century perspective. Immunol Rev. 2011;242:10–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Han Y, et al. Genome-wide analysis highlights contribution of immune system pathways to the genetic architecture of asthma. Nat Commun. 2020;11:1776.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Conroy ER, Peterson R, Phipatanakul W, Sheehan WJ. Increasing awareness regarding the relationship between environmental exposures and allergic disease. J Allergy Clin Immunol. 2024. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jaci.2024.08.008.

    Article  PubMed  Google Scholar 

  5. Liu X, et al. Bidirectional associations between asthma and types of mental disorders. J Allergy Clin Immunol Pract. 2023;11:799-808.e14.

    Article  PubMed  Google Scholar 

  6. Patel R, et al. Sex-specific associations of asthma acquisition with changes in DNA methylation during adolescence. Clin Exp Allergy. 2021;51:318–28.

    Article  PubMed  CAS  Google Scholar 

  7. Solomon O, et al. Meta-analysis of epigenome-wide association studies in newborns and children show widespread sex differences in blood DNA methylation. Mutat Res Rev Mutat Res. 2022;789: 108415.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Reese SE, et al. Epigenome-wide meta-analysis of DNA methylation and childhood asthma. J Allergy Clin Immunol. 2019;143:2062–74.

    Article  PubMed  CAS  Google Scholar 

  9. Forno E, et al. DNA methylation in nasal epithelium, atopy, and atopic asthma in children: a genome-wide study. Lancet Respir Med. 2019;7:336–46.

    Article  PubMed  CAS  Google Scholar 

  10. Edris A, den Dekker HT, Melén E, Lahousse L. Epigenome-wide association studies in asthma: a systematic review. Clin Exp Allergy. 2019;49:953–68.

    Article  PubMed  CAS  Google Scholar 

  11. den Dekker HT, et al. Newborn DNA-methylation, childhood lung function, and the risks of asthma and COPD across the life course. Eur Respir J. 2019;53:1801795.

    Article  Google Scholar 

  12. Lim RH, Kobzik L, Dahl M. Risk for asthma in offspring of asthmatic mothers versus fathers: a meta-analysis. PLoS ONE. 2010;5: e10134.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Felix JF, et al. Cohort profile: pregnancy and childhood epigenetics (PACE) consortium. Int J Epidemiol. 2018;47:22–23u.

    Article  PubMed  Google Scholar 

  14. Bakulski KM, Blostein F, London SJ. Linking prenatal environmental exposures to lifetime health with epigenome-wide association studies: State-of-the-science review and future recommendations. Environ Health Perspect. 2023;131: 126001.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Zhou W, Laird PW, Shen H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 2016. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkw967.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Everson TM, et al. Placental DNA methylation signatures of maternal smoking during pregnancy and potential impacts on fetal growth. Nat Commun. 2021;12:5095.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. David FN, Tukey JW. Exploratory data analysis. Biometrics. 1977;33:768.

    Google Scholar 

  18. Joubert BR, et al. Maternal smoking and DNA methylation in newborns: in utero effect or epigenetic inheritance? Cancer Epidemiol Biomark Prev. 2014;23:1007–17.

    Article  CAS  Google Scholar 

  19. Joubert BR, et al. DNA methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98:680–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2006;8:118–27.

    Article  PubMed  Google Scholar 

  21. Reese SE, et al. DNA methylation score as a biomarker in newborns for sustained maternal smoking during pregnancy. Environ Health Perspect. 2017;125:760–6.

    Article  PubMed  CAS  Google Scholar 

  22. Sharp GC, et al. Maternal BMI at the start of pregnancy and offspring epigenome-wide DNA methylation: findings from the pregnancy and childhood epigenetics (PACE) consortium. Hum Mol Genet. 2017;26:4067–85.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Gervin K, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenetics. 2019;11:125.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Houseman EA, Molitor J, Marsit CJ. Reference-free cell mixture adjustments in analysis of DNA methylation data. Bioinformatics. 2014;30:1431–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Salas LA, et al. Enhanced cell deconvolution of peripheral blood using DNA methylation for high-resolution immune profiling. Nat Commun. 2022;13:761.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Pedersen BS, Schwartz DA, Yang IV, Kechris KJ. Comb-p: software for combining, analyzing, grouping and correcting spatially correlated P-values. Bioinformatics. 2012;28:2986–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Mallik S, et al. An evaluation of supervised methods for identifying differentially methylated regions in Illumina methylation arrays. Brief Bioinform. 2019;20:2224–35.

    Article  PubMed  CAS  Google Scholar 

  29. Ruiz-Arenas C, et al. Identification of autosomal cis expression quantitative trait methylation (cis eQTMs) in children’s blood. Elife. 2022. https://doiorg.publicaciones.saludcastillayleon.es/10.7554/eLife.65310.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Xiong Z, et al. EWAS open platform: integrated data, knowledge and toolkit for epigenome-wide association study. Nucleic Acids Res. 2022;50:D1004–9.

    Article  PubMed  CAS  Google Scholar 

  31. Xu C-J, et al. DNA methylation in childhood asthma: an epigenome-wide meta-analysis. Lancet Respir Med. 2018;6:379–88.

    Article  PubMed  CAS  Google Scholar 

  32. Qi C, et al. Nasal DNA methylation profiling of asthma and rhinitis. J Allergy Clin Immunol. 2020;145:1655–63.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Oliva M, et al. DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits. Nat Genet. 2023;55:112–22.

    Article  PubMed  CAS  Google Scholar 

  34. Sakaue S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53:1415–24.

    Article  PubMed  CAS  Google Scholar 

  35. Barton AR, Sherman MA, Mukamel RE, Loh P-R. Whole-exome imputation within UK Biobank powers rare coding variant association and fine-mapping analyses. Nat Genet. 2021;53:1260–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Kichaev G, et al. Leveraging polygenic functional enrichment to improve GWAS power. Am J Hum Genet. 2019;104:65–75.

    Article  PubMed  CAS  Google Scholar 

  37. Clay SM, et al. Fine-mapping studies distinguish genetic risks for childhood- and adult-onset asthma in the HLA region. Genome Med. 2022;14:55.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Jeannotte L, Gotti F, Landry-Truchon K. Hoxa5: a key player in development and disease. J Dev Biol. 2016. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/jdb4020013.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Shrine N, et al. New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries. Nat Genet. 2019;51:481–93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Shrine N, et al. Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk. Nat Genet. 2023;55:410–22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Lizen B, et al. Conditional loss of Hoxa5 function early after birth impacts on expression of genes with synaptic function. Front Mol Neurosci. 2017. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fnmol.2017.00369.

    Article  PubMed  PubMed Central  Google Scholar 

  42. van Jaarsveld RH, et al. Delineation of a KDM2B-related neurodevelopmental disorder and its associated DNA methylation signature. Genet Med. 2023;25:49–62.

    Article  PubMed  Google Scholar 

  43. Yang M, et al. KCTD11 inhibits progression of lung cancer by binding to β-catenin to regulate the activity of the Wnt and Hippo pathways. J Cell Mol Med. 2021;25:9411–26.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Aberg KA, et al. A comprehensive family-based replication study of schizophrenia genes. JAMA Psychiat. 2013;70:573–81.

    Article  CAS  Google Scholar 

  45. Kurki MI, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613:508–18.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Zhu Z, et al. Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank. J Allergy Clin Immunol. 2020;145:537–49.

    Article  PubMed  CAS  Google Scholar 

  47. Kuruvilla ME, Lee FE-H, Lee GB. Understanding asthma phenotypes, endotypes, and mechanisms of disease. Clin Rev Allergy Immunol. 2019;56:219–33.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Arshad SH, et al. The effect of parental allergy on childhood allergic diseases depends on the sex of the child. J Allergy Clin Immunol. 2012;130:427-34.e6.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Koleva PT, et al. Sex-specific impact of asthma during pregnancy on infant gut microbiota. Eur Respir J. 2017;50:1700280.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Meakin AS, Saif Z, Seedat N, Clifton VL. The impact of maternal asthma during pregnancy on fetal growth and development: a review. Expert Rev Respir Med. 2020;14:1207–16.

    Article  PubMed  CAS  Google Scholar 

  51. Stevens DR, et al. Maternal asthma in relation to infant size and body composition. J Allergy Clin Immunol Glob. 2023;2: 100122.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Banno A, Reddy AT, Lakshmi SP, Reddy RC. Bidirectional interaction of airway epithelial remodeling and inflammation in asthma. Clin Sci (Lond). 2020;134:1063–79.

    Article  PubMed  CAS  Google Scholar 

  53. Lee Y-S, et al. Epigenome-scale comparison of DNA methylation between blood leukocytes and bronchial epithelial cells. Epigenomics. 2021;13:485–98.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We are extremely grateful to all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. ARIES DNA methylation data used in this analysis was generated in the Bristol Bioresource Laboratory Illumina Facility, University of Bristol. We acknowledge all participating families in the NorthPop study; the NorthPop project team for recruitment, follow-up and blood samplings of study participants; the NorthPop coordinator Richard Lundberg at the Department of Clinical Sciences, Umeå University and the personnel at Biobanken Norr,, Västerbotten county council. Methylation profiling was performed by the SNP&SEQ Technology Platform in Uppsala (www.genotyping.se). The facility is part of the National Genomics Infrastructure (NGI) Sweden and Science for Life Laboratory. The SNP&SEQ Platform is also supported by the Swedish Research Council and the Knut and Alice Wallenberg Foundation. We thank the Upstate KIDS participants and staff for their important contributions. This work utilized the computational resources of the NIH HPC Biowulf cluster (https://hpc.nih.gov/).

Funding

Open access funding provided by the National Institutes of Health. C-ETP and KB are part of COPSAC and funded by The Lundbeck Foundation (Grant no R16-A1694); The Ministry of Health (Grant no 903516); Danish Council for Strategic Research (Grant no 0603-00280B) and The Capital Region Research Foundation have provided core support to the COPSAC research center. SJL is funded by the Intramural Research Program of the NIH, National Institute of Environmental Health Sciences (ES 49019). LD received funding for projects from the European Union’s Horizon 2020 research and innovation program (LIFECYCLE, grant agreement No 733206, 2016 EUCAN-Connect grant agreement No 824989 ATHLETE, grant agreement No 874583 ENDOMIX No 101136566). AS was funded by the Lundbeck foundation and the Independent Research Fund Denmark. DS was funded by the Lundbeck foundation, The Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institute of Environmental Health Sciences, National Institute of Neurological Disorders and Stroke and the Beatrice and Samuel A. Seaver Foundation. NHS was funded by the Lundbeck foundation. HRE, RG work in the Medical Research Council Integrative Epidemiology Unit at the University of Bristol, UK, which is supported by the Medical Research Council and the University of Bristol (MC_UU_00011/5). The funding bodies had no role in the design of the study. Cohort-specific funding is mentioned in supplementary material. The researchers are independent from the funders. The study sponsors had no role in the study design, data analysis, interpretation of data or writing of this study.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

C-ETP wrote the initial draft of the paper, did the formal meta-analysis, visualized the results and administrated the project including coordinating efforts for the project pre- and post-publication. The shadow analysis was done by JJ. The concept and ideas were developed by AS, HA, NHS, HTH, JM, EH, DS, SJL, KB. The methodology and creation of models were initiated by AS, NHS, C-ETP and KB. The formal analysis was done by C-ETP, TTH, AS, YZ, HTD, LD, JFF, JS, LK, FIR, SR, HTH, SHM, JH, XZ and DS. AS, RG, HRE and GH performed experiments and or collected the data. The resources used in this study were provided by RG, HRE, HJZ, DJS, HTD, LD, JFF, MB, MC, MV, ACZ, BH, CVB, SH, MD, CW, EY. Data curation was done by HTD, LD, JFF, HTH, JH, BH, EY, DS. Writing and editing of the manuscript including critical revision, commentary in the pre- or post-publication stages were done by C-ETP, TTH, AH, HJZ, DJS, HTD, LD, JFF, MB LK, FIR HA, JWH, GH, HTH, IAM, GP, NB, BH, CVB, SH, JH, MD, CW, EY, WN, SHE, MCM, DS, SJL and KB. Visualization and creation of the data presentation was done by C-ETP. Supervision including leadership responsibility for the research, planning and execution were done by AH, EY, SJL and KB. Project administration and coordination responsibility was done C-ETP and SJL. Funding acquisition for specific cohorts was done by JWH and HA.

Corresponding authors

Correspondence to Stephanie J. London or Klaus Bønnelykke.

Ethics declarations

Competing interests

All authors declare no potential, perceived or real conflict of interest regarding the content of this manuscript. The funding agencies did not have any role in design and conduct of the study; collection, management and interpretation of the data; or preparation, review or approval of the manuscript. No pharmaceutical company was involved in the study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pedersen, CE.T., Hoang, T.T., Jin, J. et al. Maternal asthma and newborn DNA methylation. Clin Epigenet 17, 79 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13148-025-01858-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13148-025-01858-4