Asthma is characterized by airflow limitation and airway reactivity (AR). Interleukin-13 (IL-13) is involved in the pathogenesis of asthma. Two functional SNPs, rs20541 and rs1800925, of the IL-13 gene (IL13) have been frequently associated with asthma-related lung functions. However, genetic variation alone does not fully explain asthma risk. DNA-methylation (DNA-M) is an epigenetic mechanism that regulates gene expression and can be influenced by both environment and genetic variants. To explore the interplay of prenatal maternal smoking, genetic variants and DNA-M, we used a two-stage model: (1) identifying cytosine phosphate guanine (CpG) sites where DNA-M is influenced by the interaction between genetic variants and maternal smoking during pregnancy (conditional methQTL (methylation quantitative trait loci)); and (2) determining the effect of the interaction between DNA-M of CpG (from stage 1) and SNPs (modifying genetic variants; modGV) on airflow limitation and AR in 245 female participants of the Isle of Wight birth cohort. DNA-M was assessed using the Illumina Infinium HumanMethylation450 BeadChip.
Six CpG sites were analyzed in stage 1. DNA-M at cg13566430 was influenced by interaction of maternal smoking during pregnancy and rs20541. In stage 2, genotype at rs1800925 interacted with DNA-M at cg13566430 significantly affecting airflow limitation (P = 0.042) and AR (P = 0.01).
Both genetic variants and environment affect DNA-M. This study supports the proposed two-stage model (methQTL and modGV) to study genetic variants, environment and DNA-M interactions in asthma-related lung function.
Asthma genetics and epigeneticsAirway reactivityDNA methylationIL13 geneLung functionsMaternal smoking during pregnancy
Asthma is a chronic airway disease characterized by airflow limitation and airway reactivity (AR) and exhibits wide heterogeneity in disease susceptibly and phenotypes. Multiple genes, gene-gene and gene-environment interactions have been shown to play a role in determining susceptibility to asthma and associated phenotypes of lung function and AR. Maternal smoking during pregnancy is a significant risk factor for developing asthma in offspring  and can affect offspring lung function . Interleukin-13 (IL-13) is a recognized effector in airway inflammation, reactivity and remodeling. IL13 is located on chromosome 5q31, and has been consistently associated with asthma [3, 4]. Several functional genetic variants occur in IL13 including rs1800925 (−1112C/T) and rs20541 (R130Q, +2044A/G). Recent meta-analyses have shown that both SNPs are associated with asthma risk [5, 6] and with forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC) in asthmatics [7, 8]. In severe asthma, monoclonal antibody to IL-13 use is associated with improvement in lung function in humans  and airway hyper-responsiveness in a murine model .
DNA-methylation (DNA-M) represents a site of molecular interaction between the environment and genome. There is growing evidence that DNA-M plays a role in complex diseases like asthma  and can be modified by environmental exposures such as tobacco smoke  as well as by disease-associated genotypes . We have previously shown that IL13 polymorphism modifies the impact of in utero tobacco smoke exposure on childhood asthma, suggesting a role for gene-environment interaction . Recently we have also shown that genetic variants in the IL-4 receptor interact with DNA-M to determine risk of asthma . DNA-M is a potential integrator of different signals affecting disease susceptibility, with both environment and genotype influencing methylation levels. Karmaus and colleagues  proposed a two-stage model to incorporate the role of genetic variants, environment and DNA-M interactions in asthma. In stage 1, an environmental exposure and genetic variant interact to influence DNA-M at a specific site in an adjacent locus. This stage identifies the conditional methQTL (methylation quantitative trait loci) and the change in DNA-M once established can differentially regulate gene activity. In stage 2, the phenotypic effects of sequence variants of the gene (modifiable genetic variants (modGVs)) can be modified by the pre-established methylation by the conditional methQTL. This two stage model for asthma-related lung function is depicted in Figure 1. Given our previous observation of the interaction between maternal smoke exposure and genotype in determining asthma, we hypothesize that this interaction would occur through a mechanism involving methQTL and/or modGVs.
Study population and assessments
The Isle of Wight birth cohort (n = 1,456) was established in 1989 and has been followed up through childhood. At 18 years of age, spirometry and bronchial challenge tests were done according to American Thoracic Society guidelines [17, 18]. We analyzed lung function measurements and DNA-M data of 245 females aged 18 years who were randomly selected from the cohort population for epigenetic studies. Comparison of the analyzed sample with whole cohort females in terms of risk factors have been described elsewhere . A continuous measure of AR (dose response slope (DRS)) was used. DRS is the gradient of FEV1 drop from baseline with each successive incremental dose of methacholine administered. A transformation of Log10 (DRS + 10) was required to satisfy the distributional assumption of normal data. A higher positive value of this transformation signifies greater AR. A description of the 18-year follow-up of the cohort is previously reported .
Genotyping and SNP selection
DNA was isolated from peripheral blood leucocytes collected at age 18. IL13 polymorphisms were genotyped as described by Sadeghnejad and colleagues . Two known functional SNPs, rs1800925 (n = 234) and rs20541 (n = 235), frequently associated with asthma and related lung functions were used in this analysis. These SNPs were not in linkage disequilibrium in this cohort as described elsewhere .
One microgram of DNA was bisulfite-treated using the EZ 96-DNA-methylation kit (Zymo Research, Irvine, CA, USA), following the manufacturer’s protocol. Genome-wide DNA- M was assessed using the Illumina Infinium-HumanMethylation450 BeadChip (Illumina, Inc., Hayward, CA, USA) as described previously .
The pre-processed DNA-M beta (β) values, presented as the proportion of intensity of methylated (M) over the sum of methylated (M) and unmethylated (U) sites (β = M/[c + M + U] with c being a constant to prevent dividing by zero), were used to estimate the effect of DNA-M . The R-package IMA in Bioconductor (IMA is implemented in the R language and is freely available from http://www.rforge.net/IMA) was used for the pre-processing . SNP genotype-dependent methylation was analyzed using the Kruskal-Wallis test. Interaction was tested using multiple linear regressions. Statistical analyses were performed using IBM SPSS Statistics, Version-19.0 (IBM SPSS Statistics for Windows, Version 19.0. Armonk, NY, USA: IBM Corp,). The statistical significance was set at 0.05.
Six cytosine phosphate guanine (CpG) sites were identified in the promoter region of IL13. The location and description of the CpG sites and of SNPs are shown in Table 1. All the six CpG sites spanning the promoter region of IL13 were analyzed in Stage 1. The effect of maternal smoking during pregnancy interacting with both SNPs was explored independently for DNA-M at each CpG site. DNA-M at cg13566430 was influenced by the interaction of rs20541 and maternal smoking during pregnancy (P = 0.043); this remained significant after correcting for personal smoking at 18 years (P = 0.041). DNA-M at cg13566430 also showed a genotype dependent methylation for rs1800925 (Kruskall Wallis test P < 0.001; Figure 2).
Location, position and description of rs1800925, rs20541 and CpG sites on the promoter region ofIL13gene
5′ Promoter upstream
Percentiles (5%, 95%)
CpG, cytosine phosphate guanine; TSS200, 200 base pairs from transcription start site; TSS1500, 1500 base pairs from transcription start site.
In the next step, the interaction between DNA-M at cg13566430 and rs1800925 genotype on lung function (FEV1/FVC and DRS) was explored. DNA-M at cg13566430 significantly interacted with rs1800925 genotype to determine FEV1/FVC and DRS (Table 2). In subjects with the TT genotype, FEV1/FVC increases with increasing methylation of cg13566430 (coefficient 3.274, P = 0.042), for the CT genotype group the increase was smaller (coefficient 0.799, P = 0.086) and CC was the reference genotype. The graph in Figure 3 provides a visual description of the effect of DNA-M at cg13566430 and rs1800925 on FEV1/FVC. DRS decreases quickly with increasing methylation of cg13566430 in TT genotype (coefficient −27.497, P = 0.010), while the drop in DRS in CT is insignificant (coefficient −0.809, P = 0.742).
Interaction of methylation at cg13566430 and rs1800925 genotypes on FEV1/FVC and DRS
This study tested a two-stage model for integrating the interactions of maternal smoking during pregnancy, genetic variants and DNA-M for an asthma candidate gene IL13. We show that interaction of a functional IL13 SNP, rs20541, and maternal smoking during pregnancy influenced DNA-M at cg13566430. We also show that DNA-M at cg13566430 interacts with genotype of another functional SNP, rs1800925, to affect airflow limitation and AR. Michel and colleagues  examined the effect of farm exposure on DNA-M of ten asthma candidate genes and found that DNA-M at one IL13 site (spanning rs1800925; similar to the site in our study) was more methylated in the exposed group compared to the non-exposed group. They did not see significant differential methylation of IL13 between asthmatic and non-asthmatic children but asthma-related traits were not tested and also the interaction with genetic variants was not examined. Our approach identifies the effect of environmental exposure and genetic variants on DNA-M and then also combines the interaction of other modifiable genetic variants with DNA-M on the outcomes. We show genotype-dependent DNA-M in the IL13 promoter region; methylation at cg13566430 revealed variation in levels dependent on the genotype of rs1800925, as seen in other genes . This varying distribution of DNA-M across genotypes supports the plausible role for DNA-M in the pathway between genotype and phenotype.
There are some limitations to this study. While DNA-M measurement using Illumina Infinium-HumanMethylation450 arrays has been shown to have reproducibility and high validity , technical replication of the DNA-M measurements has not been undertaken. We have not measured DNA-M in airway tissue; however, the major source of IL-13 production in the airways is inflammatory cells, and methylation measurements in peripheral blood are likely to better reflect asthma-related immune mechanisms. Cell composition in the peripheral blood can influence DNA-M; however, cell composition alone cannot explain the differential methylation observed . Focusing only on female participants is a limitation; however, the results should be encouraging for further studies to replicate the model. The interaction effects seen may not imply a direct relationship between genotype and altered methylation. Genotype may also alter other epigenetic processes such as chromatin remodeling and this may then lead to altered DNA-M.
The effect of environment, genotype and DNA-M was seen for two different, but equally important and objective, characteristics of asthma: airflow limitation and AR. Similar to Michel and colleagues  we did not see an effect on asthma as an outcome (results not shown); however, the observed effect on asthma-related objective measures avoids plausible bias of reported asthma. Methylation of the promoter region can regulate gene transcription  and we have shown that promoter region methylation is dependent on genotype, and the interplay of the environment, genotype and DNA-M influence the phenotype.
While requiring replication in an independent cohort, the results show the interplay of prenatal maternal smoking, genetic variants and DNA-M of IL13 influencing asthma-related lung function. This highlights the need to consider environment, genotype and DNA-M together when seeking to understand the pathogenesis of complex disease, as DNA-M plays a role of integrator of multiple disease pathway signals.
Cytosine phosphate guanine
Dose response slope
Forced expiratory volume in 1 second
Forced vital capacity
Methylation, quantitative trait loci
Modifiable genetic variant
Single nucleotide polymorphism.
We would like to thank all the participants of the Isle of Wight birth cohort, Dr A Raza and the research team at David Hide Asthma & Allergy Research Centre, Isle of Wight, for collecting the data and Nikki Graham for technical support. We thank the High-Throughput Genomics Group at the Wellcome Trust Centre for Human Genetics (funded by Wellcome Trust grant reference 090532/Z/09/Z and MRC Hub grant G0900747 91070) for the generation of the methylation data.
The epigenetic investigation was supported by National Institutes of Health, USA (R01-AI091905, R01-AI061471, and R21-AI099367). The 18-year assessment of the 1989 Isle of Wight birth cohort was funded by grants from the National Institutes of Health, USA (R01-HL082925) and National Eczema Society/British Dermatological Nursing Group Research Awards.
David Hide Asthma and Allergy Research Centre, St Mary’s Hospital
Clinical and Experimental Sciences, Faculty of Medicine, University of Southampton
Human Development and Health, Faculty of Medicine, University of Southampton
Division of Epidemiology, Biostatistics, & Environmental Health, School of Public Health, University of Memphis
Department of Large Animal Clinical Sciences, Michigan State University
Karmaus W, Dobai AL, Ogbuanu I, Arshard SH, Matthews S, Ewart S: Long-term effects of breastfeeding, maternal smoking during pregnancy, and recurrent lower respiratory tract infections on asthma in children.J Asthma 2008, 45:688–695.PubMed CentralView ArticlePubMed
Wang L, Pinkerton KE: Detrimental effects of tobacco smoke exposure during development on postnatal lung function and asthma.Birth Defects Res C Embryo Today 2008, 84:54–60.View ArticlePubMed
Li X, Howard TD, Zheng SL, Haselkorn T, Peters SP, Meyers DA, Bleecker ER: Genome-wide association study of asthma identifies RAD50-IL13 and HLA-DR/DQ regions.J Allergy Clin Immunol 2010, 125:328–3354.PubMed CentralView ArticlePubMed
Holloway JW, Yang IA, Holgate ST: Genetics of allergic disease.J Allergy Clin Immunol 2010,125(2 Suppl 2):S81-S94.View ArticlePubMed
Nie W, Liu Y, Bian J, Li B, Xiu Q: Effects of polymorphisms -1112C/T and +2044A/G in interleukin-13 gene on asthma risk: a meta-analysis.PLoS One 2013, 8:e56065.PubMed CentralView ArticlePubMed
Cui L, Jia J, Ma CF, Li SY, Wang YP, Guo XM, Li Q, Yu HB, Liu WH, Gao LB: IL-13 polymorphisms contribute to the risk of asthma: a meta-analysis.Clin Biochem 2012, 45:285–288.View ArticlePubMed
Beghe B, Hall IP, Parker SG, Moffatt MF, Wardlaw A, Connolly MJ, Fabbri LM, Ruse C, Sayers I: Polymorphisms in IL13 pathway genes in asthma and chronic obstructive pulmonary disease.Allergy 2010, 65:474–481.View ArticlePubMed
Park H-W, Lee J-E, Kim S-H, Kim Y-K, Min K-U, Kim Y-Y, Cho S-H: Genetic variation of IL13 as a risk factor of reduced lung function in children and adolescents: a cross-sectional population-based study in Korea.Respir Med 2009, 103:284–288.View ArticlePubMed
Corren J, Lemanske RF, Hanania NA, Korenblat PE, Parsey MV, Arron JR, Harris JM, Scheerens H, Wu LC, Su Z, Mosesova S, Eisner MD, Bohen SP, Matthews JG: Lebrikizumab treatment in adults with asthma.N Engl J Med 2011, 365:1088–1098.View ArticlePubMed
Yang G, Volk A, Petley T, Emmell E, Giles-Komar J, Shang X, Li J, Das AM, Shealy D, Griswold DE, Li L: Anti-IL-13 monoclonal antibody inhibits airway hyperresponsiveness, inflammation and airway remodeling.Cytokine 2004, 28:224–232.View ArticlePubMed
Kabesch M, Michel S, Tost J: Epigenetic mechanisms and the relationship to childhood asthma.Eur Respir J 2010, 36:950–961.View ArticlePubMed
Joubert BR, Haberg SE, Nilsen RM, Wang X, Vollset SE, Murphy SK, Huang Z, Hoyo C, Midttun O, Cupul-Uicab LA, Ueland PM, Wu MC, Nystad W, Bell DA, Peddada SD, London SJ: 450 K Epigenome-wide scan identifies differential DNA methylation in newborns related to maternal smoking during pregnancy.Environ Health Perspect 2012, 120:1425–1431.PubMed CentralView ArticlePubMed
Bell CG, Finer S, Lindgren CM, Wilson GA, Rakyan VK, Teschendorff AE, Akan P, Stupka E, Down TA, Prokopenko I, Morison IM, Mill J, Pidsley R, Deloukas P, Frayling TM, Hattersley AT, McCarthy MI, Beck S, Hitman GA, International Type 2 Diabetes 1q C: Integrated genetic and epigenetic analysis identifies haplotype-specific methylation in the FTO type 2 diabetes and obesity susceptibility locus.PLoS One 2010, 5:e14040.PubMed CentralView ArticlePubMed
Sadeghnejad A, Karmaus W, Arshad SH, Kurukulaaratchy R, Huebner M, Ewart S: IL13 gene polymorphisms modify the effect of exposure to tobacco smoke on persistent wheeze and asthma in childhood, a longitudinal study.Respir Res 2008, 9:2.PubMed CentralView ArticlePubMed
Soto-Ramirez N, Arshad SH, Holloway JW, Zhang H, Schauberger E, Ewart S, Patil V, Karmaus W: The interaction of genetic variants and DNA methylation of the interleukin-4 receptor gene increase the risk of asthma at age 18 years.Clin Epigenetics 2013, 5:1.PubMed CentralView ArticlePubMed
Karmaus W, Ziyab AH, Everson T, Holloway JW: Epigenetic mechanisms and models in the origins of asthma.Curr Opin Allergy Clin Immunol 2013, 13:63–69.View ArticlePubMed
Crapo RO, Casaburi R, Coates AL, Enright PL, Hankinson JL, Irvin CG, MacIntyre NR, McKay RT, Wanger JS, Anderson SD, Cockcroft DW, Fish JE, Sterk PJ: Guidelines for methacholine and exercise challenge testing-1999. This official statement of the American Thoracic Society was adopted by the ATS board of directors, July 1999.Am J Respir Crit Care Med 2000, 161:309–329.View ArticlePubMed
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, Jensen R, Johnson DC, MacIntyre N, McKay R, Navajas D, Pedersen OF, Pellegrino R, Viegi G, Wanger J: Standardisation of spirometry.Eur Respir J 2005, 26:319–338.View ArticlePubMed
Raza A, Kurukulaaratchy RJ, Grundy JD, Clayton CB, Mitchell FA, Roberts G, Ewart S, Sadeghnejad A, Arshad SH: What does adolescent undiagnosed wheeze represent? Findings from the Isle of Wight cohort.Eur Respir J 2012, 40:580–588.View ArticlePubMed
Kuan PF, Wang S, Zhou X, Chu H: A statistical framework for Illumina DNA methylation arrays.Bioinformatics 2010, 26:2849–2855.View ArticlePubMed
Wang D, Yan L, Hu Q, Sucheston LE, Higgins MJ, Ambrosone CB, Johnson CS, Smiraglia DJ, Liu S: IMA: an R package for high-throughput analysis of Illumina’s 450K Infinium methylation data.Bioinformatics 2012, 28:729–730.View ArticlePubMed
Michel S, Busato F, Genuneit J, Pekkanen J, Dalphin JC, Riedler J, Mazaleyrat N, Weber J, Karvonen AM, Hirvonen MR, Braun-Fahrlander C, Lauener R, von Mutius E, Kabesch M, Tost J, group Ps: Farm exposure and time trends in early childhood may influence DNA methylation in genes related to asthma and allergy.Allergy 2013, 68:355–364.View ArticlePubMed
Bell JT, Tsai PC, Yang TP, Pidsley R, Nisbet J, Glass D, Mangino M, Zhai G, Zhang F, Valdes A, Shin SY, Dempster EL, Murray RM, Grundberg E, Hedman AK, Nica A, Small KS, Mu TC, Dermitzakis ET, McCarthy MI, Mill J, Spector TD, Deloukas P: Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population.PLoS Genet 2012, 8:e1002629.PubMed CentralView ArticlePubMed
Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan J-B, Shen R: High density DNA methylation array with single CpG site resolution.Genomics 2011, 98:288–295.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.