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.

Two recent metaanalyses of genome-wide association studies conducted by the CHARGE and SpiroMeta consortia identified novel loci yielding evidence of association at or near genome-wide significance (GWS) with FEV1 and FEV1/FVC. We hypothesized that a subset of these markers would also be associated with chronic obstructive pulmonary disease (COPD) susceptibility. Thirty-two single-nucleotide polymorphisms (SNPs) in or near 17 genes in 11 previously identified GWS spirometric genomic regions were tested for association with COPD status in four COPD case-control study samples (NETT/NAS, the Norway case-control study, ECLIPSE, and the first 1,000 subjects in COPDGene; total sample size, 3,456 cases and 1,906 controls). In addition to testing the 32 spirometric GWS SNPs, we tested a dense panel of imputed HapMap2 SNP markers from the 17 genes located near the 32 GWS SNPs and in a set of 21 well studied COPD candidate genes. Of the previously identified GWS spirometric genomic regions, three loci harbored SNPs associated with COPD susceptibility at a 5% false discovery rate: the 4q24 locus including FLJ20184/INTS12/GSTCD/NPNT, the 6p21 locus including AGER and PPT2, and the 5q33 locus including ADAM19. In conclusion, markers previously associated at or near GWS with spirometric measures were tested for association with COPD status in data from four COPD case-control studies, and three loci showed evidence of association with COPD susceptibility at a 5% false discovery rate.

Rationale

Traditional genome-wide association studies (GWAS) of large cohort of subjects with chronic obstructive pulmonary disease (COPD) have successfully identified novel candidate genes, but several other plausible loci do not meet strict criteria for genome-wide significance after correction for multiple testing.

Objectives

We hypothesize that by applying unbiased weights derived from unique populations we can identify additional COPD susceptibility loci.

Methods

We performed a homozygosity haplotype analysis on a group of subjects with and without COPD to identify regions of conserved homozygosity (RCHH). Weights were constructed based on the frequency of these RCHH in case vs. controls, and used to adjust the P values from a large collaborative GWAS of COPD.

Results

We identified 2,318 regions of conserved homozygosity, of which 576 were significantly (P < .05) overrepresented in cases. After applying the weights constructed from these regions to a collaborative GWAS of COPD, we identified two single nucleotide polymorphisms in a novel gene (FGF7) that gained genome-wide significance by the false discovery rate method. In a follow-up analysis, both SNPs (rs12591300 and rs4480740) were significantly associated with COPD in an independent population (combined P values of 7.9E-07 and 2.8E-06 respectively). In another independent population, increased lung tissue FGF7 expression was associated with worse measures of lung function.

Conclusion

Weights constructed from a homozygosity haplotype analysis of an isolated population successfully identify novel genetic associations from a GWAS on a separate population. This method can be used to identify promising candidate genes that fail to meet strict correction for multiple testing.

Background

Chronic Obstructive Pulmonary Disease (COPD) is defined by post-bronchodilator spirometry. Data on “normal values” come predominantly from pre-bronchodilator spirometry. The effects of this on diagnosis are unknown.

Methods

Lower limits of normal (LLN) were estimated from “normal” participants in the Burden of Obstructive Lung Disease (BOLD) programme. Values separately derived using pre- and post-bronchodilator spirometry were compared. Sensitivity and specificity of criteria derived from pre-bronchodilator spirometry and pre-bronchodilator spirometry adjusted by a constant were assessed in the remaining population. The “gold standard” was the LLN for the post-bronchodilator spirometry in the “normal population”. For FEV1/FVC, sensitivity and specificity of criteria were also assessed when a fixed value of < 70% was used rather than LLN.

Results

Of 6,600 participants with full data, 1,354 were defined as “normal”. Mean differences between pre- and post- bronchodilator measurements were small and the Bland-Altman plots showed no association between difference and mean value. Compared with using the gold standard, however, tests using pre-bronchodilator spirometry had a sensitivity and specificity of detecting a low FEV1 of 78.4% and 100%, a low FVC of 99.8% and 99.1% and a low FEV1/FVC ratio of 65% and 100%. Adjusting this by a constant improved the sensitivity without substantially altering the specificity for FEV1 (99%, 99.8%), FVC (97.4%, 99.9%) and FEV1/FVC (98.7%, 99.5%).

Conclusions

Using pre-bronchodilator spirometry to derive norms for lung function reduces sensitivity compared to a post-bronchodilator gold standard. Adjustment of these values by a constant can improve validity of the test.

Background

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.

Objective

To identify SNPs associated with lifetime average and current CPD, age at smoking initiation, and smoking cessation in COPD subjects.

Methods

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.

Results

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.

Conclusion

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.

Cachexia, whether assessed by body mass index (BMI) or fat-free mass index (FFMI), affects a significant proportion of patients with chronic obstructive pulmonary disease (COPD), and is an independent risk factor for increased mortality, increased emphysema, and more severe airflow obstruction. The variable development of cachexia among patients with COPD suggests a role for genetic susceptibility. The objective of the present study was to determine genetic susceptibility loci involved in the development of low BMI and FFMI in subjects with COPD. A genome-wide association study (GWAS) of BMI was conducted in three independent cohorts of European descent with Global Initiative for Chronic Obstructive Lung Disease stage II or higher COPD: Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-Points (ECLIPSE; n = 1,734); Norway-Bergen cohort (n = 851); and a subset of subjects from the National Emphysema Treatment Trial (NETT; n = 365). A genome-wide association of FFMI was conducted in two of the cohorts (ECLIPSE and Norway). In the combined analyses, a significant association was found between rs8050136, located in the first intron of the fat mass and obesity–associated (FTO) gene, and BMI (P = 4.97 × 10−7) and FFMI (P = 1.19 × 10−7). We replicated the association in a fourth, independent cohort consisting of 502 subjects with COPD from COPDGene (P = 6 × 10−3). Within the largest contributing cohort of our analysis, lung function, as assessed by forced expiratory volume at 1 second, varied significantly by FTO genotype. Our analysis suggests a potential role for the FTO locus in the determination of anthropomorphic measures associated with COPD.

Chronic obstructive pulmonary disease (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.

Rationale: Chronic obstructive pulmonary disease (COPD), characterized by airflow limitation, is a disorder with high phenotypic and genetic heterogeneity. Pulmonary emphysema is a major but variable component of COPD; familial data suggest that different components of COPD, such as emphysema, may be influenced by specific genetic factors.

Objectives: To identify genetic determinants of emphysema assessed through high-resolution chest computed tomography in individuals with COPD.

Methods: We performed a genome-wide association study (GWAS) of emphysema determined from chest computed tomography scans with a total of 2,380 individuals with COPD in three independent cohorts of white individuals from (1) a cohort from Bergen, Norway, (2) the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) Study, and (3) the National Emphysema Treatment Trial (NETT). We tested single-nucleotide polymorphism associations with the presence or absence of emphysema determined by radiologist assessment in two of the three cohorts and a quantitative emphysema trait (percentage of lung voxels less than –950 Hounsfield units) in all three cohorts.

Measurements and Main Results: We identified association of a single-nucleotide polymorphism in BICD1 with the presence or absence of emphysema (P = 5.2 × 10−7 with at least mild emphysema vs. control subjects; P = 4.8 × 10−8 with moderate and more severe emphysema vs. control subjects).

Conclusions: Our study suggests that genetic variants in BICD1 are associated with qualitative emphysema in COPD. Variants in BICD1 are associated with length of telomeres, which suggests that a mechanism linked to accelerated aging may be involved in the pathogenesis of emphysema.

Clinical trial registered with www.clinicaltrials.gov (NCT00292552).

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.

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.

The objective of the present study was to determine the association between CT phenotypes—emphysema by low attenuation area and bronchitis by airway wall thickness—and body composition parameters in a large cohort of subjects with and without COPD. In 452 COPD subjects and 459 subjects without COPD, CT scans were performed to determine emphysema (%LAA), airway wall thickness (AWT-Pi10), and lung mass. Muscle wasting based on FFMI was assessed by bioelectrical impedance. In both the men and women with COPD, FFMI was negatively associated with %LAA. FMI was positively associated with AWT-Pi10 in both subjects with and without COPD. Among the subjects with muscle wasting, the percentage emphysema was high, but the predictive value was moderate. In conclusion, the present study strengthens the hypothesis that the subgroup of COPD cases with muscle wasting have emphysema. Airway wall thickness is positively associated with fat mass index in both subjects with and without COPD.

Mortality statistics represent important endpoints in epidemiological studies. The diagnostic validity of cerebral stroke and ischemic heart disease recorded as the underlying cause of death in Norwegian mortality statistics was assessed by using mortality data of participants in the Bergen Clinical Blood Pressure Study in Norway and autopsy records from the Gade Institute in Bergen. In the 41 years of the study (1965–2005) 4,387 subjects had died and 1,140 (26%) had undergone a post mortem examination; 548 (12%) died from cerebral stroke and 1,120 (24%) from ischemic heart disease according to the mortality statistics, compared to 113 (10%) strokes and 323 (28%) coronary events registered in the autopsy records. The sensitivity and positive predictive value of fatal cerebral strokes in the mortality statistics were 0.75, 95% confidence interval (CI) [0.66, 0.83] and 0.86 [0.77, 0.92], respectively, whereas those of coronary deaths were 0.87 [0.84, 0.91] and 0.85 [0.81, 0.89] respectively. Cohen’s Kappa coefficients were 0.78 [0.72, 0.84] for stroke and 0.80 [0.76, 0.84] for coronary deaths. In addition to female gender and increasing age at death, cerebral stroke was a negative predictor of an autopsy being carried out (odds ratio (OR) 0.69, 95% CI [0.54, 0.87]), whereas death from coronary heart disease was not (OR 1.14, 95% CI [0.97, 1,33]), both adjusted for gender and age at death. There was substantial agreement between mortality statistics and autopsy findings for both fatal strokes and coronary deaths. Selection for post mortem examinations was associated with age, gender and cause of death.

Background

International guidelines recommend that pulmonary reference populations consist of never‐smokers without respiratory diseases or symptoms, but the diseases and symptoms are not clearly specified. The present study aimed to identify simple exclusion criteria for defining pulmonary reference populations.

Methods

Based on a random sample from a general population (the parent population), 2358 subjects aged 26–82 years performed spirometric tests. From this sample, subjects were stepwise excluded according to self‐reported obstructive lung diseases, symptoms and smoking history. Four increasingly more healthy respiratory reference populations were formed. Prediction equations for the median and lower limit of normal lung function were derived using quantile regression analysis.

Results

Subjects without self‐reported obstructive lung diseases or the cardinal respiratory symptoms of breathlessness, cough or wheeze (population B), never‐smokers without cardinal symptoms (population C) and never‐smokers without any respiratory symptoms (population D) constituted 50% (n = 1184), 23% (n = 539) and 14% (n = 331) of the parent population (population A), respectively. The largest discrepancy between prediction equations was found between the parent population and the population without cardinal respiratory symptoms (population B) (p<0.05). Minor changes in the reference equations were also seen when excluding ever‐smokers (population C). There was no additional change with exclusion of other respiratory symptoms (population D). Age‐related decline in lung function was steepest in the parent population.

Conclusions

Obstructive lung diseases, smoking history, breathlessness, cough and wheeze are optimal exclusion criteria for a pulmonary reference population. Further validation of the exclusion criteria identified in this study is recommended with identical wording in other and larger multinational populations.

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).

There is considerable variability in the susceptibility of smokers to develop chronic obstructive pulmonary disease (COPD). The only known genetic risk factor is severe deficiency of α1-antitrypsin, which is present in 1–2% of individuals with COPD. We conducted a genome-wide association study (GWAS) in a homogenous case-control cohort from Bergen, Norway (823 COPD cases and 810 smoking controls) and evaluated the top 100 single nucleotide polymorphisms (SNPs) in the family-based International COPD Genetics Network (ICGN; 1891 Caucasian individuals from 606 pedigrees) study. The polymorphisms that showed replication were further evaluated in 389 subjects from the US National Emphysema Treatment Trial (NETT) and 472 controls from the Normative Aging Study (NAS) and then in a fourth cohort of 949 individuals from 127 extended pedigrees from the Boston Early-Onset COPD population. Logistic regression models with adjustments of covariates were used to analyze the case-control populations. Family-based association analyses were conducted for a diagnosis of COPD and lung function in the family populations. Two SNPs at the α-nicotinic acetylcholine receptor (CHRNA 3/5) locus were identified in the genome-wide association study. They showed unambiguous replication in the ICGN family-based analysis and in the NETT case-control analysis with combined p-values of 1.48×10−10, (rs8034191) and 5.74×10−10 (rs1051730). Furthermore, these SNPs were significantly associated with lung function in both the ICGN and Boston Early-Onset COPD populations. The C allele of the rs8034191 SNP was estimated to have a population attributable risk for COPD of 12.2%. The association of hedgehog interacting protein (HHIP) locus on chromosome 4 was also consistently replicated, but did not reach genome-wide significance levels. Genome-wide significant association of the HHIP locus with lung function was identified in the Framingham Heart study (Wilk et al., companion article in this issue of PLoS Genetics; doi:10.1371/journal.pgen.1000429). The CHRNA 3/5 and the HHIP loci make a significant contribution to the risk of COPD. CHRNA3/5 is the same locus that has been implicated in the risk of lung cancer.

Author Summary

There is considerable variability in the susceptibility of smokers to develop chronic obstructive pulmonary disease (COPD), which is a heritable multi-factorial trait. Identifying the genetic determinants of COPD risk will have tremendous public health importance. This study describes the first genome-wide association study (GWAS) in COPD. We conducted a GWAS in a homogenous case-control cohort from Norway and evaluated the top 100 single nucleotide polymorphisms in the family-based International COPD Genetics Network. The polymorphisms that showed replication were further evaluated in subjects from the US National Emphysema Treatment Trial and controls from the Normative Aging Study and then in a fourth cohort of extended pedigrees from the Boston Early-Onset COPD population. Two polymorphisms in the α-nicotinic acetylcholine receptor 3/5 locus on chromosome 15 showed unambiguous evidence of association with COPD. This locus has previously been implicated in both smoking behavior and risk of lung cancer, suggesting the possibility of multiple functional polymorphisms in the region or a single polymorphism with wide phenotypic consequences. The hedgehog interacting protein (HHIP) locus on chromosome 4, which is associated with COPD, is also a significant risk locus for COPD.

Background

To our knowledge, no studies of the possible association of early life environment with snoring in adulthood have been published. We aimed to investigate whether early life environment is associated with snoring later in life.

Methods

A questionnaire including snoring frequency in adulthood and environmental factors in early life was obtained from 16,190 randomly selected men and women, aged 25–54 years, in Sweden, Norway, Iceland, Denmark and Estonia (response rate 74%).

Results

A total of 15,556 subjects answered the questions on snoring. Habitual snoring, defined as loud and disturbing snoring at least 3 nights a week, was reported by 18%. Being hospitalized for a respiratory infection before the age of two years (adjusted odds ratio (OR) = 1.27; 95% confidence interval (CI) 1.01–1.59), suffering from recurrent otitis as a child (OR = 1.18; 95%CI 1.05–1.33), growing up in a large family (OR = 1.04; 95%CI 1.002–1.07) and being exposed to a dog at home as a newborn (OR = 1.26; 95%CI 1.12–1.42) were independently related to snoring later in life and independent of a number of possible confounders in adulthood. The same childhood environmental factors except household size were also related with snoring and daytime sleepiness combined.

Conclusion

The predisposition for adult snoring may be partly established early in life. Having had severe airway infections or recurrent otitis in childhood, being exposed to a dog as a newborn and growing up in a large family are environmental factors associated with snoring in adulthood.