Rationale: Genome-wide association studies have shown significant associations between variants near hedgehog interacting protein HHIP, FAM13A, and cholinergic nicotinic acetylcholine receptor CHRNA3/5 with increased risk of chronic obstructive pulmonary disease (COPD) in smokers; however, the disease mechanisms behind these associations are not well understood.
Objectives: To identify the association between replicated loci and COPD-related phenotypes in well-characterized patient populations.
Methods: The relationship between these three loci and COPD-related phenotypes was assessed in the Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-point (ECLIPSE) cohort. The results were validated in the family-based International COPD Genetics Network (ICGN).
Measurements and Main Results: The CHRNA3/5 locus was significantly associated with pack-years of smoking (P = 0.002 and 3 × 10−4), emphysema assessed by a radiologist using high-resolution computed tomography (P = 2 × 10−4 and 4.8 × 10−5), and airflow obstruction (P = 0.004 and 1.8 × 10−5) in the ECLIPSE and ICGN populations, respectively. However, variants in the IREB2 gene were only significantly associated with FEV1. The HHIP locus was not associated with smoking intensity but was associated with FEV1/FVC (P = 1.9 × 10−4 and 0.004 in the ECLIPSE and ICGN populations). The HHIP locus was also associated with fat-free body mass (P = 0.007) and with both retrospectively (P = 0.015) and prospectively (P = 0.024) collected COPD exacerbations in the ECLIPSE cohort. Single-nucleotide polymorphisms in the FAM13A locus were associated with lung function.
Conclusions: The CHRNA3/5 locus was associated with increased smoking intensity and emphysema in individuals with COPD, whereas the HHIP and FAM13A loci were not associated with smoking intensity. The HHIP locus was associated with the systemic components of COPD and with the frequency of COPD exacerbations. FAM13A locus was associated with lung function.
COPD exacerbations; nicotine addiction; high-resolution CT; genetic association analysis; emphysema
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.
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
Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.
Superoxide dismutase-3 (SOD3) is a major extracellular antioxidant enzyme, and previous studies have indicated a possible role of this gene in chronic obstructive pulmonary disease (COPD). We hypothesized that polymorphisms in the SOD3 gene would be associated with COPD and COPD-related phenotypes.
We genotyped three SOD3 polymorphisms (rs8192287 (E1), rs8192288 (I1) and rs1799895 (R213G)) in a case-control cohort, with severe COPD cases from the National Emphysema Treatment Trial (NETT, n=389) and smoking controls from the Normative Aging Study (NAS, n=472). We examined whether the SNPs were associated with COPD status, lung function variables, and quantitative CT measurements of emphysema and airway wall thickness. Further, we tried to replicate our initial findings in two family-based studies, the International COPD Genetics Network (ICGN, n=3061) and the Boston Early-Onset COPD Study (EOCOPD, n=949).
In NETT COPD cases, the minor alleles of SNPs E1 and I1 were associated with a higher percentage of emphysema (%LAA950) on chest CT scan (p=0.029 and p=0.0058). The association with E1 was replicated in the ICGN family study, where the minor allele was associated with more emphysema (p=0.048). Airway wall thickness was positively associated with the E1 SNP in ICGN; however, this finding was not confirmed in NETT. Quantitative CT data were not available in EOCOPD. The SNPs were not associated with lung function variables or COPD status in any of the populations.
In conclusion, polymorphisms in the SOD3 gene were associated with CT emphysema but not COPD susceptibility, highlighting the importance of phenotype definition in COPD genetics studies.
To identify COPD associated gene susceptibility and lung function in a longitudinal cohort including COPD and subjects who were at risk for developing COPD, and to replicate this in two cross-sectional and longitudinal populations in Chinese Han population.
Three cohorts were recruited in this study, including an 18-year follow-up population (306 COPD and 743 control subjects) in one village in 1992 and it changed to 409 COPD and 611 controls in 2010, a 2 year follow-up study in another village (374 COPD and 377 controls) and another 2 year follow-up one in a city (541 COPD and 560 controls) in 2010. Sixteen candidate single nucleotide polymorphisms (SNPs) were selected for genotyping. Among them, 5SNPs in or near HHIP, 1SNP in IREB2 and 1SNP in FAM13A were previously reported to be associated with COPD susceptibility or lung function decline. And another 9SNPs were selected from HapMap website as HHIP tags. In 2010, totaling 1,324 COPD patients and 1,548 healthy controls were finally included in our genetic susceptibility analyses.
We identified two new regions showing an association with COPD susceptibility in the Human Hedgehog interacting protein (HHIP) rs11100865 and rs7654947, and we confirmed that the family with sequence similarity 13 member A gene (FAM13A) rs7671167 was associated with the development of COPD in Chinese Han population. And the HHIP rs7654947 and FAM13A rs7671167 were associated with lung function decline, and this result was replicated in other two populations.
These results suggest an important role of the HHIP and FAM13A regions as genetic risk factors for COPD development and lung function decline in Chinese Han population. Future research on these genes should focus on the molecular mechanisms of these genes on developing COPD and creating therapies to alleviate reduced lung function.
Electronic supplementary material
The online version of this article (doi:10.1186/s12931-015-0209-3) contains supplementary material, which is available to authorized users.
Chronic obstructive pulmonary disease; Cross-sectional; Longitudinal; Gene susceptibility; Lung function
We examined the association between single-nucleotide polymorphisms (SNPs) previously associated with chronic obstructive pulmonary disease (COPD) and/or lung function with COPD and COPD-related phenotypes in a novel cohort of patients with severe to very severe COPD. We examined 315 cases of COPD and 330 Caucasian control smokers from Poland. We included three SNPs previously associated with COPD: rs7671167 (FAM13A), rs13180 (IREB2), and rs8034191 (CHRNA 3/5), and four SNPs associated with lung function in a genome-wide association study of general population samples: rs2070600 (AGER), rs11134242 (ADCY2), rs4316710 (THSD4), and rs17096090 (INTS12). We tested for associations with severe COPD and COPD-related phenotypes, including lung function, smoking behavior, and body mass index. Subjects with COPD were older (average age 62 versus 58 years, P < 0.01), with more pack-years of smoking (45 versus 33 pack-years, P < 0.01). CHRNA3/5 (odds ratio [OR], 1.89; 95% confidence interval [CI], 1.5–2.4; P = 7.4 × 10−7), IREB2 (OR, 0.69; 95% CI, 0.5–0.9; P = 3.4 × 10−3), and ADCY2 (OR, 1.35; 95% CI, 1.1–1.7; P = 0.01) demonstrated significant associations with COPD. FAM13A (OR, 0.8; 95% CI, 0.7–1.0; P = 0.11) approached statistical significance. FAM13A and ADCY2 also demonstrated a significant association with lung function. Thus, in severe to very severe COPD, we demonstrate a replication of association between two SNPs previously associated with COPD (CHRNA3/5 and IREB2), as well as an association with COPD of one locus initially associated with lung function (ADCY2).
chronic obstructive pulmonary disease; genetic association analysis; lung function; smoking; nicotine addiction
Lung cancer is a leading cause of cancer death and disease burden in many countries. Understanding of the biological pathways involved in lung cancer aetiology is required to identify key biomolecules that could be of significant clinical value, either as predictive, prognostic or diagnostic markers, or as targets for the development of novel therapies to treat this disease, in addition to smoking avoidance strategies. Genome-wide association studies (GWAS) have enabled significant progress in the past 5 years in investigating genetic susceptibility to lung cancer. Large scale, multi-cohort GWAS of mainly Caucasian, smoking, populations have identified strong associations for lung cancer mapped to chromosomal regions 15q [nicotinic acetylcholine receptor (nAChR) subunits: CHRNA3, CHRNA5], 5p (TERT-CLPTM1L locus) and 6p (BAT3-MSH5). Some studies in Asian populations of smokers have found similar risk loci, whereas GWAS in never smoking Asian females have identified associations in other chromosomal regions, e.g., 3q (TP63), that are distinct from smoking-related lung cancer risk loci. GWAS of smoking behaviour have identified risk loci for smoking quantity at 15q (similar genes to lung cancer susceptibility: CHRNA3, CHRNA5) and 19q (CYP2A6). Other genes have been mapped for smoking initiation and smoking cessation. In chronic obstructive pulmonary disease (COPD), which is a known risk factor for lung cancer, GWAS in large cohorts have also found CHRNA3 and CHRNA5 single nucleotide polymorphisms (SNPs) mapping at 15q as risk loci, as well as other regions at 4q31 (HHIP), 4q24 (FAM13A) and 5q (HTR4). The overlap in risk loci between lung cancer, smoking behaviour and COPD may be due to the effects of nicotine addiction; however, more work needs to be undertaken to explore the potential direct effects of nicotine and its metabolites in gene-environment interaction in these phenotypes. Goals of future genetic susceptibility studies of lung cancer should focus on refining the strongest risk loci in a wide range of populations with lung cancer, and integrating other clinical and biomarker information, in order to achieve the aim of personalised therapy for lung cancer.
Lung cancer; genetics; pulmonary disease; chronic obstructive; genome-wide association study (GWAS)
Cigarette smoking is the major environmental risk factor for chronic obstructive pulmonary disease (COPD). Genome-wide association studies have provided compelling associations for three loci with COPD. In this study, we aimed to estimate direct, i.e., independent from smoking, and indirect effects of those loci on COPD development using mediation analysis. We included a total of 3,424 COPD cases and 1,872 unaffected controls with data on two smoking-related phenotypes: lifetime average smoking intensity and cumulative exposure to tobacco smoke (pack years). Our analysis revealed that effects of two linked variants (rs1051730 and rs8034191) in the AGPHD1/CHRNA3 cluster on COPD development are significantly, yet not entirely, mediated by the smoking-related phenotypes. Approximately 30 % of the total effect of variants in the AGPHD1/CHRNA3 cluster on COPD development was mediated by pack years. Simultaneous analysis of modestly (r2 = 0.21) linked markers in CHRNA3 and IREB2 revealed that an even larger (~42 %) proportion of the total effect of the CHRNA3 locus on COPD was mediated by pack years after adjustment for an IREB2 single nucleotide polymorphism. This study confirms the existence of direct effects of the AGPHD1/CHRNA3, IREB2, FAM13A and HHIP loci on COPD development. While the association of the AGPHD1/CHRNA3 locus with COPD is significantly mediated by smoking-related phenotypes, IREB2 appears to affect COPD independently of smoking.
The HHIP gene, encoding Hedgehog interacting protein, has been implicated in chronic obstructive pulmonary disease (COPD) by genome-wide association studies (GWAS), and our subsequent studies identified a functional upstream genetic variant that decreased HHIP transcription. However, little is known about how HHIP contributes to COPD pathogenesis.
We exposed Hhip haploinsufficient mice (Hhip+/-) to cigarette smoke (CS) for 6 months to model the biological consequences caused by CS in human COPD risk-allele carriers at the HHIP locus. Gene expression profiling in murine lungs was performed followed by an integrative network inference analysis, PANDA (Passing Attributes between Networks for Data Assimilation) analysis.
We detected more severe airspace enlargement in Hhip+/- mice vs. wild-type littermates (Hhip+/+) exposed to CS. Gene expression profiling in murine lungs suggested enhanced lymphocyte activation pathways in CS-exposed Hhip+/- vs. Hhip+/+ mice, which was supported by increased numbers of lymphoid aggregates and enhanced activation of CD8+ T cells after CS-exposure in the lungs of Hhip+/-mice compared to Hhip+/+ mice. Mechanistically, results from PANDA network analysis suggested a rewired and dampened Klf4 signaling network in Hhip+/- mice after CS exposure.
In summary, HHIP haploinsufficiency exaggerated CS-induced airspace enlargement, which models CS-induced emphysema in human smokers carrying COPD risk alleles at the HHIP locus. Network modeling suggested rewired lymphocyte activation signaling circuits in the HHIP haploinsufficiency state.
Electronic supplementary material
The online version of this article (doi:10.1186/s13073-015-0137-3) contains supplementary material, which is available to authorized users.
Rationale: Genomic loci are associated with FEV1 or the ratio of FEV1 to FVC in population samples, but their association with chronic obstructive pulmonary disease (COPD) has not yet been proven, nor have their combined effects on lung function and COPD been studied.
Objectives: To test association with COPD of variants at five loci (TNS1, GSTCD, HTR4, AGER, and THSD4) and to evaluate joint effects on lung function and COPD of these single-nucleotide polymorphisms (SNPs), and variants at the previously reported locus near HHIP.
Methods: By sampling from 12 population-based studies (n = 31,422), we obtained genotype data on 3,284 COPD case subjects and 17,538 control subjects for sentinel SNPs in TNS1, GSTCD, HTR4, AGER, and THSD4. In 24,648 individuals (including 2,890 COPD case subjects and 13,862 control subjects), we additionally obtained genotypes for rs12504628 near HHIP. Each allele associated with lung function decline at these six SNPs contributed to a risk score. We studied the association of the risk score to lung function and COPD.
Measurements and Main Results: Association with COPD was significant for three loci (TNS1, GSTCD, and HTR4) and the previously reported HHIP locus, and suggestive and directionally consistent for AGER and TSHD4. Compared with the baseline group (7 risk alleles), carrying 10–12 risk alleles was associated with a reduction in FEV1 (β = –72.21 ml, P = 3.90 × 10−4) and FEV1/FVC (β = –1.53%, P = 6.35 × 10−6), and with COPD (odds ratio = 1.63, P = 1.46 × 10−5).
Conclusions: Variants in TNS1, GSTCD, and HTR4 are associated with COPD. Our highest risk score category was associated with a 1.6-fold higher COPD risk than the population average score.
FEV1; FVC; genome-wide association study; modeling risk
Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.
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
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: A genome-wide association study (GWAS) for circulating chronic obstructive pulmonary disease (COPD) biomarkers could identify genetic determinants of biomarker levels and COPD susceptibility.
Objectives: To identify genetic variants of circulating protein biomarkers and novel genetic determinants of COPD.
Methods: GWAS was performed for two pneumoproteins, Clara cell secretory protein (CC16) and surfactant protein D (SP-D), and five systemic inflammatory markers (C-reactive protein, fibrinogen, IL-6, IL-8, and tumor necrosis factor-α) in 1,951 subjects with COPD. For genome-wide significant single nucleotide polymorphisms (SNPs) (P < 1 × 10−8), association with COPD susceptibility was tested in 2,939 cases with COPD and 1,380 smoking control subjects. The association of candidate SNPs with mRNA expression in induced sputum was also elucidated.
Measurements and Main Results: Genome-wide significant susceptibility loci affecting biomarker levels were found only for the two pneumoproteins. Two discrete loci affecting CC16, one region near the CC16 coding gene (SCGB1A1) on chromosome 11 and another locus approximately 25 Mb away from SCGB1A1, were identified, whereas multiple SNPs on chromosomes 6 and 16, in addition to SNPs near SFTPD, had genome-wide significant associations with SP-D levels. Several SNPs affecting circulating CC16 levels were significantly associated with sputum mRNA expression of SCGB1A1 (P = 0.009–0.03). Several SNPs highly associated with CC16 or SP-D levels were nominally associated with COPD in a collaborative GWAS (P = 0.001–0.049), although these COPD associations were not replicated in two additional cohorts.
Conclusions: Distant genetic loci and biomarker-coding genes affect circulating levels of COPD-related pneumoproteins. A subset of these protein quantitative trait loci may influence their gene expression in the lung and/or COPD susceptibility.
Clinical trial registered with www.clinicaltrials.gov (NCT 00292552).
biomarker; chronic obstructive pulmonary disease; genome-wide association study
Chronic obstructive pulmonary disease (COPD) is a heterogeneous syndrome, including emphysema and airway disease. Phenotypes defined on the basis of chest computed tomography (CT) may decrease disease heterogeneity and aid in the identification of candidate genes for COPD subtypes. To identify these genes, we performed genome-wide linkage analysis in extended pedigrees from the Boston Early-Onset COPD Study, stratified by emphysema status (defined by chest CT scans) of the probands, followed by genetic association analysis of positional candidate genes. A region on chromosome 1p showed strong evidence of linkage to lung function traits in families of emphysema-predominant probands in the stratified analysis (LOD score = 2.99 in families of emphysema-predominant probands versus 1.98 in all families). Association analysis in 949 individuals from 127 early-onset COPD pedigrees revealed association for COPD-related traits with an intronic single-nucleotide polymorphism (SNP) in transforming growth factor-β receptor-3 (TGFBR3) (P = 0.005). This SNP was significantly associated with COPD affection status comparing 389 cases from the National Emphysema Treatment Trial to 472 control smokers (P = 0.04), and with FEV1 (P = 0.004) and CT emphysema (P = 0.05) in 3,117 subjects from the International COPD Genetics Network. Gene-level replication of association with lung function was seen in 427 patients with COPD from the Lung Health Study. In conclusion, stratified linkage analysis followed by association testing identified TGFBR3 (betaglycan) as a potential susceptibility gene for COPD. Published human microarray and murine linkage studies have also demonstrated the importance of TGFBR3 in emphysema and lung function, and our group and others have previously found association of COPD-related traits with TGFB1, a ligand for TGFBR3.
betaglycan; chronic obstructive pulmonary disease; computed tomography; linkage; single nucleotide polymorphism
Several single nucleotide polymorphisms (SNPs) in an α-neuronal nicotinic acetylcholine receptor subunit (CHRNA3/5) were identified to be associated with chronic obstructive pulmonary disease (COPD) in a study based on a Norwegian population. However, results from subsequent studies have been controversial, particularly in studies recruiting Asians. In the present study, we conducted a comprehensive search and meta-analyses to identify susceptibility SNPs for COPD in the CHRNA3/5 locus.
A comprehensive literature search was conducted to find studies that have reported an association between SNPs in the CHRNA3/5 locus and COPD risk. Pooled odds ratios (ORs) with 95% confidence intervals (CIs) for each SNP were calculated with the major allele or genotype as the reference group. The influence of individual studies on pooled measures was assessed, in addition to publication bias.
A total of 12 articles with 14 eligible studies were included in this analysis. Association between 4 SNPs in the CHRNA3/5 locus and COPD was evaluated and included rs1051730, rs8034191, rs6495309, and rs16969968. Significant associations between the 4 SNPs and COPD were identified under allele (rs1051730: OR = 1.14, 95%CI = 1.10–1.18; rs8034191: OR = 1.29, 95%CI = 1.18–1.41; rs6495309: OR = 1.26, 95%CI = 1.09–1.45; rs16969968: OR = 1.27, 95%CI = 1.17–1.39) and genotype models. Subgroup analysis conducted for rs1051730 showed a significant association between this SNP and COPD risk in non-Asians (OR = 1.14, 95%CI = 1.10–1.18), but not Asians (OR = 1.23, 95%CI = 0.91–1.67). Rs1051730 and rs6495309 were also significantly associated with COPD after adjusting for multiple variables, including age and smoking status.
Our results indicate that 4 SNPs in the CHRNA3/5 locus are associated with COPD risk. Rs1051730 was particularly associated with COPD in non-Asians, but its role in Asians still needs to be verified. Additional studies will be necessary to assess the effect of rs6495309 on COPD. Although rs1051730 and rs6495309 were shown to be independent risk factors for COPD, validation studies should be performed.
Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.
The marked variability in individual susceptibility to the detrimental effects of smoking on lung function and findings suggest a significant genetic contribution to COPD, which has been demonstrated in several studies. The only known genetic risk factor for COPD, severe alpha 1 antitrypsin (AAT) deficiency, explains only 1–2% of cases of this disease. Screening for severe AAT should be conducted in all cases of COPD. Intravenous augmentation therapy should be combined with currently recommended treatment modalities for COPD when treating patients with COPD due to severe AAT deficiency. There is considerable interest in identifying susceptibility genes for COPD unrelated to severe AAT deficiency, as this could greatly enhance current efforts to prevent, diagnose and treat this disease by yielding novel insights into its pathogenesis. Genome-wide association studies (GWAS) of COPD and its intermediate phenotypes (e.g., lung function measures) have identified novel susceptibility loci for COPD. Some of these susceptibility loci may also influence lung function in the general population (e.g., HHIP and FAM13A), while others may affect not only COPD but other diseases related to smoking behavior (e.g., CHRNA3/CHRNA5). Although much work remains to be done, recent advances and the implementation of novel approaches to study COPD genetics (e.g., sequencing) and epigenetics are promising, and could have a profound impact on COPD management.
COPD; genetics; alpha 1 antitrypsin deficiency
Cigarette smoking is a major risk factor for COPD and COPD severity. Previous genome-wide association studies (GWAS) have identified numerous single nucleotide polymorphisms (SNPs) associated with the number of cigarettes smoked per day (CPD) and a Dopamine Beta-Hydroxylase (DBH) locus associated with smoking cessation in multiple populations.
To identify SNPs associated with lifetime average and current CPD, age at smoking initiation, and smoking cessation in COPD subjects.
GWAS were conducted in 4 independent cohorts encompassing 3,441 ever-smoking COPD subjects (GOLD stage II or higher). Untyped SNPs were imputed using HapMap (phase II) panel. Results from all cohorts were meta-analyzed.
Several SNPs near the HLA region on chromosome 6p21 and in an intergenic region on chromosome 2q21 showed associations with age at smoking initiation, both with the lowest p=2×10−7. No SNPs were associated with lifetime average CPD, current CPD or smoking cessation with p<10−6. Nominally significant associations with candidate SNPs within alpha-nicotinic acetylcholine receptors 3/5 (CHRNA3/CHRNA5; e.g. p=0.00011 for SNP rs1051730) and Cytochrome P450 2A6 (CYP2A6; e.g. p=2.78×10−5 for a nonsynonymous SNP rs1801272) regions were observed for lifetime average CPD, however only CYP2A6 showed evidence of significant association with current CPD. A candidate SNP (rs3025343) in the DBH was significantly (p=0.015) associated with smoking cessation.
We identified two candidate regions associated with age at smoking initiation in COPD subjects. Associations of CHRNA3/CHRNA5 and CYP2A6 loci with CPD and DBH with smoking cessation are also likely of importance in the smoking behaviors of COPD patients.
Chronic Obstructive Pulmonary Disease (COPD); Genome Wide Association study (GWAS); smoking behaviors; Single Nucleotide Polymorphism (SNP)
The ratio of forced expiratory volume in one second to forced vital capacity (FEV1/FVC) is a measure used to diagnose airflow obstruction and is highly heritable. We performed a genome-wide association study in 7,691 Framingham Heart Study participants to identify single-nucleotide polymorphisms (SNPs) associated with the FEV1/FVC ratio, analyzed as a percent of the predicted value. Identified SNPs were examined in an independent set of 835 Family Heart Study participants enriched for airflow obstruction. Four SNPs in tight linkage disequilibrium on chromosome 4q31 were associated with the percent predicted FEV1/FVC ratio with p-values of genome-wide significance in the Framingham sample (best p-value = 3.6e-09). One of the four chromosome 4q31 SNPs (rs13147758; p-value 2.3e-08 in Framingham) was genotyped in the Family Heart Study and produced evidence of association with the same phenotype, percent predicted FEV1/FVC (p-value = 2.0e-04). The effect estimates for association in the Framingham and Family Heart studies were in the same direction, with the minor allele (G) associated with higher FEV1/FVC ratio levels. Results from the Family Heart Study demonstrated that the association extended to FEV1 and dichotomous airflow obstruction phenotypes, particularly among smokers. The SNP rs13147758 was associated with the percent predicted FEV1/FVC ratio in independent samples from the Framingham and Family Heart Studies producing a combined p-value of 8.3e-11, and this region of chromosome 4 around 145.68 megabases was associated with COPD in three additional populations reported in the accompanying manuscript. The associated SNPs do not lie within a gene transcript but are near the hedgehog-interacting protein (HHIP) gene and several expressed sequence tags cloned from fetal lung. Though it is unclear what gene or regulatory effect explains the association, the region warrants further investigation.
Cigarette smoking is the primary risk factor for impaired lung function, yet only 20% of smokers develop chronic obstructive pulmonary disease (COPD). This observation, along with family studies of lung function and COPD, suggests that genetic factors influence susceptibility to cigarette smoke. We examined the relationship between common genetic variants and measures of lung function in a sample of 7,691 participants from the Framingham Heart Study and confirmed our observations in 835 participants from the Family Heart Study selected to include cases of airflow obstruction. We identified a variant on chromosome 4 that was strongly associated with FEV1/FVC in the Framingham Study and confirmed the association in the Family Heart Study. The accompanying manuscript identified the same region to be associated with COPD. Several interesting genes are present in the region that we identified, including a gene (HHIP) interacting with a biological pathway involved in lung development, but it is not yet clear which gene in the region explains the association. Our results identified a region of chromosome 4 that warrants further study to understand the genetic effects influencing lung function.
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: Genome-wide association studies (GWAS) have identified loci influencing lung function, but fewer genes influencing chronic obstructive pulmonary disease (COPD) are known.
Objectives: Perform meta-analyses of GWAS for airflow obstruction, a key pathophysiologic characteristic of COPD assessed by spirometry, in population-based cohorts examining all participants, ever smokers, never smokers, asthma-free participants, and more severe cases.
Methods: Fifteen cohorts were studied for discovery (3,368 affected; 29,507 unaffected), and a population-based family study and a meta-analysis of case-control studies were used for replication and regional follow-up (3,837 cases; 4,479 control subjects). Airflow obstruction was defined as FEV1 and its ratio to FVC (FEV1/FVC) both less than their respective lower limits of normal as determined by published reference equations.
Measurements and Main Results: The discovery meta-analyses identified one region on chromosome 15q25.1 meeting genome-wide significance in ever smokers that includes AGPHD1, IREB2, and CHRNA5/CHRNA3 genes. The region was also modestly associated among never smokers. Gene expression studies confirmed the presence of CHRNA5/3 in lung, airway smooth muscle, and bronchial epithelial cells. A single-nucleotide polymorphism in HTR4, a gene previously related to FEV1/FVC, achieved genome-wide statistical significance in combined meta-analysis. Top single-nucleotide polymorphisms in ADAM19, RARB, PPAP2B, and ADAMTS19 were nominally replicated in the COPD meta-analysis.
Conclusions: These results suggest an important role for the CHRNA5/3 region as a genetic risk factor for airflow obstruction that may be independent of smoking and implicate the HTR4 gene in the etiology of airflow obstruction.
chronic obstructive pulmonary disease; single-nucleotide polymorphism; genes
Chronic bronchitis (CB) is one of the classic phenotypes of COPD. The aims of our study were to investigate genetic variants associated with COPD subjects with CB relative to smokers with normal spirometry, and to assess for genetic differences between subjects with CB and without CB within the COPD population.
We analyzed data from current and former smokers from three cohorts: the COPDGene Study; GenKOLS (Bergen, Norway); and the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE). CB was defined as having a cough productive of phlegm on most days for at least 3 consecutive months per year for at least 2 consecutive years. CB COPD cases were defined as having both CB and at least moderate COPD based on spirometry. Our primary analysis used smokers with normal spirometry as controls; secondary analysis was performed using COPD subjects without CB as controls. Genotyping was performed on Illumina platforms; results were summarized using fixed-effect meta-analysis.
For CB COPD relative to smoking controls, we identified a new genome-wide significant locus on chromosome 11p15.5 (rs34391416, OR = 1.93, P = 4.99 × 10-8) as well as significant associations of known COPD SNPs within FAM13A. In addition, a GWAS of CB relative to those without CB within COPD subjects showed suggestive evidence for association on 1q23.3 (rs114931935, OR = 1.88, P = 4.99 × 10-7).
We found genome-wide significant associations with CB COPD on 4q22.1 (FAM13A) and 11p15.5 (EFCAB4A, CHID1 and AP2A2), and a locus associated with CB within COPD subjects on 1q23.3 (RPL31P11 and ATF6). This study provides further evidence that genetic variants may contribute to phenotypic heterogeneity of COPD.
ClinicalTrials.gov NCT00608764, NCT00292552
Electronic supplementary material
The online version of this article (doi:10.1186/s12931-014-0113-2) contains supplementary material, which is available to authorized users.
Pulmonary disease; Chronic obstructive; Chronic bronchitis; Genome-wide association study
Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02–1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72–1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.