Integrative pathway genomics of lung function and airflow obstruction

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signalling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analysed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10’s role in influencing lung’s susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unravelled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung disease

The Ser82 RAGE Variant Affects Lung Function and Serum RAGE in Smokers and sRAGE Production In Vitro

Introduction Genome-Wide Association Studies have identified associations between lung function measures and Chronic Obstructive Pulmonary Disease (COPD) and chromosome region 6p21 containing the gene for the Advanced Glycation End Product Receptor (AGER, encoding RAGE). We aimed to (i) characterise RAGE expression in the lung, (ii) identify AGER transcripts, (iii) ascertain if SNP rs2070600 (Gly82Ser C/T) is associated with lung function and serum sRAGE levels and (iv) identify whether the Gly82Ser variant is functionally important in altering sRAGE levels in an airway epithelial cell model. Methods Immunohistochemistry was used to identify RAGE protein expression in 26 human tissues and qPCR was used to quantify AGER mRNA in lung cells. Gene expression array data was used to identify AGER expression during lung development in 38 fetal lung samples. RNA-Seq was used to identify AGER transcripts in lung cells. sRAGE levels were assessed in cells and patient serum by ELISA. BEAS2B-R1 cells were transfected to overexpress RAGE protein with either the Gly82 or Ser82 variant and sRAGE levels identified. Results Immunohistochemical assessment of 6 adult lung samples identified high RAGE expression in the alveoli of healthy adults and individuals with COPD. AGER/RAGE expression increased across developmental stages in human fetal lung at both the mRNA (38 samples) and protein levels (20 samples). Extensive AGER splicing was identified. The rs2070600T (Ser82) allele is associated with higher FEV1, FEV1/FVC and lower serum sRAGE levels in UK smokers. Using an airway epithelium model overexpressing the Gly82 or Ser82 variants we found that HMGB1 activation of the RAGE-Ser82 receptor results in lower sRAGE production. Conclusions This study provides new information regarding the expression profile and potential role of RAGE in the human lung and shows a functional role of the Gly82Ser variant. These findings advance our understanding of the potential mechanisms underlying COPD particularly for carriers of this AGER polymorphism

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 characterised 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 ratios per standard deviation of the risk score (~6 alleles) (95% confidence interval) 1.24 (1.20-1.27), P=5.05x10^-49) and we observed a 3.7 fold difference in COPD risk between highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in development, elastic fibres 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