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
The development of COPD in subjects with alpha-1 antitrypsin (AAT) deficiency is likely to be influenced by modifier genes. Genome-wide association studies and integrative genomics approaches in COPD have demonstrated significant associations with SNPs in the chromosome 15q region that includes CHRNA3 (cholinergic nicotine receptor alpha3) and IREB2 (iron regulatory binding protein 2).
We investigated whether SNPs in the chromosome 15q region would be modifiers for lung function and COPD in AAT deficiency.
The current analysis included 378 PIZZ subjects in the AAT Genetic Modifiers Study and a replication cohort of 458 subjects from the UK AAT Deficiency National Registry. Nine SNPs in LOC123688, CHRNA3 and IREB2 were selected for genotyping. FEV1 percent of predicted and FEV1/FVC ratio were analyzed as quantitative phenotypes. Family-based association analysis was performed in the AAT Genetic Modifiers Study. In the replication set, general linear models were used for quantitative phenotypes and logistic regression models were used for the presence/absence of emphysema or COPD.
Three SNPs (rs2568494 in IREB2, rs8034191 in LOC123688, and rs1051730 in CHRNA3) were associated with pre-bronchodilator FEV1 percent of predicted in the AAT Genetic Modifiers Study. Two SNPs (rs2568494 and rs1051730) were associated with the post-bronchodilator FEV1 percent of predicted and pre-bronchodilator FEV1/FVC ratio; SNP-by-gender interactions were observed. In the UK National Registry dataset, rs2568494 was significantly associated with emphysema in the male subgroup; significant SNP-by-smoking interactions were observed.
IREB2 and CHRNA3 are potential genetic modifiers of COPD phenotypes in individuals with severe AAT deficiency and may be sex-specific in their impact.
CHRNA3; Chronic obstructive pulmonary disease; Genetic association analysis; Genetic modifiers; IREB2
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
To identify non-invasive gene expression markers for chronic obstructive pulmonary disease (COPD), we performed genome-wide expression profiling of peripheral blood samples from 12 subjects with significant airflow obstruction and an equal number of non-obstructed controls. RNA was isolated from Peripheral Blood Mononuclear Cells (PBMCs) and gene expression was assessed using Affymetrix U133 Plus 2.0 arrays.
Tests for gene expression changes that discriminate between COPD cases (FEV1< 70% predicted, FEV1/FVC < 0.7) and controls (FEV1> 80% predicted, FEV1/FVC > 0.7) were performed using Significance Analysis of Microarrays (SAM) and Bayesian Analysis of Differential Gene Expression (BADGE). Using either test at high stringency (SAM median FDR = 0 or BADGE p < 0.01) we identified differential expression for 45 known genes. Correlation of gene expression with lung function measurements (FEV1 & FEV1/FVC), using both Pearson and Spearman correlation coefficients (p < 0.05), identified a set of 86 genes. A total of 16 markers showed evidence of significant correlation (p < 0.05) with quantitative traits and differential expression between cases and controls. We further compared our peripheral gene expression markers with those we previously identified from lung tissue of the same cohort. Two genes, RP9and NAPE-PLD, were identified as decreased in COPD cases compared to controls in both lung tissue and blood. These results contribute to our understanding of gene expression changes in the peripheral blood of patients with COPD and may provide insight into potential mechanisms involved in the disease.
Microarray; Biomarkers; PBMC
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)
Airflow limitation in COPD patients is not fully reversible. However, there may be large variability in bronchodilator responsiveness (BDR) among COPD patients, and familial aggregation of BDR suggests a genetic component. Therefore we investigated the association between six candidate genes and BDR in subjects with severe COPD. A total of 389 subjects from the National Emphysema Treatment Trial (NETT) were analyzed. Bronchodilator responsiveness to albuterol was expressed in three ways: absolute change in FEV1, change in FEV1 as a percent of baseline FEV1, and change in FEV1 as a percent of predicted FEV1. Genotyping was completed for 122 single nucleotide polymorphisms (SNPs) in six candidate genes (EPHX1, SFTPB, TGFB1, SERPINE2, GSTP1, ADRB2). Associations between BDR phenotypes and SNP genotypes were tested using linear regression, adjusting for age, sex, pack-years of smoking, and height. Genes associated with BDR phenotypes in the NETT subjects were assessed for replication in 127 pedigrees from the Boston Early-Onset COPD (EOCOPD) Study. Three SNPs in EPHX1 (p = 0.009 – 0.04), three SNPs in SERPINE2 (p = 0.004 – 0.05) and two SNPs in ADRB2 (0.04 – 0.05) were significantly associated with BDR phenotypes in NETT subjects. BDR. One SNP in EPHX1 (rs1009668, p = 0.04) was significantly replicated in EOCOPD subjects. SNPs in SFTPB, TGFB1, and GSTP1 genes were not associated with BDR. In conclusion, a polymorphism of EPHX1 was associated with bronchodilator responsiveness phenotypes in subjects with severe COPD.
bronchodilator responsiveness; chronic obstructive pulmonary disease; genetics; association analysis
Chronic obstructive pulmonary disease (COPD) is an inflammatory lung disorder with complex pathological features and largely unknown etiology. The identification of biomarkers for this disease could aid the development of methods to facilitate earlier diagnosis, the classification of disease subtypes, and provide a means to define therapeutic response. To identify gene expression biomarkers, we completed expression profiling of RNA derived from the lung tissue of 56 subjects with varying degrees of airflow obstruction using the Affymetrix U133 Plus 2.0 array. We applied multiple, independent analytical methods to define biomarkers for either discrete or quantitative disease phenotypes. Analysis of differential expression between cases (n = 15) and controls (n = 18) identified a set of 65 discrete biomarkers. Correlation of gene expression with quantitative measures of airflow obstruction (FEV1%predicted or FEV1/FVC) identified a set of 220 biomarkers. Biomarker genes were enriched in functions related to DNA binding and regulation of transcription. We used this group of biomarkers to predict disease in an unrelated data set, generated from patients with severe emphysema, with 97% accuracy. Our data contribute to the understanding of gene expression changes occurring in the lung tissue of patients with obstructive lung disease and provide additional insight into potential mechanisms involved in the disease process. Furthermore, we present the first gene expression biomarker for COPD validated in an independent data set.
microarray; gene expression; emphysema; lung function
COPD exacerbations reduce quality of life and increase mortality. Genetic variation may explain the substantial variability seen in exacerbation frequency among COPD subjects with similar lung function. We analyzed whether polymorphisms in five candidate genes previously associated with COPD susceptibility also demonstrate association with COPD exacerbations.
Eighty-eight single nucleotide polymorphisms in microsomal epoxide hydrolase (EPHX1), transforming growth factor beta 1 (TGFB1), SERPINE2, glutathione S-transferase pi (GSTP1), and surfactant protein B (SFTPB) were genotyped in 389 non-Hispanic white participants in the National Emphysema Treatment Trial. Exacerbations were defined as COPD-related emergency room visits or hospitalizations using Centers for Medicare and Medicaid Services claims data.
Measurements and Main Results
216 subjects (56%) experienced one or more exacerbations during the study period. An SFTPB promoter polymorphism, rs3024791, was associated with COPD exacerbations (p=0.008). Logistic regression models confirmed the association with rs3024791 (p = 0.007). Poisson regression models demonstrated association of multiple SFTPB SNPs with exacerbation rates: rs2118177 (p = 0.006), rs2304566 (p = 0.002), rs1130866 (p = 0.04), and rs3024791 (p = 0.002). Polymorphisms in EPHX1, GSTP1, TGFB1, and SERPINE2 did not demonstrate association with COPD exacerbations.
Variants in SFTPB are associated with COPD susceptibility and COPD exacerbation frequency.
association analysis; COPD; exacerbations; genetics; surfactant protein B; single nucleotide polymorphisms
Although a hereditary contribution to emphysema has been long suspected, severe α1-antitrypsin deficiency remains the only conclusively proven genetic risk factor for chronic obstructive pulmonary disease (COPD). Recently, genome-wide linkage analysis has led to the identification of two promising candidate genes for COPD: TGFB1 and SERPINE2. Like multiple other COPD candidate gene associations, even these positionally identified genes have not been universally replicated across all studies. Differences in phenotype definition may contribute to nonreplication in genetic studies of heterogeneous disorders such as COPD. The use of precisely measured phenotypes, including emphysema quantification on high-resolution chest computed tomography scans, has aided in the discovery of additional genes for clinically relevant COPD-related traits. The use of computed tomography scans to assess emphysema and airway disease as well as newer genetic technologies, including gene expression microarrays and genome-wide association studies, has great potential to detect novel genes affecting COPD susceptibility, severity, and response to treatment.
α1-antitrypsin deficiency; chronic obstructive pulmonary disease; genetic linkage; single-nucleotide polymorphism
In the National Emphysema Treatment Trial (NETT), marked variability in response to lung volume reduction surgery (LVRS) was observed. We sought to identify genetic differences which may explain some of this variability.
In 203 subjects from the NETT Genetics Ancillary Study, four outcome measures were used to define response to LVRS at six months: modified BODE index, post-bronchodilator FEV1, maximum work achieved on a cardiopulmonary exercise test, and University of California, San Diego shortness of breath questionnaire. Sixty-four single nucleotide polymorphisms (SNPs) were genotyped in five genes previously shown to be associated with chronic obstructive pulmonary disease susceptibility, exercise capacity, or emphysema distribution.
A SNP upstream from glutathione S-transferase pi (GSTP1; p = 0.003) and a coding SNP in microsomal epoxide hydrolase (EPHX1; p = 0.02) were each associated with change in BODE score. These effects appeared to be strongest in patients in the non-upper lobe predominant, low exercise subgroup. A promoter SNP in EPHX1 was associated with change in BODE score (p = 0.008), with the strongest effects in patients with upper lobe predominant emphysema and low exercise capacity. One additional SNP in GSTP1 and three additional SNPs in EPHX1 were associated (p < 0.05) with additional LVRS outcomes. None of these SNP effects were seen in 166 patients randomized to medical therapy.
Genetic variants in GSTP1 and EPHX1, two genes encoding xenobiotic metabolizing enzymes, were predictive of response to LVRS. These polymorphisms may identify patients most likely to benefit from LVRS.
Rationale: Patients with severe chronic obstructive pulmonary disease (COPD) may have varying levels of disability despite similar levels of lung function. This variation may reflect different COPD subtypes, which may have different genetic predispositions.
Objectives: To identify genetic associations for COPD-related phenotypes, including measures of exercise capacity, pulmonary function, and respiratory symptoms.
Methods: In 304 subjects from the National Emphysema Treatment Trial, we genotyped 80 markers in 22 positional and/or biologically plausible candidate genes. Regression models were used to test for association, using a test–replication approach to guard against false-positive results. For significant associations, effect estimates were recalculated using the entire cohort. Positive associations with dyspnea were confirmed in families from the Boston Early-Onset COPD Study.
Results: The test–replication approach identified four genes—microsomal epoxide hydrolase (EPHX1), latent transforming growth factor-β binding protein-4 (LTBP4), surfactant protein B (SFTPB), and transforming growth factor-β1 (TGFB1)—that were associated with COPD-related phenotypes. In all subjects, single-nucleotide polymorphisms (SNPs) in EPHX1 (p ⩽ 0.03) and in LTBP4 (p ⩽ 0.03) were associated with maximal output on cardiopulmonary exercise testing. Markers in LTBP4 (p ⩽ 0.05) and SFTPB (p = 0.005) were associated with 6-min walk test distance. SNPs in EPHX1 were associated with carbon monoxide diffusing capacity (p ⩽ 0.04). Three SNPs in TGFB1 were associated with dyspnea (p ⩽ 0.002), one of which replicated in the family study (p = 0.02).
Conclusions: Polymorphisms in several genes seem to be associated with COPD-related traits other than FEV1. These associations may identify genes in pathways important for COPD pathogenesis.
dyspnea; emphysema; exercise tolerance; genetic association; pulmonary function tests
Case-control studies have successfully identified many significant genetic associations for complex diseases, but lack of replication has been a criticism of case-control genetic association studies in general. We selected 12 candidate genes with reported associations to chronic obstructive pulmonary disease (COPD) and genotyped 29 polymorphisms in a family-based study and in a case-control study. In the Boston Early-Onset COPD Study families, significant associations with quantitative and/or qualitative COPD-related phenotypes were found for the tumor necrosis factor (TNF)-α −308G>A promoter polymorphism (P < 0.02), a coding variant in surfactant protein B (SFTPB Thr131Ile) (P = 0.03), and the (GT)31 allele of the heme oxygenase (HMOX1) promoter short tandem repeat (P = 0.02). In the case-control study, the SFTPB Thr131Ile polymorphism was associated with COPD, but only in the presence of a gene-by-environment interaction term (P = 0.01 for both main effect and interaction). The 30-repeat, but not the 31-repeat, allele of HMOX1 was associated (P = 0.04). The TNF −308G>A polymorphism was not significant. In addition, the microsomal epoxide hydrolase “fast” allele (EPHX1 His139Arg) was significantly associated in the case-control study (P = 0.03). Although some evidence for replication was found for SFTPB and HMOX1, none of the previously published COPD genetic associations was convincingly replicated across both study designs.
association studies; case-control studies; emphysema; genetics; single nucleotide polymorphism