Identifying COPD subtypes using multi-trait genetics.

Chronic Obstructive Pulmonary Disease (COPD) has a simple physiological diagnostic criterion but a wide range of clinical characteristics. The mechanisms underlying this variability in COPD phenotypes are unclear. To investigate the potential contribution of genetic variants to phenotypic heterogeneity, we examined the association of genome-wide associated lung function, COPD, and asthma variants with other phenotypes using phenome-wide association results derived in the UK Biobank. Our clustering analysis of the variants-phenotypes association matrix identified three clusters of genetic variants with different effects on white blood cell counts, height, and body mass index (BMI). To assess the potential clinical and molecular effects of these groups of variants, we investigated the association between cluster-specific genetic risk scores and phenotypes in the COPDGene cohort. We observed differences in steroid use, BMI, lymphocyte counts, chronic bronchitis, and differential gene and protein expression across the three genetic risk scores. Our results suggest that multi-phenotype analysis of obstructive lung disease-related risk variants may identify genetically driven phenotypic patterns in COPD.

Comprehensive drug response profiling and pan-omic analysis identified therapeutic candidates and prognostic biomarkers for Asian cholangiocarcinoma.

Cholangiocarcinoma (CCA) is rare cancer with the highest incidence in Eastern and Southeast Asian countries. Advanced CCA patients rely on chemotherapeutic regimens that offer unsatisfied clinical outcomes. We developed a comprehensive drug response profiling to investigate potential new drugs using CCA cell lines from Thai and Japanese patients against 100 approved anti-cancer drugs. We identified two major CCA subgroups that displayed unique molecular pathways from our integrative pan-omic and ligand-induced pathway activation analyses. MEK and Src inhibitors specifically killed the CCA1 subgroup without causing cytotoxicity to the normal cholangiocyte. Next, we developed the CCA45 signature to classify CCA patients based on their transcriptomic data. Our CCA45 signature could accurately predict prognosis, especially for Asian CCA patients. Our study provides a comprehensive public resource for drug repurposing in CCA and introduces analytical strategies for prioritizing cancer therapeutic agents for other rare cancer.

A genome-wide association study of bronchodilator response in participants of European and African ancestry from six independent cohorts.

Introduction: Bronchodilator response (BDR) is a measurement of acute bronchodilation in response to short-acting ß2-agonists, with a heritability between 10 and 40%. Identifying genetic variants associated with BDR may lead to a better understanding of its complex pathophysiology. Methods: We performed a genome-wide association study (GWAS) of BDR in six adult cohorts with participants of European ancestry (EA) and African ancestry (AA) including community cohorts and cohorts ascertained on the basis of obstructive pulmonary disease. Validation analysis was carried out in two paediatric asthma cohorts. Results: A total of 10 623 EA and 3597 AA participants were included in the analyses. No single nucleotide polymorphism (SNP) was associated with BDR at the conventional genome-wide significance threshold (p<5×10-8). Performing fine mapping and using a threshold of p<5×10-6 to identify suggestive variants of interest, we identified three SNPs with possible biological relevance: rs35870000 (within FREM1), which may be involved in IgE- and IL5-induced changes in airway smooth muscle cell responsiveness; rs10426116 (within ZNF284), a zinc finger protein, which has been implicated in asthma and BDR previously; and rs4782614 (near ATP2C2), involved in calcium transmembrane transport. Validation in paediatric cohorts yielded no significant SNPs, possibly due to age-genotype interaction effects. Conclusion: Ancestry-stratified and ancestry-combined GWAS meta-analyses of over 14 000 participants did not identify genetic variants associated with BDR at the genome-wide significance threshold, although a less stringent threshold identified three variants showing suggestive evidence of association. A common definition and protocol for measuring BDR in research may improve future efforts to identify variants associated with BDR.

Interaction of Cigarette Smoking and Polygenic Risk Score on Reduced Lung Function.

Importance: The risk of airflow limitation and chronic obstructive pulmonary disease (COPD) is influenced by combinations of cigarette smoking and genetic susceptibility, yet it remains unclear whether gene-by-smoking interactions are associated with quantitative measures of lung function. Objective: To assess the interaction of cigarette smoking and polygenic risk score in association with reduced lung function. Design, Setting, and Participants: This UK Biobank cohort study included UK citizens of European ancestry aged 40 to 69 years with genetic and spirometry data passing quality control metrics. Data was analyzed from July 2020 to March 2021. Exposures: PRS of combined forced expiratory volume in 1 second (FEV1) and percent of forced vital capacity exhaled in the first second (FEV1/FVC), self-reported pack-years of smoking, ever- vs never-smoking status, and current- vs former- or never-smoking status. Main Outcomes and Measures: FEV1/FVC was the primary outcome. Models were used to test for interactions with models, including the main effects of PRS, different smoking variables, and their cross-product terms. The association between pack-years of smoking and FEV1/FVC were compared for those in the highest vs lowest decile of estimated genetic risk for low lung function. Results: We included 319 730 individuals, of whom 24 915 (8%) had moderate-to-severe COPD cases, and 44.4% were men. Participants had a mean (SD) age 56.5 of (8.02) years. The PRS and pack-years were significantly associated with lower FEV1/FVC (PRS: ß, -0.03; 95% CI, -0.031 to -0.03; pack-years: ß, -0.0064; 95% CI, -0.0064 to -0.0063) and the interaction term (ß, -0.0028; 95% CI, -0.0029 to -0.0026). A stepwise increment in estimated effect sizes for these interaction terms was observed per 10 pack-years of smoking exposure. The interaction of PRS with 11 to 20, 31 to 40, and more than 50 pack-years categories were ß (interaction) -0.0038 (95% CI, -0.0046 to -0.0031); -0.013 (95% CI, -0.014 to -0.012); and -0.017 (95% CI, -0.019 to -0.016), respectively. There was evidence of significant interaction between PRS with ever- or never- smoking status (ß, interaction; -0.0064; 95% CI, -0.0068 to -0.0060) and current or not-current smoking (ß, interaction; -0.0091; 95% CI, -0.0097 to -0.0084). For any given level of pack-years of smoking exposure, FEV1/FVC was significantly lower for individuals in the tenth decile (ie, highest risk) than the first decile (ie, lowest risk) of genetic risk. For every 20 pack-years of smoking, those in the tenth decile compared with the first decile of genetic risk showed nearly a 2-fold reduction in FEV1/FVC. Conclusions and Relevance: COPD is characterized by diminished lung function, and our analyses suggest there is substantial interaction between genome-wide PRS and smoking exposures. While smoking was associated with decreased lung function across all genetic risk categories, the associations were strongest in individuals with higher estimated genetic risk.

A polygenic risk score for asthma in a large racially diverse population.

BACKGROUND: Polygenic risk scores (PRSs) will have important utility for asthma and other chronic diseases as a tool for predicting disease incidence and subphenotypes. OBJECTIVE: We utilized findings from a large multiancestry GWAS of asthma to compute a PRS for asthma with relevance for racially diverse populations. METHODS: We derived two PRSs for asthma using a standard approach (based on genome-wide significant variants) and a lasso sum regression approach (allowing all genetic variants to potentially contribute). We used data from the racially diverse Kaiser Permanente GERA cohort (68 638 non-Hispanic Whites, 5874 Hispanics, 6870 Asians and 2760 Blacks). Race was self-reported by questionnaire. RESULTS: For the standard PRS, non-Hispanic Whites showed the highest odds ratio for a standard deviation increase in PRS for asthma (OR = 1.16 (95% CI 1.14-1.18)). The standard PRS was also associated with asthma in Hispanic (OR = 1.12 (95% CI 1.05-1.19)) and Asian (OR = 1.10 (95% CI 1.04-1.17)) subjects, with a trend towards increased risk in Blacks (OR = 1.05 (95% CI 0.97-1.15)). We detected an interaction by sex, with men showing a higher risk of asthma with an increase in PRS as compared to women. The lasso sum regression-derived PRS showed stronger associations with asthma in non-Hispanic White subjects (OR = 1.20 (95% CI 1.18-1.23)), Hispanics (OR = 1.17 (95% 1.10-1.26)), Asians (OR = 1.18 (95% CI 1.10-1.27)) and Blacks (OR = 1.10 (95% CI 0.99-1.22)). CONCLUSION: Polygenic risk scores across multiple racial/ethnic groups were associated with increased asthma risk, suggesting that PRSs have potential as a tool for predicting disease development.

A systematic analysis of protein-altering exonic variants in chronic obstructive pulmonary disease.

Genome-wide association studies (GWASs) have identified regions associated with chronic obstructive pulmonary disease (COPD). GWASs of other diseases have shown an approximately 10-fold overrepresentation of nonsynonymous variants, despite limited exonic coverage on genotyping arrays. We hypothesized that a large-scale analysis of coding variants could discover novel genetic associations with COPD, including rare variants with large effect sizes. We performed a meta-analysis of exome arrays from 218,399 controls and 33,851 moderate-to-severe COPD cases. All exome-wide significant associations were present in regions previously identified by GWAS. We did not identify any novel rare coding variants with large effect sizes. Within GWAS regions on chromosomes 5q, 6p, and 15q, four coding variants were conditionally significant (P < 0.00015) when adjusting for lead GWAS single-nucleotide polymorphisms A common gasdermin B (GSDMB) splice variant (rs11078928) previously associated with a decreased risk for asthma was nominally associated with a decreased risk for COPD [minor allele frequency (MAF) = 0.46, P = 1.8e-4]. Two stop variants in coiled-coil α-helical rod protein 1 (CCHCR1), a gene involved in regulating cell proliferation, were associated with COPD (both P < 0.0001). The SERPINA1 Z allele was associated with a random-effects odds ratio of 1.43 for COPD (95% confidence interval = 1.17-1.74), though with marked heterogeneity across studies. Overall, COPD-associated exonic variants were identified in genes involved in DNA methylation, cell-matrix interactions, cell proliferation, and cell death. In conclusion, we performed the largest exome array meta-analysis of COPD to date and identified potential functional coding variants. Future studies are needed to identify rarer variants and further define the role of coding variants in COPD pathogenesis.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease.

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008-2010) and SpiroMeta Consortium (multiple countries, 1947-2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic ß4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

Relative contributions of family history and a polygenic risk score on COPD and related outcomes: COPDGene and ECLIPSE studies.

INTRODUCTION: Family history is a risk factor for chronic obstructive pulmonary disease (COPD). We previously developed a COPD risk score from genome-wide genetic markers (Polygenic Risk Score, PRS). Whether the PRS and family history provide complementary or redundant information for predicting COPD and related outcomes is unknown. METHODS: We assessed the predictive capacity of family history and PRS on COPD and COPD-related outcomes in non-Hispanic white (NHW) and African American (AA) subjects from COPDGene and ECLIPSE studies. We also performed interaction and mediation analyses. RESULTS: In COPDGene, family history and PRS were significantly associated with COPD in a single model (PFamHx <0.0001; PPRS<0.0001). Similar trends were seen in ECLIPSE. The area under the receiver operator characteristic curve for a model containing family history and PRS was significantly higher than a model with PRS (p=0.00035) in NHWs and a model with family history (p<0.0001) alone in NHWs and AAs. Both family history and PRS were significantly associated with measures of quantitative emphysema and airway thickness. There was a weakly positive interaction between family history and the PRS under the additive, but not multiplicative scale in NHWs (relative excess risk due to interaction=0.48, p=0.04). Mediation analyses found that a significant proportion of the effect of family history on COPD was mediated through PRS in NHWs (16.5%, 95% CI 9.4% to 24.3%), but not AAs. CONCLUSION: Family history and the PRS provide complementary information for predicting COPD and related outcomes. Future studies can address the impact of obtaining both measures in clinical practice.

Whole genome sequence analysis of pulmonary function and COPD in 19,996 multi-ethnic participants.

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry.

Chronic obstructive pulmonary disease and related phenotypes: polygenic risk scores in population-based and case-control cohorts.

BACKGROUND: Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes. METHODS: We constructed a polygenic risk score using a genome-wide association study of lung function (FEV1 and FEV1/forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV1/FVC <0·7 and FEV1 <80% of predicted). Associations were tested using logistic regression models, adjusting for age, sex, height, smoking pack-years, and principal components of genetic ancestry. We assessed predictive performance of models by area under the curve. In a subset of studies, we also studied quantitative and qualitative CT imaging phenotypes that reflect parenchymal and airway pathology, and patterns of reduced lung growth. FINDINGS: The polygenic risk score was associated with COPD in European (odds ratio [OR] per SD 1·81 [95% CI 1·74-1·88] and non-European (1·42 [1·34-1·51]) populations. Compared with the first decile, the tenth decile of the polygenic risk score was associated with COPD, with an OR of 7·99 (6·56-9·72) in European ancestry and 4·83 (3·45-6·77) in non-European ancestry cohorts. The polygenic risk score was superior to previously described genetic risk scores and, when combined with clinical risk factors (ie, age, sex, and smoking pack-years), showed improved prediction for COPD compared with a model comprising clinical risk factors alone (AUC 0·80 [0·79-0·81] vs 0·76 [0·75-0·76]). The polygenic risk score was associated with CT imaging phenotypes, including wall area percent, quantitative and qualitative measures of emphysema, local histogram emphysema patterns, and destructive emphysema subtypes. The polygenic risk score was associated with a reduced lung growth pattern. INTERPRETATION: A risk score comprised of genetic variants can identify a small subset of individuals at markedly increased risk for moderate-to-severe COPD, emphysema subtypes associated with cigarette smoking, and patterns of reduced lung growth. FUNDING: US National Institutes of Health, Wellcome Trust.

Genome-Wide Association Study of Susceptibility to Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a complex lung disease characterized by scarring of the lung that is believed to result from an atypical response to injury of the epithelium. Genome-wide association studies have reported signals of association implicating multiple pathways including host defense, telomere maintenance, signaling, and cell-cell adhesion.Objectives: To improve our understanding of factors that increase IPF susceptibility by identifying previously unreported genetic associations.Methods: We conducted genome-wide analyses across three independent studies and meta-analyzed these results to generate the largest genome-wide association study of IPF to date (2,668 IPF cases and 8,591 controls). We performed replication in two independent studies (1,456 IPF cases and 11,874 controls) and functional analyses (including statistical fine-mapping, investigations into gene expression, and testing for enrichment of IPF susceptibility signals in regulatory regions) to determine putatively causal genes. Polygenic risk scores were used to assess the collective effect of variants not reported as associated with IPF.Measurements and Main Results: We identified and replicated three new genome-wide significant (P < 5 × 10-8) signals of association with IPF susceptibility (associated with altered gene expression of KIF15, MAD1L1, and DEPTOR) and confirmed associations at 11 previously reported loci. Polygenic risk score analyses showed that the combined effect of many thousands of as yet unreported IPF susceptibility variants contribute to IPF susceptibility.Conclusions: The observation that decreased DEPTOR expression associates with increased susceptibility to IPF supports recent studies demonstrating the importance of mTOR signaling in lung fibrosis. New signals of association implicating KIF15 and MAD1L1 suggest a possible role of mitotic spindle-assembly genes in IPF susceptibility.

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene.

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38.

DSP variants may be associated with longitudinal change in quantitative emphysema.

BACKGROUND: Emphysema, characterized by lung destruction, is a key component of Chronic Obstructive Pulmonary Disease (COPD) and is associated with increased morbidity and mortality. Genome-wide association studies (GWAS) have identified multiple genetic factors associated with cross-sectional measures of quantitative emphysema, but the genetic determinants of longitudinal change in quantitative measures of emphysema remain largely unknown. Our study aims to identify genetic variants associated with longitudinal change in quantitative emphysema measured by computed tomography (CT) imaging. METHODS: We included current and ex-smokers from two longitudinal cohorts: COPDGene, a study of Non-Hispanic Whites (NHW) and African Americans (AA), and the Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points (ECLIPSE). We calculated annual change in two quantitative measures of emphysema based on chest CT imaging: percent low attenuation area (≤ - 950HU) (%LAA-950) and adjusted lung density (ALD). We conducted GWAS, separately in 3030 NHW and 1158 AA from COPDGene and 1397 Whites from ECLIPSE. We further explored effects of 360 previously reported variants and a lung function based polygenic risk score on annual change in quantitative emphysema. RESULTS: In the genome-wide association analysis, no variants achieved genome-wide significance (P < 5e-08). However, in the candidate region analysis, rs2076295 in the DSP gene, previously associated with COPD, lung function and idiopathic pulmonary fibrosis, was associated with change in %LAA-950 (ß (SE) = 0.09 (0.02), P = 3.79e-05) and in ALD (ß (SE) = - 0.06 (0.02), P = 2.88e-03). A lung function based polygenic risk score was associated with annual change in %LAA-950 (P = 4.03e-02) and with baseline measures of quantitative emphysema (P < 1e-03) and showed a trend toward association with annual change in ALD (P = 7.31e-02). CONCLUSIONS: DSP variants may be associated with longitudinal change in quantitative emphysema. Additional investigation of the DSP gene are likely to provide further insights into the disease progression in emphysema and COPD. TRIAL REGISTRATION: Clinicaltrials.gov Identifier: NCT00608764 , NCT00292552 .

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries.

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function-associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

Genetic landscape of chronic obstructive pulmonary disease identifies heterogeneous cell-type and phenotype associations.

Chronic obstructive pulmonary disease (COPD) is the leading cause of respiratory mortality worldwide. Genetic risk loci provide new insights into disease pathogenesis. We performed a genome-wide association study in 35,735 cases and 222,076 controls from the UK Biobank and additional studies from the International COPD Genetics Consortium. We identified 82 loci associated with P < 5 × 10-8; 47 of these were previously described in association with either COPD or population-based measures of lung function. Of the remaining 35 new loci, 13 were associated with lung function in 79,055 individuals from the SpiroMeta consortium. Using gene expression and regulation data, we identified functional enrichment of COPD risk loci in lung tissue, smooth muscle, and several lung cell types. We found 14 COPD loci shared with either asthma or pulmonary fibrosis. COPD genetic risk loci clustered into groups based on associations with quantitative imaging features and comorbidities. Our analyses provide further support for the genetic susceptibility and heterogeneity of COPD.

Genome-Wide Association Analysis of Single-Breath DlCO.

DlCO is a widely used pulmonary function test in clinical practice and a particularly useful measure for assessing patients with chronic obstructive pulmonary disease (COPD). We hypothesized that elucidating genetic determinants of DlCO could lead to better understanding of the genetic architecture of COPD. We estimated the heritability of DlCO using common genetic variants and performed genome-wide association analyses in four cohorts enriched for subjects with COPD (COPDGene [Genetic Epidemiology of COPD], NETT [National Emphysema Treatment Trial], GenKOLS [Genetics of Chronic Obstructive Lung Disease study], and TESRA [Treatment of Emphysema With a Gamma-Selective Retinoid Agonist study]) using a combined European ancestry white dataset and a COPDGene African American dataset. We assessed our genome-wide significant and suggestive associations for DlCO in previously reported genome-wide association studies of COPD and related traits. We also characterized associations of known COPD-associated variants and DlCO. We estimated the SNP-based heritability of DlCO in the European ancestry white population to be 22% (P = 0.0004). We identified three genome-wide significant associations with DlCO: variants near TGFB2, CHRNA3, and PDE11A loci (P < 5 × 10-8). In addition, 12 loci were suggestively associated with DlCO in European ancestry white (P < 1 × 10-5 in the combined analysis and P < 0.05 in both COPDGene and GenKOLS), including variants near NEGR1, CADM2, PCDH7, RETREG1, DACT2, NRG1, ANKRD18A, KRT86, NTN4, ARHGAP28, INSR, and PCBP3. Some DlCO-associated variants were also associated with COPD, emphysema, and/or spirometric values. Among 25 previously reported COPD loci, TGFB2, CHRNA3/CHRNA5, FAM13A, DSP, and CYP2A6 were associated with DlCO (P < 0.001). We identified several genetic loci that were significantly associated with DlCO and characterized effects of known COPD-associated loci on DlCO. These results could lead to better understanding of the heterogeneous nature of COPD.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution.

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4+, CD8+, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-ß signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

GWAS and systems biology analysis of depressive symptoms among smokers from the COPDGene cohort.

BACKGROUND: Large sample GWAS is needed to identify genetic factors associated with depression. This study used genome-wide genotypic and phenotypic data from the COPDGene study to identify genetic risk factors for depression. METHODS: Data were from 9716 COPDGene subjects with ≥10 pack-year history. Depression was defined as antidepressant use and/or a HADS depression subscale score ≥8. Non-Hispanic White (6576) and African-American (3140) subsets were analyzed. A GWAS pipeline identified SNPs associated with depression in each group. Network analysis software analyzed gene interactions through common biological pathways, genetic interactions, and tissue-specific gene expression. RESULTS: The mean age was 59.4 years (SD 9.0) with 46.5% female subjects. Depression was in 24.7% of the NHW group (1622) and 12.5% of the AA group (391). No SNPs had genome-wide significance. One of the top SNPs, rs12036147 (p = 1.28 × 10-6), is near CHRM3. Another SNP was near MDGA2 (rs17118176, p = 3.52 × 10-6). Top genes formed networks for synaptic transmission with a statistically significant level of more co-expression in brain than other tissues, particularly in the basal ganglia (p = 1.00 × 10-4). LIMITATIONS: Limitations included a depression definition based on antidepressant use and a limited HADS score subgroup, which could increase false negatives in depressed patients not on antidepressants. Antidepressants used for smoking cessation in non-depressed patients could lead to false positives. CONCLUSIONS: Systems biology analysis identified statistically significant pathways whereby multiple genes influence depression. The gene set pathway analysis and COPDGene data can help investigate depression in future studies.