- Research article
- Open Access
Genome-wide association study of lung function and clinical implication in heavy smokers
BMC Medical Genetics volume 19, Article number: 134 (2018)
The aim of this study is to identify genetic loci associated with post-bronchodilator FEV1/FVC and FEV1, and develop a multi-gene predictive model for lung function in COPD.
Genome-wide association study (GWAS) of post-bronchodilator FEV1/FVC and FEV1 was performed in 1645 non-Hispanic White European descent smokers.
A functional rare variant in SERPINA1 (rs28929474: Glu342Lys) was significantly associated with post-bronchodilator FEV1/FVC (p = 1.2 × 10− 8) and FEV1 (p = 2.1 × 10− 9). In addition, this variant was associated with COPD (OR = 2.3; p = 7.8 × 10− 4) and severity (OR = 4.1; p = 0.0036). Heterozygous subjects (CT genotype) had significantly lower lung function and higher percentage of COPD and more severe COPD than subjects with the CC genotype. 8.6% of the variance of post-bronchodilator FEV1/FVC can be explained by SNPs in 10 genes with age, sex, and pack-years of cigarette smoking (P < 2.2 × 10− 16).
This study is the first to show genome-wide significant association of rs28929474 in SERPINA1 with lung function. Of clinical importance, heterozygotes of rs28929474 (4.7% of subjects) have significantly reduced pulmonary function, demonstrating a major impact in smokers. The multi-gene model is significantly associated with CT-based emphysema and clinical outcome measures of severity. Combining genetic information with demographic and environmental factors will further increase the predictive power for assessing reduced lung function and COPD severity.
Chronic obstructive pulmonary disease (COPD) is a common respiratory disease caused by the interaction of genetic susceptibility with environmental influences, primarily tobacco exposure. COPD is defined as a reduced ratio of post-bronchodilator forced expiratory volume in 1 s (FEV1) to forced vital capacity (FVC) (post-bronchodilator FEV1/FVC < 0.70) . COPD severity is measured by the reduction in post-bronchodilator percent predicted FEV1, i.e., GOLD stages 1–4 (mild, moderate, severe, and very severe COPD) have post-bronchodilator percent predicted FEV1 ≥ 80%, ≥50%, ≥30%, or < 30%, respectively .
Twenty-eight genomic loci associated with baseline FEV1/FVC or FEV1 have been identified by meta-analyses of genome-wide association studies (GWAS) in general populations of European descent [2,3,4]. A recent GWAS comparing extremes of high and low baseline FEV1 in subjects of European ancestry from the UK Biobank has identified five loci (KANSL1, HLA-DQ, NPNT, TET2, and TSEN54) in never smokers and RBM19-TBX5 in heavy smokers . HHIP, FAM13A1, CHRNA3, RIN3, MMP12, and TGFB2 have been associated with COPD at genome-wide significant levels . To our knowledge, no GWAS study has been performed on post-bronchodilator FEV1/FVC and FEV1 in smokers, which defines a diagnosis of COPD and determines COPD severity, respectively.
GWAS of post-bronchodilator FEV1/FVC and percent predicted FEV1 were performed in non-Hispanic White smokers (n = 1645, GOLD stage 0–4, smoking≥20 packs/year) from the NHLBI-sponsored SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS). In addition to evaluating previously reported loci associated with baseline lung function in general populations, we aimed to identify novel genes associated with abnormal post-bronchodilator lung function in smokers enriched for COPD and develop a model to predict lung function using multiple genes and demographic/environmental factors.
SPIROMICS is a prospective cohort study that enrolled 2981 participants with the goals of identifying new COPD subgroups and intermediate markers of disease progression [7, 8]. SPIROMICS is a well-characterized longitudinal cohort with comprehensive phenotyping including measurements of lung function and quantitative CT scan. Spirometry was performed before and after four inhalations with 90 μg albuterol and 18 μg ipratropium per inhalation according to ATS recommendations. Non-Hispanic White smokers (ever or current smoking≥20 packs/year) with genotyping information available were included in this analysis. Smokers with COPD were defined as smokers (smoking≥20 packs/year) with post-bronchodilator FEV1/FVC < 0.7 (GOLD stage 1–4) and ‘healthy’ smoking controls were defined as smokers (smoking≥20 packs/year) with post-bronchodilator FEV1/FVC ≥ 0.7 (GOLD stage 0). DNA was isolated using standard protocols, and SNP genotyping performed using Illumina HumanOmniExpressExome BeadChip and BeadStudio (Illumina, Inc., San Diego, CA).
Participants were recruited at each center through physician referral, advertisement in clinical areas or self-referral using the SPIROMICS study website (www.spiromics.com). The research protocol was approved by the institutional review boards of all participating institutions with written informed consent from all participants.
For quality control, subjects were removed if they 1) had genotyping call rates< 95%, 2) were discrepant for genetic sex, 3) failed the check for family relatedness, or 4) were detected as an outlier. After subjects meeting these criteria were excluded, SNPs were removed if 1) call rates< 95%, 2) inconsistent with Hardy-Weinberg Equilibrium (HWE) (p < 10− 6), or 3) minor allele frequency (MAF) < 0.01.
A linear additive model was used for analysis of pre−/post-bronchodilator FEV1/FVC, percent predicted FEV1, FVC, and % change in FEV1 bronchodilator response using PLINK software (URL: zzz.bwh.harvard.edu/plink/) , adjusted for age, sex, current smoking status, pack-years of cigarette smoking, and the first two principal components from the multidimensional scaling analysis of genotypes on the chip. Association analyses of Pre-/Post-bronchodilator FEV1 and FVC in ml were performed using linear regression adjusted for sex, age, age2, height, height2, weight, current smoking status, pack-years of cigarette smoking, and the first two principal components. Association analyses of COPD and COPD severity were performed using logistic regression adjusted for age, sex, current smoking status, pack-years of cigarette smoking, and the first two principal components. P values≤5 × 10− 8 were considered genome-wide significant. P values ≤0.05 were considered significant for SNP-level evaluation of previously reported candidate SNPs associated with baseline lung function. SNAP software (URL: http://www.broad.mit.edu/mpg/snap/) was used to generate the association plots .
Joint analysis of 10 confirmed candidate SNPs was performed, in which eight subjects with homozygous TT genotype of rs28929474 in SERPINA1 (PiZZ genotype) were not included in joint analysis to avoid bias. Genetic scores were defined by the number of risk alleles presented in these 10 SNPs. A linear model was used for analysis of post-bronchodilator FEV1/FVC and percent predicted FEV1 with genetic scores in 1632 current or former smokers. Joint analysis of these 10 candidate SNPs was also performed for post-bronchodilator percent predicted FEV1 and percentage of subjects with severe COPD (GOLD stage 3–4) in 1077 smokers with COPD.
GWAS of post-bronchodilator pulmonary function
After quality control analysis, 1645 non-Hispanic White subjects (1086 subjects with COPD and 559 current and former smokers with preserved lung function ) were included in the analysis (Table 1). GWAS of post-bronchodilator FEV1/FVC and percent predicted FEV1 were performed for 635,970 single nucleotide polymorphisms (SNPs) with MAF ≥ 0.01 in 1645 non-Hispanic White smokers with age, sex, current smoking status, pack-years of cigarette smoking, and the first two principal components as covariates in the linear additive model. Genomic inflation factors are 1.013 and 1.017 for GWAS of post-bronchodilator FEV1/FVC and percent predicted FEV1, respectively, indicating limited genomic inflation.
SNPs in nine genes previously identified for baseline FEV1/FVC or FEV1 in general populations [2,3,4], extremes of high and low baseline FEV1  or COPD  were also associated (p < 0.05) with post-bronchodilator FEV1/FVC or FEV1 (Table 2). SNPs in RARB, HDAC4, CHRNA3, and RIN3 were associated with post-bronchodilator FEV1/FVC and FEV1. SNPs in HHIP, AGER, FAM13A1, and PID1 were only associated with post-bronchodilator FEV1/FVC. A SNP in MMP12 was only associated with post-bronchodilator FEV1. The associations were significant at the SNP level with same effect direction as previous findings [2,3,4, 6].
rs28929474 (Glu342Lys) in alpha-1 antitrypsin member 1 (SERPINA1) was associated with post-bronchodilator FEV1/FVC (β = − 0.087, p = 1.2 × 10− 8) and percent predicted FEV1 (β = − 13.6, p = 3.5 × 10− 8) at a genome-wide significant level (Table 2 and Additional file 1: Tables S1-S2). No other SNPs in the SERPINA1 region were in strong linkage disequilibrium (LD) with rs28929474 or strongly associated with post-bronchodilator lung function (Figs. 1 and 2). rs4537555 in hedgehog acyltransferase (HHAT) and rs8079868 in myosin heavy chain 3 (MYH3) were strongly associated with post-bronchodilator FEV1/FVC (p = 2.1 × 10− 7) and percent predicted FEV1 (p = 5.9 × 10− 7), respectively (Table 2 and Additional file 1: Tables S1-S2).
Association of SERPINA1 with lung function, COPD, and COPD severity
Pre−/post-bronchodilator lung function was stratified by genotypes of rs28929474 (Table 3). rs28929474 was associated in a stepwise fashion with pre−/post-bronchodilator FEV1/FVC ratio (0.39, 0.54, and 0.61 for genotype TT, CT, and CC, respectively; p = 1.2 × 10− 8). rs28929474 was associated with pre−/post-bronchodilator FEV1 (33.5, 61.3, and 72.5 or 1210, 1841, and 2115 ml for genotype TT, CT, and CC, respectively; p = 2.1 × 10− 9). Pre−/post-bronchodilator lung function was significantly different between CT and CC or TT and CC genotype groups, however differences between TT and CT genotype groups were not as marked. rs28929474 was associated with post-bronchodilator FEV1/FVC (β = − 0.060, p = 1.1 × 10− 5) and percent predicted FEV1 (β = − 8.73; p = 2.6 × 10− 4) in subjects with COPD (GOLD stage 1–4), but not in subjects without COPD (GOLD stage 0; data not shown). Thus, the association of rs28929474 with lung function was driven by subjects with COPD.
Additional COPD-related phenotypes were analyzed for association with rs28929474 (Table 3). rs28929474 was also associated with COPD status (odds ratio = 2.3, p = 7.8 × 10− 4) and COPD severity (odds ratio = 4.1, p = 0.0036) (Table 3). The percentage of subjects with COPD or severe COPD was significantly higher in subjects with CT genotype than CC genotype. rs28929474 was a less common SNP with minor allele frequency (MAF) of 0.029 in SPIROMICS (Additional file 1: Table S3). Homozygous risk genotype TT was present only in subjects (n = 8) with severe COPD (GOLD stage 3–4).
Prediction of post-bronchodilator pulmonary function
Joint analysis of the most consistently associated 10 SNPs, based on our analyses and previous findings was performed. Genetic scores (the number of risk alleles) and pack-years of cigarette smoking were significantly associated with post-bronchodilator FEV1/FVC and percent predicted FEV1 (Table 4). Age at enrollment and sex were significantly associated with post-bronchodilator FEV1/FVC but not associated with percent predicted FEV1. In 1632 SPIROMICS non-Hispanic White smokers (GOLD stage 0–4), genetic score, age, sex, and pack-years of cigarette smoking explained 3.6, 1.5, 1.9, 3.0%, and together 8.6% of the variance of post-bronchodilator FEV1/FVC (Table 4). Genetic score and pack-years of cigarette smoking explained 3.0, 2.9%, and together 5.8% of the variance of post-bronchodilator percent predicted FEV1 (Table 4). In 1077 SPIROMICS non-Hispanic White smokers with COPD (GOLD stage 1–4), post-bronchodilator FEV1 decreased significantly with the increase in the number of risk alleles, from 65.4 to 54.0 (p = 1.2 × 10− 5) and the percentage of subjects with severe COPD (GOLD stage 3–4) increased significantly from 25.6 to 48.3% (p = 5.5 × 10− 5) (Fig. 3).
Joint analysis of the top 10 SNPs associated with post-bronchodilator % predicted FEV1 in this study was also performed (Additional file 1: Table S2). In 1634 SPIROMICS non-Hispanic White smokers (GOLD stage 0–4), genetic score, age, sex, and pack-years of cigarette smoking explained 7.5, 1.4, 1.9, 3.1%, and together 12.8% of the variance of post-bronchodilator FEV1/FVC (Additional file 1: Table S4). Genetic score and pack-years of cigarette smoking explained 9.9, 3.0%, and together 12.9% of the variance of post-bronchodilator percent predicted FEV1 (Additional file 1: Table S4). Increase in the number of risk alleles from 4 to 6 to 11–13 was associated with significant decrease in post-bronchodilator FEV1 from 69.4 to 45.6 (p < 2.2 × 10− 16) and with a significant increase in the percentage of subjects with severe COPD (GOLD stage 3–4) from 21.4 to 57.9% (p = 2.2 × 10− 12) (Fig. S1).
Joint analysis of 10 SNPs with emphysema, clinical symptoms, and exacerbation
Joint analysis of 10 candidate SNPs was further performed on quantitative Computed Tomography (CT) evidence of emphysema (TLC % area < − 950 HU) and airtrapping (RV % area < − 856 HU), BODE index, COPD Assessment Test (CAT) score, St. George’s Respiratory Questionnaire (SGRQ) total score, 6-Minute Walk Distance (6MWD), and exacerbations requiring ED visit or hospitalization in last 12 month (Table 5). In general, with the increase of genetic scores, emphysema (p < 0.0001) and airtrapping (p < 0.0001) increased, BODE index (p < 0.0001) and SGRQ total score (p = 0.0044) increased, 6MWD (p = 0.0086) decreased, and the percentage of subjects with exacerbations (p = 0.001) increased. Two extreme genetic score groups (8 to 11 risk alleles vs. 16 to 18 risk alleles) showed statistical and clinical difference for emphysema (5.54 vs. 12.5 of TLC % area < − 950 HU), airtrapping (21.8 vs. 33.9 RV % area < − 856 HU), BODE index (1.15 vs. 2.21), SGRQ total score (30.4 vs. 35.5), 6MWD (416 m vs. 390 m), and percentage of exacerbations (7.5% vs. 14%).
In this study, we performed GWAS of post-bronchodilator FEV1/FVC and percent predicted FEV1, and identified rs28929474 in SERPINA1. In 1963, Laurell and Eriksson identified the connection between alpha 1-antitrypsin (A1AT) deficiency and degenerative pulmonary disease . The SERPINA1 gene on chr14q32 encodes A1AT protein. The most common variant of SERPINA1 causing A1AT deficiency is the Z allele (rs28929474: Glu342Lys), which is a missense mutation of glutamic acid to lysine at position 342 of A1AT protein. The homozygous TT genotype of rs28929474 (PiZZ genotype) is consistently associated with emphysema, decreased lung function, and COPD [12, 13].
Previous GWAS of COPD, emphysema, and lung function did not identify rs28929474 in SERPINA1 [2,3,4,5,6, 14]. There are several potential reasons for missing this association. rs28929474 is relatively rare in the general population, for example, approximately 2 and 0.01% of the population in the United States are heterozygous or homozygous for the T allele, respectively . The largest meta-analyses of GWAS of baseline lung function in general populations of European descent [2,3,4,5] have included tens of thousands subjects, however very few subjects may have been homozygous for the T allele and more importantly these studies did not ascertain subjects with a significant history of cigarette smoking, a necessary environmental exposure. Thus, these studies in general populations have limited power to identify the association between rs28929474 and lung function. In this study, we performed GWAS of post-bronchodilator lung function in heavy smokers enriched for COPD. As expected the number of subjects with homozygous TT genotype was rare (n = 8 in 1645 or 0.49%) but the heterozygous CT genotype was more common (n = 78 or 4.74%). In addition, rs28929474 is not included in the previously designed GWAS chips nor are there other SNPs in strong LD (r2 > 0.5) with rs28929474, preventing the identification of association with COPD and emphysema [6, 14]. The Illumina OmniExpressExome BeadChip used in this study includes exonic markers identified from exome and whole genome sequencing projects. rs28919474 (exm1124179) was directly genotyped. This study found rs28929474 in SERPINA1 to be associated with pre- and post-bronchodilator FEV1/FVC and FEV1 at a genome-wide significant level (Table 3).
Although the function of homozygous TT has been known for a long while, the effect of heterozygous CT is more controversial and has been questioned in candidate-gene studies in the past [16,17,18]. For example, in a general population (n = 4600), baseline FEV1/FVC and FEV1 were not significantly different between PiMM and PiMZ . In a case-control study (834 COPD cases and 835 controls), post-bronchodilator FEV1/FVC and FEV1 were not significantly different between PiMM and PiMZ . In a small study composed of mainly healthy subjects, post-bronchodilator FEV1/FVC (0.77 or 0.71 for PiMM or PiMZ) and percent predicted FEV1 (96.4 or 84.6 for PiMM or PiMZ) were significantly different in ever-smokers but not in never-smokers . In a recent candidate-gene study (5518 non-Hispanic Whites and 2753 African Americans with ≥10 pack-years of smoking), subjects with PiMZ had significant lower lung function than subjects with PiMM in both Whites and African Americans . In the current study, subjects with CT genotype had intermediate values for lung function between subjects with TT and CC genotype (Table 3). Subjects with CT genotypes had significantly lower post-bronchodilator FEV1/FVC and percent predicted FEV1, and higher percentage of COPD and more severe COPD than subjects with CC genotype. Thus, SERPINA1 CT heterozygosity has important functional effects on COPD and lung function. All subjects included in our study had a history of tobacco smoking with at least a 20-pack-years. Association results were unaffected by the number of pack-years of cigarette smoking in our study. Compared with results from COPDGene study , this study included heavier smokers, and thus had lower lung function. More importantly, this study is a hypothesis-free GWAS study, which identified association of rs28929474 with lung function at genome-wide significant level for the first time. More than a hundred common and rare variants exist in the SERPINA1 gene. Thun et al. have identified synthetic association between common variants in SERPINA1 and serum A1AT levels, suggesting A1AT levels are causally determined by rare variants such as Z allele and S allele (rs17580) . Cho et al. have identified rs45505795 in SERPINA10 with MAF of 0.04 (not in strong LD with rs28929474: r2 = 0.295) associated with emphysema . We found no SNP other than rs28929474 in SERPINA1 region to be strongly associated with lung function (Figs. 1 and 2).
To develop a multi-gene predictive model for lung function, genes associated with lung function and COPD in previous published studies were evaluated. We identified the association of HHIP, FAM13A1, AGER, PID1, HDAC4, RARB, CHRNA3, RIN3, and MMP12 with post-bronchodilator lung function at the SNP level (Table 2). In a previous study, we have showed that HHIP, FAM13A1, AGER and RARB associated with pre-bronchodilator lung function in subjects with asthma . The lung expression quantitative trait locus (eQTL) analysis has identified cis-eQTL SNPs in HHIP, FAM13A1, and AGER . All the evidence indicates rs1980057 in HHIP, rs2869967 in FAM13A1, and rs2070600 in AGER are functionally relevant SNPs important for lung function in the general population and in subjects with COPD or asthma. rs4537555 in HHAT was strongly associated with post-bronchodilator FEV1/FVC (Table 2). HHAT is a hedgehog acyltransferase which catalyzes N-terminal palmitoylation of sonic hedgehog (SHH). Hedgehog interacting protein (HHIP) and patched homolog 1 (PTCH1) are the other two genes involved in hedgehog signaling pathway and associated with lung function [2,3,4, 21], indicating the importance of this pathway in lung development and function. Independent replication and functional study of HHAT are warranted.
Since each of these variants alone had smaller effects, we performed a joint analysis of 10 confirmed candidate SNPs. This analysis explained 3.63 and 2.96% variance of post-bronchodilator FEV1/FVC and percent predicted FEV1, respectively (Table 4). In contrast, pack-years of cigarette smoking explained 3.04 and 2.94% variance of post-bronchodilator FEV1/FVC and percent predicted FEV1. A genetic score using these 10 candidate SNPs, age, sex, and pack-years of cigarette smoking together explained 8.59 and 5.83% variance of post-bronchodilator FEV1/FVC and percent predicted FEV1. In addition, joint analysis of 10 confirmed candidate SNPs (with Z allele homozygotes removed) was performed on CT evidence of emphysema and airtrapping, BODE index, COPD, CAT score, SGRQ total score, 6MWD, and exacerbations (Table 5) in all heavy smokers (Gold stage 0–4). Statistical and clinical significant difference was shown between two extreme genetic score groups (8–11 vs. 16–18) for emphysema, airtrapping, BODE index, SGRQ total score, 6MWD, and exacerbations, indicating the potential usefulness of genetic information to distinguish clinical subgroups of heavy smokers. It will be important to evaluate the power of this model to predict decline in lung function and progression of COPD severity longitudinally in clinical settings.
In summary, rs28929474 in SERPINA1 is clearly associated with post-bronchodilator FEV1/FVC and FEV1 among heavy smokers. This study is the first to show genome-wide significant association of rs28929474 with lung function. In addition, rs28929474 is associated with COPD and COPD severity. While well-established rare ZZ homozygotes have severe COPD and emphysema, this study establishes that more common heterozygotes (4.7% of subjects) at this locus lead to pulmonary abnormality in smokers and COPD. Thus, in future clinical studies, this largely ignored heterozygotes group should be carefully examined. A joint genetic model combined with environmental factors is associated with reduced lung function, emphysema, exacerbation, and clinical symptoms. The models should be tested in other populations as well as longitudinally to evaluate potential value of predicting COPD progression and severity.
Chronic obstructive pulmonary disease
- FEV1 :
Forced expiratory volume in 1 s
Genome-wide association studies
The SubPopulations and InteRmediate Outcome Measures In COPD Study
Vestbo J, Hurd SS, Agustí AG, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. 2013;187:347–65.
Soler Artigas M, Loth DW, Wain LV, et al. Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat Genet. 2011;43:1082–90.
Hancock DB, Eijgelsheim M, Wilk JB, et al. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet. 2010;42:45–52.
Repapi E, Sayers I, Wain LV, et al. Genome-wide association study identifies five loci associated with lung function. Nat Genet. 2010;42:36–44.
Wain LV, Shrine N, Miller S, et al. Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK biobank. Lancet Respir Med. 2015;3:769–81.
Cho MH, McDonald ML, Zhou X, et al. Risk loci for chronic obstructive pulmonary disease: a genome-wide association study and meta-analysis. Lancet Respir Med. 2014;2:214–25.
Couper D, LaVange LM, Han M, et al. Design of the Subpopulations and Intermediate Outcomes in COPD study (SPIROMICS). Thorax. 2014;69:491–4.
Woodruff PG, Barr RG, Bleecker E, et al. Clinical significance of symptoms in smokers with preserved pulmonary function. N Engl J Med. 2016;374:1811–21.
Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Johnson AD, Handsaker RE, Pulit SL, et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–9.
Laurell CB, Eriksson S. The electrophoretic alpha 1 globulin pattern of serum in alpha 1 antitrypsin deficiency. Scand J Clin Lab Invest. 1963;15:132–40.
DeMeo DL, Silverman EK. Alpha1-antitrypsin deficiency. 2: genetic aspects of alpha(1)-antitrypsin deficiency: phenotypes and genetic modifiers of emphysema risk. Thorax. 2004;59:259–64.
Brebner JA, Stockley RA. Recent advances in α-1-antitrypsin deficiency-related lung disease. Expert Rev Respir Med. 2013;7:213–29.
Cho MH, Castaldi PJ, Hersh CP, et al. A genome-wide association study of emphysema and airway quantitative imaging phenotypes. Am J Respir Crit Care Med. 2015;192:559–69.
de Serres FJ, Blanco I. Prevalence of α1-antitrypsin deficiency alleles PI*S and PI*Z worldwide and effective screening for each of the five phenotypic classes PI*MS, PI*MZ, PI*SS, PI*SZ, and PI*ZZ: a comprehensive review. Ther Adv Respir Dis. 2012;6:277–95.
Sørheim IC, Bakke P, Gulsvik A, et al. α1-Antitrypsin protease inhibitor MZ heterozygosity is associated with airflow obstruction in two large cohorts. Chest. 2010;138:1125–32.
Thun GA, Ferrarotti I, Imboden M, et al. SERPINA1 PiZ and PiS heterozygotes and lung function decline in the SAPALDIA cohort. PLoS One. 2012;7:e42728.
Molloy K, Hersh CP, Morris VB, et al. Clarification of the risk of chronic obstructive pulmonary disease in α1-antitrypsin deficiency PiMZ heterozygotes. Am J Respir Crit Care Med. 2014;189:419–27.
Foreman MG, Wilson C, DeMeo DL, et al. Alpha-1 antitrypsin PI MZ genotype is associated with COPD in two racial groups. Ann Am Thorac Soc 2017. https://doi.org/10.1513/AnnalsATS.201611-838OC. Epub ahead of print.
Thun GA, Imboden M, Ferrarotti I, et al. Causal and synthetic associations of variants in the SERPINA gene cluster with alpha1-antitrypsin serum levels. PLoS Genet. 2013;9:e1003585.
Li X, Hawkins GA, Ampleford EJ, et al. Genome-wide association study identifies TH1 pathway genes associated with lung function in asthmatic patients. J Allergy Clin Immunol. 2013;132:313–20.
Obeidat M, Hao K, Bossé Y, et al. Molecular mechanisms underlying variations in lung function: a systems genetics analysis. Lancet Respir Med. 2015;3:782–95.
The authors thank the SPIROMICS participants and participating physicians, investigators and staff for making this research possible. More information about the study and how to access SPIROMICS data is at www.spiromics.org. We would like to acknowledge the following current and former investigators of the SPIROMICS sites and reading centers: Neil E Alexis, PhD; Wayne H Anderson, PhD; R Graham Barr, MD, DrPH; Eugene R Bleecker, MD; Richard C Boucher, MD; Russell P Bowler, MD, PhD; Elizabeth E Carretta, MPH; Stephanie A Christenson, MD; Alejandro P Comellas, MD; Christopher B Cooper, MD, PhD; David J Couper, PhD; Gerard J Criner, MD; Ronald G Crystal, MD; Jeffrey L Curtis, MD; Claire M Doerschuk, MD; Mark T Dransfield, MD; Christine M Freeman, PhD; MeiLan K Han, MD, MS; Nadia N Hansel, MD, MPH; Annette T Hastie, PhD; Eric A Hoffman, PhD; Robert J Kaner, MD; Richard E Kanner, MD; Eric C Kleerup, MD; Jerry A Krishnan, MD, PhD; Lisa M LaVange, PhD; Stephen C Lazarus, MD; Fernando J Martinez, MD, MS; Deborah A Meyers, PhD; John D Newell Jr., MD; Elizabeth C Oelsner, MD, MPH; Wanda K O’Neal, PhD; Robert Paine, III, MD; Nirupama Putcha, MD, MHS; Stephen I. Rennard, MD; Donald P Tashkin, MD; Mary Beth Scholand, MD; J Michael Wells, MD; Robert A Wise, MD; and Prescott G Woodruff, MD, MPH. The project officers from the Lung Division of the National Heart, Lung, and Blood Institute were Lisa Postow, PhD, and Thomas Croxton, PhD, MD.
SPIROMICS was supported by contracts from the NIH/NHLBI (HHSN268200900013C, HHSN268200900014C, HHSN268200900015C, HHSN268200900016C, HHSN268200900017C, HHSN268200900018C, HHSN268200900019C, HHSN268200900020C), which were supplemented by contributions made through the Foundation for the NIH from AstraZeneca; Bellerophon Pharmaceuticals; Boehringer-Ingelheim Pharmaceuticals, Inc.; Chiesi Farmaceutici SpA; Forest Research Institute, Inc.; GSK; Grifols Therapeutics, Inc.; Ikaria, Inc.; Nycomed GmbH; Takeda Pharmaceutical Company; Novartis Pharmaceuticals Corporation; Regeneron Pharmaceuticals, Inc.; and Sanofi. The funders had no role in the study design, data collection, data analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files. Raw genotype data may be accessible by contacting SPIROMICS (https://www.spiromics.org).
Ethics approval and consent to participate
Participants were recruited at each center through physician referral, advertisement in clinical areas or self-referral using the SPIROMICS study website (www.spiromics.com). The research protocol was approved by the institutional review boards of all participating institutions (Wake Forest School of Medicine, Columbia University, University of California at San Francisco, University of California at Los Angeles, University of North Carolina at Chapel Hill, University of Alabama at Birmingham, University of Michigan, Johns Hopkins University School of Medicine, University of Iowa, University of Utah, Weill Cornell Medical College of Cornell University) with written informed consent from all participants.
Consent for publication
X.L.: Associate Editor of BMC Medical Genetics.
V.E.O.: funding from the Foundation for the NIH NHLBI in the form of a K08 training award; consultancy fees from CSL Behring.
E.J.A.: no conflicts of interest to disclose.
R.G.B.: grants from the Foundation for the NIH, Alpha1 Foundation and personal fees from UpToDate and the COPD Foundation all outside of the submitted work.
S.A.C.: no conflicts of interest to disclose.
C.B.C.: grants from the Foundation for the NIH and NIH NHLBI; part-time employment by the Global Respiratory Franchise in the GlaxoSmithKline.
D.C.: grants from the Foundation for the NIH and NIH NHLBI.
M.T.D.: grants from the NIH, the Department of Defense, and the American Heart.
Association; consultancy fees from Boehringer Ingelheim, Boston Scientific, and GlaxoSmithKline and contracted clinical trials from Boehringer Ingelheim, Boston Scientific, GlaxoSmithKline, Pearl, Pulmonx, PneumRx, AstraZeneca, Novartis, and Yungjin.
M.H.: grants from the NIH NHLBI and the Foundation for the NIH; consultancy fees from GlaxoSmithKline, Boehringer-Ingelheim, Novartis, and AstraZeneca.
N.N.H.: grants from the Foundation for the NIH and NIH NHLBI.
E.A.H.: grants from the Foundation for the NIH and NIH NHLBI; founder and shareholder of VIDA Diagnostics.
R.E.K.: grants from the Foundation for the NIH and NIH NHLBI.
E.K.: grants from the Foundation for the NIH and NIH NHLBI; grants from Boehringer-Ingelheim, Novartis, Pearl, AstraZeneca, and Sunovion outside of the submitted work.
F.J.M.: grants from National Institutes of Health, Clarion, Continuing Education, Potomac, Afferent, and Adept; personal fees from Forest, Janssen, GlaxoSmithKline, Nycomed/Takeda, Amgen, Astra Zeneca, Boehringer-Ingelheim, Ikaria/Bellerophon, Genentech, Janssen, Johnson & Johnson, Novartis, Pearl, Pfizer, Roche, Sunovion, Theravance, Axon Communication, CME Incite, California Society for Allergy and Immunology, Annenberg, Integritas, InThought, Miller Medical, National Association for Continuing Education, Paradigm, Peer Voice, UpToDate, Haymarket Communications, Western Society of Allergy and Immunology, Bioscale, Unity Biotechnology, ConCert, Lucid, Methodist Hospital, Prime, WebMD, Mereo, Kadmon, Pfizer, Veracyte, American Thoracic Society, Academic CME, Falco, and the National Association for Continuing Education.
R.P.: grants from the Foundation for the NIH and NIH NHLBI.
P.G.W.: grants from Medimmune and consultancy fees from Genentech/Roche, Astra Zeneca, Novartis, Neostem, Janssen outside the submitted work; a patent with Asthma diagnostics pending.
G.A.H.: no conflicts of interest to disclose.
E.R.B.: grants from the NIH NHLBI for the Severe Asthma Research Program, AsthmaNet, SPIROMICS, and the Foundation for the NIH; consultancy fees from Amgen, AstraZeneca-MedImmune, Boehringer-Ingelheim, Genentech/Roche, GlaxoSmithKline, Knopp, Novartis, and Sanofi/Regeneron; funds for clinical trials administered through the Wake Forest School of Medicine from Amgen, AstraZeneca-MedImmune, Boehringer-Ingelheim, Genentech/Roche, GlaxoSmithKline, Janssen/Johnson & Johnson, Novartis, Pzifer, Sanofi-Regeneron, and Teva.
D.A.M.: no conflicts of interest to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Association Results of the Top SNPs (P < 10− 4) with Post-bronchodilator FEV1/FVC. Table S2. Association Results of the Top SNPs (P < 10− 4) with Post-bronchodilator % Predicted FEV1.Table S3. Genotype Frequency of rs28929474 in SERPINA1 Stratified by GOLD Stages. Table S4. Prediction Models for Post-bronchodilator Lung Function Using Top 10 SNPs for Post-bronchodilator % Predicted FEV1.Figure S1. Joint analysis of the top10 SNPs for post-bronchodilator % predicted FEV1 in 1075 SPIROMICS non-Hispanic White smokers with COPD. (DOCX 141 kb)
About this article
Cite this article
Li, X., Ortega, V.E., Ampleford, E.J. et al. Genome-wide association study of lung function and clinical implication in heavy smokers. BMC Med Genet 19, 134 (2018). https://doi.org/10.1186/s12881-018-0656-z
- Lung function