Refining susceptibility loci of chronic obstructive pulmonary disease with lung eqtls

Breathing new life into the study of COPD with genes identified from genome-wide association studies

COPD is a major cause of morbidity and mortality globally. While the significance of environmental exposures in disease pathogenesis is well established, the functional contribution of genetic factors has only in recent years drawn attention. Notably, many genes associated with COPD risk are also linked with lung function. Because reduced lung function precedes COPD onset, this association is consistent with the possibility that derangements leading to COPD could arise during lung development. In this review, we summarise the role of leading genes (HHIP, FAM13A, DSP, AGER and TGFB2) identified by genome-wide association studies in lung development and COPD. Because many COPD genome-wide association study genes are enriched in lung epithelial cells, we focus on the role of these genes in the lung epithelium in development, homeostasis and injury.

AKT Phosphorylates FAM13A and Promotes Its Degradation via CUL4A/DDB1/DCAF1 E3 Complex

SNPs within FAM13A (family with sequence similarity 13 member A) gene are significantly associated with chronic obstructive pulmonary disease and lung function in genome-wide association studies (GWAS). However, how FAM13A protein is regulated under physiological and pathological conditions remains largely elusive. Herein, we report that FAM13A is phosphorylated at the serine 312 residue by AKT kinase after cigarette smoke extract treatment and thereby recognized by the CULLIN4A/DCAF1 (DDB1 and CUL4 associated factor 1) E3 ligase complex, rendering the ubiquitination-mediated degradation of FAM13A. More broadly, downregulation of FAM13A protein upon AKT activation, as a general cellular response to acute stress, was also detected in influenza- or naphthalene-injured lungs in mice. Functionally, reduced protein levels of FAM13A lead to accelerated epithelial cell proliferation in murine lungs during the recovery phase after injury. In summary, we characterized a novel molecular mechanism that regulates the stability of FAM13A protein, which enables the fine-tuning of lung epithelial repair after injury. These significant findings will expand our molecular understanding of the regulation of protein stability, which may modulate lung epithelial repair implicated in the development of chronic obstructive pulmonary disease and other lung diseases.

Assessing the contribution of rare genetic variants to phenotypes of chronic obstructive pulmonary disease using whole-genome sequence data

RATIONALE

Genetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear.

OBJECTIVES

We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.

METHODS

Using the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants.

MEASUREMENTS AND MAIN RESULTS

In European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability.

CONCLUSIONS

Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

FAM13A as potential therapeutic target in modulating TGF-β-induced airway tissue remodeling in COPD

Genome-wide association studies have shown that a gene variant in the Family with sequence similarity 13, member A (FAM13A) is strongly associated with reduced lung function and the appearance of respiratory symptoms in patients with chronic obstructive pulmonary disease (COPD). A key player in smoking-induced tissue injury and airway remodeling is the transforming growth factor-β1 (TGF-β1). To determine the role of FAM13A in TGF-β1 signaling, FAM13A^-/- airway epithelial cells were generated using CRISPR-Cas9, whereas overexpression of FAM13A was achieved using lipid nanoparticles. Wild-type (WT) and FAM13A^-/- cells were treated with TGF-β1, followed by gene and/or protein expression analyses. FAM13A^-/- cells augmented TGF-β1-induced increase in collagen type 1 (COL1A1), matrix metalloproteinase 2 (MMP2), expression compared with WT cells. This effect was mediated by an increase in β-catenin (CTNNB1) expression in FAM13A^-/- cells compared with WT cells after TGF-β1 treatment. FAM13A overexpression was partially protective from TGF-β1-induced COL1A1 expression. Finally, we showed that airway epithelial-specific FAM13A protein expression is significantly increased in patients with severe COPD compared with control nonsmokers, and negatively correlated with lung function. In contrast, β-catenin (CTNNB1), which has previously been linked to be regulated by FAM13A, is decreased in the airway epithelium of smokers with COPD compared with non-COPD subjects. Together, our data showed that FAM13A may be protective from TGF-β1-induced fibrotic response in the airway epithelium via sequestering CTNNB1 from its regulation on downstream targets. Therapeutic increase in FAM13A expression in the airway epithelium of smokers at risk for COPD, and those with mild COPD, may reduce the extent of airway tissue remodeling.

Variants associated with HHIP expression have sex-differential effects on lung function

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

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.

The SNP rs13147758 in the HHIP Gene Is Associated With COPD Susceptibility, Serum, and Sputum Protein Levels in Smokers

BACKGROUND

Genetic association studies have identified single nucleotide polymorphisms (SNPs) related to chronic obstructive pulmonary disease (COPD) susceptibility. The aim of this study was to identify HHIP genetic variants associated with COPD, pulmonary function, and serum and sputum HHIP protein levels in Mexican mestizo smokers.

MATERIALS AND METHODS

Association analysis was performed by carrying out a case-control study in Mexican mestizo smokers comprised of two groups: tobacco-smoking subjects with COPD (COPD-TS, n = 222) and smokers without COPD (SWOC, n = 333). We evaluated three SNPs (rs13147758, rs1828591, and rs13118928) in the HHIP gene. Allele discrimination was accomplished by qPCR using TaqMan probes, and determination of protein levels in the serum and sputum supernatants (SS) was performed using ELISA.

RESULTS

Statistically significant differences were observed in the rs13147758 GG genotype (adjusted p = 0.014, OR = 1.95) and the rs13147758-rs1828591 GA haplotype (p = 6.6E-06, OR = 2.65) in the case-control comparison. HHIP protein levels were elevated in SS samples from the COPD-TS group compared to those from the SWOC group (p = 0.03). Based on genotype analysis, HHIP protein levels were lower in the serum samples of rs13147758 GG genotype carriers in the COPD-TS group than in the serum samples of rs13147758 GG genotype carriers from the SWOC group (p < 0.05), but there were no differences in the sputum samples.

CONCLUSION

The rs13147758 GG genotype and the rs13147758-rs1828591 GA haplotype are associated with susceptibility to COPD. Furthermore, an association in protein levels was observed between the HHIP rs13147758 genotype and COPD in Mexican mestizo smokers.

Variants associated with HHIP expression have sex-differential effects on lung function

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

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.

Why is Disease Penetration So Variable? Role of Genetic Modifiers of Lung Function in Alpha-1 Antitrypsin Deficiency

Individuals with alpha-1 antitrypsin deficiency (AATD) have marked heterogeneity in lung function, suspected to be related to a combination of both environmental (e.g., cigarette smoking) and genetic factors. Lung function is heritable in the general population and in persons with severe AATD. Several genetic modifiers of lung function in persons with AATD have been described; however, replication is lacking. A genome-wide association study (GWAS) of lung function in persons with AATD has yet to be performed and may inform whether genetic determinants of lung function are overlapping in persons with AATD and in the general population. As GWASs require large sample sizes for adequate power, genetic risk scores offer an alternate approach to assess the overlap of genetic determinants of lung function in the general population in persons with AATD. Where GWASs are limited to common genetic variant discovery, whole genome sequencing (for rare variant discovery) and integrative genomic studies (examining the influence of genetic variants on gene, protein, and metabolite levels) offer potential for an expanded discovery of genetic modifiers of lung function in AATD. In the following review we examine past descriptions of genetic modifiers of lung function in AATD and describe a path forward to further investigate and define the likely genetic modifiers of lung function in AATD.

Genetic effects of welding fumes on the development of respiratory system diseases

BACKGROUND

The welding process releases potentially hazardous gases and fumes, mainly composed of metallic oxides, fluorides and silicates. Long term welding fume (WF) inhalation is a recognized health issue that carries a risk of developing chronic health problems, particularly respiratory system diseases (RSDs). Aside from general airway irritation, WF exposure may drive direct cellular responses in the respiratory system which increase risk of RSD, but these are not well understood.

METHODS

We developed a quantitative framework to identify gene expression effects of WF exposure that may affect RSD development. We analyzed gene expression microarray data from WF-exposed tissues and RSD-affected tissues, including chronic bronchitis (CB), asthma (AS), pulmonary edema (PE), lung cancer (LC) datasets. We built disease-gene (diseasome) association networks and identified dysregulated signaling and ontological pathways, and protein-protein interaction sub-network using neighborhood-based benchmarking and multilayer network topology.

RESULTS

We observed many genes with altered expression in WF-exposed tissues were also among differentially expressed genes (DEGs) in RSD tissues; for CB, AS, PE and LC there were 34, 27, 50 and 26 genes respectively. DEG analysis, using disease association networks, pathways, ontological analysis and protein-protein interaction sub-network suggest significant links between WF exposure and the development of CB, AS, PE and LC.

CONCLUSIONS

Our network-based analysis and investigation of the genetic links of WFs and RSDs confirm a number of genes and gene products are plausible participants in RSD development. Our results are a significant resource to identify causal influences on the development of RSDs, particularly in the context of WF exposure.

Mixed graphical models for integrative causal analysis with application to chronic lung disease diagnosis and prognosis

MOTIVATION

Integration of data from different modalities is a necessary step for multi-scale data analysis in many fields, including biomedical research and systems biology. Directed graphical models offer an attractive tool for this problem because they can represent both the complex, multivariate probability distributions and the causal pathways influencing the system. Graphical models learned from biomedical data can be used for classification, biomarker selection and functional analysis, while revealing the underlying network structure and thus allowing for arbitrary likelihood queries over the data.

RESULTS

In this paper, we present and test new methods for finding directed graphs over mixed data types (continuous and discrete variables). We used this new algorithm, CausalMGM, to identify variables directly linked to disease diagnosis and progression in various multi-modal datasets, including clinical datasets from chronic obstructive pulmonary disease (COPD). COPD is the third leading cause of death and a major cause of disability and thus determining the factors that cause longitudinal lung function decline is very important. Applied on a COPD dataset, mixed graphical models were able to confirm and extend previously described causal effects and provide new insights on the factors that potentially affect the longitudinal lung function decline of COPD patients.

AVAILABILITY AND IMPLEMENTATION

The CausalMGM package is available on http://www.causalmgm.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Hedgehog signaling in the airway epithelium of patients with chronic obstructive pulmonary disease

Genome-wide association studies have linked gene variants of the receptor patched homolog 1 (PTCH1) with chronic obstructive pulmonary disease (COPD). However, its biological role in the disease is unclear. Our objective was to determine the expression pattern and biological role of PTCH1 in the lungs of patients with COPD. Airway epithelial-specific PTCH1 protein expression and epithelial morphology were assessed in lung tissues of control and COPD patients. PTCH1 mRNA expression was measured in bronchial epithelial cells obtained from individuals with and without COPD. The effects of PTCH1 siRNA knockdown on epithelial repair and mucous expression were evaluated using human epithelial cell lines. Ptch1^+/− mice were used to assess the effect of decreased PTCH1 on mucous expression and airway epithelial phenotypes. Airway epithelial-specific PTCH1 protein expression was significantly increased in subjects with COPD compared to controls, and its expression was associated with total airway epithelial cell count and thickness. PTCH1 knockdown attenuated wound closure and mucous expression in airway epithelial cell lines. Ptch1^+/− mice had reduced mucous expression compared to wildtype mice following mucous induction. PTCH1 protein is up-regulated in COPD airway epithelium and may upregulate mucous expression. PTCH1 provides a novel target to reduce chronic bronchitis in COPD patients.

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.

Leveraging lung tissue transcriptome to uncover candidate causal genes in COPD genetic associations

Causal genes of chronic obstructive pulmonary disease (COPD) remain elusive. The current study aims at integrating genome-wide association studies (GWAS) and lung expression quantitative trait loci (eQTL) data to map COPD candidate causal genes and gain biological insights into the recently discovered COPD susceptibility loci. Two complementary genomic datasets on COPD were studied. First, the lung eQTL dataset which included whole-genome gene expression and genotyping data from 1038 individuals. Second, the largest COPD GWAS to date from the International COPD Genetics Consortium (ICGC) with 13 710 cases and 38 062 controls. Methods that integrated GWAS with eQTL signals including transcriptome-wide association study (TWAS), colocalization and Mendelian randomization-based (SMR) approaches were used to map causality genes, i.e. genes with the strongest evidence of being the functional effector at specific loci. These methods were applied at the genome-wide level and at COPD risk loci derived from the GWAS literature. Replication was performed using lung data from GTEx. We collated 129 non-overlapping risk loci for COPD from the GWAS literature. At the genome-wide scale, 12 new COPD candidate genes/loci were revealed and six replicated in GTEx including CAMK2A, DMPK, MYO15A, TNFRSF10A, BTN3A2 and TRBV30. In addition, we mapped candidate causal genes for 60 out of the 129 GWAS-nominated loci and 23 of them were replicated in GTEx. Mapping candidate causal genes in lung tissue represents an important contribution to the genetics of COPD, enriches our biological interpretation of GWAS findings, and brings us closer to clinical translation of genetic associations.

Trade-offs in aging lung diseases: a review on shared but opposite genetic risk variants in idiopathic pulmonary fibrosis, lung cancer and chronic obstructive pulmonary disease

PURPOSE OF REVIEW

The process of aging involves biological changes that increases susceptibility for disease. In the aging lung disease IPF, GWAS studies identified genes associated with risk for disease. Recently, several of these genes were also found to be involved in risk for COPD or lung cancer. This review describes GWAS-derived risk genes for IPF that overlap with risk genes for lung cancer or COPD.

RECENT FINDINGS

Risk genes that overlap between aging lung diseases, include FAM13A, DSP and TERT. Most interestingly, disease predisposing alleles for IPF are opposite to those for COPD or lung cancer. Studies show that the alleles are associated with differential gene expression and with physiological traits in the general population. The opposite allelic effect sizes suggest the presence of trade-offs in the aging lung. For TERT, the trade-off involves cellular senescence versus proliferation and repair. For FAM13A and DSP, trade-offs may involve protection from noxious gases or tissue integrity.

SUMMARY

The overlap in risk genes in aging lung diseases provides evidence that processes associated with FAM13A, DSP and TERT are important for healthy aging. The opposite effect size of the disease risk alleles may represent trade-offs, for which a model involving an apicobasal gene expression gradient is presented.

Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study

BACKGROUND

Few genetic studies that focus on moderate-to-severe asthma exist. We aimed to identity novel genetic variants associated with moderate-to-severe asthma, see whether previously identified genetic variants for all types of asthma contribute to moderate-to-severe asthma, and provide novel mechanistic insights using expression analyses in patients with asthma.

METHODS

In this genome-wide association study, we used a two-stage case-control design. In stage 1, we genotyped patient-level data from two UK cohorts (the Genetics of Asthma Severity and Phenotypes [GASP] initiative and the Unbiased BIOmarkers in PREDiction of respiratory disease outcomes [U-BIOPRED] project) and used data from the UK Biobank to collect patient-level genomic data for cases and controls of European ancestry in a 1:5 ratio. Cases were defined as having moderate-to-severe asthma if they were taking appropriate medication or had been diagnosed by a doctor. Controls were defined as not having asthma, rhinitis, eczema, allergy, emphysema, or chronic bronchitis as diagnosed by a doctor. For stage 2, an independent cohort of cases and controls (1:5) was selected from the UK Biobank only, with no overlap with stage 1 samples. In stage 1 we undertook a genome-wide association study of moderate-to-severe asthma, and in stage 2 we followed up independent variants that reached the significance threshold of p less than 1 × 10^-6 in stage 1. We set genome-wide significance at p less than 5 × 10^-8. For novel signals, we investigated their effect on all types of asthma (mild, moderate, and severe). For all signals meeting genome-wide significance, we investigated their effect on gene expression in patients with asthma and controls.

FINDINGS

We included 5135 cases and 25 675 controls for stage 1, and 5414 cases and 21 471 controls for stage 2. We identified 24 genome-wide significant signals of association with moderate-to-severe asthma, including several signals in innate or adaptive immune-response genes. Three novel signals were identified: rs10905284 in GATA3 (coded allele A, odds ratio [OR] 0·90, 95% CI 0·88-0·93; p=1·76 × 10^-10), rs11603634 in the MUC5AC region (coded allele G, OR 1·09, 1·06-1·12; p=2·32 × 10^-8), and rs560026225 near KIAA1109 (coded allele GATT, OR 1·12, 1·08-1·16; p=3·06 × 10^-9). The MUC5AC signal was not associated with asthma when analyses included mild asthma. The rs11603634 G allele was associated with increased expression of MUC5AC mRNA in bronchial epithelial brush samples via proxy SNP rs11602802; (p=2·50 × 10^-5) and MUC5AC mRNA was increased in bronchial epithelial samples from patients with severe asthma (in two independent analyses, p=0·039 and p=0·022).

INTERPRETATION

We found substantial shared genetic architecture between mild and moderate-to-severe asthma. We also report for the first time genetic variants associated with the risk of developing moderate-to-severe asthma that regulate mucin production. Finally, we identify candidate causal genes in these loci and provide increased insight into this difficult to treat population.

FUNDING

Asthma UK, AirPROM, U-BIOPRED, UK Medical Research Council, and Rosetrees Trust.

The Overlap of Lung Tissue Transcriptome of Smoke Exposed Mice with Human Smoking and COPD

Genome-wide mRNA profiling in lung tissue from human and animal models can provide novel insights into the pathogenesis of chronic obstructive pulmonary disease (COPD). While 6 months of smoke exposure are widely used, shorter durations were also reported. The overlap of short term and long-term smoke exposure in mice is currently not well understood, and their representation of the human condition is uncertain. Lung tissue gene expression profiles of six murine smoking experiments (n = 48) were obtained from the Gene Expression Omnibus (GEO) and analyzed to identify the murine smoking signature. The "human smoking" gene signature containing 386 genes was previously published in the lung eQTL study (n = 1,111). A signature of mild COPD containing 7 genes was also identified in the same study. The lung tissue gene signature of "severe COPD" (n = 70) contained 4,071 genes and was previously published. We detected 3,723 differentially expressed genes in the 6 month-exposure mice datasets (FDR <0.1). Of those, 184 genes (representing 48% of human smoking) and 1,003 (representing 27% of human COPD) were shared with the human smoking-related genes and the COPD severity-related genes, respectively. There was 4-fold over-representation of human and murine smoking-related genes (P = 6.7 × 10^-26) and a 1.4 fold in the severe COPD -related genes (P = 2.3 × 10^-12). There was no significant enrichment of the mice and human smoking-related genes in mild COPD signature. These data suggest that murine smoke models are strongly representative of molecular processes of human smoking but less of COPD.

Correlation between lung cancer and the HHIP polymorphisms of chronic obstructive pulmonary disease (COPD) in the Chinese Han population

To further investigate the relationship between lung cancer and hedgehog interacting protein (HHIP) polymorphisms of chronic obstructive pulmonary disease (COPD) patients, we conducted a case-control study in a Chinese Han population. Six HHIP SNPs with minor allele frequencies >5% (rs1489758, rs1489759, rs10519717, rs13131837, rs1492820, and rs7689420) were analyzed in 1,017 COPD patients (767 males and 246 females) and 430 non-COPD patients. Using logistic regression analysis, we found that rs7689420 was significantly associated with lung cancer in COPD patients in the Chinese Han population (P < 0.001). The recessive allele of rs7689420 was associated with the occurrence of lung cancer in all COPD patients (odds ratios [OR] of 0.609 and 0.424 for the CT and TT genotypes, respectively) as well as in serious COPD patients (OR of 0.403 and 0.305 for CT and TT, respectively). Additionally, rs1489759 and rs3131837 were associated with lung cancer in various genetic models. rs1489758, rs1489759, and rs10519717 were also associated with lung cancer in serious COPD patients. However, none of the SNPs were significantly associated with lung cancer in mild COPD patients or healthy subjects. Therefore, the HHIP SNPs of COPD patients likely play a role in lung cancer pathology in the Chinese Han population.

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV ₁), forced vital capacity (FVC) and the ratio of FEV ₁ to FVC (FEV ₁/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 ^-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

What do polymorphisms tell us about the mechanisms of COPD?

COPD (chronic obstructive pulmonary disease) is characterized by irreversible lung airflow obstruction. Cigarette smoke is the major risk factor for COPD development. However, only a minority number of smokers develop COPD, and there are substantial variations in lung function among smokers, suggesting that genetic determinants in COPD susceptibility. During the past decade, genome-wide association studies and exome sequencing have been instrumental to identify the genetic determinants of complex traits, including COPD. Focused studies have revealed mechanisms by which genetic variants contribute to COPD and have led to novel insights in COPD pathogenesis. Through functional investigations of causal variants in COPD, from the proteinase-antiproteinase theory to emerging roles of developmental pathways (such as Hedgehog and Wnt pathways) in COPD, we have greatly expanded our understanding on this complex pulmonary disease. In this review, we critically review functional investigations on roles of genetic polymorphisms in COPD, and discuss future challenges and opportunities in discovering novel mechanisms of functional variants.

Genetic variants associated with susceptibility to idiopathic pulmonary fibrosis in people of European ancestry: a genome-wide association study

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a chronic progressive lung disease with high mortality, uncertain cause, and few treatment options. Studies have identified a significant genetic risk associated with the development of IPF; however, mechanisms by which genetic risk factors promote IPF remain unclear. We aimed to identify genetic variants associated with IPF susceptibility and provide mechanistic insight using gene and protein expression analyses.

METHODS

We used a two-stage approach: a genome-wide association study in patients with IPF of European ancestry recruited from nine different centres in the UK and controls selected from UK Biobank (stage 1) matched for age, sex, and smoking status; and a follow-up of associated genetic variants in independent datasets of patients with IPF and controls from two independent US samples from the Chicago consortium and the Colorado consortium (stage 2). We investigated the effect of novel signals on gene expression in large transcriptomic and genomic data resources, and examined expression using lung tissue samples from patients with IPF and controls.

FINDINGS

602 patients with IPF and 3366 controls were selected for stage 1. For stage 2, 2158 patients with IPF and 5195 controls were selected. We identified a novel genome-wide significant signal of association with IPF susceptibility near A-kinase anchoring protein 13 (AKAP13; rs62025270, odds ratio [OR] 1·27 [95% CI 1·18-1·37], p=1·32 × 10-9) and confirmed previously reported signals, including in mucin 5B (MUC5B; rs35705950, OR 2·89 [2·56-3·26], p=1·12 × 10-66) and desmoplakin (DSP; rs2076295, OR 1·44 [1·35-1·54], p=7·81 × 10-28). For rs62025270, the allele A associated with increased susceptibility to IPF was also associated with increased expression of AKAP13 mRNA in lung tissue from patients who had lung resection procedures (n=1111). We showed that AKAP13 is expressed in the alveolar epithelium and lymphoid follicles from patients with IPF, and AKAP13 mRNA expression was 1·42-times higher in lung tissue from patients with IPF (n=46) than that in lung tissue from controls (n=51).

INTERPRETATION

AKAP13 is a Rho guanine nucleotide exchange factor regulating activation of RhoA, which is known to be involved in profibrotic signalling pathways. The identification of AKAP13 as a susceptibility gene for IPF increases the prospect of successfully targeting RhoA pathway inhibitors in patients with IPF.

FUNDING

UK Medical Research Council, National Heart, Lung, and Blood Institute of the US National Institutes of Health, Agencia Canaria de Investigación, Innovación y Sociedad de la Información, Spain, UK National Institute for Health Research, and the British Lung Foundation.

Integrative Genomics of Emphysema-Associated Genes Reveals Potential Disease Biomarkers

Chronic obstructive pulmonary disease is the third leading cause of death worldwide. Gene expression profiling across multiple regions of the same lung identified genes significantly related to emphysema. We sought to determine whether the lung and epithelial expression of 127 emphysema-related genes was also related to lung function in independent cohorts, and whether any of these genes could be used as biomarkers in the peripheral blood of patients with chronic obstructive pulmonary disease. To that end, we examined whether the expression levels of these genes were under genetic control in lung tissue (n = 1,111). We then determined whether the mRNA levels of these genes in lung tissue (n = 727), small airway epithelial cells (n = 238), and peripheral blood (n = 620) were significantly related to lung function measurements. The expression of 63 of the 127 genes (50%) was under genetic control in lung tissue. The lung and epithelial mRNA expression of a subset of the emphysema-associated genes, including ASRGL1, LPHN2, and EDNRB, was strongly associated with lung function. In peripheral blood, the expression of 40 genes was significantly associated with lung function. Twenty-nine of these genes (73%) were also associated with lung function in lung tissue, but with the opposite direction of effect for 24 of the 29 genes, including those involved in hypoxia and B cell-related responses. The integrative genomics approach uncovered a significant overlap of emphysema genes associations with lung function between lung and blood with opposite directions between the two. These results support the use of peripheral blood to detect disease biomarkers.

Comparison of the role of HHIP SNPs in susceptibility to chronic obstructive pulmonary disease between Chinese Han and Mongolian populations

The contribution of risk alleles to chronic obstructive pulmonary disease (COPD) may vary between populations. This study aimed to investigate the role of single nucleotide polymorphisms (SNPs) in HHIP on COPD susceptibility in the Chinese Mongolian and Han populations from Inner Mongolia autonomous region. In this case-control study, five SNPs in HHIP gene were detected in 700 COPD patients (350 Chinese Han and 350 Chinese Mongolian) and 700 healthy controls (350 Chinese Han and 350 Chinese Mongolian). The genotype, the association with COPD, and the differences between Chinese Han and Mongolian populations were evaluated using the chi-squared (χ²) test, genetic models, and logistic analysis. The minor allele C in SNP rs10519717 was associated with COPD in Mongolian (Odds ration (OR)=1.401, 95% confidence interval (CI): 1.110-1.769, P=0.044); however, not in the Han. The CC genotype in SNP rs10519717 was a risk factor for COPD in Mongolian (OR=2.667, 95% CI: 1.479-4.809, P=0.044); however, the TC genotype in Han played the same role (OR=1.396, 95% CI: 1.018-1.915, P=0.044). The GG genotype in SNP rs13147758 was protective in the Han (OR=0.546, 95% CI: 0.332-0.897, P=0.017). The homozygote of the minor alleles in SNPs rs12504628, rs1828591, and rs13118928 had a protective role in the both of COPD populations. For the minor allele distribution, the differences between the Han and Mongolian were presented only in the case group for rs12504628 (P=0.003), rs13147758 (P=0.002), rs1828591 (P=0.001), and rs13118928 (P=0.002); for the genotypes, differences was presented at the frequency of the minor allele homozygote in rs13147758 (P=0.048), rs10519717 (P=0.027), rs1828591 (P=0.041), and rs13118928 (P=0.044) in Mongolian. Our findings suggested that HHIP rs10519717 might be associated with susceptibility of Mongolian COPD. For the other SNPs, the differences between the two populations were represented by minor allele distribution and frequency of the minor allele homozygotes.

An independent component analysis confounding factor correction framework for identifying broad impact expression quantitative trait loci

Genome-wide expression Quantitative Trait Loci (eQTL) studies in humans have provided numerous insights into the genetics of both gene expression and complex diseases. While the majority of eQTL identified in genome-wide analyses impact a single gene, eQTL that impact many genes are particularly valuable for network modeling and disease analysis. To enable the identification of such broad impact eQTL, we introduce CONFETI: Confounding Factor Estimation Through Independent component analysis. CONFETI is designed to address two conflicting issues when searching for broad impact eQTL: the need to account for non-genetic confounding factors that can lower the power of the analysis or produce broad impact eQTL false positives, and the tendency of methods that account for confounding factors to model broad impact eQTL as non-genetic variation. The key advance of the CONFETI framework is the use of Independent Component Analysis (ICA) to identify variation likely caused by broad impact eQTL when constructing the sample covariance matrix used for the random effect in a mixed model. We show that CONFETI has better performance than other mixed model confounding factor methods when considering broad impact eQTL recovery from synthetic data. We also used the CONFETI framework and these same confounding factor methods to identify eQTL that replicate between matched twin pair datasets in the Multiple Tissue Human Expression Resource (MuTHER), the Depression Genes Networks study (DGN), the Netherlands Study of Depression and Anxiety (NESDA), and multiple tissue types in the Genotype-Tissue Expression (GTEx) consortium. These analyses identified both cis-eQTL and trans-eQTL impacting individual genes, and CONFETI had better or comparable performance to other mixed model confounding factor analysis methods when identifying such eQTL. In these analyses, we were able to identify and replicate a few broad impact eQTL although the overall number was small even when applying CONFETI. In light of these results, we discuss the broad impact eQTL that have been previously reported from the analysis of human data and suggest that considerable caution should be exercised when making biological inferences based on these reported eQTL.

Genome-wide association analyses for lung function and chronic obstructive pulmonary disease identify new loci and potential druggable targets

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10-49), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

Genetic loci associated with chronic obstructive pulmonary disease overlap with loci for lung function and pulmonary fibrosis

Chronic obstructive pulmonary disease (COPD) is a leading cause of mortality worldwide. We performed a genetic association study in 15,256 cases and 47,936 controls, with replication of select top results (P < 5 × 10-6) in 9,498 cases and 9,748 controls. In the combined meta-analysis, we identified 22 loci associated at genome-wide significance, including 13 new associations with COPD. Nine of these 13 loci have been associated with lung function in general population samples, while 4 (EEFSEC, DSP, MTCL1, and SFTPD) are new. We noted two loci shared with pulmonary fibrosis (FAM13A and DSP) but that had opposite risk alleles for COPD. None of our loci overlapped with genome-wide associations for asthma, although one locus has been implicated in joint susceptibility to asthma and obesity. We also identified genetic correlation between COPD and asthma. Our findings highlight new loci associated with COPD, demonstrate the importance of specific loci associated with lung function to COPD, and identify potential regions of genetic overlap between COPD and other respiratory diseases.

Targeted high-throughput sequencing of candidate genes for chronic obstructive pulmonary disease

Background

Reduced lung function in patients with chronic obstructive pulmonary disease (COPD) is likely due to both environmental and genetic factors. We report here a targeted high-throughput DNA sequencing approach to identify new and previously known genetic variants in a set of candidate genes for COPD.

Methods

Exons in 22 genes implicated in lung development as well as 61 genes and 10 genomic regions previously associated with COPD were sequenced using individual DNA samples from 68 cases with moderate or severe COPD and 66 controls matched for age, gender and smoking. Cases and controls were selected from the Obstructive Lung Disease in Northern Sweden (OLIN) studies.

Results

In total, 37 genetic variants showed association with COPD (p < 0.05, uncorrected). Several variants previously discovered to be associated with COPD from genetic genome-wide analysis studies were replicated using our sample. Two high-risk variants were followed-up for functional characterization in a large eQTL mapping study of 1,111 human lung specimens. The C allele of a synonymous variant, rs8040868, predicting a p.(S45=) in the gene for cholinergic receptor nicotinic alpha 3 (CHRNA3) was associated with COPD (p = 8.8 x 10⁻³). This association remained (p = 0.003 and OR = 1.4, 95 % CI 1.1-1.7) when analysing all available cases and controls in OLIN (n = 1,534). The rs8040868 variant is in linkage disequilibrium with rs16969968 previously associated with COPD and altered expression of the CHRNA5 gene. A follow-up analysis for detection of expression quantitative trait loci revealed that rs8040868-C was found to be significantly associated with a decreased expression of the nearby gene cholinergic receptor, nicotinic, alpha 5 (CHRNA5) in lung tissue.

Conclusion

Our data replicate previous result suggesting CHRNA5 as a candidate gene for COPD and rs8040868 as a risk variant for the development of COPD in the Swedish population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12890-016-0309-y) contains supplementary material, which is available to authorized users.

A Chronic Obstructive Pulmonary Disease Susceptibility Gene, FAM13A, Regulates Protein Stability of β-Catenin

RATIONALE

A genetic locus within the FAM13A gene has been consistently associated with chronic obstructive pulmonary disease (COPD) in genome-wide association studies. However, the mechanisms by which FAM13A contributes to COPD susceptibility are unknown.

OBJECTIVES

To determine the biologic function of FAM13A in human COPD and murine COPD models and discover the molecular mechanism by which FAM13A influences COPD susceptibility.

METHODS

Fam13a null mice (Fam13a(-/-)) were generated and exposed to cigarette smoke. The lung inflammatory response and airspace size were assessed in Fam13a(-/-) and Fam13a(+/+) littermate control mice. Cellular localization of FAM13A protein and mRNA levels of FAM13A in COPD lungs were assessed using immunofluorescence, Western blotting, and reverse transcriptase-polymerase chain reaction, respectively. Immunoprecipitation followed by mass spectrometry identified cellular proteins that interact with FAM13A to reveal insights on FAM13A's function.

MEASUREMENTS AND MAIN RESULTS

In murine and human lungs, FAM13A is expressed in airway and alveolar type II epithelial cells and macrophages. Fam13a null mice (Fam13a(-/-)) were resistant to chronic cigarette smoke-induced emphysema compared with Fam13a(+/+) mice. In vitro, FAM13A interacts with protein phosphatase 2A and recruits protein phosphatase 2A with glycogen synthase kinase 3β and β-catenin, inducing β-catenin degradation. Fam13a(-/-) mice were also resistant to elastase-induced emphysema, and this resistance was reversed by coadministration of a β-catenin inhibitor, suggesting that FAM13A could increase the susceptibility of mice to emphysema development by inhibiting β-catenin signaling. Moreover, human COPD lungs had decreased protein levels of β-catenin and increased protein levels of FAM13A.

CONCLUSIONS

We show that FAM13A may influence COPD susceptibility by promoting β-catenin degradation.

Airway Epithelial Expression Quantitative Trait Loci Reveal Genes Underlying Asthma and Other Airway Diseases

Genome-wide association studies (GWASs) have identified loci that are robustly associated with asthma and related phenotypes; however, the molecular mechanisms underlying these associations need to be explored. The most relevant tissues to study the functional consequences of asthma are the airways. We used publically available data to derive expression quantitative trait loci (eQTLs) for human epithelial cells from small and large airways and applied the eQTLs in the interpretation of GWAS results of asthma and related phenotypes. For the small airways (n = 105), we discovered 660 eQTLs at a 10% false discovery rate (FDR), among which 315 eQTLs were not previously reported in a large-scale eQTL study of whole lung tissue. A large fraction of the identified eQTLs is supported by data from Encyclopedia of DNA Elements (ENCODE) showing that the eQTLs reside in regulatory elements (57.5 and 67.6% of cis- and trans-eQTLs, respectively). Published pulmonary GWAS hits were enriched as airway epithelial eQTLs (9.2-fold). Further, genes regulated by asthma GWAS loci in epithelium are significantly enriched in immune response pathways, such as IL-4 signaling (FDR, 5.2 × 10(-4)). The airway epithelial eQTLs described in this study are complementary to previously reported lung eQTLs and represent a powerful resource to link GWAS-associated variants to their regulatory function and thus elucidate the molecular mechanisms underlying asthma and airway-related conditions.

Sixteen new lung function signals identified through 1000 Genomes Project reference panel imputation

Lung function measures are used in the diagnosis of chronic obstructive pulmonary disease. In 38,199 European ancestry individuals, we studied genome-wide association of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC with 1000 Genomes Project (phase 1)-imputed genotypes and followed up top associations in 54,550 Europeans. We identify 14 novel loci (P<5 × 10(-8)) in or near ENSA, RNU5F-1, KCNS3, AK097794, ASTN2, LHX3, CCDC91, TBX3, TRIP11, RIN3, TEKT5, LTBP4, MN1 and AP1S2, and two novel signals at known loci NPNT and GPR126, providing a basis for new understanding of the genetic determinants of these traits and pulmonary diseases in which they are altered.

Candidate genes for COPD: current evidence and research

COPD is a common complex disease characterized by progressive airflow limitation. Several genome-wide association studies (GWASs) have discovered genes that are associated with COPD. Recently, candidate genes for COPD identified by GWASs include CHRNA3/5 (cholinergic nicotine receptor alpha 3/5), IREB2 (iron regulatory binding protein 2), HHIP (hedgehog-interacting protein), FAM13A (family with sequence similarity 13, member A), and AGER (advanced glycosylation end product–specific receptor). Their association with COPD susceptibility has been replicated in multiple populations. Since these candidate genes have not been considered in COPD, their pathological roles are still largely unknown. Herein, we review some evidences that they can be effective drug targets or serve as biomarkers for diagnosis or subtyping. However, more study is required to understand the functional roles of these candidate genes. Future research is needed to characterize the effect of genetic variants, validate gene function in humans and model systems, and elucidate the genes’ transcriptional and posttranscriptional regulatory mechanisms.

Publisher’s Note:Abstraction for data integration:Fusing mammalian molecular, cellular and phenotype big datasets for better knowledge extraction

With advances in genomics, transcriptomics, metabolomics and proteomics, and more expansive electronic clinical record monitoring, as well as advances in computation, we have entered the Big Data era in biomedical research. Data gathering is growing rapidly while only a small fraction of this data is converted to useful knowledge or reused in future studies. To improve this, an important concept that is often overlooked is data abstraction. To fuse and reuse biomedical datasets from diverse resources, data abstraction is frequently required. Here we summarize some of the major Big Data biomedical research resources for genomics, proteomics and phenotype data, collected from mammalian cells, tissues and organisms. We then suggest simple data abstraction methods for fusing this diverse but related data. Finally, we demonstrate examples of the potential utility of such data integration efforts, while warning about the inherit biases that exist within such data.

Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank

BACKGROUND

Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health.

METHODS

We sampled individuals of European ancestry from UK Biobank, from the middle and extremes of the forced expiratory volume in 1 s (FEV1) distribution among heavy smokers (mean 35 pack-years) and never smokers. We developed a custom array for UK Biobank to provide optimum genome-wide coverage of common and low-frequency variants, dense coverage of genomic regions already implicated in lung health and disease, and to assay rare coding variants relevant to the UK population. We investigated whether there were shared genetic causes between different phenotypes defined by extremes of FEV1. We also looked for novel variants associated with extremes of FEV1 and smoking behaviour and assessed regions of the genome that had already shown evidence for a role in lung health and disease. We set genome-wide significance at p<5 × 10(-8).

FINDINGS

UK Biobank participants were recruited from March 15, 2006, to July 7, 2010. Sample selection for the UK BiLEVE study started on Nov 22, 2012, and was completed on Dec 20, 2012. We selected 50,008 unique samples: 10,002 individuals with low FEV1, 10,000 with average FEV1, and 5002 with high FEV1 from each of the heavy smoker and never smoker groups. We noted a substantial sharing of genetic causes of low FEV1 between heavy smokers and never smokers (p=2.29 × 10(-16)) and between individuals with and without doctor-diagnosed asthma (p=6.06 × 10(-11)). We discovered six novel genome-wide significant signals of association with extremes of FEV1, including signals at four novel loci (KANSL1, TSEN54, TET2, and RBM19/TBX5) and independent signals at two previously reported loci (NPNT and HLA-DQB1/HLA-DQA2). These variants also showed association with COPD, including in individuals with no history of smoking. The number of copies of a 150 kb region containing the 5' end of KANSL1, a gene that is important for epigenetic gene regulation, was associated with extremes of FEV1. We also discovered five new genome-wide significant signals for smoking behaviour, including a variant in NCAM1 (chromosome 11) and a variant on chromosome 2 (between TEX41 and PABPC1P2) that has a trans effect on expression of NCAM1 in brain tissue.

INTERPRETATION

By sampling from the extremes of the lung function distribution in UK Biobank, we identified novel genetic causes of lung function and smoking behaviour. These results provide new insight into the specific mechanisms underlying airflow obstruction, COPD, and tobacco addiction, and show substantial shared genetic architecture underlying airflow obstruction across individuals, irrespective of smoking behaviour and other airway disease.

FUNDING

Medical Research Council.

Molecular mechanisms underlying variations in lung function: a systems genetics analysis

BACKGROUND

Lung function measures reflect the physiological state of the lung, and are essential to the diagnosis of chronic obstructive pulmonary disease (COPD). The SpiroMeta-CHARGE consortium undertook the largest genome-wide association study (GWAS) so far (n=48,201) for forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC) in the general population. The lung expression quantitative trait loci (eQTLs) study mapped the genetic architecture of gene expression in lung tissue from 1111 individuals. We used a systems genetics approach to identify single nucleotide polymorphisms (SNPs) associated with lung function that act as eQTLs and change the level of expression of their target genes in lung tissue; termed eSNPs.

METHODS

The SpiroMeta-CHARGE GWAS results were integrated with lung eQTLs to map eSNPs and the genes and pathways underlying the associations in lung tissue. For comparison, a similar analysis was done in peripheral blood. The lung mRNA expression levels of the eSNP-regulated genes were tested for associations with lung function measures in 727 individuals. Additional analyses identified the pleiotropic effects of eSNPs from the published GWAS catalogue, and mapped enrichment in regulatory regions from the ENCODE project. Finally, the Connectivity Map database was used to identify potential therapeutics in silico that could reverse the COPD lung tissue gene signature.

FINDINGS

SNPs associated with lung function measures were more likely to be eQTLs and vice versa. The integration mapped the specific genes underlying the GWAS signals in lung tissue. The eSNP-regulated genes were enriched for developmental and inflammatory pathways; by comparison, SNPs associated with lung function that were eQTLs in blood, but not in lung, were only involved in inflammatory pathways. Lung function eSNPs were enriched for regulatory elements and were over-represented among genes showing differential expression during fetal lung development. An mRNA gene expression signature for COPD was identified in lung tissue and compared with the Connectivity Map. This in-silico drug repurposing approach suggested several compounds that reverse the COPD gene expression signature, including a nicotine receptor antagonist. These findings represent novel therapeutic pathways for COPD.

INTERPRETATION

The system genetics approach identified lung tissue genes driving the variation in lung function and susceptibility to COPD. The identification of these genes and the pathways in which they are enriched is essential to understand the pathophysiology of airway obstruction and to identify novel therapeutic targets and biomarkers for COPD, including drugs that reverse the COPD gene signature in silico.

FUNDING

The research reported in this article was not specifically funded by any agency. See Acknowledgments for a full list of funders of the lung eQTL study and the Spiro-Meta CHARGE GWAS.

Biomarkers for precision medicine in airways disease

Chronic obstructive pulmonary disease (COPD) is a complex clinical entity. In contrast to previously limited diagnostic definitions, it is now apparent that COPD is a clinically and biologically heterogeneous disease process, overlapping with other airways diseases like chronic asthma. As such, symptomatic response to current standard treatment practices is variable. New clinical guidelines have been altered to reflect this, with the inclusion of symptoms and risk factors in diagnostic and management algorithms. However, as our understanding of COPD pathophysiology deepens, many novel physiological, cellular, proteomic, and genetic markers have been identified. Several have been observed to be independently predictive of distinct clinical disease patterns, which at present are not illustrated by conventional measurements of lung impairment. The potential use of these predictive biomarkers to stratify this diverse patient population could transform the care we offer. We should aim for precision medicine to optimize diagnosis and treatment choices and to monitor and improve clinical outcomes in this disease.

Gene susceptibility identification in a longitudinal study confirms new loci in the development of chronic obstructive pulmonary disease and influences lung function decline

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Informed genome-wide association analysis with family history as a secondary phenotype identifies novel loci of lung cancer

Lung cancer is the leading cause of cancer death worldwide. Although several genetic variants associated with lung cancer have been identified in the past, stringent selection criteria of genome-wide association studies (GWAS) can lead to missed variants. The objective of this study was to uncover missed variants by using the known association between lung cancer and first-degree family history of lung cancer to enrich the variant prioritization for lung cancer susceptibility regions. In this two-stage GWAS study, we first selected a list of variants associated with both lung cancer and family history of lung cancer in four GWAS (3,953 cases, 4,730 controls), then replicated our findings for 30 variants in a meta-analysis of four additional studies (7,510 cases, 7,476 controls). The top ranked genetic variant rs12415204 in chr10q23.33 encoding FFAR4 in the Discovery set was validated in the Replication set with an overall OR of 1.09 (95% CI=1.04, 1.14, P=1.63×10(-4)). When combining the two stages of the study, the strongest association was found in rs1158970 at Ch4p15.2 encoding KCNIP4 with an OR of 0.89 (95% CI=0.85, 0.94, P=9.64×10(-6)). We performed a stratified analysis of rs12415204 and rs1158970 across all eight studies by age, gender, smoking status, and histology, and found consistent results across strata. Four of the 30 replicated variants act as expression quantitative trait loci (eQTL) sites in 1,111 nontumor lung tissues and meet the genome-wide 10% FDR threshold.

Secreted phosphoprotein 1 is a determinant of lung function development in mice

Secreted phosphoprotein 1 (Spp1) is located within quantitative trait loci associated with lung function that was previously identified by contrasting C3H/HeJ and JF1/Msf mouse strains that have extremely divergent lung function. JF1/Msf mice with diminished lung function had reduced lung SPP1 transcript and protein during the peak stage of alveologenesis (postnatal day [P]14-P28) as compared with C3H/HeJ mice. In addition to a previously identified genetic variant that altered runt-related transcription factor 2 (RUNX2) binding in the Spp1 promoter, we identified another promoter variant in a putative RUNX2 binding site that increased the DNA protein binding. SPP1 induced dose-dependent mouse lung epithelial-15 cell proliferation. Spp1((-/-)) mice have decreased specific total lung capacity/body weight, higher specific compliance, and increased mean airspace chord length (Lm) compared with Spp1((+/+)) mice. Microarray analysis revealed enriched gene ontogeny categories, with numerous genes associated with lung development and/or respiratory disease. Insulin-like growth factor 1, Hedgehog-interacting protein, wingless-related mouse mammary tumor virus integration site 5A, and NOTCH1 transcripts decreased in the lung of P14 Spp1((-/-)) mice as determined by quantitative RT-PCR analysis. SPP1 promotes pneumocyte growth, and mice lacking SPP1 have smaller, more compliant lungs with enlarged airspace (i.e., increased Lm). Microarray analysis suggests a dysregulation of key lung developmental transcripts in gene-targeted Spp1((-/-)) mice, particularly during the peak phase of alveologenesis. In addition to its known roles in lung disease, this study supports SPP1 as a determinant of lung development in mice.

Genetic control of gene expression at novel and established chronic obstructive pulmonary disease loci

Genetic risk loci have been identified for a wide range of diseases through genome-wide association studies (GWAS), but the relevant functional mechanisms have been identified for only a small proportion of these GWAS-identified loci. By integrating results from the largest current GWAS of chronic obstructive disease (COPD) with expression quantitative trait locus (eQTL) analysis in whole blood and sputum from 121 subjects with COPD from the ECLIPSE Study, this analysis identifies loci that are simultaneously associated with COPD and the expression of nearby genes (COPD eQTLs). After integrative analysis, 19 COPD eQTLs were identified, including all four previously identified genome-wide significant loci near HHIP, FAM13A, and the 15q25 and 19q13 loci. For each COPD eQTL, fine mapping and colocalization analysis to identify causal shared eQTL and GWAS variants identified a subset of sites with moderate-to-strong evidence of harboring at least one shared variant responsible for both the eQTL and GWAS signals. Transcription factor binding site (TFBS) analysis confirms that multiple COPD eQTL lead SNPs disrupt TFBS, and enhancer enrichment analysis for loci with the strongest colocalization signals showed enrichment for blood-related cell types (CD3 and CD4+ T cells, lymphoblastoid cell lines). In summary, integrative eQTL and GWAS analysis confirms that genetic control of gene expression plays a key role in the genetic architecture of COPD and identifies specific blood-related cell types as likely participants in the functional pathway from GWAS-associated variant to disease phenotype.

Integrative genomics of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a complex disease with both environmental and genetic determinants, the most important of which is cigarette smoking. There is marked heterogeneity in the development of COPD among persons with similar cigarette smoking histories, which is likely partially explained by genetic variation. Genomic approaches such as genomewide association studies and gene expression studies have been used to discover genes and molecular pathways involved in COPD pathogenesis; however, these "first generation" omics studies have limitations. Integrative genomic studies are emerging which can combine genomic datasets to further examine the molecular underpinnings of COPD. Future research in COPD genetics will likely use network-based approaches to integrate multiple genomic data types in order to model the complex molecular interactions involved in COPD pathogenesis. This article reviews the genomic research to date and offers a vision for the future of integrative genomic research in COPD.

Moving beyond genetics: is FAM13A a major biological contributor in lung physiology and chronic lung diseases?

Variants in FAM13A have been found in genome-wide association studies (GWAS) to associate with lung function in the general population as well as in several common chronic lung diseases (CLD) such as chronic obstructive pulmonary disease (COPD), asthma, as well as in idiopathic interstitial pneumonias (IIP). The gene was cloned in 2004, yet the encoded protein has not been characterised and its function is unknown. The redundancy of its genetic contribution in CLD suggests a major function of this gene both in lung physiology and CLD. This review provides a brief summary of the current knowledge of FAM13A, and demonstrates the necessity to resolve its biological function besides its well accepted genetic contribution. Further interpretations of FAM13A variants may help in the understanding of CLD mechanisms and reveal opportunity for intervention.