Genome-wide association study of lung function and clinical implication in heavy smokers

Multi-trait and multi-ancestry genetic analysis of comorbid lung diseases and traits improves genetic discovery and polygenic risk prediction

While respiratory diseases such as COPD and asthma share many risk factors, most studies investigate them in insolation and in predominantly European ancestry populations. Here, we conducted the most powerful multi-trait and -ancestry genetic analysis of respiratory diseases and auxiliary traits to date. Our approach improves the power of genetic discovery across traits and ancestries, identifying 44 novel loci associated with lung function in individuals of East Asian ancestry. Using these results, we developed PRSxtra (cross TRait and Ancestry), a multi-trait and -ancestry polygenic risk score approach that leverages shared components of heritable risk via pleiotropic effects. PRSxtra significantly improved the prediction of asthma, COPD, and lung cancer compared to trait- and ancestry-matched PRS in a multi-ancestry cohort from the All of Us Research Program, especially in diverse populations. PRSxtra identified individuals in the top decile with over four-fold odds of asthma and COPD compared to the first decile. Our results present a new framework for multi-trait and -ancestry studies of respiratory diseases to improve genetic discovery and polygenic prediction.

Air Pollution Exposure and Interstitial Lung Features in SPIROMICS Participants with Chronic Obstructive Pulmonary Disease

Rationale: It is unknown whether air pollution is associated with radiographic features of interstitial lung disease in individuals with chronic obstructive pulmonary disease (COPD). Objectives: To determine whether air pollution increases the prevalence of interstitial lung abnormalities (ILA) or percent high-attenuation areas (HAA) on computed tomography (CT) in individuals with a heavy smoking history and COPD. Methods: We performed a cross-sectional study of SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study), focused on current or former smokers with COPD. Ten-year exposure to particulate matter ⩽2.5 μm in aerodynamic diameter (PM2.5), nitrogen oxides (NOx), nitrogen dioxide (NO2), and ozone before enrollment CT (completed between 2010 and 2015) were estimated with validated spatiotemporal models at residential addresses. We applied adjusted multivariable modified Poisson regression and linear regression to investigate associations between pollution exposure and relative risk (RR) of ILA or increased percent HAA (between -600 and -250 Hounsfield units), respectively. We assessed for effect modification by MUC5B-promoter polymorphism (variant allele carriers GT or TT vs. GG at rs3705950), smoking status, sex, and percent emphysema. Results: Among 1,272 participants with COPD assessed for HAA, 424 were current smokers, and 249 were carriers of the variant MUC5B allele. A total of 519 participants were assessed for ILA. We found no association between pollution exposure and ILA or HAA. Associations between pollutant exposures and risk of ILA were modified by the presence of MUC5B polymorphism (P value interaction term for NOx = 0.04 and PM2.5 = 0.05) and smoking status (P value interaction term for NOx = 0.05; NO2 = 0.01; and ozone = 0.05). With higher exposure to NOx and PM2.5, MUC5B variant carriers had an increased risk of ILA (RR per 26 ppb NOx, 2.41; 95% confidence interval [CI], 0.97-6.0; and RR per 4 μg ⋅ m-3 PM2.5, 1.43; 95% CI, 0.93-2.2, respectively). With higher exposure to NO2, former smokers had an increased risk of ILA (RR per 10 ppb, 1.64; 95% CI, 1.0-2.7). Conclusions: Exposure to ambient air pollution was not associated with interstitial features on CT in this population of heavy smokers with COPD. MUC5B modified the association between pollution and ILA, suggesting that gene-environment interactions may influence prevalence of interstitial lung features in COPD.

Asthma and Chronic Obstructive Pulmonary Disease

In the United States, asthma and chronic obstructive pulmonary disease (COPD) disproportionately affect African Americans, Puerto Ricans, and other minority groups. Compared with non-Hispanic whites, minorities have been marginalized and more frequently exposed to environmental risk factors such as tobacco smoke and outdoor and indoor pollutants. Such divergent environmental exposures, alone or interacting with heredity, lead to disparities in the prevalence, morbidity, and mortality of asthma and COPD, which are worsened by lack of access to health care. In this article, we review the burden and risk factors for racial or ethnic disparities in asthma and COPD and discuss future directions in this field.

Large scale proteomic studies create novel privacy considerations

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

[Alpha 1-antitrypsin deficiency]

INTRODUCTION

Pulmonary emphysema and liver disease are the clinical expressions of alpha 1-antitrypsin deficiency, an autosomal recessive genetic disease.

STATE OF THE ART

Alpha 1-antitrypsin deficiency is usually associated with the homozygous Z variant of the SERPINA1 gene. Its clinical expression always consists in a substantial reduction of alpha 1-antitrypsin serum concentration and its variants are analyzed by isoelectric focalization or molecular techniques. Assessed by CO transfer alteration and CT scan, risk of pulmonary emphysema is increased by tobacco consumption. Assessed by transient elastography and liver ultrasound, risk of liver disease is increased by alcohol consumption or obesity. Treatment of COPD-associated alpha 1-antitrypsin deficiency does not differ from that of other forms of COPD. In patients presenting with severe deficiency, augmentation therapy with plasma-derived alpha 1-antitrypsin reduces the progression of emphysema, as shown in terms of CT-based lung density metrics. Patients with alpha 1-antitrypsin deficiency with a ZZ genotype should refrain from alcohol or tobacco consumption, and watch their weight; so should their close relatives.

PERSPECTIVES

Modulation of alpha 1-antitrypsin liver production offers an interesting new therapeutic perspective.

CONCLUSION

Homozygous (Z) variants of the SERPINA1 gene confer an increased risk of pulmonary emphysema and liver disease, particularly among smokers, drinkers and obese persons.

Alpha-1 Antitrypsin MZ Heterozygosity Is an Endotype of Chronic Obstructive Pulmonary Disease

Rationale: Multiple studies have demonstrated an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygous carriers of the AAT (alpha-1 antitrypsin) Z allele. However, it is not known if MZ subjects with COPD are phenotypically different from noncarriers (MM genotype) with COPD. Objectives: To assess if MZ subjects with COPD have different clinical features compared with MM subjects with COPD. Methods: Genotypes of SERPINA1 were ascertained by using whole-genome sequencing data in three independent studies. We compared outcomes between MM subjects with COPD and MZ subjects with COPD in each study and combined the results in a meta-analysis. We performed longitudinal and survival analyses to compare outcomes in MM and MZ subjects with COPD over time. Measurements and Main Results: We included 290 MZ subjects with COPD and 6,184 MM subjects with COPD across the three studies. MZ subjects had a lower FEV1% predicted and greater quantitative emphysema on chest computed tomography scans compared with MM subjects. In a meta-analysis, the FEV1 was 3.9% lower (95% confidence interval [CI], -6.55% to -1.26%) and emphysema (the percentage of lung attenuation areas <-950 HU) was 4.14% greater (95% CI, 1.44% to 6.84%) in MZ subjects. We found one gene, PGF (placental growth factor), to be differentially expressed in lung tissue from one study between MZ subjects and MM subjects. Conclusions: Carriers of the AAT Z allele (those who were MZ heterozygous) with COPD had lower lung function and more emphysema than MM subjects with COPD. Taken with the subtle differences in gene expression between the two groups, our findings suggest that MZ subjects represent an endotype of COPD.

A large-scale genome-wide association analysis of lung function in the Chinese population identifies novel loci and highlights shared genetic aetiology with obesity

BACKGROUND

Lung function is a heritable complex phenotype with obesity being one of its important risk factors. However, knowledge of their shared genetic basis is limited. Most genome-wide association studies (GWASs) for lung function have been based on European populations, limiting the generalisability across populations. Large-scale lung function GWASs in other populations are lacking.

METHODS

We included 100 285 subjects from the China Kadoorie Biobank (CKB). To identify novel loci for lung function, single-trait GWAS analyses were performed on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC in the CKB. We then performed genome-wide cross-trait analysis between lung function and obesity traits (body mass index (BMI), BMI-adjusted waist-to-hip ratio and BMI-adjusted waist circumference) to investigate the shared genetic effects in the CKB. Finally, polygenic risk scores (PRSs) of lung function were developed in the CKB and their interaction with BMI's association on lung function were examined. We also conducted cross-trait analysis in parallel with the CKB using up to 457 756 subjects from the UK Biobank (UKB) for replication and investigation of ancestry-specific effects.

RESULTS

We identified nine genome-wide significant novel loci for FEV1, six for FVC and three for FEV1/FVC in the CKB. FEV1 and FVC showed significant negative genetic correlation with obesity traits in both the CKB and UKB. Genetic loci shared between lung function and obesity traits highlighted important biological pathways, including cell proliferation, embryo, skeletal and tissue development, and regulation of gene expression. Mendelian randomisation analysis suggested significant negative causal effects of BMI on FEV1 and on FVC in both the CKB and UKB. Lung function PRSs significantly modified the effect of change in BMI on change in lung function during an average follow-up of 8 years.

CONCLUSION

This large-scale GWAS of lung function identified novel loci and shared genetic aetiology between lung function and obesity. Change in BMI might affect change in lung function differently according to a subject's polygenic background. These findings may open new avenues for the development of molecular-targeted therapies for obesity and lung function improvement.

LncRNA CASC2 is involved in the development of chronic obstructive pulmonary disease via targeting miR-18a-5p/IGF1 axis

AIMS

Chronic obstructive pulmonary disease (COPD) is a systemic disease. Several long non-coding RNAs (lncRNAs) have been identified to be aberrantly expressed in COPD patients. This study investigated the role of lncRNA cancer susceptibility candidate 2 (CASC2) in COPD, as well as its potential mechanism.

METHODS

Fifty smokers with COPD and another 50 smokers without COPD were recruited. Receiver operating characteristic curve was constructed to assess the diagnostic value of CASC2 in COPD patients. 16HBE cells were treated with cigarette smoke extract (CSE) to establish a cell model. qRT-PCR was used for the measurement of mRNA levels. The cell viability and apoptosis were detected by using Cell Counting Kit-8 and flow cytometry assay. Enzyme-linked immunosorbent assay was performed to detect the levels of proinflammatory cytokines. Luciferase reporter assay was performed for the target gene analysis.

RESULTS

Serum CASC2 was dramatically decreased in COPD patients compared with smokers without COPD, and was positively associated with FEV1 (forced expiratory volume in one second). Serum CASC2 was overexpressed in severe COPD patients, and had the diagnostic accuracy to distinguish COPD patients from smokers. CASC2 overexpression alleviated CSE-induced apoptosis and inflammation in 16HBE cells. CASC2 functions as a ceRNA of miR-18a-5p. Upregulation of miR-18a-5p reversed the influence of CASC2 on cell apoptosis and inflammation in 16HBE cells. IGF1 was the target gene of miR-18a-5p.

CONCLUSION

CASC2 was downregulated in COPD patients and it might be a promising biomarker for the disease diagnosis. Overexpression of CASC2 might inhibit the bronchial epithelial cell apoptosis and inflammation via targeting miR-18a-5p/IGF1 axis.The reviews of this paper are available via the supplemental material section.

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.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (β=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (β=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (β=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Plasma Metabolomic Signatures of Chronic Obstructive Pulmonary Disease and the Impact of Genetic Variants on Phenotype-Driven Modules

Background: Small studies have recently suggested that there are specific plasma metabolic signatures in chronic obstructive pulmonary disease (COPD), but there have been no large comprehensive study of metabolomic signatures in COPD that also integrate genetic variants. Materials and Methods: Fresh frozen plasma from 957 non-Hispanic white subjects in COPDGene was used to quantify 995 metabolites with Metabolon's global metabolomics platform. Metabolite associations with five COPD phenotypes (chronic bronchitis, exacerbation frequency, percent emphysema, post-bronchodilator forced expiratory volume at one second [FEV1]/forced vital capacity [FVC], and FEV1 percent predicted) were assessed. A metabolome-wide association study was performed to find genetic associations with metabolite levels. Significantly associated single-nucleotide polymorphisms were tested for replication with independent metabolomic platforms and independent cohorts. COPD phenotype-driven modules were identified in network analysis integrated with genetic associations to assess gene-metabolite-phenotype interactions. Results: Of metabolites tested, 147 (14.8%) were significantly associated with at least 1 COPD phenotype. Associations with airflow obstruction were enriched for diacylglycerols and branched chain amino acids. Genetic associations were observed with 109 (11%) metabolites, 72 (66%) of which replicated in an independent cohort. For 20 metabolites, more than 20% of variance was explained by genetics. A sparse network of COPD phenotype-driven modules was identified, often containing metabolites missed in previous testing. Of the 26 COPD phenotype-driven modules, 6 contained metabolites with significant met-QTLs, although little module variance was explained by genetics. Conclusion: A dysregulation of systemic metabolism was predominantly found in COPD phenotypes characterized by airflow obstruction, where we identified robust heritable effects on individual metabolite abundances. However, network analysis, which increased the statistical power to detect associations missed previously in classic regression analyses, revealed that the genetic influence on COPD phenotype-driven metabolomic modules was modest when compared with clinical and environmental factors.

Genetic regulation of gene expression of MIF family members in lung tissue

Macrophage migration inhibitory factor (MIF) is a cytokine found to be associated with chronic obstructive pulmonary disease (COPD). However, there is no consensus on how MIF levels differ in COPD compared to control conditions and there are no reports on MIF expression in lung tissue. Here we studied gene expression of members of the MIF family MIF, D-Dopachrome Tautomerase (DDT) and DDT-like (DDTL) in a lung tissue dataset with 1087 subjects and identified single nucleotide polymorphisms (SNPs) regulating their gene expression. We found higher MIF and DDT expression in COPD patients compared to non-COPD subjects and found 71 SNPs significantly influencing gene expression of MIF and DDTL. Furthermore, the platform used to measure MIF (microarray or RNAseq) was found to influence the splice variants detected and subsequently the direction of the SNP effects on MIF expression. Among the SNPs found to regulate MIF expression, the major LD block identified was linked to rs5844572, a SNP previously found to be associated with lower diffusion capacity in COPD. This suggests that MIF may be contributing to the pathogenesis of COPD, as SNPs that influence MIF expression are also associated with symptoms of COPD. Our study shows that MIF levels are affected not only by disease but also by genetic diversity (i.e. SNPs). Since none of our significant eSNPs for MIF or DDTL have been described in GWAS for COPD or lung function, MIF expression in COPD patients is more likely a consequence of disease-related factors rather than a cause of the disease.

Genetic analyses identify GSDMB associated with asthma severity, exacerbations, and antiviral pathways

BACKGROUND

The Chr17q12-21.2 region is the strongest and most consistently associated region with asthma susceptibility. The functional genes or single nucleotide polymorphisms (SNPs) are not obvious due to linkage disequilibrium.

OBJECTIVES

We sought to comprehensively investigate whole-genome sequence and RNA sequence from human bronchial epithelial cells to dissect functional genes/SNPs for asthma severity in the Severe Asthma Research Program.

METHODS

Expression quantitative trait loci analysis (n = 114), correlation analysis (n = 156) of gene expression and asthma phenotypes, and pathway analysis were performed in bronchial epithelial cells and replicated. Genetic association for asthma severity (426 severe vs 531 nonsevere asthma) and longitudinal asthma exacerbations (n = 273) was performed.

RESULTS

Multiple SNPs in gasdermin B (GSDMB) associated with asthma severity (odds ratio, >1.25) and longitudinal asthma exacerbations (P < .05). Expression quantitative trait loci analyses identified multiple SNPs associated with expression levels of post-GPI attachment to proteins 3, GSDMB, or gasdermin A (3.1 × 10^-9 <P < 1.8 × 10^-4). Higher expression levels of GSDMB correlated with asthma and greater number of exacerbations (P < .05). Expression levels of GSDMB correlated with genes involved in IFN signaling, MHC class I antigen presentation, and immune system pathways (false-discovery rate-adjusted P < .05). rs1031458 and rs3902920 in GSDMB colocalized with IFN regulatory factor binding sites and associated with GSDMB expression, asthma severity, and asthma exacerbations (P < .05).

CONCLUSIONS

By using a unique set of gene expression data from lung cells obtained using bronchoscopy from comprehensively characterized subjects with asthma, we show that SNPs in GSDMB associated with asthma severity, exacerbations, and GSDMB expression levels. Furthermore, its expression levels correlated with asthma exacerbations and antiviral pathways. Thus, GSDMB is a functional gene for both asthma susceptibility and severity.

Power and Sample Size Calculations for Genetic Association Studies in the Presence of Genetic Model Misspecification

INTRODUCTION

When analyzing data from large-scale genetic association studies, such as targeted or genome-wide resequencing studies, it is common to assume a single genetic model, such as dominant or additive, for all tests of association between a given genetic variant and the phenotype. However, for many variants, the chosen model will result in poor model fit and may lack statistical power due to model misspecification.

OBJECTIVE

We develop power and sample size calculations for tests of gene and gene × environment interaction, allowing for misspecification of the true mode of genetic susceptibility.

METHODS

The power calculations are based on a likelihood ratio test framework and are implemented in an open-source R package ("genpwr").

RESULTS

We use these methods to develop an analysis plan for a resequencing study in idiopathic pulmonary fibrosis and show that using a 2-degree of freedom test can increase power to detect recessive genetic effects while maintaining power to detect dominant and additive effects.

CONCLUSIONS

Understanding the impact of model misspecification can aid in study design and developing analysis plans that maximize power to detect a range of true underlying genetic effects. In particular, these calculations help identify when a multiple degree of freedom test or other robust test of association may be advantageous.

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.

The Effects of Rare SERPINA1 Variants on Lung Function and Emphysema in SPIROMICS

Rationale: The role of PI (protease inhibitor) type Z heterozygotes and additional rare variant genotypes in the gene encoding alpha-1 antitrypsin, SERPINA1 (serpin peptidase inhibitor, clade A, member 1), in determining chronic obstructive pulmonary disease risk and severity is controversial.Objectives: To comprehensively evaluate the effects of rare SERPINA1 variants on lung function and emphysema phenotypes in subjects with significant tobacco smoke exposure using deep gene resequencing and alpha-1 antitrypsin concentrations.Methods: DNA samples from 1,693 non-Hispanic white individuals, 385 African Americans, and 90 Hispanics with ≥20 pack-years smoking were resequenced for the identification of rare variants (allele frequency < 0.05) in 16.9 kB of SERPINA1.Measurements and Main Results: White PI Z heterozygotes confirmed by sequencing (MZ; n = 74) had lower post-bronchodilator FEV1 (P = 0.007), FEV1/FVC (P = 0.003), and greater computed tomography-based emphysema (P = 0.02) compared with 1,411 white individuals without PI Z, S, or additional rare variants denoted as VR. PI Z-containing compound heterozygotes (ZS/ZVR; n = 7) had lower FEV1/FVC (P = 0.02) and forced expiratory flow, midexpiratory phase (P = 0.009). Nineteen white heterozygotes for five non-S/Z coding variants associated with lower alpha-1 antitrypsin had greater computed tomography-based emphysema compared with those without rare variants. In African Americans, a 5' untranslated region insertion (rs568223361) was associated with lower alpha-1 antitrypsin and functional small airway disease (P = 0.007).Conclusions: In this integrative deep sequencing study of SERPINA1 with alpha-1 antitrypsin concentrations in a heavy smoker and chronic obstructive pulmonary disease cohort, we confirmed the effects of PI Z heterozygote and compound heterozygote genotypes. We demonstrate the cumulative effects of multiple SERPINA1 variants on alpha-1 antitrypsin deficiency, lung function, and emphysema, thus significantly increasing the frequency of SERPINA1 variation associated with respiratory disease in at-risk smokers.

Integrated transcriptomic correlation network analysis identifies COPD molecular determinants

Chronic obstructive pulmonary disease (COPD) is a complex and heterogeneous syndrome. Network-based analysis implemented by SWIM software can be exploited to identify key molecular switches - called "switch genes" - for the disease. Genes contributing to common biological processes or defining given cell types are usually co-regulated and co-expressed, forming expression network modules. Consistently, we found that the COPD correlation network built by SWIM consists of three well-characterized modules: one populated by switch genes, all up-regulated in COPD cases and related to the regulation of immune response, inflammatory response, and hypoxia (like TIMP1, HIF1A, SYK, LY96, BLNK and PRDX4); one populated by well-recognized immune signature genes, all up-regulated in COPD cases; one where the GWAS genes AGER and CAVIN1 are the most representative module genes, both down-regulated in COPD cases. Interestingly, 70% of AGER negative interactors are switch genes including PRDX4, whose activation strongly correlates with the activation of known COPD GWAS interactors SERPINE2, CD79A, and POUF2AF1. These results suggest that SWIM analysis can identify key network modules related to complex diseases like COPD.

Genetics of COPD

Although chronic obstructive pulmonary disease (COPD) risk is strongly influenced by cigarette smoking, genetic factors are also important determinants of COPD. In addition to Mendelian syndromes such as alpha-1 antitrypsin deficiency, many genomic regions that influence COPD susceptibility have been identified in genome-wide association studies. Similarly, multiple genomic regions associated with COPD-related phenotypes, such as quantitative emphysema measures, have been found. Identifying the functional variants and key genes within these association regions remains a major challenge. However, newly identified COPD susceptibility genes are already providing novel insights into COPD pathogenesis. Network-based approaches that leverage these genetic discoveries have the potential to assist in decoding the complex genetic architecture of COPD.

Genetic underpinnings of lung function and COPD

Spirometry based measurement of lung function is a global initiative for chronic obstructive lung disease (GOLD) standard to diagnose chronic obstructive pulmonary disease (COPD), one of the leading causes of mortality worldwide. The environmental and behavioural risk factors for COPD includes tobacco smoking, air pollutants and biomass fuel exposure, which can induce one or more abnormal lung function patterns. While smoking remains the primary risk factor, only 15-20% smokers develop COPD, indicating that the genetic factors are also likely to play a role. According to the study of Global Burden of Disease 2015, ∼174 million people across the world have COPD. From a comprehensive literature search conducted using the 'PubMed' and 'GWAS Catalogue' databases, and reviewing the literature available, only a limited number of studies were identified which had attempted to investigate the genetics of COPD and lung volumes, implying a huge research gap. With the advent of genomewide association studies several genetic variants linked to lung function and COPD, like HHIP, HTR4, ADAM19 and GSTCD etc., have been found and validated in different population groups, suggesting their potential role in determining lung volume and risk for COPD. This article aims at reviewing the present knowledge of the genetics of lung function and COPD.

Defining a role for lung function associated gene GSTCD in cell homeostasis

Genome wide association (GWA) studies have reproducibly identified signals on chromosome 4q24 associated with lung function and COPD. GSTCD (Glutathione S-transferase C-terminal domain containing) represents a candidate causal gene in this locus, however little is currently known about the function of this protein. We set out to further our understanding of the role of GSTCD in cell functions and homeostasis using multiple molecular and cellular approaches in airway relevant cells. Recombinant expression of human GSTCD in conjunction with a GST activity assay did not identify any enzymatic activity for two GSTCD isoforms questioning the assignment of this protein to this family of enzymes. Protein structure analyses identified a potential methyltransferase domain contained within GSTCD, with these enzymes linked to cell viability and apoptosis. Targeted knockdown (siRNA) of GSTCD in bronchial epithelial cells identified a role for GSTCD in cell viability as proliferation rates were not altered. To provide greater insight we completed transcriptomic analyses on cells with GSTCD expression knocked down and identified several differentially expressed genes including those implicated in airway biology; fibrosis e.g. TGFBR1 and inflammation e.g. IL6R. Pathway based transcriptomic analyses identified an over-representation of genes related to adipogenesis which may suggest additional functions for GSTCD. These findings identify potential additional functions for GSTCD in the context of airway biology beyond the hypothesised GST activity and warrant further investigation.

Identification of an emphysema-associated genetic variant near TGFB2 with regulatory effects in lung fibroblasts

Murine studies have linked TGF-β signaling to emphysema, and human genome-wide association studies (GWAS) studies of lung function and COPD have identified associated regions near genes in the TGF-β superfamily. However, the functional regulatory mechanisms at these loci have not been identified. We performed the largest GWAS of emphysema patterns to date, identifying 10 GWAS loci including an association peak spanning a 200 kb region downstream from TGFB2. Integrative analysis of publicly available eQTL, DNaseI, and chromatin conformation data identified a putative functional variant, rs1690789, that may regulate TGFB2 expression in human fibroblasts. Using chromatin conformation capture, we confirmed that the region containing rs1690789 contacts the TGFB2 promoter in fibroblasts, and CRISPR/Cas-9 targeted deletion of a ~ 100 bp region containing rs1690789 resulted in decreased TGFB2 expression in primary human lung fibroblasts. These data provide novel mechanistic evidence linking genetic variation affecting the TGF-β pathway to emphysema in humans.

The etiologic origins for chronic obstructive pulmonary disease

COPD, characterized by long-term poorly irreversible airway limitation and persistent respiratory symptoms, has resulted in enormous challenges to human health worldwide, with increasing rates of prevalence, death, and disability. Although its origin was thought to be in the interactions of genetic with environmental factors, the effects of environmental factors on the disease during different life stages remain little known. Without clear mechanisms and radical cure for it, early screening and prevention of COPD seem to be important. In this review, we will discuss the etiologic origins for poor lung function and COPD caused by specific adverse effects during corresponding life stages, as well as try to find new insights and potential prevention strategies for this disease.