Chronic obstructive pulmonary disease (COPD) is linked to both cigarette smoking and genetic determinants. We have previously identified iron-responsive element binding protein 2 (IRP2) as an important COPD susceptibility gene, with IRP2 protein increased in the lungs of individuals with COPD. Here we demonstrate that mice deficient in Irp2 were protected from cigarette smoke (CS)-induced experimental COPD. By integrating RIP-Seq, RNA-Seq, gene expression and functional enrichment clustering analysis, we identified IRP2 as a regulator of mitochondrial function in the lung. IRP2 increased mitochondrial iron loading and cytochrome c oxidase (COX), which led to mitochondrial dysfunction and subsequent experimental COPD. Frataxin-deficient mice with higher mitochondrial iron loading had impaired airway mucociliary clearance (MCC) and higher pulmonary inflammation at baseline, whereas synthesis of cytochrome c oxidase (Sco2)-deficient mice with reduced COX were protected from CS-induced pulmonary inflammation and impairment of MCC. Mice treated with a mitochondrial iron chelator or mice fed a low-iron diet were protected from CS-induced COPD. Mitochondrial iron chelation also alleviated CS-impairment of MCC, CS-induced pulmonary inflammation and CS-associated lung injury in mice with established COPD, suggesting a critical functional role and potential therapeutic intervention for the mitochondrial-iron axis in COPD.
COPD; Mitochondria; Iron; IRP2; Lung
Klinefelter syndrome (KS) (47 XXY) is a common sex-chromosome aneuploidy with an estimated prevalence of 1 in every 660 male births. Investigations into the associations between DNA methylation and the highly variable clinical manifestations of KS have largely focused on the supernumerary X chromosome; systematic investigations of the epigenome have been limited. We obtained genome-wide DNA methylation data from peripheral blood using the Illumina HumanMethylation450K platform in 5 KS (47 XXY), 102 male (46 XY), and 113 female (46 XX) control subjects participating in the chronic obstructive pulmonary disease (COPD) Gene Study. Empirical Bayes-mediated models were used to test for differential methylation by KS status. CpG sites with a false-discovery rate <0.05 from the first-generation HumanMethylation27K platform were further examined in an independent replication cohort of 2 KS subjects, 590 male, and 495 female controls drawn from the International COPD Genetics Network (ICGN). Differential methylation at sites throughout the genome were identified, including 86 CpG sites that were differentially methylated in KS subjects relative to both male and female controls. CpG sites annotated to the HEN1 methyltransferase homolog 1 (HENMT1), calcyclin-binding protein (CACYBP), and GTPase-activating protein (SH3 domain)-binding protein 1 (G3BP1) genes were among the “KS-specific” loci that were replicated in ICGN. We therefore conclude that site-specific differential methylation exists throughout the genome in KS. The functional impact and clinical relevance of these differentially methylated loci should be explored in future studies.
[MeSH]: Klinefelter syndrome; DNA methylation; epigenomics; XXY syndrome
Rationale: Former smoking history and chronic obstructive pulmonary disease (COPD) are potential risk factors for osteoporosis and fractures. Under existing guidelines for osteoporosis screening, women are included but men are not, and only current smoking is considered.
Objectives: To demonstrate the impact of COPD and smoking history on the risk of osteoporosis and vertebral fracture in men and women.
Methods: Characteristics of participants with low volumetric bone mineral density (vBMD) were identified and related to COPD and other risk factors. We tested associations of sex and COPD with both vBMD and fractures adjusting for age, race, body mass index (BMI), smoking, and glucocorticoid use.
Measurements and Main Results: vBMD by calibrated quantitative computed tomography (QCT), visually scored vertebral fractures, and severity of lung disease were determined from chest CT scans of 3,321 current and ex-smokers in the COPDGene study. Low vBMD as a surrogate for osteoporosis was calculated from young adult normal values. Male smokers had a small but significantly greater risk of low vBMD (2.5 SD below young adult mean by calibrated QCT) and more fractures than female smokers. Low vBMD was present in 58% of all subjects, was more frequent in those with worse COPD, and rose to 84% among subjects with very severe COPD. Vertebral fractures were present in 37% of all subjects and were associated with lower vBMD at each Global Initiative for Chronic Obstructive Lung Disease stage of severity. Vertebral fractures were most common in the midthoracic region. COPD and especially emphysema were associated with both low vBMD and vertebral fractures after adjustment for steroid use, age, pack-years of smoking, current smoking, and exacerbations. Airway disease was associated with higher bone density after adjustment for other variables. Calibrated QCT identified more subjects with abnormal values than the standard dual-energy X-ray absorptiometry in a subset of subjects and correlated well with prevalent fractures.
Conclusions: Male smokers, with or without COPD, have a significant risk of low vBMD and vertebral fractures. COPD was associated with low vBMD after adjusting for race, sex, BMI, smoking, steroid use, exacerbations, and age. Screening for low vBMD by using QCT in men and women who are smokers will increase opportunities to identify and treat osteoporosis in this at-risk population.
low bone density; COPD; vertebral fractures; quantitative computed tomography; smoking
Airflow obstruction on spirometry is universally used to define chronic obstructive pulmonary disease (COPD), and current or former smokers without airflow obstruction may assume that they are disease free.
To identify clinical and radiologic evidence of smoking-related disease in a cohort of current and former smokers who did not meet spirometric criteria for COPD, for whom we adopted the discarded label of Global Initiative for Obstructive Lung Disease (GOLD) 0.
DESIGN, SETTING, AND PARTICIPANTS
Individuals from the Genetic Epidemiology of COPD (COPDGene) cross-sectional observational study completed spirometry, chest computed tomography (CT) scans, a 6-minute walk, and questionnaires. Participants were recruited from local communities at 21 sites across the United States. The GOLD 0 group (n = 4388) (ratio of forced expiratory volume in the first second of expiration [FEV1] to forced vital capacity >0.7 and FEV1 ≥80% predicted) from the COPDGene study was compared with a GOLD 1 group (n = 794), COPD groups (n = 3690), and a group of never smokers (n = 108). Recruitment began in January 2008 and ended in July 2011.
MAIN OUTCOMES AND MEASURES
Physical function impairments, respiratory symptoms, CT abnormalities, use of respiratory medications, and reduced respiratory-specific quality of life.
One or more respiratory-related impairments were found in 54.1% (2375 of 4388) of the GOLD 0 group. The GOLD 0 group had worse quality of life (mean [SD] St George’s Respiratory Questionnaire total score, 17.0 [18.0] vs 3.8 [6.8] for the never smokers; P < .001) and a lower 6-minute walk distance, and 42.3% (127 of 300) of the GOLD 0 group had CT evidence of emphysema or airway thickening. The FEV1 percent predicted distribution and mean for the GOLD 0 group were lower but still within the normal range for the population. Current smoking was associated with more respiratory symptoms, but former smokers had greater emphysema and gas trapping. Advancing age was associated with smoking cessation and with more CT findings of disease. Individuals with respiratory impairments were more likely to use respiratory medications, and the use of these medications was associated with worse disease.
CONCLUSIONS AND RELEVANCE
Lung disease and impairments were common in smokers without spirometric COPD. Based on these results, we project that there are 35 million current and former smokers older than 55 years in the United States who may have unrecognized disease or impairment. The effect of chronic smoking on the lungs and the individual is substantially underestimated when using spirometry alone.
Pulmonary function decline is a major contributor to morbidity and mortality among smokers. Post bronchodilator FEV1 and FEV1/FVC ratio are considered the standard assessment of airflow obstruction. We performed a genome-wide association study (GWAS) in 9919 current and former smokers in the COPDGene study (6659 non-Hispanic Whites [NHW] and 3260 African Americans [AA]) to identify associations with spirometric measures (post-bronchodilator FEV1 and FEV1/FVC). We also conducted meta-analysis of FEV1 and FEV1/FVC GWAS in the COPDGene, ECLIPSE, and GenKOLS cohorts (total n = 13,532).
Among NHW in the COPDGene cohort, both measures of pulmonary function were significantly associated with SNPs at the 15q25 locus [containing CHRNA3/5, AGPHD1, IREB2, CHRNB4] (lowest p-value = 2.17 × 10−11), and FEV1/FVC was associated with a genomic region on chromosome 4 [upstream of HHIP] (lowest p-value = 5.94 × 10−10); both regions have been previously associated with COPD. For the meta-analysis, in addition to confirming associations to the regions near CHRNA3/5 and HHIP, genome-wide significant associations were identified for FEV1 on chromosome 1 [TGFB2] (p-value = 8.99 × 10−9), 9 [DBH] (p-value = 9.69 × 10−9) and 19 [CYP2A6/7] (p-value = 3.49 × 10−8) and for FEV1/FVC on chromosome 1 [TGFB2] (p-value = 8.99 × 10−9), 4 [FAM13A] (p-value = 3.88 × 10−12), 11 [MMP3/12] (p-value = 3.29 × 10−10) and 14 [RIN3] (p-value = 5.64 × 10−9).
In a large genome-wide association study of lung function in smokers, we found genome-wide significant associations at several previously described loci with lung function or COPD. We additionally identified a novel genome-wide significant locus with FEV1 on chromosome 9 [DBH] in a meta-analysis of three study populations.
Electronic supplementary material
The online version of this article (doi:10.1186/s12863-015-0299-4) contains supplementary material, which is available to authorized users.
Chronic obstructive pulmonary disease; DBH; FEV1; FEV1/FVC; Genome-wide association study; Spirometry
Rationale: The characterization of young adults who develop late-onset diseases may augment the detection of novel genes and promote new pathogenic insights.
Methods: We analyzed data from 2,500 individuals of African and European ancestry in the COPDGene Study. Subjects with severe, early-onset chronic obstructive pulmonary disease (COPD) (n = 70, age < 55 yr, FEV1 < 50% predicted) were compared with older subjects with COPD (n = 306, age > 64 yr, FEV1 < 50% predicted).
Measurements and Main Results: Subjects with severe, early-onset COPD were predominantly females (66%), P = 0.0004. Proportionally, early-onset COPD was seen in 42% (25 of 59) of African Americans versus 14% (45 of 317) of non-Hispanic whites, P < 0.0001. Other risk factors included current smoking (56 vs. 17%, P < 0.0001) and self-report of asthma (39 vs. 25%, P = 0.008). Maternal smoking (70 vs. 44%, P = 0.0001) and maternal COPD (23 vs. 12%, P = 0.03) were reported more commonly in subjects with early-onset COPD. Multivariable regression analysis found association with African American race, odds ratio (OR), 7.5 (95% confidence interval [CI], 2.3–24; P = 0.0007); maternal COPD, OR, 4.7 (95% CI, 1.3–17; P = 0.02); female sex, OR, 3.1 (95% CI, 1.1–8.7; P = 0.03); and each pack-year of smoking, OR, 0.98 (95% CI, 0.96–1.0; P = 0.03).
Conclusions: These observations support the hypothesis that severe, early-onset COPD is prevalent in females and is influenced by maternal factors. Future genetic studies should evaluate (1) gene-by-sex interactions to address sex-specific genetic contributions and (2) gene-by-race interactions.
chronic obstructive pulmonary disease; female; African Americans
Cachexia, whether assessed by body mass index (BMI) or fat-free mass index (FFMI), affects a significant proportion of patients with chronic obstructive pulmonary disease (COPD), and is an independent risk factor for increased mortality, increased emphysema, and more severe airflow obstruction. The variable development of cachexia among patients with COPD suggests a role for genetic susceptibility. The objective of the present study was to determine genetic susceptibility loci involved in the development of low BMI and FFMI in subjects with COPD. A genome-wide association study (GWAS) of BMI was conducted in three independent cohorts of European descent with Global Initiative for Chronic Obstructive Lung Disease stage II or higher COPD: Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-Points (ECLIPSE; n = 1,734); Norway-Bergen cohort (n = 851); and a subset of subjects from the National Emphysema Treatment Trial (NETT; n = 365). A genome-wide association of FFMI was conducted in two of the cohorts (ECLIPSE and Norway). In the combined analyses, a significant association was found between rs8050136, located in the first intron of the fat mass and obesity–associated (FTO) gene, and BMI (P = 4.97 × 10−7) and FFMI (P = 1.19 × 10−7). We replicated the association in a fourth, independent cohort consisting of 502 subjects with COPD from COPDGene (P = 6 × 10−3). Within the largest contributing cohort of our analysis, lung function, as assessed by forced expiratory volume at 1 second, varied significantly by FTO genotype. Our analysis suggests a potential role for the FTO locus in the determination of anthropomorphic measures associated with COPD.
chronic obstructive pulmonary disease genetics; chronic obstructive pulmonary disease epidemiology; chronic obstructive pulmonary disease metabolism; genome-wide association study
Chronic obstructive pulmonary disease (COPD) is characterized by marked phenotypic heterogeneity. Most previous studies have focused on COPD subjects with FEV1 < 80% predicted. We investigated the clinical and genetic heterogeneity in subjects with mild airflow limitation in spirometry grade 1 defined by the Global Initiative for chronic Obstructive Lung Disease (GOLD 1).
Data from current and former smokers participating in the COPDGene Study (NCT00608764) were analyzed. K-means clustering was performed to explore subtypes within 794 GOLD 1 subjects. For all subjects with GOLD 1 and with each cluster, a genome-wide association study and candidate gene testing were performed using smokers with normal lung function as a control group. Combinations of COPD genome-wide significant single nucleotide polymorphisms (SNPs) were tested for association with FEV1 (% predicted) in GOLD 1 and in a combined group of GOLD1 and smoking control subjects.
K-means clustering of GOLD 1 subjects identified putative “near-normal”, “airway-predominant”, “emphysema-predominant” and “lowest FEV1 % predicted” subtypes. In non-Hispanic whites, the only SNP nominally associated with GOLD 1 status relative to smoking controls was rs7671167 (FAM13A) in logistic regression models with adjustment for age, sex, pack-years of smoking, and genetic ancestry. The emphysema-predominant GOLD 1 cluster was nominally associated with rs7671167 (FAM13A) and rs161976 (BICD1). The lowest FEV1 % predicted cluster was nominally associated with rs1980057 (HHIP) and rs1051730 (CHRNA3). Combinations of COPD genome-wide significant SNPs were associated with FEV1 (% predicted) in a combined group of GOLD 1 and smoking control subjects.
Our results indicate that GOLD 1 subjects show substantial clinical heterogeneity, which is at least partially related to genetic heterogeneity.
pulmonary disease; chronic obstructive; population characteristics; cluster analysis; genetic association
Severe α1-antitrypsin (AAT) deficiency is a proven genetic risk factor for chronic obstructive pulmonary disease (COPD), especially in individuals who smoke. There is marked variability in the development of lung disease in individuals homozygous (PI ZZ) for this autosomal recessive condition, suggesting that modifier genes could be important. We hypothesized that genetic determinants of obstructive lung disease may be modifiers of airflow obstruction in individuals with severe AAT deficiency. To identify modifier genes, we performed family-based association analyses for 10 genes previously associated with asthma and/or COPD, including IL10, TNF, GSTP1, NOS1, NOS3, SERPINA3, SERPINE2, SFTPB, TGFB1, and EPHX1. All analyses were performed in a cohort of 378 PI ZZ individuals from 167 families. Quantitative spirometric phenotypes included forced expiratory volume in one second (FEV1) and the ratio of FEV1/forced vital capacity (FVC). A qualitative phenotype of moderate-to-severe COPD was defined for individuals with FEV1 ⩽ 50 percent predicted. Six of 11 single-nucleotide polymorphisms (SNPs) in IL10 (P = 0.0005–0.05) and 3 of 5 SNPs in TNF (P = 0.01–0.05) were associated with FEV1 and/or FEV1/FVC. IL10 SNPs also demonstrated association with the qualitative COPD phenotype. When phenotypes of individuals with a physician's diagnosis of asthma were excluded, IL10 SNPs remained significantly associated, suggesting that the association with airflow obstruction was independent of an association with asthma. Haplotype analysis of IL10 SNPs suggested the strongest association with IL10 promoter SNPs. IL10 is likely an important modifier gene for the development of COPD in individuals with severe AAT deficiency.
chronic obstructive pulmonary disease; genetic modifiers; interleukin 10; family-based association analysis
Rationale: Computed tomography (CT) scanning of the lung may reduce phenotypic heterogeneity in defining subjects with chronic obstructive pulmonary disease (COPD), and allow identification of genetic determinants of emphysema severity and distribution.
Objectives: We sought to identify genes associated with CT scan distribution of emphysema in individuals without α1-antitrypsin deficiency but with severe COPD.
Methods: We evaluated baseline CT densitometry phenotypes in 282 individuals with emphysema enrolled in the Genetics Ancillary Study of the National Emphysema Treatment Trial, and used regression models to identify genetic variants associated with emphysema distribution.
Measurements and Main Results: Emphysema distribution was assessed by two methods—assessment by radiologists and by computerized density mask quantitation, using a threshold of −950 Hounsfield units. A total of 77 polymorphisms in 20 candidate genes were analyzed for association with distribution of emphysema. GSTP1, EPHX1, and MMP1 polymorphisms were associated with the densitometric, apical-predominant distribution of emphysema (p value range = 0.001–0.050). When an apical-predominant phenotype was defined by the radiologist scoring method, GSTP1 and EPHX1 single-nucleotide polymorphisms were found to be significantly associated. In a case–control analysis of COPD susceptibility limited to cases with densitometric upper-lobe–predominant cases, the EPHX1 His139Arg single-nucleotide polymorphism was associated with COPD (p = 0.005).
Conclusions: Apical and basal emphysematous destruction appears to be influenced by different genes. Polymorphisms in the xenobiotic enzymes, GSTP1 and EPHX1, are associated with apical-predominant emphysema. Altered detoxification of cigarette smoke metabolites may contribute to emphysema distribution, and these findings may lead to further insight into genetic determinants of emphysema.
COPD; genetics; association analysis; computed tomography; emphysema
The genetic risk factors for susceptibility to chronic obstructive
pulmonary disease (COPD) are still largely unknown. Additional genetic
variants are likely to be identified by genome-wide association studies in
larger cohorts or specific subgroups.
Genome-wide association analysis in COPDGene (non-Hispanic whites and
African-Americans) was combined with existing data from the ECLIPSE,
NETT/NAS, and GenKOLS (Norway) studies. Analyses were performed both using
all moderate-to-severe cases and the subset of severe cases. Top loci not
previously described as genome-wide significant were genotyped in the ICGN
study, and results combined in a joint meta-analysis.
Analysis of a total of 6,633 moderate-to-severe cases and 5,704
controls confirmed association at three known loci:
CHRNA3/CHRNA5/IREB2, FAM13A, and HHIP
(10−12 < P < 10−14),
and also showed significant evidence of association at a novel locus near
RIN3 (overall P, including ICGN =
5•4×10−9). In the severe COPD analysis
(n=3,497), the effects at two of three previously described loci were
significantly stronger; we also identified two additional loci previously
reported to affect gene expression of MMP12 and
TGFB2 (overall P = 2•6x10−9
and 8•3×10−9). RIN3 and
TGFB2 expression levels were reduced in a set of Lung
Tissue Research Consortium COPD lung tissue samples compared with
In a genome-wide study of COPD, we confirmed associations at three
known loci and found additional genome-wide significant associations with
moderate-to-severe COPD near RIN3 and with severe COPD near
MMP12 and TGFB2. Genetic variants,
apart from alpha-1 antitrypsin deficiency, increase the risk of COPD. Our
analysis of severe COPD suggests additional genetic variants may be
identified by focusing on this subgroup.
National Heart, Lung, and Blood Institute; the COPD Foundation
through contributions from AstraZeneca, Boehringer Ingelheim, Novartis, and
Sepracor; GlaxoSmithKline; Centers for Medicare and Medicaid Services;
Agency for Healthcare Research and Quality; US Department of Veterans
Genetic epidemiology studies of end-stage lung disease are potentially hindered by low numbers of participants due to early death of patients from the underlying disease, or due to exclusion from studies after patients have had lung transplants, because of concern about bias of genotype data due to chimerism. The number of participants enrolled in genetic studies of end-stage lung disease could be increased by including those individuals who have undergone lung transplant. We hypothesized that individuals who have had lung transplants can be included in genetic epidemiology studies that use single nucleotide polymorphism and short tandem repeat marker data, without confounding due to chimerism. Ten probands with severe, early-onset chronic obstructive pulmonary disease were included in this analysis. Pre– and post–lung transplant DNA samples were used in the investigation of concordance of genotype results for 12 short tandem repeat markers and 23 single nucleotide polymorphisms. Concordance was observed for all genotypes before and after lung transplant. We conclude that the risk of biasing genetic epidemiology studies due to donor lung–related DNA microchimerism is low, and that the inclusion of post–lung transplantation participants will allow for larger genetic epidemiology studies of individuals with end-stage lung disease.
genetic epidemiology; lung; chimerism; transplantation
Whether African Americans (AA) are more susceptible to COPD than non-Hispanic Whites (NHW) and whether racial differences in disease phenotype exist is controversial. The objective is to determine racial differences in the extent of emphysema and airway remodeling in COPD.
First, 2,500 subjects enrolled in the COPDGene study were used to evaluate racial differences in quantitative CT (QCT) parameters of % emphysema, air trapping and airway wall thickness. Independent variables studied included race, age, gender, education, BMI, pack-years, smoking status, age at smoking initiation, asthma, previous work in dusty job, CT scanner and center of recruitment.
Of the 1,063 subjects with GOLD Stage II-IV COPD, 200 self-reported as AA. AAs had a lower mean % emphysema (13.1 % vs. 16.1%, p = 0.005) than NHW and proportionately less emphysema in the lower lung zones. After adjustment for covariates, there was no statistical difference by race in air trapping or airway wall thickness. Measured QCT parameters were more predictive of poor functional status in NHWs compared to AAs.
AAs have less emphysema than NHWs but the same degree of airway disease. Additional factors not easily assessed by current QCT techniques may account for the poor functional status in AAs.
Airway wall thickness; Air trapping; Chronic obstructive pulmonary disease; Emphysema; Quantitative CT; Race
Comorbidities are common in COPD, but quantifying their burden is difficult. Currently there is a COPD-specific comorbidity index to predict mortality and another to predict general quality of life. We sought to develop and validate a COPD-specific comorbidity score that reflects comorbidity burden on patient-centered outcomes.
Materials and Methods
Using the COPDGene study (GOLD II-IV COPD), we developed comorbidity scores to describe patient-centered outcomes employing three techniques: 1) simple count, 2) weighted score, and 3) weighted score based upon statistical selection procedure. We tested associations, area under the Curve (AUC) and calibration statistics to validate scores internally with outcomes of respiratory disease-specific quality of life (St. George's Respiratory Questionnaire, SGRQ), six minute walk distance (6MWD), modified Medical Research Council (mMRC) dyspnea score and exacerbation risk, ultimately choosing one score for external validation in SPIROMICS.
Associations between comorbidities and all outcomes were comparable across the three scores. All scores added predictive ability to models including age, gender, race, current smoking status, pack-years smoked and FEV1 (p<0.001 for all comparisons). Area under the curve (AUC) was similar between all three scores across outcomes: SGRQ (range 0·7624–0·7676), MMRC (0·7590–0·7644), 6MWD (0·7531–0·7560) and exacerbation risk (0·6831–0·6919). Because of similar performance, the comorbidity count was used for external validation. In the SPIROMICS cohort, the comorbidity count performed well to predict SGRQ (AUC 0·7891), MMRC (AUC 0·7611), 6MWD (AUC 0·7086), and exacerbation risk (AUC 0·7341).
Quantifying comorbidity provides a more thorough understanding of the risk for patient-centered outcomes in COPD. A comorbidity count performs well to quantify comorbidity in a diverse population with COPD.
We propose a method for testing gene–environment (G × E) interactions on a complex trait in family-based studies in which a phenotypic ascertainment criterion has been imposed. This novel approach employs G-estimation, a semiparametric estimation technique from the causal inference literature, to avoid modeling of the association between the environmental exposure and the phenotype, to gain robustness against unmeasured confounding due to population substructure, and to acknowledge the ascertainment conditions. The proposed test allows for incomplete parental genotypes. It is compared by simulation studies to an analogous conditional likelihood–based approach and to the QBAT-I test, which also invokes the G-estimation principle but ignores ascertainment. We apply our approach to a study of chronic obstructive pulmonary disorder.
Causal inference; COPD; Family-based association; G-estimation; Gene–environment interaction
The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10−9). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV1 (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.
Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.
It is useful to have robust gene-environment interaction tests that can utilize a variety of family structures in an efficient way. This paper focuses on tests for gene-environment interaction in the presence of main genetic and environmental effects. The objective is to develop powerful tests that can combine trio data with parental genotypes and discordant sibships when parents genotypes are missing. We first make a modest improvement on a method for discordant sibs (discordant on phenotype), but the approach does not allow one to use families when all offspring are affected, e.g. trios. We then make a modest improvement on a Mendelian transmission-based approach that is inefficient when discordant sibs are available, but can be applied to any nuclear family. Finally, we propose a hybrid approach that utilizes the most efficient method for a specific family type, then combines over families. We utilize this hybrid approach to analyze a chronic obstructive pulmonary disorder dataset to test for gene-environment interaction in the Serpine2 gene with smoking. The methods are freely available in the R package fbati.
Gene-Environment Interaction; Family-Based Association Tests; Candidate Gene Analysis; Binary Trait; COPD; Serpine2
Acute exacerbations are a significant source of morbidity and mortality associated with chronic obstructive pulmonary disease. Among patients with COPD, some patients suffer an inordinate number of exacerbations while others remain relatively protected. We undertook a study to determine the clinical factors associated with "frequent exacerbator" status within a population of subjects with severe COPD.
Case-control cohort recruited from two Boston-area practices. All subjects had GOLD stage 3 or 4 (FEV1 ≤50% predicted) COPD. "Frequent exacerbators" (n=192) had an average of ≥2 moderate-to-severe exacerbations per year while "non-exacerbators" (n=153) had no exacerbations in the preceding 12 months. Multivariate logistic regression was performed to determine the significant clinical predictors of "frequent exacerbator" status.
Physician-diagnosed asthma was a significant predictor of frequent exacerbations. Within a subset of our cohort, the modified Medical Research Council dyspnea score and FEF25–75 % predicted were also significant clinical predictors of frequent exacerbator status (p<0.05). Differences in exacerbation frequency were not found to be due to increased current tobacco use or decreased rates of maintenance medication use.
Within our severe COPD cohort, a history of physician-diagnosed asthma was found to be a significant clinical predictor of frequent exacerbations. Although traditional risk factors such as decreased FEV1% predicted were not significantly associated with frequent exacerbator status, lower mid-expiratory flow rates, as assessed by FEF 25–75 % predicted, were significantly associated with frequent exacerbations in a subset of our cohort.
Chronic obstructive pulmonary disease (COPD) is characterized by alveolar destruction and abnormal inflammatory responses to noxious stimuli. Surfactant protein–D (SFTPD) is immunomodulatory and essential to host defense. We hypothesized that polymorphisms in SFTPD could influence the susceptibility to COPD. We genotyped six single-nucleotide polymorphisms (SNPs) in surfactant protein D in 389 patients with COPD in the National Emphysema Treatment Trial (NETT) and 472 smoking control subjects from the Normative Aging Study (NAS). Case-control association analysis was performed using Cochran–Armitage trend tests and multivariate logistic regression. The replication of significant associations was attempted in the Boston Early-Onset COPD Study, the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) Study, and the Bergen Cohort. We also correlated SFTPD genotypes with serum concentrations of surfactant protein–D (SP-D) in the ECLIPSE Study. In the NETT–NAS case-control analysis, four SFTPD SNPs were associated with susceptibility to COPD: rs2245121 (P = 0.01), rs911887 (P = 0.006), rs6413520 (P = 0.004), and rs721917 (P = 0.006). In the family-based analysis of the Boston Early-Onset COPD Study, rs911887 was associated with prebronchodilator and postbronchodilator FEV1 (P = 0.003 and P = 0.02, respectively). An intronic SNP in SFTPD, rs7078012, was associated with COPD in the ECLIPSE Study and the Bergen Cohort. Multiple SFTPD SNPs were associated with serum SP-D concentrations in the ECLIPSE Study. We demonstrated an association of polymorphisms in SFTPD with COPD in multiple populations. We demonstrated a correlation between SFTPD SNPs and SP-D protein concentrations. The SNPs associated with COPD and SP-D concentrations differed, suggesting distinct genetic influences on susceptibility to COPD and SP-D concentrations.
COPD; surfactant protein–D; single-nucleotide polymorphisms; genetics
Rationale: Several family-based studies have identified genetic linkage for lung function and airflow obstruction to chromosome 2q.
Objectives: We hypothesized that merging results of high-resolution single nucleotide polymorphism (SNP) mapping in four separate populations would lead to the identification of chronic obstructive pulmonary disease (COPD) susceptibility genes on chromosome 2q.
Methods: Within the chromosome 2q linkage region, 2,843 SNPs were genotyped in 806 COPD cases and 779 control subjects from Norway, and 2,484 SNPs were genotyped in 309 patients with severe COPD from the National Emphysema Treatment Trial and 330 community control subjects. Significant associations from the combined results across the two case-control studies were followed up in 1,839 individuals from 603 families from the International COPD Genetics Network (ICGN) and in 949 individuals from 127 families in the Boston Early-Onset COPD Study.
Measurements and Main Results: Merging the results of the two case-control analyses, 14 of the 790 overlapping SNPs had a combined P < 0.01. Two of these 14 SNPs were consistently associated with COPD in the ICGN families. The association with one SNP, located in the gene XRCC5, was replicated in the Boston Early-Onset COPD Study, with a combined P = 2.51 × 10−5 across the four studies, which remains significant when adjusted for multiple testing (P = 0.02). Genotype imputation confirmed the association with SNPs in XRCC5.
Conclusions: By combining data from COPD genetic association studies conducted in four independent patient samples, we have identified XRCC5, an ATP-dependent DNA helicase, as a potential COPD susceptibility gene.
emphysema; genetic linkage; metaanalysis; single nucleotide polymorphism
We introduce a stepwise approach for family-based designs for selecting a set of markers in a gene that are independently associated with the disease. The approach is based on testing the effect of a set of markers conditional on another set of markers. Several likelihood-based approaches have been proposed for special cases, but no model-free based tests have been proposed. We propose two types of tests in a family-based framework that are applicable to arbitrary family structures and completely robust to population stratification. We propose methods for ascertained dichotomous traits and unascertained quantitative traits. We first propose a completely model-free extension of the FBAT main genetic effect test. Then, for power issues, we introduce two model-based tests, one for dichotomous traits and one for continuous traits. Lastly, we utilize these tests to analyze a continuous lung function phenotype as a proxy for asthma in the Childhood Asthma Management Program. The methods are implemented in the free R package fbati.
Binary trait; Candidate gene analysis; Family-based association tests; FBAT-C; Linkage disequilibrium (LD); Model-based test; Model-free test; Nuclear families; Quantitative trait
Substantial evidence suggests that there is genetic susceptibility to chronic obstructive pulmonary disease (COPD). To identify common genetic risk variants, we performed a genome-wide association study in 2940 cases and 1380 smoking controls with normal lung function. We demonstrate a novel susceptibility locus at 4q22.1 in FAM13A (rs7671167, OR=0.76, P=8.6×10−8) and provide evidence of replication in one case-control and two family-based cohorts (for all studies, combined P=1.2×10−11).
Genetic variants influencing lung function in children and adults may ultimately lead to the development of chronic obstructive pulmonary disease (COPD), particularly in high-risk groups.
We tested for an association between single-nucleotide polymorphisms (SNPs) in the gene encoding matrix metalloproteinase 12 (MMP12) and a measure of lung function (prebronchodilator forced expiratory volume in 1 second [FEV1]) in more than 8300 subjects in seven cohorts that included children and adults. Within the Normative Aging Study (NAS), a cohort of initially healthy adult men, we tested for an association between SNPs that were associated with FEV1 and the time to the onset of COPD. We then examined the relationship between MMP12 SNPs and COPD in two cohorts of adults with COPD or at risk for COPD.
The minor allele (G) of a functional variant in the promoter region of MMP12 (rs2276109 [−82A→G]) was positively associated with FEV1 in a combined analysis of children with asthma and adult former and current smokers in all cohorts (P=2×10−6). This allele was also associated with a reduced risk of the onset of COPD in the NAS cohort (hazard ratio, 0.65; 95% confidence interval [CI], 0.46 to 0.92; P = 0.02) and with a reduced risk of COPD in a cohort of smokers (odds ratio, 0.63; 95% CI, 0.45 to 0.88; P = 0.005) and among participants in a family-based study of early-onset COPD (P = 0.006).
The minor allele of a SNP in MMP12 (rs2276109) is associated with a positive effect on lung function in children with asthma and in adults who smoke. This allele is also associated with a reduced risk of COPD in adult smokers.
Interleukin-6 (IL6) is a pleiotropic pro-inflammatory and immunomodulatory cytokine which likely plays an important role in the pathogenesis of COPD. There is a functional single nucleotide polymorphism (SNP), −174G/C, in the promoter region of IL6. We hypothesized that IL6 SNPs influence susceptibility for impaired lung function and COPD in smokers.
Seven and 5 SNPs in IL6 were genotyped in two nested case-control samples derived from the Lung Health Study (LHS) based on phenotypes of rate of decline of forced expiratory volume in one second (FEV1) over 5 years and baseline FEV1 at the beginning of the LHS. Serum IL6 concentrations were measured for all subjects. A partially overlapping panel of 9 IL6 SNPs was genotyped in 389 COPD cases from the National Emphysema Treatment Trial (NETT) and 420 controls from the Normative Aging Study (NAS).
In the LHS, three IL6 SNPs were associated with FEV1 decline (0.023 ≤ P ≤ 0.041 in additive models). Among them the IL6_−174C allele was associated with rapid decline of lung function. The association was more significant in a genotype-based analysis (P = 0.006). In the NETT-NAS study, IL6_−174G/C and four other IL6 SNPs, all of which are in linkage disequilibrium with IL6_−174G/C, were associated with susceptibility to COPD (0.01 ≤ P ≤ 0.04 in additive genetic models).
Our results suggest that the IL6_−174G/C SNP is associated with rapid decline of FEV1 and susceptibility to COPD in smokers.
genetic polymorphism; IL6; forced expiratory volume in one second (FEV1); lung function; chronic obstructive pulmonary disease (COPD)