Severe α1‐antitrypsin (AAT) deficiency is an autosomal recessive genetic condition associated with an increased but variable risk for chronic obstructive pulmonary disease (COPD). A study was undertaken to assess the impact of chronic bronchitis, pneumonia, asthma and sex on the development of COPD in individuals with severe AAT deficiency.
The AAT Genetic Modifier Study is a multicentre family‐based cohort study designed to study the genetic and epidemiological determinants of COPD in AAT deficiency. 378 individuals (age range 33–80 years), confirmed to be homozygous for the SERPINA1 Z mutation, were included in the analyses. The primary outcomes of interest were a quantitative outcome, forced expiratory volume in 1 s (FEV1) percentage predicted, and a qualitative outcome, severe airflow obstruction (FEV1 <50% predicted).
In multivariate analysis of the overall cohort, cigarette smoking, sex, asthma, chronic bronchitis and pneumonia were risk factors for reduced FEV1 percentage predicted and severe airflow obstruction (p<0.01). Index cases had lower FEV1 values, higher smoking histories and more reports of adult asthma, pneumonia and asthma before age 16 than non‐index cases (p<0.01). Men had lower pre‐ and post‐bronchodilator FEV1 percentage predicted than women (p<0.0001); the lowest FEV1 values were observed in men reporting a history of childhood asthma (26.9%). This trend for more severe obstruction in men remained when index and non‐index groups were examined separately, with men representing the majority of non‐index individuals with airflow obstruction (71%). Chronic bronchitis (OR 3.8, CI 1.8 to 12.0) and a physician's report of asthma (OR 4.2, CI 1.4 to 13.1) were predictors of severe airflow obstruction in multivariate analysis of non‐index men but not women.
In individuals with severe AAT deficiency, sex, asthma, chronic bronchitis and pneumonia are risk factors for severe COPD, in addition to cigarette smoking. These results suggest that, in subjects severely deficient in AAT, men, individuals with symptoms of chronic bronchitis and/or a past diagnosis of asthma or pneumonia may benefit from closer monitoring and potentially earlier treatment.
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 impact of cigarette smoking can persist for extended periods following smoking cessation and may involve epigenetic reprogramming. Changes in DNA methylation associated with smoking may help to identify molecular pathways that contribute to the latency between exposure and disease onset. Cross-sectional cohort data from subjects in the International COPD Genetics Network (n = 1085) and the Boston Early-Onset COPD study (n = 369) were analyzed as the discovery and replication cohorts, respectively. Genome-wide methylation data on 27 578 CpG sites in 14 475 genes were obtained on DNA from peripheral blood leukocytes using the Illumina HumanMethylation27K Beadchip in both cohorts. We identified 15 sites significantly associated with current smoking, 2 sites associated with cumulative smoke exposure, and, within the subset of former smokers, 3 sites associated with time since quitting cigarettes. Two loci, factor II receptor-like 3 (F2RL3) and G-protein-coupled receptor 15 (GPR15), were significantly associated in all three analyses and were validated by pyrosequencing. These findings (i) identify a novel locus (GPR15) associated with cigarette smoking and (ii) suggest the existence of dynamic, site-specific methylation changes in response to smoking which may contribute to the extended risks associated with cigarette smoking that persist after cessation.
Alpha-1 antitrypsin (AAT) deficiency and tobacco smoking are confirmed risk factors for Chronic Obstructive Pulmonary Disease. We hypothesized that variable DNA methylation would be associated with smoking and inflammation, as reflected by the level of C-Reactive Protein (CRP) in AAT-deficient subjects. Methylation levels of 1,411 autosomal CpG sites from the Illumina GoldenGate Methylation Cancer Panel I were analyzed in 316 subjects. Associations of five smoking behaviors and CRP levels with individual CpG sites and average methylation levels were assessed using non-parametric testing, linear regression and linear mixed effect models, with and without adjustment for age and gender. Univariate linear regression analysis revealed that methylation levels of 16 CpG sites significantly associated with ever-smoking status. A CpG site in the TGFBI gene was the only site associated with ever-smoking after adjustment for age and gender. No highly significant associations existed between age at smoking initiation, pack-years smoked, duration of smoking, and time since quitting smoking as predictors of individual CpG site methylation levels. However, ever-smoking and younger age at smoking initiation associated with lower methylation level averaged across all sites. DNA methylation at CpG sites in the RUNX3, JAK3 and KRT1 genes associated with CRP levels. The most significantly associated CpG sites with gender and age mapped to the CASP6 and FZD9 genes, respectively. In summary, this study identified multiple potential candidate CpG sites associated with ever-smoking and CRP level in AAT-deficient subjects. Phenotypic variability in Mendelian diseases may be due to epigenetic factors.
68kDa (TGFBI); C-Reactive Protein (CRP); Chronic Obstructive Pulmonary Disease (COPD); Illumina GoldenGate Methylation Cancer Panel I; alpha-1 antitrypsin (AAT) deficiency; beta-induced; methylation; smoking behaviors; transforming growth factor
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
Rationale: Chronic obstructive pulmonary disease (COPD) is associated with local (lung) and systemic (blood) inflammation and manifestations. DNA methylation is an important regulator of gene transcription, and global and specific gene methylation marks may vary with cigarette smoke exposure.
Objectives: To perform a comprehensive assessment of methylation marks in DNA from subjects well phenotyped for nonneoplastic lung disease.
Methods: We conducted array-based methylation screens, using a test-replication approach, in two family-based cohorts (n = 1,085 and 369 subjects).
Measurements and Main Results: We observed 349 CpG sites significantly associated with the presence and severity of COPD in both cohorts. Seventy percent of the associated CpG sites were outside of CpG islands, with the majority of CpG sites relatively hypomethylated. Gene ontology analysis based on these 349 CpGs (330 genes) suggested the involvement of a number of genes responsible for immune and inflammatory system pathways, responses to stress and external stimuli, as well as wound healing and coagulation cascades. Interestingly, our observations include significant, replicable associations between SERPINA1 hypomethylation and COPD and lower average lung function phenotypes (combined P values: COPD, 1.5 × 10−23; FEV1/FVC, 1.5 × 10−35; FEV1, 2.2 × 10−40).
Conclusions: Genetic and epigenetic pathways may both contribute to COPD. Many of the top associations between COPD and DNA methylation occur in biologically plausible pathways. This large-scale analysis suggests that DNA methylation may be a biomarker of COPD and may highlight new pathways of COPD pathogenesis.
chronic obstructive pulmonary disease; epigenetics; DNA methylation; smoking
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.
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
To investigate the association between methylation of transposable elements Alu and long-interspersed nuclear elements (LINE-1) and lung function.
Outpatient Veterans Administration facilities in greater Boston, Massachusetts, USA.
Individuals from the Veterans Administration Normative Aging Study, a longitudinal study of aging in men, evaluated between 1999 and 2007. The majority (97%) were white.
Primary and secondary outcome measures
Primary predictor was methylation, assessed using PCR-pyrosequencing after bisulphite treatment. Primary outcome was lung function as assessed by spirometry, performed according to American Thoracic Society/European Respiratory Society guidelines at the same visit as the blood draws.
In multivariable models adjusted for age, height, body mass index (BMI), pack-years of smoking, current smoking and race, Alu hypomethylation was associated with lower forced expiratory volume in 1 s (FEV1) (β=28 ml per 1% change in Alu methylation, p=0.017) and showed a trend towards association with a lower forced vital capacity (FVC) (β=27 ml, p=0.06) and lower FEV1/FVC (β=0.3%, p=0.058). In multivariable models adjusted for age, height, BMI, pack-years of smoking, current smoking, per cent lymphocytes, race and baseline lung function, LINE-1 hypomethylation was associated with more rapid decline of FEV1 (β=6.9 ml/year per 1% change in LINE-1 methylation, p=0.005) and of FVC (β=9.6 ml/year, p=0.002).
In multiple regression analysis, Alu hypomethylation was associated with lower lung function, and LINE-1 hypomethylation was associated with more rapid lung function decline in a cohort of older and primarily white men from North America. Future studies should aim to replicate these findings and determine if Alu or LINE-1 hypomethylation may be due to specific and modifiable environmental exposures.
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
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
COPD; Genetics; Association analysis; Consortium
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).
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
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
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.