SERPINA1 methylation and lung function in tobacco-smoke exposed European children and adults: a meta-analysis of ALEC population-based cohorts
Respiratory Research volume 19, Article number: 156 (2018)
The pathophysiological role of SERPINA1 in respiratory health may be more strongly determined by the regulation of its expression than by common genetic variants. A family based study of predominantly smoking adults found methylation at two Cytosine-phosphate-Guanine sites (CpGs) in SERPINA1 gene to be associated with chronic obstructive pulmonary disease risk. The objective of this study was to confirm the association of lung function with SERPINA1 methylation in general population samples by testing a comprehensive set of CpGs in the SERPINA gene cluster. We considered lung function level and decline in adult smokers from three European population-based cohorts and lung function level and growth in tobacco-smoke exposed children from a birth cohort.
DNA methylation using Illumina Infinium Human Methylation 450 k and EPIC beadchips and lung function were measured at two time points in 1076 SAPALDIA, ECRHS and NFBC adult cohort participants and 259 ALSPAC children. Associations of methylation at 119 CpG sites in the SERPINA gene cluster (PP4R4-SERPINA13P) with lung functions and circulating alpha-1-antitripsin (AAT) were assessed using multivariable cross-sectional and longitudinal regression models.
Methylation at cg08257009 in the SERPINA gene cluster, located 32 kb downstream of SERPINA1, not annotated to a gene, was associated with FEV1/FVC at the Bonferroni corrected level in adults, but not in children. None of the methylation signals in the SERPINA1 gene showed associations with lung function after correcting for multiple testing.
The results do not support a role of SERPINA1 gene methylation as determinant of lung function across the life course in the tobacco smoke exposed general population exposed.
Tobacco smoke represents the most important COPD risk factor. Yet, not all smokers develop the disease during their lifetime . The interaction between rare genetically determined alpha-1 antitrypsin (AAT) deficiencies and smoking on emphysema risk illustrates the relevance of genetic susceptibility and gene-environment interactions. The protease inhibitor AAT, encoded by SERPINA1, prevents the extracellular matrix degradation and destruction of alveoli by neutrophil. Both, tobacco smoke and AAT deficiency cause neutrophil recruitment in the lung. Pro-inflammatory and oxidative processes further diminish the anti-proteolytic AAT activity [2, 3].
Genome-wide association studies (GWAS) initially suggested that common SERPINA1 variants might influence COPD risk and associated lung function phenotypes. However, reported associations with common SERPINA1 single nucleotide polymorphisms (SNPs) were shown to reflect their linkage with more penetrant rare variants [4, 5]. In a recent GWAS, restricted to smokers, SERPINA1 variants were not associated with the level of lung function . In the largest GWAS meta-analysis to date, a genetic risk score including 95SNPs, independently associated with lung function or COPD, did not contain SNPs from SERPINA1. Nevertheless, an over-representation of genetic variants related to elastic-fibre pathway was observed .
SERPINA1 gene is complex with eleven splicing isoforms differing in the SERPINA1 5’-UTR, tissue expression and secondary structure. This suggests that the gene’s pathophysiological role may be more strongly determined by differences in expression, regulation or posttranscriptional modification than by genetic variation . A recent family-based study of predominantly smoking adults without severe AAT deficiency investigated SERPINA1 methylation. It was associated with COPD risk, forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC)  at two CpG sites.
Our candidate gene study is the first to investigate the influence of SERPINA1 methylation on lung function in a general population. Specifically, we tested the association of a comprehensive set of methylation signals in the SERPINA gene cluster with lung function levels; with 10 to 15-year lung function decline in adult smokers from three population-based European cohorts; and with lung function growth in tobacco-smoke exposed children from the ALSPAC birth cohort.
Cross-sectional and longitudinal analyses used data from three European adult cohorts and one European birth cohort with longitudinal data on lung function and DNA methylation (DNAm).
Ethics approval and consent to participate
All studies were approved by the local ethics committees and participants or their guardians provided written consent prior to taking part in the study.
Adult cohorts and participants
Participants came from three population-based studies: Swiss Study on Air Pollution and Lung and Heart Disease in Adults (SAPALDIA) [10, 11], European Community Respiratory Health Survey (ECRHS)  and Northern Finland Birth Cohort 1966 . Participation across the studies included structured questionnaires, pre-bronchodilation spirometry, and blood sampling for DNA extraction and analysis. SAPALDIA and ECRHS share a harmonized respiratory health protocol. The study sample was restricted to ever smokers, aged ≥25 years, with data on valid lung function, relevant covariates, and DNA samples from two follow-ups subjected to methylome typing in the context of the Aging Lungs in European Cohorts (ALEC) project. The final sample size was 1076 (n = 561 SAPALDIA, n = 267 ECRHS, and n = 248 NFBC).
Children and adolescent cohort and participants
The Avon Longitudinal Study of Parents and Children (ALSPAC) consisted of 68 follow-up assessments between birth and 18 years . Spirometry was performed at ages 8.5 and 15 years. Participants were restricted to children and adolescents with DNA methylome and valid lung function data from two time points, as well as information on relevant covariates. The study sample included 259 children exposed to tobacco-smoke (mother smoked during pregnancy and/or lived with a smoker and/or reported smoking ≥twice in their lifetime).
Pre-bronchodilation spirometry was performed by trained personnel according to the ATS/ERS recommendations . FEV1, FVC, and FEV1/FVC were the lung function parameters considered. In SAPALDIA, parameters were derived from 2001 and 2010 measurements and corrected for change in spirometers from SensorMedics to ndd EasyOne  in ERCHS they were derived from 1998 and 2008 measurements and were corrected from several spirometers to ndd EasyOne. In NFBC they were measures by Vitalograph P-model in 1997 and MasterScreen Pneumo spirometer in 2012. ALSPAC parameters were derived from 2000 and 2016 spirometries obtained from the same brand spirometer .
DNA methylation measurement
DNAm of autosomes obtained at two time points of lung function measurements was the predictor of interest. DNA was extracted from peripheral whole blood in all cohorts. In SAPALDIA and ALSPAC, DNAm was measured at both time points using the Infinium HumanMethylation450K BeadChip (Illumina, Inc.), in NFBC using the Infinium HumanMethylation450K BeadChip at the first time point and the EPIC BeadChip at the second time point and in ECRHS, at both time points using the EPIC BeadChip.
For ECRHS and SAPALDIA, randomized distribution of samples for bisulphite conversion was applied and for methylome typing batches, the samples from each time point from the same person were placed next to each other on the array. In NFBC1966 DNAm data were recorded in two batches following the clinical assessments of participants aged 31 and 46.
The methylation level (β value) was derived from raw intensities after pre-processing using R package using minfi  followed by beta-mixture quantile normalization (BMIQ)  in SAPALDIA or RnBeads  followed by quantile normalization (QN) in ECRHS. In NFBC, CPACOR pipeline  was used to pre-process and prepare β values. In all adult cohorts, normalized beta scores were regressed on principle components derived from the array control probes reflecting technical bias. The resulting residuals were used as predictors in the associations with lung function.
In ALSPAC, analogous standard quality control procedures were applied, in addition, genotype probes on the HumanMethylation450K were compared between samples from the same individual and against SNP-chip data to identify and remove any sample mismatches. Data were pre-processed in R (version 3.0.1) with the WateRmelon package according to the subset quantile normalization approach  to reduce the non-biological differences between probes. Technical batch effect for each methylation time-point was adjusted for by including ten surrogate variables into the models.
Methylation signals considered in SERPINA gene cluster
The human serine protease inhibitor (serpin) gene cluster is located at 14q32. It consists of eleven functionally diverse serpin genes within a region of approximately 400 kb in length. Gene sub-clusters consist of four, three, and four genes each. The best characterized and proximal sub-cluster with a length of about 107 kb includes SERPINA1 as well as an antitrypsin-related pseudogene (ATR, SERPINA2;∼13 kb downstream), the corticosteroid-binding globulin gene (CBG, SERPINA6;∼68 kb downstream), and the protein Z inhibitor gene (ZPI, SERPINA10;∼100 kb downstream).
Because the tissue-specific expression of different genes in the serpin cluster is regulated by chromosomal elements and chromatin structure, and in the absence of knowledge about the relevance of more distal methylation signals on gene expression, we included in the analysis all 119 CpGs located 99 kb downstream (PPP4R4) and 376 kb upstream (GSC) from the SERPINA1 gene. The CpGs were allocated to 12 genes: PPP4R4, SERPINA10, SERPINA6, SERPINA1, SERPINA11, SERPINA9, SERPINA12, SERPINA4, SERPINA5, SERPINA3, SERPINA13, and GSC (Fig. 1: location of CpG sites considered; Fig. 2: correlation of methylation at different CpGs at both time points).
Information on participant’s age, sex, education, height, and smoking status was derived from questionnaires administered during the two time-points. Cell proportions in the respective blood samples were estimated using the Houseman method  implemented in the minfi package .
For sensitivity analysis, SAPALDIA provided information on concentrations of high-sensitive C-reactive protein (CRP) and AAT in blood samples collected at the first time point. CRP and AAT concentrations were measured by latex-enhanced immunoturbidimetric assays (Roche Cobas Integra analyzer; Roche Diagnostics). Inter-assay coefficients of variation were below 5% and lower detection rate were 0.1 mg/l (CRP) and 0.21 (AAT) . In addition, SAPALDIA provided genetic information on rare AAT deficiency variants PiZ and PiS .
Cohort-specific analyses on lung function and its decline in adults
Within each adult cohort, four statistical models were run: a) cross-sectional associations of DNAm with lung function at T1, b) cross-sectional associations of DNAm with lung function at T2, c) repeat cross-sectional associations of DNAm and lung function at both time points, and d) predictive associations of DNAm at T1 on annual lung function decline between T1 and T2, calculated as (lung function at T2 - lung function at T1 divided by the time of follow-up in years).
The absolute levels of lung functions (FEV1, FVC, FEV1/FVC) were regressed on residuals of methylation by fitting linear regression models and a mixed-model with a random effect for the subject (repeat cross-sectional model), respectively. Sample size was kept constant across all models and analyses. Models were a priori adjusted for the following covariates associated with lung function level at P-value < 0.05: study center, age, age2, education, height, (height-mean(height))2, sex, sex*age, sex*age2, sex*height, sex*(height-mean(height))2, and cell composition (CD8T; CD4T; NK; Bcell; Mono and Eos).
Meta-analysis of adult cohort results
The fixed-effect meta-analysis weighted by the inverse of variance was completed using METAL on the 119CpG sites common to all three cohorts . Associations with P-values< 0.05 were considered as nominally significant. Since methylation of CpGs in the SERPINA cluster did not correlate throughout the whole genetic region (Fig. 2), the total number of CpGs was used as correction for multiple testing. A Bonferroni corrected P-value< 4.2 × 10− 4 was considered statistically significant. Given the generally medium-to-high correlation between lung function parameters, associations of CpGs with different lung function parameters were not considered independent .
Cohort-specific analysis on lung function and its growth in children and adolescents
Four regression models equivalent to those assessed in adult cohorts were run to investigate cross-sectional and predictive associations of methylation at 119 CpGs in the SERPINA gene cluster with the absolute level and increase in FEV1, FVC, and FEV1/FVC in children and adolescents. The models were a priori adjusted for covariates associated with lung function at P-value< 0.05: study center, age, mother’s education, height, (height-mean(height))2, sex, sex*age, sex*height, and estimated cell composition (CD8T; CD4T; NK; Bcell; Mono; Eos) .
Seven sensitivity analyses were conducted in SAPALDIA:(1) 119 CpGs from T1 were regressed on circulating AAT measured at T1  by including the same set of covariates as in cross-sectional models on lung function at T1;(2) cross-sectional models at T1 were additionally adjusted for CRP and AAT;(3) for all phenotypes, all four statistical models were additionally adjusted for the presence of PiS and PiZ: 0, 1 or 2 alleles;(4) for neutrophils;(5) stratified by gender;(6) stratified by obesity (BMI 30 < vs. ≥ 30 kg/m2;(7) and stratified by self-reported asthma (“Have you ever had asthma?”).
Characteristics of adult cohort participants
Characteristics of the adult participants at the first time point (T1) and second time point (T2) are presented in Table 1. SAPALDIA participants were oldest, reported the highest number of pack-years, had lung function levels similar to ECHRS, but lower than NFBC, and exhibited the highest prevalence of COPD based on pre-bronchodilation spirometry and defined as FEV1/FVC below the lower limit of normal: in SAPALDIA 15.3% (T1) and 17.1% (T2), respectively, compared to 10.5% (T1) and 12.7% (T2) in ECHRS and 2.2% (T1) and 11.4% (T2) in NFBC. The prevalence of self-reported doctor’s diagnosed asthma was between 10 and 20%, increasing with aging in all three cohorts.
Association of methylation at 119 CpGs in the SERPINA1 cluster with lung function in adult ever smokers
Of the 119 CpGs, 17 were located in the SERPINA1 gene. DNAm at these 17 sites was not associated with any of the three lung function parameters at a Bonferroni-corrected P-value, irrespective of the model considered (Tables 2 and 3; Additional file 1: Table ES1). However, for FEV1, meta-analysis revealed three nominally significant associations at T1 (cg09968361, cg25968219, cg04179148) (Table 2). The positive association of cg09968361 with levels of FEV1, nominally significant at T1 (β:1.43, P-value = 0.04) was consistent in direction for T2 and the repeat cross-sectional analysis (β:0.72, P-value = 0.02). However for change in lung function, an increase in methylation at this CpG site was associated with accelerated FEV1 decline (β:-0.10, P-value = 0.011). For FEV1/FVC (Table 3), meta-analysis also revealed three nominally significant associations at T1 (cg25968219, cg24621042, cg04179148), which in part overlapped with those associated with FEV1, also in terms of direction of association. No association of methylation at SERPINA1 CpGs and FVC was observed (Additional file 1: Table ES1).
SERPINA gene cluster
Results from the meta-analysis on the association of methylation at the 119 CpGs with cross-sectional lung function and lung function decline for FEV1, FVC and FEV1/FVC are presented in Additional file 1: Tables ES2-S4. A single CpG at cg08257009, located 32 kb downstream of SERPINA1, withstood Bonferroni-correction for multiple testing. Methylation at this site was positively associated with FEV1/FVC in the repeat cross-sectional analysis (β:0.11; P-value = 2.6 × 10− 4) (Additional file 1: Table ES3). The associations of this signal for the two cross-sectional time points were comparable (β:0.10; P-value = 0.01 at T1; β:0.14; P-value = 2.2 × 10− 3 at T2). Consistent with the observation that higher methylation at this site was associated with higher level of FEV1/FVC cross-sectionally, it also predicted attenuated decline of FEV1/FVC (β:0.0089; P-value = 9.1 × 10− 3).
Association of methylation at 119 CpGs with circulating alpha-1 antitrypsin in adult ever smokers in SAPALDIA
None of the associations between methylation and circulating AAT in SAPALDIA participants reached multiple-testing-corrected statistical significance, but nominal statistical significance was observed at four CpG sites in the SERPINA1 gene (Tables 2 and 3) and at two additional sites outside SERPINA1 (Additional file 1: Tables ES2-S4), one of which was cg08257009, the only lung function associated signal withstanding multiple testing. In all instances, higher methylation was associated with lower AAT concentrations. These inverse associations with circulating AAT did not translate into statistically significant and inverse associations with any of the lung function parameters measured.
Association of methylation at 119 CpGs with lung function in tobacco-smoke exposed children and adolescents
Characteristics of ALSPAC children and adolescents from the first (T1) and second (T2) time points are presented in Table 4. Height and weight increased during follow-up from a mean of 133 to 170 cm and from 31 to 63 kg, respectively. Twenty-two and 29% of children and adolescents had asthma and at the second time point, the majority (86%) of adolescents reported intake of asthma medication in the last 12 months. While FEV1 and FVC increased during follow-up, on average, FEV1/FVC remained stable.
Focusing on the methylation sites in SERPINA1 (Additional file 1: Tables ES5-S7), again none of the signals reached Bonferroni-corrected significance in association with any lung function parameter or model. There were several nominally statistically significant associations, particularly for lung function growth (FEV1: cg26938334; FVC: cg26938334, cg13826459, cg 24,621,042; FEV1/FVC: cg10070185). Methylation at cg26938334 was inversely associated with growth in FEV1 and FVC, and methylation at cg24621042 and cg10070185 was inversely associated with growth in FVC and in FEV1/FVC, respectively. The associations of cg24621042 with FEV1/FVC change were inconsistent with those observed for FEV1/FVC decline in adults.
SERPINA gene cluster
Results from meta-analysis on the association of methylation at all 119 CpGs with cross-sectional lung function and lung function growth of FEV1, FVC and FEV1/FVC are presented in Additional file 1: Tables ES5-S7. No methylation signals remained statistically significant after Bonferroni correction. Methylation at cg08257009, 32 kb downstream of SERPINA1 gene, which showed statistically significant positive association with FEV1/FVC cross-sectionally in the adult cohort, was not associated with FEV1 or FEV1/FVC in children. However, there was evidence for an inverse, rather than a positive association with FVC at T1 (β:-0.91, P-value = 4.4 × 10− 3, Additional file 1: Table ES6). Overall, none of the SERPINA1 CpGs showed consistent associations across cross-sectional and longitudinal models.
Replication of previously reported COPD signals
Relative hypomethylation of cg02181506 and cg24621042 were previously associated with COPD . Table 5 summarizes the association of both signals with circulating AAT in SAPALDIA and with lung function and its change in adults as well as children. Contrary to expectation, relative hypermethylation, not hypomethylation at these two sites was associated with lower circulating AAT. For cg02181506, irrespective of age, no association of methylation with either FEV1/FVC or FEV1 was observed. For cg24621042, in adults, hypermethylation was positively associated with FEV1/FVC at T1 only, but not in the repeat cross-sectional analysis, and the inverse association of methylation at cg24621042 and FVC change was only observed in children. No lung function associations reached statistical significance at the nominal P-value< 0.05.
Sensitivity analyses in the SAPALDIA cohort
In the absence of consistent CpG and lung function associations, we restricted sensitivity analyses to the FEV1/FVC association between the CpG withstanding multiple testing (cg08257009) and the two CpGs (cg02181506 and cg24621042) previously associated with COPD  (Table 6, Additional file 1: Tables ES8-S9). Sensitivity analyses did not reveal a more consistent pattern of associations between methylation at these three sites and lung function.
This is the first study to investigate the association of DNAm in the SERPINA gene cluster with lung function and its longitudinal change in ever-smoking adults from three European population-based adult cohort studies and in tobacco-smoke exposed children and adolescents from England. We observed methylation at cg08257009 not annotated to a gene in the SERPINA gene cluster, located 32 kb downstream of the SERPINA1 gene to be significantly associated with FEV1/FVC at the Bonferroni-corrected level. No methylation signals in the SERPINA1 gene showed associations with lung function level or change over time after correcting for multiple testing.
Few studies have reported associations between DNAm and lung function or COPD. Our results obtained in general population samples contradict the previous finding of hypomethylation at two CpG sites in SERPINA1 and COPD risk in two family-based studies. The first of these studies was restricted to smokers and the second consisted of participants with and without a history of smoking . The functional relevance of the two CpG sites remains unclear, given that the associations of hypomethylation with COPD risk and circulating AAT are inconsistent in direction. Consistent with the inverse association observed with AAT in the blood, hypomethylation of the AAT gene has been associated with increased gene expression in rat models .
Our results are consistent with the results from a study comparing gene methylation in lung tissue of former smoker COPD patients and controls with normal lung function [27, 28]. Methylation in SERPINA1 was not associated with COPD risk. Similarly no cross-sectional association with post-bronchodilation lung function was observed for SERPINA1 gene methylation measured in blood samples from a small rural Korean COPD case-control study . Furthermore, in a study of middle-aged monozygotic twins, DNAm in SERPINA1 was not associated with intra-pair differences in lung function decline .
Smoking exerts strong adverse effects on lung function and increases COPD risk. In a small proportion of smokers, it interacts with genetically determined rare AAT deficiency to cause COPD. The adverse respiratory effects on the lung tissue may therefore be mediated by altering DNAm in SERPINA1. Yet, the largest epigenome-wide association study for smoking did not identify a relevant association between epigenetic signatures in SERPINA1 and smoking . Interestingly, the cg08257009 which reached Bonferroni corrected significance in the repeat cross-sectional analysis was previously reported as a smoking-related CpG in buccal cells in epithelial cancer . In addition, cg08257009 and cg02181506 were reported to be smoking related signals in blood-derived DNA samples from participants in 16 cohorts, with cg02181506 showing weaker associations .
The use of blood instead of lung tissue to assess lung disease related methylation may be considered a limitation of the study. However as previously discussed, COPD is a systemic disease related to low grade inflammation, which supports studying the peripheral blood methylome . Teschendorff and others demonstrated the correlation between smoking related DNAm in blood and lung tissue as well as between normal and malignant lung tissue. Their results also pointed to the prognostic value of these signals in lung cancer patients [32,33,34].
Additional limitations of the study include the lack of cross–omics data to investigate the biological network related to AAT more comprehensively. This is important, given the fact that the most recent lung function GWAS suggests an important role of elastic-fibre pathways . Furthermore, misclassification of smoking exposure could have biased the observed associations, most likely towards the null.
The advantages of the current study are several fold. First, the studies providing data for this investigation are well established respiratory population-based cohorts. They are known for high quality testing of lung function and large sample sizes. They share similarities in study protocols. Methylome analyses followed stringent quality control. Second, the study is based on longitudinal lung function data and prospectively and repeatedly measured DNAm. Repeated assessment of lung function and predictors may have helped improve the statistical power of the cross-sectional analysis. Prospective assessment of DNAm with change in lung function decreases the problem of reverse causation. Third, the integration of lung function from both adults and children offered the opportunity to study SERPINA1 methylation and lung function over the life course and investigate whether relevant time windows for genome-environment interactions may exist.
In conclusion, this first comprehensive study on DNAm in the SERPINA gene cluster provides weak evidence of an association with lung functions and its change across the life course. Larger studies based on post-bronchodilation lung function need to be followed by investigating associations of SERPINA1 and elastin-related pathways and networks through cross–omics approaches.
Salvi SS, Barnes PJ. Chronic obstructive pulmonary disease in non-smokers. Lancet. 2009;374:733–43.
Hazari YM, Bashir A, Habib M, Bashir S, Habib H, Qasim MA, Shah NN, Haq E, Teckman J, Fazili KM. Alpha-1-antitrypsin deficiency: genetic variations, clinical manifestations and therapeutic interventions. Mutat Res. 2017;773:14–25.
Sandhaus RA, Stoller JK. Introduction to the 50th anniversary of the description of alpha-1 antitrypsin deficiency. COPD. 2013;10(Suppl 1):1–2.
Thun GA, Imboden M, Ferrarotti I, Kumar A, Obeidat M, Zorzetto M, Haun M, Curjuric I, Couto Alves A, Jackson VE, et al. Causal and synthetic associations of variants in the SERPINA gene cluster with alpha1-antitrypsin serum levels. PLoS Genet. 2013;9:e1003585.
Busch R, Hobbs BD, Zhou J, Castaldi PJ, McGeachie MJ, Hardin ME, Hawrylkiewicz I, Sliwinski P, Yim JJ, Kim WJ, et al. Genetic association and risk scores in a chronic obstructive pulmonary disease Meta-analysis of 16,707 subjects. Am J Respir Cell Mol Biol. 2017;57:35–46.
Lutz SM, Cho MH, Young K, Hersh CP, Castaldi PJ, McDonald ML, Regan E, Mattheisen M, DeMeo DL, Parker M, et al. A genome-wide association study identifies risk loci for spirometric measures among smokers of European and African ancestry. BMC Genet. 2015;16:138.
Wain LV, Shrine N, Artigas MS, Erzurumluoglu AM, Noyvert B, Bossini-Castillo L, Obeidat M, Henry AP, Portelli MA, Hall RJ, et al. Genome-wide association analyses for lung function and chronic obstructive pulmonary disease identify new loci and potential druggable targets. Nat Genet. 2017;49:416–25.
Corley M, Solem A, Phillips G, Lackey L, Ziehr B, Vincent HA, Mustoe AM, Ramos SBV, Weeks KM, Moorman NJ, Laederach A. An RNA structure-mediated, posttranscriptional model of human alpha-1-antitrypsin expression. Proc Natl Acad Sci U S A. 2017;114:E10244–53.
Qiu W, Baccarelli A, Carey VJ, Boutaoui N, Bacherman H, Klanderman B, Rennard S, Agusti A, Anderson W, Lomas DA, DeMeo DL. Variable DNA methylation is associated with chronic obstructive pulmonary disease and lung function. Am J Respir Crit Care Med. 2012;185:373–81.
Martin BW, Ackermann-Liebrich U, Leuenberger P, Kunzli N, Stutz EZ, Keller R, Zellweger JP, Wuthrich B, Monn C, Blaser K, et al. SAPALDIA: methods and participation in the cross-sectional part of the Swiss study on air pollution and lung diseases in adults. Soz Praventivmed. 1997;42:67–84.
Ackermann-Liebrich U, Kuna-Dibbert B, Probst-Hensch NM, Schindler C, Felber Dietrich D, Stutz EZ, Bayer-Oglesby L, Baum F, Brandli O, Brutsche M, et al. Follow-up of the Swiss Cohort Study on Air Pollution and Lung Diseases in Adults (SAPALDIA 2) 1991-2003: methods and characterization of participants. Soz Praventivmed. 2005;50:245–63.
Janson C, Anto J, Burney P, Chinn S, de Marco R, Heinrich J, Jarvis D, Kuenzli N, Leynaert B, Luczynska C, et al. The European Community respiratory health survey: what are the main results so far? European Community respiratory health survey II. Eur Respir J. 2001;18:598–611.
University of Oulu. Northern Finland Cohorts. 31-year follow-up study. 2018. http://www.oulu.fi/nfbc/node/18097. Accessed 18 Mar 2018.
Boyd A, Golding J, Macleod J, Lawlor DA, Fraser A, Henderson J, Molloy L, Ness A, Ring S, Davey Smith G. Cohort profile: the ‘children of the 90s’--the index offspring of the Avon longitudinal study of parents and children. Int J Epidemiol. 2013;42:111–27.
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38.
Bridevaux PO, Dupuis-Lozeron E, Schindler C, Keidel D, Gerbase MW, Probst-Hensch NM, Bettschart R, Burdet L, Pons M, Rothe T, et al. Spirometer replacement and serial lung function measurements in population studies: results from the SAPALDIA study. Am J Epidemiol. 2015;181:752–61.
Kotecha SJ, Watkins WJ, Lowe J, Henderson AJ, Kotecha S. Effect of early-term birth on respiratory symptoms and lung function in childhood and adolescence. Pediatr Pulmonol. 2016;51:1212–21.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.
Assenov Y, Muller F, Lutsik P, Walter J, Lengauer T, Bock C. Comprehensive analysis of DNA methylation data with RnBeads. Nat Methods. 2014;11:1138–40.
Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan ST, Afzal U, Scott J, Jarvelin MR, Elliott P, et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015;16:37.
Touleimat N, Tost J. Complete pipeline for Infinium((R)) human methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012;4:325–41.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, Wiencke JK, Kelsey KT. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.
Senn O, Russi EW, Schindler C, Imboden M, von Eckardstein A, Brandli O, Zemp E, Ackermann-Liebrich U, Berger W, Rochat T, et al. Circulating alpha1-antitrypsin in the general population: determinants and association with lung function. Respir Res. 2008;9:35.
Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–1.
Barton DE, Francke U. Activation of human alpha 1-antitrypsin genes in rat hepatoma x human fibroblast hybrid cell lines is correlated with demethylation. Somat Cell Mol Genet. 1987;13:635–44.
Morrow JD, Glass K, Cho MH, Hersh CP, Pinto-Plata V, Celli B, Marchetti N, Criner G, Bueno R, Washko G, et al. Human lung DNA methylation quantitative trait loci Colocalize with COPD genome-wide association loci. Am J Respir Crit Care Med. 2018;197:1275–84.
Morrow JD, Cho MH, Hersh CP, Pinto-Plata V, Celli B, Marchetti N, Criner G, Bueno R, Washko G, Glass K, et al. DNA methylation profiling in human lung tissue identifies genes associated with COPD. Epigenetics. 2016;11:730–9.
Lee MK, Hong Y, Kim SY, Kim WJ, London SJ. Epigenome-wide association study of chronic obstructive pulmonary disease and lung function in Koreans. Epigenomics. 2017;9:971–84.
Bolund ACS, Starnawska A, Miller MR, Schlunssen V, Backer V, Borglum AD, Christensen K, Tan Q, Christiansen L, Sigsgaard T. Lung function discordance in monozygotic twins and associated differences in blood DNA methylation. Clin Epigenetics. 2017;9:132.
Joehanes R, Just AC, Marioni RE, Pilling LC, Reynolds LM, Mandaviya PR, Guan W, Xu T, Elks CE, Aslibekyan S, et al. Epigenetic signatures of cigarette smoking. Circ Cardiovasc Genet. 2016;9:436–47.
Teschendorff AE, Yang Z, Wong A, Pipinikas CP, Jiao Y, Jones A, Anjum S, Hardy R, Salvesen HB, Thirlwell C, et al. Correlation of smoking-associated DNA methylation changes in Buccal cells with DNA methylation changes in epithelial Cancer. JAMA Oncol. 2015;1:476–85.
Stueve TR, Li WQ, Shi J, Marconett CN, Zhang T, Yang C, Mullen D, Yan C, Wheeler W, Hua X, et al. Epigenome-wide analysis of DNA methylation in lung tissue shows concordance with blood studies and identifies tobacco smoke-inducible enhancers. Hum Mol Genet. 2017;26:3014–27.
Russo AL, Thiagalingam A, Pan H, Califano J, Cheng KH, Ponte JF, Chinnappan D, Nemani P, Sidransky D, Thiagalingam S. Differential DNA hypermethylation of critical genes mediates the stage-specific tobacco smoke-induced neoplastic progression of lung cancer. Clin Cancer Res. 2005;11:2466–70.
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) for the generation of Methylation data in SAPALDIA and the National Public Health Institute, Biomedicum Helsinki in Finland to have performed DNA extractions, sample quality controls, biobank up-keeping and aliquotting in NFBC. We would like to acknowledge the scientific teams, the field workers and the administrative staff of the SAPALDIA, ECRHS, NFBC, and ALSPAC teams. We are most grateful to all study participants and their families for having donated time, data, and biospecimens to these studies. We thank the late Professor Paula Rantakallio for the launch of NFBC1966.
The ALEC Study is funded by the European Union’s Horizon 2020 Research and Innovation programme under grant agreement No. 633212. The ALEC study leader is Deborah Jarvis. The manuscript was done under ALEC Work-package 5 led by Nicole Probst-Hensch. Other Work-package leaders in ALEC are Cecilie Svanes, John Henderson, Judith Garcia-Aymerich, and Cosetta Minelli. The ALEC International Scientific Advisory Board is: Marike Boezen (University Medical Center Groningen, University of Groningen, Groningen, The Netherlands); Bernice Elger (Institute for Biomedical Ethics, University of Basel, Basel, Switzerland); Bo Alexander Gleditsch (The Norwegian Asthma and Allergy Association, Norway); Bas Heijmans (Department of Medical Statistics and Bioinformatics, Leiden University Medical Center, Leiden, The Netherlands); Isabelle Romieu (National Institute of Public Health, Cuernavaca, Mexico; and Emory University, Atlanta, US); John Thompson (Department of Health Sciences, University of Leicester, Leicester, UK).
The SAPALDIA Cohort is funded by the Swiss National Science Foundation (grants no 33CS30–148470/1&2, 33CSCO-134276/1, 33CSCO-108796,, 324730_135673, 3247BO-104283, 3247BO-104288, 3247BO-104284, 3247–065896, 3100–059302, 3200–052720, 3200–042532, 4026–028099, PMPDP3_129021/1, PMPDP3_141671/1), the Federal Office for the Environment, the Federal Office of Public Health, the Federal Office of Roads and Transport, the canton’s government of Aargau, Basel-Stadt, Basel-Land, Geneva, Luzern, Ticino, Valais, and Zürich, the Swiss Lung League, the canton’s Lung League of Basel Stadt/ Basel Landschaft, Geneva, Ticino, Valais, Graubünden and Zurich, Stiftung ehemals Bündner Heilstätten, SUVA, Freiwillige Akademische Gesellschaft, UBS Wealth Foundation, Talecris Biotherapeutics GmbH, Abbott Diagnostics, Klinik Barmelweid, Hirslanden Klinik Aarau, European Commission 018996 (GABRIEL), Wellcome Trust WT 084703MA, Exposomics EC FP7 grant(Grant agreement No: 308610).
The ECRHS cohort was supported by France: Ministère de la Santé, Glaxo France, Insitut Pneumologique d’Aquitaine, Contrat de Plan Etat-Région Languedoc-Rousillon, CNMATS, CNMRT (90MR/10, 91AF/6), Ministre Délégué de la Santé, RNSP, GSF, and Programme Hospitalier de Recherche Clinique National 2010. France (Grenoble): Direction de la Recherche Clinique de Grenoble 2000 (no. 2610), Ministère de l’Emploi et de la Solidarité, Direction Générale de la Sante, CHU Grenoble, Comite des Maladies Respiratoires de l’Isère, and Comité Scientifique AGIRadom 2011. France (Paris): Ministère de l’Emploi et de la Solidarité, Direction Générale de la Sante, Union Chimique Belge-Pharma, Aventis, Glaxo France, Agence Nationale de la Santé, Région Ile de France, and Domaine d’intérêt majeur. Germany: Bundesminister für Forschung und Technologie. Germany (Erfurt): DFG—German Research Foundation (FR1526/1–1, HE 3294/10–1). Norway (Bergen): Norwegian Research Council (no. 101422/310, no. 214123), Norwegian Asthma and Allergy Association, Glaxo Wellcome AS, Norway Research Fund, Western Norway Regional Health Authorities (no. 911631), and the Bergen Medical Research Foundation. Spain: Ministerio de Sanidad y Consumo FIS (no. 91/0016060/00E-05E, no. 93/0393, no. 97/0035–01, no. 99/0034–01, no. 99/0034–02). Spain (Galdakao): Basque Health Department and FIS (no. 09/01511). Sweden (Gothenburg, Umeå, and Uppsala): the Swedish Medical Research Council, Swedish Heart-Lung Foundation, Swedish Association against Asthma and Allergy, Swedish Cancer and Allergy Foundation, and Swedish Council for Working Life and Social Research. European Commission (no. 018996–GABRIEL); and the Wellcome Trust (WT084703MA). United Kingdom: Asthma UK (formerly known as National Asthma Campaign), Department of Health, South Thames Regional Health Authority, and the Medical Research Council (G0901214/1). The coordination of the European Community Respiratory Health Survey II was supported by the European Commission. The coordination of ECRHS III was supported by the Medical Research Council (G0901214/1).
NFBC1966 received financial support related to the present work from the Academy of Finland (project grants 104781, 120315, 129269, 1114194, 24300796, 85547, Center of Excellence in Complex Disease Genetics), University Hospital Oulu, Biocenter, University of Oulu, Finland (75617), NHLBI grant 5R01HL087679–02 through the STAMPEED program (1RL1MH083268–01), NIH/NIMH (5R01MH63706:02), the EU FP5 EURO-BLCS, QLG1-CT-2000-01643, ENGAGE project and grant agreement HEALTH-F4–2007-201413, EU FP7 EurHEALTHAgeing 277849, the Medical Research Council (MRC), UK (G0500539, G0600705, G1002319, PrevMetSyn/SALVE) and ERDF European Regional Development Fund Grant no. 539/2010 A31592. The program is currently being funded by the EU H2020--PHC-2014 DynaHEALTH action (grant agreements No. 633595), EU H2020-HCO-2004 iHEALTH Action, EU H2020-PHC-2014 ALEC Action (grant agreement No. 633212), EU H2020-SC1–2016-2017 LIFECYCLE Action, EU H2020-MSCA-ITN-2016 CAPICE Action, Academy of Finland EGEA-project (285547) and MRC Grant nro MR/M013138/1. In NFBC, DNA extractions, sample quality controls, biobank up-keeping and aliquotting was supported financially by the Academy of Finland and Biocentrum Helsinki.
The ALSPAC work was supported by the MRC Integrative Epidemiology Unit, The UK Medical Research Council and Wellcome (Grant ref.: 102215/2/13/2) and the University of Bristol (MC_UU_12013_2, MC_UU_12013_5 and MC_UU_12013_8).
The Accessible Resource for Integrated Epigenomics Studies (ARIES) which generated large scale methylation data was funded by the UK Biotechnology and Biological Sciences Research Council (BB/I025751/1 and BB/I025263/1). Additional epigenetic profiling on the ALSPAC cohort was supported by the UK Medical Research Council Integrative Epidemiology Unit and the University of Bristol (MC_UU_12013_1, MC_UU_12013_2, MC_UU_12013_5 and MC_UU_12013_8), the Wellcome Trust (WT088806) and the United States National Institute of Diabetes and Digestive and Kidney Diseases (R01 DK10324).
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request and with permission of principal investigators of each cohort.
Ethics approval and consent to participate
All studies were approved by the local ethics committees and participants provided written consent prior to taking part in the study.
Consent for publication
The following authors report no competing interests: ABB, MI, FR, MW, AFSA, AJ, ES, JA, SS, VK, RB, AT, MP, FK, RA, GCS, CR, MRJ, NMPH.
DS reports consultancy at Roche AG, Novartis AG, Astra-Zeneca AG, GSK AG; grants and pending grants at Astra-Zeneca AG, Pan Gas AG, Weimann AG, Curetis AG, Swiss National foundation (PP00-P3–128412/1), Astra Zeneca AG and payments for lectures including speakers bureaus at Novartis AG.
AJH reports grants from Medical Research Council and grants from Wellcome during the conduct of the study.
DJ and JWH report grants from European Union during the conduct of the study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Anna Beckmeyer-Borowko, Medea Imboden, Faisal I. Rezwan, Matthias Wielscher, Andre F.S. Amaral are equal first author.
Marjo-Riitta Jarvelin, Deborah Jarvis, John W. Holloway, Nicole M. Probst-Hensch are equal senior author.
Table ES1. Meta-Analysis of the association of methylation at CpGs in the SERPINA* Gene Cluster with FVC level and decline in adult ever smokers from SAPALDIA, ECRHS and NFBC (n = 1076) and with circulating AAT concentrations in SAPALDIA (n = 561). Table ES2. Meta-Analysis of the association of methylation at 119 CpGs in the SERPINA* cluster with FEV1 level and decline in adult ever smokers from SAPALDIA, ECRHS and NFBC (n = 1076) and with circulating AAT concentrations in SAPALDIA (n = 561). Table ES3. Meta-analysis of the association of methylation at 119 CpGs in the SERPINA* gene cluster with FVC level and decline in adult ever smokers from SAPALDIA, ECRHS and NFBC (n = 1076) and with circulating AAT concentrations in SAPALDIA (n = 561). Table ES4. Meta-analysis of the association of methylation at 119 CpGs in the SERPINA* gene cluster with FEV1/FVC level and decline in adult ever smokers from SAPALDIA, ECRHS and NFBC (n = 1076) and with circulating AAT concentrations in SAPALDIA (n = 561). Table ES5. Association of methylation at 119 CpGs in the SERPINA* gene cluster with FEV1 level and decline in ALSPAC children exposed to tobacco smoke (n = 259). Table ES6. Association of methylation at 119 CpGs in the SERPINA* gene cluster with FVC level and decline in in ALSPAC children exposed to tobacco smoke (n = 259). Table ES7. Association of methylation at 119 CpGs in the SERPINA* gene cluster with FEV1/FVC level and decline in ALSPAC children exposed to tobacco smoke (n = 259). Table ES8. Association of methylation at 119 CpGs in the SERPINA* gene cluster with FEV1 level and decline in adult smokers from SAPALDIA, basic adjustment and adjustment for PIS and PIZ genotypes. Table ES9. Association of methylation at 119 CpGs in the SERPINA* gene cluster with FEV1/FVC level and decline in adult smokers from SAPALDIA, basic adjustment and adjustment for PIS and PIZ genotypes. (DOCX 349 kb)
About this article
Cite this article
Beckmeyer-Borowko, A., Imboden, M., Rezwan, F.I. et al. SERPINA1 methylation and lung function in tobacco-smoke exposed European children and adults: a meta-analysis of ALEC population-based cohorts. Respir Res 19, 156 (2018). https://doi.org/10.1186/s12931-018-0850-8