Chromatin Landscapes of Human Lung Cells Predict Potentially Functional Chronic Obstructive Pulmonary Disease Genome-Wide Association Study Variants

A novel application of data-consistent inversion to overcome spurious inference in genome-wide association studies

The genome-wide association studies (GWAS) typically use linear or logistic regression models to identify associations between phenotypes (traits) and genotypes (genetic variants) of interest. However, the use of regression with the additive assumption has potential limitations. First, the normality assumption of residuals is the one that is rarely seen in practice, and deviation from normality increases the Type-I error rate. Second, building a model based on such an assumption ignores genetic structures, like, dominant, recessive, and protective-risk cases. Ignoring genetic variants may result in spurious conclusions about the associations between a variant and a trait. We propose an assumption-free model built upon data-consistent inversion (DCI), which is a recently developed measure-theoretic framework utilized for uncertainty quantification. This proposed DCI-derived model builds a nonparametric distribution on model inputs that propagates to the distribution of observed data without the required normality assumption of residuals in the regression model. This characteristic enables the proposed DCI-derived model to cover all genetic variants without emphasizing on additivity of the classic-GWAS model. Simulations and a replication GWAS with data from the COPDGene demonstrate the ability of this model to control the Type-I error rate at least as well as the classic-GWAS (additive linear model) approach while having similar or greater power to discover variants in different genetic modes of transmission.

Identification of a distal enhancer regulating hedgehog interacting protein gene in human lung epithelial cells

BACKGROUND

An intergenic region at chromosome 4q31 is one of the most significant regions associated with COPD susceptibility and lung function in GWAS. In this region, the implicated causal gene HHIP has a unique epithelial expression pattern in adult human lungs, in contrast to dominant expression in fibroblasts in murine lungs. However, the mechanism underlying the species-dependent cell type-specific regulation of HHIP remains largely unknown.

METHODS

We employed snATAC-seq analysis to identify open chromatin regions within the COPD GWAS region in various human lung cell types. ChIP-quantitative PCR, reporter assays, chromatin conformation capture assays and Hi-C assays were conducted to characterize the regulatory element in this region. CRISPR/Cas9-editing was performed in BEAS-2B cells to generate single colonies with stable knockout of the regulatory element. RT-PCR and Western blot assays were used to evaluate expression of HHIP and epithelial-mesenchymal transition (EMT)-related marker genes.

FINDINGS

We identified a distal enhancer within the COPD 4q31 GWAS locus that regulates HHIP transcription at baseline and after TGFβ treatment in a SMAD3-dependent, but Hedgehog-independent manner in human bronchial epithelial cells. The distal enhancer also maintains chromatin topological domains near 4q31 locus and HHIP gene. Reduced HHIP expression led to increased EMT induced by TGFβ in human bronchial epithelial cells.

INTERPRETATION

A distal enhancer regulates HHIP expression both under homeostatic condition and upon TGFβ treatment in human bronchial epithelial cells. The interaction between HHIP and TGFβ signalling possibly contributes to COPD pathogenesis.

FUNDING

Supported by NIH grants R01HL127200, R01HL148667 and R01HL162783 (to X. Z).

Precision Approaches to Chronic Obstructive Pulmonary Disease Management

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide. COPD heterogeneity has hampered progress in developing pharmacotherapies that affect disease progression. This issue can be addressed by precision medicine approaches, which focus on understanding an individual's disease risk, and tailoring management based on pathobiology, environmental exposures, and psychosocial issues. There is an urgent need to identify COPD patients at high risk for poor outcomes and to understand at a mechanistic level why certain individuals are at high risk. Genetics, omics, and network analytic techniques have started to dissect COPD heterogeneity and identify patients with specific pathobiology. Drug repurposing approaches based on biomarkers of specific inflammatory processes (i.e., type 2 inflammation) are promising. As larger data sets, additional omics, and new analytical approaches become available, there will be enormous opportunities to identify high-risk individuals and treat COPD patients based on their specific pathophysiological derangements. These approaches show great promise for risk stratification, early intervention, drug repurposing, and developing novel therapeutic approaches for COPD.

Integrating Genetics, Transcriptomics, and Proteomics in Lung Tissue to Investigate Chronic Obstructive Pulmonary Disease

The integration of transcriptomic and proteomic data from lung tissue with chronic obstructive pulmonary disease (COPD)-associated genetic variants could provide insight into the biological mechanisms of COPD. Here, we assessed associations between lung transcriptomics and proteomics with COPD in 98 subjects from the Lung Tissue Research Consortium. Low correlations between transcriptomics and proteomics were generally observed, but higher correlations were found for COPD-associated proteins. We integrated COPD risk SNPs or SNPs near COPD-associated proteins with lung transcripts and proteins to identify regulatory cis-quantitative trait loci (QTLs). Significant expression QTLs (eQTLs) and protein QTLs (pQTLs) were found regulating multiple COPD-associated biomarkers. We investigated mediated associations from significant pQTLs through transcripts to protein levels of COPD-associated proteins. We also attempted to identify colocalized effects between COPD genome-wide association studies and eQTL and pQTL signals. Evidence was found for colocalization between COPD genome-wide association study signals and a pQTL for RHOB and an eQTL for DSP. We applied weighted gene co-expression network analysis to find consensus COPD-associated network modules. Two network modules generated by consensus weighted gene co-expression network analysis were associated with COPD with a false discovery rate lower than 0.05. One network module is related to the catenin complex, and the other module is related to plasma membrane components. In summary, multiple cis-acting determinants of transcripts and proteins associated with COPD were identified. Colocalization analysis, mediation analysis, and correlation-based network analysis of multiple omics data may identify key genes and proteins that work together to influence COPD pathogenesis.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Association of three single nucleotide polymorphisms in interleukin 6 gene with risk of chronic obstructive pulmonary disease

BACKGROUNDS

Experimental studies have implicated the involvement of interleukin 6 (IL6) in chronic obstructive pulmonary disease (COPD), supporting the possible association between genetic alterations in IL6 gene and risk of COPD.

OBJECTIVES

We aimed to examine the association of three polymorphisms in IL6 gene with COPD by integrating the summary results of available reports, and summarize the changes circulating IL6 concentrations across genotypes.

METHODS

Two investigators independently searched medical literature, selected qualified publications, assessed study quality, and abstracted detailed information. STATA software v14.1 was applied for statistical analyses.

RESULTS

Ten articles, including 18 independent studies and 4 815 participants, were comprehensively pooled. When all studies were analyzed together, the association with COPD was not significant for rs1800795 and rs1800797, yet statistical significance was attained for rs1818879 (P < 0.05) across four genetic modes of inheritance. When analyses were restricted to preset subgroups, race, matched condition, and COPD diagnosis were identified as the factors attributed to heterogeneity between studies. Moreover, changes in circulating IL6 concentrations remained nonsignificant across the genotypes of rs1800795 and rs1818879 (P > 0.05). Funnel plots and related tests revealed a low likelihood of publication bias for three studied polymorphisms (P > 0.1).

CONCLUSIONS

Our findings indicated that IL6 gene rs1818879 might be a COPD-susceptibility locus. Further validations in longitudinal studies are necessary.

Genetics of chronic obstructive pulmonary disease: understanding the pathobiology and heterogeneity of a complex disorder

Chronic obstructive pulmonary disease (COPD) is a deadly and highly morbid disease. Susceptibility to and heterogeneity of COPD are incompletely explained by environmental factors such as cigarette smoking. Family-based and population-based studies have shown that a substantial proportion of COPD risk is related to genetic variation. Genetic association studies have identified hundreds of genetic variants that affect risk for COPD, decreased lung function, and other COPD-related traits. These genetic variants are associated with other pulmonary and non-pulmonary traits, demonstrate a genetic basis for at least part of COPD heterogeneity, have a substantial effect on COPD risk in aggregate, implicate early-life events in COPD pathogenesis, and often involve genes not previously suspected to have a role in COPD. Additional progress will require larger genetic studies with more ancestral diversity, improved profiling of rare variants, and better statistical methods. Through integration of genetic data with other omics data and comprehensive COPD phenotypes, as well as functional description of causal mechanisms for genetic risk variants, COPD genetics will continue to inform novel approaches to understanding the pathobiology of COPD and developing new strategies for management and treatment.