Gene expression profiling of human lung tissue from smokers with severe emphysema

Pathophysiology and genomics of bronchiectasis

Bronchiectasis is a complex and heterogeneous inflammatory chronic respiratory disease with an unknown cause in around 30-40% of patients. The presence of airway infection together with chronic inflammation, airway mucociliary dysfunction and lung damage are key components of the vicious vortex model that better describes its pathophysiology. Although bronchiectasis research has significantly increased over the past years and different endotypes have been identified, there are still major gaps in the understanding of the pathophysiology. Genomic approaches may help to identify new endotypes, as has been shown in other chronic airway diseases, such as COPD.Different studies have started to work in this direction, and significant contributions to the understanding of the microbiome and proteome diversity have been made in bronchiectasis in recent years. However, the systematic application of omics approaches to identify new molecular insights into the pathophysiology of bronchiectasis (endotypes) is still limited compared with other respiratory diseases.Given the complexity and diversity of these technologies, this review describes the key components of the pathophysiology of bronchiectasis and how genomics can be applied to increase our knowledge, including the study of new techniques such as proteomics, metabolomics and epigenomics. Furthermore, we propose that the novel concept of trained innate immunity, which is driven by microbiome exposures leading to epigenetic modifications, can complement our current understanding of the vicious vortex. Finally, we discuss the challenges, opportunities and implications of genomics application in clinical practice for better patient stratification into new therapies.

Thresholding Gini variable importance with a single-trained random forest: An empirical Bayes approach

Random forests (RFs) are a widely used modelling tool capable of feature selection via a variable importance measure (VIM), however, a threshold is needed to control for false positives. In the absence of a good understanding of the characteristics of VIMs, many current approaches attempt to select features associated to the response by training multiple RFs to generate statistical power via a permutation null, by employing recursive feature elimination, or through a combination of both. However, for high-dimensional datasets these approaches become computationally infeasible. In this paper, we present RFlocalfdr, a statistical approach, built on the empirical Bayes argument of Efron, for thresholding mean decrease in impurity (MDI) importances. It identifies features significantly associated with the response while controlling the false positive rate. Using synthetic data and real-world data in health, we demonstrate that RFlocalfdr has equivalent accuracy to currently published approaches, while being orders of magnitude faster. We show that RFlocalfdr can successfully threshold a dataset of 106 datapoints, establishing its usability for large-scale datasets, like genomics. Furthermore, RFlocalfdr is compatible with any RF implementation that returns a VIM and counts, making it a versatile feature selection tool that reduces false discoveries.

CEACAM6 as a Novel Therapeutic Target to Boost HO-1-mediated Antioxidant Defense in COPD

Rationale: Tobacco smoking and air pollution are primary causes of chronic obstructive pulmonary disease (COPD). However, only a minority of smokers develop COPD. The mechanisms underlying the defense against nitrosative/oxidative stress in nonsusceptible smokers to COPD remain largely unresolved. Objectives: To investigate the defense mechanisms against nitrosative/oxidative stress that possibly prevent COPD development or progression. Methods: Four cohorts were investigated: 1) sputum samples (healthy, n = 4; COPD, n = 37), 2) lung tissue samples (healthy, n = 13; smokers without COPD, n = 10; smoker+COPD, n = 17), 3) pulmonary lobectomy tissue samples (no/mild emphysema, n = 6), and 4) blood samples (healthy, n = 6; COPD, n = 18). We screened 3-nitrotyrosine (3-NT) levels, as indication of nitrosative/oxidative stress, in human samples. We established a novel in vitro model of a cigarette smoke extract (CSE)-resistant cell line and studied 3-NT formation, antioxidant capacity, and transcriptomic profiles. Results were validated in lung tissue, isolated primary cells, and an ex vivo model using adeno-associated virus-mediated gene transduction and human precision-cut lung slices. Measurements and Main Results: 3-NT levels correlate with COPD severity of patients. In CSE-resistant cells, nitrosative/oxidative stress upon CSE treatment was attenuated, paralleled by profound upregulation of heme oxygenase-1 (HO-1). We identified carcinoembryonic antigen cell adhesion molecule 6 (CEACAM6) as a negative regulator of HO-1-mediated nitrosative/oxidative stress defense in human alveolar type 2 epithelial cells (hAEC2s). Consistently, inhibition of HO-1 activity in hAEC2s increased the susceptibility toward CSE-induced damage. Epithelium-specific CEACAM6 overexpression increased nitrosative/oxidative stress and cell death in human precision-cut lung slices on CSE treatment. Conclusions: CEACAM6 expression determines the hAEC2 sensitivity to nitrosative/oxidative stress triggering emphysema development/progression in susceptible smokers.

Severity of Lung Function Impairment Drives Transcriptional Phenotypes of COPD and Relates to Immune and Metabolic Processes

Purpose

This study sought to characterize transcriptional phenotypes of COPD through unsupervised clustering of sputum gene expression profiles, and further investigate mechanisms underlying the characteristics of these clusters.

Patients and methods

Induced sputum samples were collected from patients with stable COPD (n = 72) and healthy controls (n = 15). Induced sputum was collected for inflammatory cell counts, and RNA extracted. Transcriptional profiles were generated (Illumina Humanref-8 V2) and analyzed by GeneSpring GX14.9.1. Unsupervised hierarchical clustering and differential gene expression analysis were performed, and gene alterations validated in the ECLIPSE dataset (GSE22148).

Results

We identified 2 main clusters (Cluster 1 [n = 35] and Cluster 2 [n = 37]), which further divided into 4 sub-clusters (Sub-clusters 1.1 [n = 14], 1.2 [n = 21], 2.1 [n = 20] and 2.2 [n = 17]). Compared with Cluster 1, Cluster 2 was associated with significantly lower lung function (p = 0.014), more severe disease (p = 0.009) and breathlessness (p = 0.035), and increased sputum neutrophils (p = 0.031). Sub-cluster 1.1 had significantly higher proportion of people with comorbid cardiovascular disease compared to the other 3 sub-clusters (92.5% vs 57.1%, 50% and 52.9%, p < 0.013). Through supervised analysis we determined that degree of airflow limitation (GOLD stage) was the predominant factor driving gene expression differences in our transcriptional clusters. There were 452 genes (adjusted p < 0.05 and ≥2 fold) altered in GOLD stage 3 and 4 versus 1 and 2, of which 281 (62%) were also found to be significantly expressed between these GOLD stages in the ECLIPSE data set (GSE22148). Differentially expressed genes were largely downregulated in GOLD stages 3 and 4 and connected in 5 networks relating to lipoprotein and cholesterol metabolism; metabolic processes in oxidation/reduction and mitochondrial function; antigen processing and presentation; regulation of complement activation and innate immune responses; and immune and metabolic processes.

Conclusion

Severity of lung function drives 2 distinct transcriptional phenotypes of COPD and relates to immune and metabolic processes.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Dysregulated Polycomb Repressive Complex 2 contributes to chronic obstructive pulmonary disease by rewiring stem cell fate

Aberrant lung cell differentiation is a hallmark of many lung diseases including chronic obstructive pulmonary disease (COPD). The EZH2-containing Polycomb Repressive Complex 2 (PRC2) regulates embryonic lung stem cell fate, but its role in adult lung is obscure. Histological analysis of patient tissues revealed that loss of PRC2 activity was correlated with aberrant bronchiolar cell differentiation in COPD lung. Histological and single-cell RNA-sequencing analyses showed that loss of EZH2 in mouse lung organoids led to lowered self-renewal capability, increased squamous morphological development, and marked shifts in progenitor cell populations. Evaluation of in vivo models revealed that heterozygosity of Ezh2 in mice with ovalbumin-induced lung inflammation led to epithelial cell differentiation patterns similar to those in COPD lung. We also identified cystathionine-β-synthase as a possible upstream factor for PRC2 destabilization. Our findings suggest that PRC2 is integral to facilitating proper lung stem cell differentiation in humans and mice.

Clustering-based COPD subtypes have distinct longitudinal outcomes and multi-omics biomarkers

Introduction

Chronic obstructive pulmonary disease (COPD) can progress across several domains, complicating the identification of the determinants of disease progression. In our previous work, we applied k-means clustering to spirometric and chest radiological measures to identify four COPD-related subtypes: ‘relatively resistant smokers (RRS)’, ‘mild upper lobe-predominant emphysema (ULE)’, ‘airway-predominant disease (AD)’ and ‘severe emphysema (SE)’. In the current study, we examined the associations of these subtypes to longitudinal COPD-related health measures as well as blood transcriptomic and plasma proteomic biomarkers.

Methods

We included 8266 non-Hispanic white and African-American smokers from the COPDGene study. We used linear regression to investigate cluster associations to 5-year prospective changes in spirometric and radiological measures and to gene expression and protein levels. We used Cox-proportional hazard test to test for cluster associations to prospective exacerbations, comorbidities and mortality.

Results

The RRS, ULE, AD and SE clusters represented 39%, 15%, 26% and 20% of the studied cohort at baseline, respectively. The SE cluster had the greatest 5-year FEV₁ (forced expiratory volume in 1 s) and emphysema progression, and the highest risks of exacerbations, cardiovascular disease and mortality. The AD cluster had the highest diabetes risk. After adjustments, only the SE cluster had an elevated respiratory mortality risk, while the ULE, AD and SE clusters had elevated all-cause mortality risks. These clusters also demonstrated differential protein and gene expression biomarker associations, mostly related to inflammatory and immune processes.

Conclusion

COPD k-means subtypes demonstrate varying rates of disease progression, prospective comorbidities, mortality and associations to transcriptomic and proteomic biomarkers. These findings emphasise the clinical and biological relevance of these subtypes, which call for more study for translation into clinical practice.

Trail registration number

NCT00608764.

Association of childhood BMI trajectory with post-adolescent and adult lung function is mediated by pre-adolescent DNA methylation

Background

Body mass index (BMI) has been shown to be associated with lung function. Recent findings showed that DNA methylation (DNAm) variation is likely to be a consequence of changes in BMI. However, whether DNAm mediates the association of BMI with lung function is unknown. We examined the mediating role of DNAm on the association of pre-adolescent BMI trajectories with post-adolescent and adulthood lung function (forced expiratory volume (FEV₁), forced vital capacity (FVC), and FEV₁/FVC).

Methods

Analyses were undertaken in the Isle of Wight birth cohort (IOWBC). Group-based trajectory modelling was applied to infer latent BMI trajectories from age 1 to 10 years. An R package, ttscreening, was applied to identify CpGs at 10 years potentially associated with BMI trajectories for each sex. Linear regressions were implemented to further screen CpGs for their association with lung function at 18 years. Path analysis, stratified by sex, was applied to each screened CpG to assess its role of mediation. Internal validation was applied to further examine the mediation consistency of the detected CpGs based on lung function at 26 years. Mendelian randomization (MR-base) was used to test possible causal effects of the identified CpGs.

Results

Two BMI trajectories (high vs. low) were identified. Of the 442,475 CpG sites, 18 CpGs in males and 33 in females passed screening. Eight CpGs in males and 16 CpGs in females (none overlapping) were identified as mediators. For subjects with high BMI trajectory, high DNAm at all CpGs in males were associated with decreased lung function, while 8 CpGs in females were associated with increased lung function at 18 years. At 26 years, 6 CpGs in males and 14 CpGs in females showed the same direction of indirect effects as those at 18 years. DNAm at CpGs cg19088553 (GRIK2) and cg00612625 (HPSE2) showed a potential causal effect on FEV₁.

Conclusions

The effects of BMI trajectory in early childhood on post-adolescence lung function were likely to be mediated by pre-adolescence DNAm in both males and females, but such mediation effects were likely to diminish over time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-022-02089-4.

Mitochondrial ribosomal stress in lung diseases

Mitochondria are involved in a variety of critical cellular functions, and their impairment drives cell injury. The mitochondrial ribosome (mitoribosome) is responsible for the protein synthesis of mitochondrial DNA encoded genes. These proteins are involved in oxidative phosphorylation, respiration, and ATP production required in the cell. Mitoribosome components originate from both mitochondrial and nuclear genomes. Their dysfunction can be caused by impaired mitochondrial protein synthesis or mitoribosome misassembly, leading to a decline in mitochondrial translation. This decrease can trigger mitochondrial ribosomal stress and contribute to pulmonary cell injury, death, and diseases. This review focuses on the contribution of the impaired mitoribosome structural components and function to respiratory disease pathophysiology. We present recent findings in the fields of lung cancer, chronic obstructive pulmonary disease, interstitial lung disease, and asthma. We also include reports on the mitoribosome dysfunction in pulmonary hypertension, high altitude pulmonary edema, bacterial and viral infections. Studies of the mitoribosome alterations in respiratory diseases can lead to novel therapeutic targets.

Pulmonary MicroRNA Changes Alter Angiogenesis in Chronic Obstructive Pulmonary Disease and Lung Cancer

The pulmonary endothelium is dysfunctional in chronic obstructive pulmonary disease (COPD), a known risk factor for lung cancer. The pulmonary endothelium is altered in emphysema, which is disproportionately affected by cancers. Gene and microRNA expression differs between COPD and non-COPD lung. We hypothesised that the alteration in microRNA expression in the pulmonary endothelium contributes to its dysfunction. A total of 28 patients undergoing pulmonary resection were recruited and endothelial cells were isolated from healthy lung and tumour. MicroRNA expression was compared between COPD and non-COPD patients. Positive findings were confirmed by quantitative polymerase chain reaction (qPCR). Assays assessing angiogenesis and cellular migration were conducted in Human Umbilical Vein Endothelial Cells (n = 3-4) transfected with microRNA mimics and compared to cells transfected with negative control RNA. Expression of miR-181b-3p, miR-429 and miR-23c (all p < 0.05) was increased in COPD. Over-expression of miR-181b-3p was associated with reduced endothelial sprouting (p < 0.05). miR-429 was overexpressed in lung cancer as well and exhibited a reduction in tubular formation. MicroRNA-driven changes in the pulmonary endothelium thus represent a novel mechanism driving emphysema. These processes warrant further study to determine if they may be therapeutic targets in COPD and lung cancer.

Alveolar epithelial inter-alpha-trypsin inhibitor heavy chain 4 deficiency associated with senescence-regulated apoptosis by air pollution

Inter-alpha-trypsin inhibitor heavy chain 4 (ITIH4) is a type II acute-phase protein; however, the role of pulmonary ITIH4 after exposure to air pollution remains unclear. In this study, we investigated the role of ITIH4 in the lungs in response to air pollution. ITIH4 expression in bronchoalveolar lavage fluid (BAL) of 47 healthy human subjects and of Sprague-Dawley rats whole-body exposed to air pollution was determined, and the underlying antiapoptotic and matrix-stabilizing pathways in alveolar epithelial A549 cells induced by diesel exhaust particles (DEPs) as well as ITIH4-knockdown were investigated. We found that an interquartile range (IQR) increase in PM_2.5 was associated with a decrease of 2.673 ng/mL in ITIH4, an increase of 1.104 pg/mL of 8-isoprostane, and an increase of 6.918 pg/mL of interleukin (IL)-6 in human BAL. In rats, increases in 8-isoprostane, IL-6, and p53 and a decrease in sirtuin-1 (Sirt1) in the lungs and decreases in ITIH4 in the BAL, lungs, and serum were observed after PM_2.5 and gaseous exposure. ITIH4 levels in lung lysates were correlated with levels in BAL samples (r = 0.377, p < 0.01), whereas ITIH4 levels in BAL were correlated with IL-6 levels (r = -0.420, p < 0.01). ITIH4 expression was significantly reduced in alveolar epithelial A549 cells by DEP in a dose-dependent manner. A decrease in Sirt1 and increases in phosphorylated extracellular signal-regulated kinase (p-ERK) and caspase-3 were observed after DEP exposure and ITIH4-knockdown. In conclusion, air pollution decreased ITIH4 expression in the lungs, which was associated with alveolar epithelial cell senescence and apoptosis. ITIH4 could be a vital protein in regulating alveolar cell destruction and its inhibition after exposure to air pollution.

Single cell RNA sequencing identifies IGFBP5 and QKI as ciliated epithelial cell genes associated with severe COPD

Background

Whole lung tissue transcriptomic profiling studies in chronic obstructive pulmonary disease (COPD) have led to the identification of several genes associated with the severity of airflow limitation and/or the presence of emphysema, however, the cell types driving these gene expression signatures remain unidentified.

Methods

To determine cell specific transcriptomic changes in severe COPD, we conducted single-cell RNA sequencing (scRNA seq) on n = 29,961 cells from the peripheral lung parenchymal tissue of nonsmoking subjects without underlying lung disease (n = 3) and patients with severe COPD (n = 3). The cell type composition and cell specific gene expression signature was assessed. Gene set enrichment analysis (GSEA) was used to identify the specific cell types contributing to the previously reported transcriptomic signatures.

Results

T-distributed stochastic neighbor embedding and clustering of scRNA seq data revealed a total of 17 distinct populations. Among them, the populations with more differentially expressed genes in cases vs. controls (log fold change >|0.4| and FDR = 0.05) were: monocytes (n = 1499); macrophages (n = 868) and ciliated epithelial cells (n = 590), respectively. Using GSEA, we found that only ciliated and cytotoxic T cells manifested a trend towards enrichment of the previously reported 127 regional emphysema gene signatures (normalized enrichment score [NES] = 1.28 and = 1.33, FDR = 0.085 and = 0.092 respectively). Among the significantly altered genes present in ciliated epithelial cells of the COPD lungs, QKI and IGFBP5 protein levels were also found to be altered in the COPD lungs.

Conclusions

scRNA seq is useful for identifying transcriptional changes and possibly individual protein levels that may contribute to the development of emphysema in a cell-type specific manner.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-021-01675-2.

Exposure to Polyhexamethylene Guanidine Exacerbates Bronchial Hyperresponsiveness and Lung Inflammation in a Mouse Model of Ovalbumin-Induced Asthma

Humidifier disinfectants (HDs) exposure has now been associated with acute lung injury and pulmonary fibrosis; polyhexamethylene guanidine (PHMG) has been confirmed to cause severe lung inflammation and fibrosis in mice. Recent evidence also indicates that HDs exposure increases the asthma risk in children, but the underlying mechanisms remain unclear. We aimed to investigate the effects of PHMG exposure on asthma in mice and the potential underlying mechanisms. BALB/c mice were intranasally administered PHMG (0.1 mg/kg/day; 5 days per week) during 2 episodes of ovalbumin (OVA) sensitization and were then challenged with 1% OVA by inhalation. Bronchial hyperresponsiveness (BHR), inflammatory cell influx into bronchoalveolar lavage (BAL) fluid, serum total and OVA-specific immunoglobulin (Ig) E levels, and histopathological changes in the lung were analyzed. The levels of asthma-related cytokines and chemokines were assayed in the lung tissues to evaluate possible mechanisms. Exposure to PHMG following OVA sensitization and challenge significantly enhanced BHR, inflammatory cell counts in BAL fluid, airway inflammation, and total serum IgE levels in the asthma mouse model. In addition, the levels of chemokine ligand (CCL) 11 and serpine F1/pigment epithelium-derived factor (SERPINF1) were significantly elevated in the lungs of these mice compared to those in the control and OVA-treated only groups. Our findings suggest that PHMG can enhance the development of allergic responses and lung inflammation via CCL11- and SERPINF1-induced signaling in a mouse model of asthma.

Multiple Gene Expression Dataset Analysis Reveals Toll-Like Receptor Signaling Pathway is Strongly Associated With Chronic Obstructive Pulmonary Disease Pathogenesis

Chronic obstructive pulmonary disease is a complex pulmonary disease that causes airflow obstruction in humans. To identify the core genes in COPD pathogenesis, seven diverse microarray datasets (GSE475, GSE1122, GSE1650, GSE3212, GSE8823, GSE37768, and GSE22148) were downloaded from the gene expression omnibus database. All the datasets were analyzed independently with the R/Bioconductor package to screen the differentially expressed genes (DEGs). The gene ontology and pathway enrichment analysis were performed for the acquired DEGs using DAVID (Database for Annotation, Visualization, and Integrated Discovery). Further protein-protein interaction network was constructed for the DEGs and their potential hub genes and sub-networks were identified using Cytoscape software. From the selected seven datasets, 188 overlapped DEGs were perceived eventually based on considering the repetitive genes between at-least one dataset. Gene Ontology analysis reveals that most of the DEGs were significantly enriched in immune response, inflammatory response, extracellular region, lipid binding, cytokine, and chemokine activity. Moreover, genes from the sub-network analysis were again submitted to the DAVID server to validate the results which uncover the Toll-like receptor signaling pathway was significantly enriched and all the genes present in this pathway were likewise detected as hub genes from Cytoscape software. CXCL9, CXCL10, CXCL11, CCL4, TLR7, and SPP1 hub genes in the toll-like receptor signaling pathway were explored in this study as potential biomarker genes associated with COPD pathogenesis.

Pulmonary Artery Thrombosis: A Diagnosis That Strives for Its Independence

According to a widespread theory, thrombotic masses are not formed in the pulmonary artery (PA) but result from migration of blood clots from the venous system. This concept has prevailed in clinical practice for more than a century. However, a new technologic era has brought forth more diagnostic possibilities, and it has been shown that thrombotic masses in the PA could, in many cases, be found without any obvious source of emboli. Chronic obstructive pulmonary disease, asthma, sickle cell anemia, emergency and elective surgery, viral pneumonia, and other conditions could be complicated by PA thrombosis development without concomitant deep vein thrombosis (DVT). Different pathologies have different causes for local PA thrombotic process. As evidenced by experimental results and clinical observations, endothelial and platelet activation are the crucial mechanisms of this process. Endothelial dysfunction can impair antithrombotic function of the arterial wall through downregulation of endothelial nitric oxide synthase (eNOS) or via stimulation of adhesion receptor expression. Hypoxia, proinflammatory cytokines, or genetic mutations may underlie the procoagulant phenotype of the PA endothelium. Both endotheliocytes and platelets could be activated by protease mediated receptor (PAR)- and receptors for advanced glycation end (RAGE)-dependent mechanisms. Hypoxia, in particular induced by high altitudes, could play a role in thrombotic complications as a trigger of platelet activity. In this review, we discuss potential mechanisms of PA thrombosis in situ.

Gene expression profiling of bronchial brushes is associated with the level of emphysema measured by computed tomography-based parametric response mapping

Parametric response mapping (PRM) is a computed tomography (CT)-based method to phenotype patients with chronic obstructive pulmonary disease (COPD). It is capable of differentiating emphysema-related air trapping with nonemphysematous air trapping (small airway disease), which helps to identify the extent and localization of the disease. Most studies evaluating the gene expression in smokers and COPD patients related this to spirometric measurements, but none have investigated the relationship with CT-based measurements of lung structure. The current study aimed to examine gene expression profiles of brushed bronchial epithelial cells in association with the PRM-defined CT-based measurements of emphysema (PRM^Emph) and small airway disease (PRM^fSAD). Using the Top Institute Pharma (TIP) study cohort (COPD = 12 and asymptomatic smokers = 32), we identified a gene expression signature of bronchial brushings, which was associated with PRM^Emph in the lungs. One hundred thirty-three genes were identified to be associated with PRM^Emph. Among the most significantly associated genes, CXCL11 is a potent chemokine involved with CD8⁺ T cell activation during inflammation in COPD, indicating that it may play an essential role in the development of emphysema. The PRM^Emph signature was then replicated in two independent data sets. Pathway analysis showed that the PRM^Emph signature is associated with proinflammatory and notch signaling pathways. Together these findings indicate that airway epithelium may play a role in the development of emphysema and/or may act as a biomarker for the presence of emphysema. In contrast, its role in relation to functional small airways disease is less clear.

A protective polymorphism in MMP16, improved blood gas levels, and chronic obstructive pulmonary diseases: Family and two population-based studies

The aberrant expression of matrix metalloproteinases (MMPs) is known to contribute to the pathogenesis of airway remodeling and alveolar disruption in chronic obstructive pulmonary disease (COPD). In the discovery stage, 11 COPD from five families were subjected to whole-genome sequencing, and 21 common polymorphisms in MMPs and TIMPs were identified. These polymorphisms were genotyped in two subsequent verification studies. Of these polymorphisms, c.2392G>A (rs2664370T>C) and c.4158C>A (rs2664369T>G) in MMP16 remained significantly different. Functionally, we found that MMP16 expression was significantly increased in peripheral blood monocytes (PBMCs) from COPD and in cigarette smoke extract-treated 16HBE cells compared with controls. This was also shown by bioinformatics analysis. COPD carrying rs2664370CC showed decreased levels of MMP16 in the plasma and in PBMCs compared with those carrying CT and TT. Treatment with hsa-miR-576-5p mimics led to a greater reduction in luciferase reporter activity in cells transfected with rs2664370CC. Moreover, blood levels of base excess, PCO₂ , and PO₂ in COPD with rs2664370CC were significantly lower than those with rs2664370CT+TT. Taken together, these results demonstrate that the rs2664370T>C polymorphism in MMP16 protects against the risk of COPD, likely by favoring interaction with hsa-miR-576-5p, leading to reduced MMP16 expression and improved blood gas levels.

Joint Transcriptomic Analysis of Lung Cancer and Other Lung Diseases

Background: Epidemiological and clinical evidence points cancer comorbidity with pulmonary chronic disease. The acquisition of some hallmarks of cancer by cells affected with lung pathologies as a cell adaptive mechanism to a shear stress, suggests that could be associated with the establishment of tumoral processes. Objective: To propose a bioinformatic pipeline for the identification of all deregulated genes and the transcriptional regulators (TFs) that are coexpressed during lung cancer establishment, and therefore could be important for the acquisition of the hallmarks of cancer. Methods: Ten microarray datasets (six of lung cancer, four of lung diseases) comparing normal and diseases-related lung tissue were selected to identify hub differentiated expressed genes (DEGs) in common between lung pathologies and lung cancer, along with transcriptional regulators through the utilization of specialized libraries from R language. DAVID bioinformatics tool for gene enrichment analyses was used to identify genes with experimental evidence associated to tumoral processes and signaling pathways. Coexpression networks of DEGs and TFs in lung cancer establishment were created with Coexnet library, and a survival analysis of the main hub genes was made. Results: Two hundred ten DEGs were identified in common between lung cancer and other lung diseases related to the acquisition of tumoral characteristics, which are coexpressed in a lung cancer network with TFs, suggesting that could be related to the establishment of the tumoral pathology in lung. The comparison of the coexpression networks of lung cancer and other lung diseases allowed the identification of common connectivity patterns (CCPs) with DEGs and TFs correlated to important tumoral processes and signaling pathways, that haven´t been studied to experimentally validate their role in the early stages of lung cancer. Some of the TFs identified showed a correlation between its expression levels and the survival of lung cancer patients. Conclusion: Our findings indicate that lung diseases share genes with lung cancer which are coexpressed in lung cancer, and might be able to explain the epidemiological observations that point to direct and inverse comorbid associations between some chronic lung diseases and lung cancer and represent a complex transcriptomic scenario.

Protein phosphatase 2A (PP2A): a key phosphatase in the progression of chronic obstructive pulmonary disease (COPD) to lung cancer

Lung cancer (LC) has the highest relative risk of development as a comorbidity of chronic obstructive pulmonary disease (COPD). The molecular mechanisms that mediate chronic inflammation and lung function impairment in COPD have been identified in LC. This suggests the two diseases are more linked than once thought. Emerging data in relation to a key phosphatase, protein phosphatase 2A (PP2A), and its regulatory role in inflammatory and tumour suppression in both disease settings suggests that it may be critical in the progression of COPD to LC. In this review, we uncover the importance of the functional and active PP2A holoenzyme in the context of both diseases. We describe PP2A inactivation via direct and indirect means and explore the actions of two key PP2A endogenous inhibitors, cancerous inhibitor of PP2A (CIP2A) and inhibitor 2 of PP2A (SET), and the role they play in COPD and LC. We explain how dysregulation of PP2A in COPD creates a favourable inflammatory micro-environment and promotes the initiation and progression of tumour pathogenesis. Finally, we highlight PP2A as a druggable target in the treatment of COPD and LC and demonstrate the potential of PP2A re-activation as a strategy to halt COPD disease progression to LC. Although further studies are required to elucidate if PP2A activity in COPD is a causal link for LC progression, studies focused on the potential of PP2A reactivating agents to reduce the risk of LC formation in COPD patients will be pivotal in improving clinical outcomes for both COPD and LC patients in the future.

Roundabout signaling pathway involved in the pathogenesis of COPD by integrative bioinformatics analysis

Purpose

To explore the potential mechanism underpinning the development of chronic obstructive pulmonary disease (COPD) and to investigate the role of the Roundabout signaling pathway in COPD.

Methods

Three microarray datasets (GSE1650, GSE38974 and GSE76925) including 139 cases of severe COPD and 52 cases of normal smokers without carcinoma, were integrated to screen differentially expressed genes (DEGs) using bioinformatics methods. Gene ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of the DEGs were performed by a DAVID online tool. Finally, a cigarette smoke (CS)- induced emphysema mice model was established, the lung mRNA expression levels of genes associated with Slit guidance ligand 2 (SLIT2) -Roundabout (ROBO) signaling pathway were detected by reverse transcription quantitative polymerase chain reaction (RT-qPCR), and the protein level of SLIT2 was examined by immunohistochemistry staining.

Results

A total of 315 DEGs were identified in three databases. GO and KEGG pathway analyses suggested that the inflammatory response, extracellular matrix disassembly, immune response, the apoptotic signaling pathway, ubiquitination and the Roundabout signaling pathway all together were involved in the development of COPD. The genes SLIT2 and ROBO2 were decreased in patients with COPD and these decreases were significantly negatively correlated with the disease stages of COPD. Consistently, the mRNA expression levels of SLIT2, ROBO1 and ROBO2, and the protein level of SLIT2 were revealed to be lower in the lungs of CS-induced emphysema mice compared with the air-exposed control mice. In addition, the SLIT2 protein level was negatively associated with alveolar mean linear intercept.

Conclusion

Integrated bioinformatics analysis may provide novel insights into the complicated pathogenesis of COPD, and to the best of our knowledge, this study is the first to provide evidence to suggest that the Roundabout signaling pathway may be involved in the pathogenesis of COPD.

Combined Pulmonary Fibrosis and Emphysema: Pulmonary Function Testing and a Pathophysiology Perspective

Combined pulmonary fibrosis and emphysema (CPFE) has been increasingly recognized over the past 10–15 years as a clinical entity characterized by rather severe imaging and gas exchange abnormalities, but often only mild impairment in spirometric and lung volume indices. In this review, we explore the gas exchange and mechanical pathophysiologic abnormalities of pulmonary emphysema, pulmonary fibrosis, and combined emphysema and fibrosis with the goal of understanding how individual pathophysiologic observations in emphysema and fibrosis alone may impact clinical observations on pulmonary function testing (PFT) patterns in patients with CPFE. Lung elastance and lung compliance in patients with CPFE are likely intermediate between those of patients with emphysema and fibrosis alone, suggesting a counter-balancing effect of each individual process. The outcome of combined emphysema and fibrosis results in higher lung volumes overall on PFTs compared to patients with pulmonary fibrosis alone, and the forced expiratory volume in one second (FEV₁)/forced vital capacity (FVC) ratio in CPFE patients is generally preserved despite the presence of emphysema on chest computed tomography (CT) imaging. Conversely, there appears to be an additive deleterious effect on gas exchange properties of the lungs, reflecting a loss of normally functioning alveolar capillary units and effective surface area available for gas exchange, and manifested by a uniformly observed severe reduction in the diffusing capacity for carbon monoxide (D_LCO). Despite normal or only mildly impaired spirometric and lung volume indices, patients with CPFE are often severely functionally impaired with an overall rather poor prognosis. As chest CT imaging continues to be a frequent imaging modality in patients with cardiopulmonary disease, we expect that patients with a combination of pulmonary emphysema and pulmonary fibrosis will continue to be observed. Understanding the pathophysiology of this combined process and the abnormalities that manifest on PFT testing will likely be helpful to clinicians involved with the care of patients with CPFE.

PD-L1, FGFR1, PIK3CA, PTEN, and p16 expression in pulmonary emphysema and chronic obstructive pulmonary disease with resected lung squamous cell carcinoma

BACKGROUND

Emphysema and chronic obstructive pulmonary disease (COPD) are well known independent risk factors for lung cancer. However, the developmental mechanisms between emphysema/COPD and lung cancer remain unknown. The purpose of this study was to evaluate PD-L1, FGFR1, PIK3CA, PTEN, and p16 expression in squamous cell carcinoma (SCC) associated with emphysema/COPD.

METHODS

A total of 59 patients with squamous cell lung carcinoma (SCC) resected between 2008 and 2012 were retrospectively reviewed. Emphysema was assessed according to the Goddard score. Total severity was divided into none-mild (0-7), moderate (8-15), and severe (≥ 16). Local severity around the existing tumor was divided into no emphysema (0) and presence of emphysema (1-4). COPD severity was based on the Global Initiative for Chronic Obstructive Lung Disease (GOLD) criteria. PD-L1, FGFR1, PIK3CA, PTEN, and p16 expression were evaluated by immunohistochemistry (IHC). Expression level was classified as tumor cells (TC) 3 (≥ 50%), TC2 (5-49%), TC1 (1-4%), or TC0 (< 1%), and as tumor-infiltrating immune cells (IC) 3 (≥ 50%), IC2 (5-49%), IC1 (1-4%), or IC0 (< 1%) for PD-L1. Expression level was compared between none-mild/moderate-severe total emphysema, no/presence of local emphysema, no COPD/COPD, and GOLD 1/GOLD 2, 3.

RESULTS

PD-L1 expression was significantly correlated with severity of emphysema in TC0, 1, 2 vs. TC3 (P = 0.012). PD-L1 was significantly higher inversely in none-mild emphysema compared to moderate-severe (95% CI, 0.061-5.852, P = 0.045). There were no other significant associations between PD-L1, FGFR1, PIK3CA, PTEN, and p16 expression and total/local severity of emphysema or presence of COPD/GOLD stage.

CONCLUSIONS

PD-L1 expression in SCC was correlated with severity of emphysema in TC0, 1, 2 vs. TC3 and more frequent in none-mild emphysema than moderate-severe emphysema.

RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD

BACKGROUND

Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants.

METHODS

We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance.

RESULTS

The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies.

CONCLUSIONS

Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.

Identification of novel candidate genes involved in the progression of emphysema by bioinformatic methods

Purpose

By reanalyzing the gene expression profile GSE76925 in the Gene Expression Omnibus database using bioinformatic methods, we attempted to identify novel candidate genes promoting the development of emphysema in patients with COPD.

Patients and methods

According to the Quantitative CT data in GSE76925, patients were divided into mild emphysema group (%LAA-950<20%, n=12) and severe emphysema group (%LAA-950>50%, n=11). Differentially expressed genes (DEGs) were identified using Agilent GeneSpring GX v11.5 (corrected P-value <0.05 and |Fold Change|>1.3). Known driver genes of COPD were acquired by mining literatures and retrieving databases. Direct protein–protein interaction network (PPi) of DEGs and known driver genes was constructed by STRING.org to screen the DEGs directly interacting with driver genes. In addition, we used STRING.org to obtain the first-layer proteins interacting with DEGs’ products and constructed the indirect PPi of these interaction proteins. By merging the indirect PPi with driver genes’ PPi using Cytoscape v3.6.1, we attempted to discover potential pathways promoting emphysema’s development.

Results

All the patients had COPD with severe airflow limitation (age=62±8, FEV₁%=28±12). A total of 57 DEGs (including 12 pseudogenes) and 135 known driving genes were identified. Direct PPi suggested that GPR65, GNB4, P2RY13, NPSR1, BCR, BAG4, and IMPDH2 were potential pathogenic genes. GPR65 could regulate the response of immune cells to the acidic microenvironment, and NPSR1’s expression on eosinophils was associated with asthma’s severity and IgE level. Indirect merging PPi demonstrated that the interacting network of TP53, IL8, CCR2, HSPA1A, ELANE, PIK3CA was associated with the development of emphysema. IL8, ELANE, and PIK3CA were molecules involved in the pathological mechanisms of emphysema, which also in return proved the role of TP53 in emphysema.

Conclusion

Candidate genes such as GPR65, NPSR1, and TP53 may be involved in the progression of emphysema.

A candidate gene identification strategy utilizing mouse to human big-data mining: "3R-tenet" in COPD genetic research

Background

Early life impairments leading to lower lung function by adulthood are considered as risk factors for chronic obstructive pulmonary disease (COPD). Recently, we compared the lung transcriptomic profile between two mouse strains with extreme total lung capacities to identify plausible pulmonary function determining genes using microarray analysis (GSE80078). Advancement of high-throughput techniques like deep sequencing (eg. RNA-seq) and microarray have resulted in an explosion of genomic data in the online public repositories which however remains under-exploited. Strategic curation of publicly available genomic data with a mouse-human translational approach can effectively implement “3R- Tenet” by reducing screening experiments with animals and performing mechanistic studies using physiologically relevant in vitro model systems. Therefore, we sought to analyze the association of functional variations within human orthologs of mouse lung function candidate genes in a publicly available COPD lung RNA-seq data-set.

Methods

Association of missense single nucleotide polymorphisms, insertions, deletions, and splice junction variants were analyzed for susceptibility to COPD using RNA-seq data of a Korean population (GSE57148). Expression of the associated genes were studied using the Gene Paint (mouse embryo) and Human Protein Atlas (normal adult human lung) databases. The genes were also assessed for replication of the associations and expression in COPD−/mouse cigarette smoke exposed lung tissues using other datasets.

Results

Significant association (p <  0.05) of variations in 20 genes to higher COPD susceptibility have been detected within the investigated cohort. Association of HJURP, MCRS1 and TLR8 are novel in relation to COPD. The associated ADAM19 and KIT loci have been reported earlier. The remaining 15 genes have also been previously associated to COPD. Differential transcript expression levels of the associated genes in COPD- and/ or mouse emphysematous lung tissues have been detected.

Conclusion

Our findings suggest strategic mouse-human datamining approaches can identify novel COPD candidate genes using existing datasets in the online repositories. The candidates can be further evaluated for mechanistic role through in vitro studies using appropriate primary cells/cell lines. Functional studies can be limited to transgenic animal models of only well supported candidate genes. This approach will lead to a significant reduction of animal experimentation in respiratory research.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0795-y) contains supplementary material, which is available to authorized users.

COPD as an endothelial disorder: endothelial injury linking lesions in the lungs and other organs? (2017 Grover Conference Series)

Chronic obstructive pulmonary disease (COPD) is characterized by chronic expiratory airflow obstruction that is not fully reversible. COPD patients develop varying degrees of emphysema, small and large airway disease, and various co-morbidities. It has not been clear whether these co-morbidities share common underlying pathogenic processes with the pulmonary lesions. Early research into the pathogenesis of COPD focused on the contributions of injury to the extracellular matrix and pulmonary epithelial cells. More recently, cigarette smoke-induced endothelial dysfunction/injury have been linked to the pulmonary lesions in COPD (especially emphysema) and systemic co-morbidities including atherosclerosis, pulmonary hypertension, and chronic renal injury. Herein, we review the evidence linking endothelial injury to COPD, and the pathways underlying endothelial injury and the "vascular COPD phenotype" including: (1) direct toxic effects of cigarette smoke on endothelial cells; (2) generation of auto-antibodies directed against endothelial cells; (3) vascular inflammation; (4) increased oxidative stress levels in vessels inducing increases in lipid peroxidation and increased activation of the receptor for advanced glycation end-products (RAGE); (5) reduced activation of the anti-oxidant pathways in endothelial cells; (6) increased endothelial cell release of mediators with vasoconstrictor, pro-inflammatory, and remodeling activities (endothelin-1) and reduced endothelial cell expression of mediators that promote vasodilation and homeostasis of endothelial cells (nitric oxide synthase and prostacyclin); and (7) increased endoplasmic reticular stress and the unfolded protein response in endothelial cells. We also review the literature on studies of drugs that inhibit RAGE signaling in other diseases (angiotensin-converting enzyme inhibitors and angiotensin receptor blockers), or vasodilators developed for idiopathic pulmonary arterial hypertension that have been tested on cell culture systems, animal models of COPD, and/or smokers and COPD patients.

Meta-analysis of peripheral blood gene expression modules for COPD phenotypes

Chronic obstructive pulmonary disease (COPD) occurs typically in current or former smokers, but only a minority of people with smoking history develops the disease. Besides environmental factors, genetics is an important risk factor for COPD. However, the relationship between genetics, environment and phenotypes is not well understood. Sample sizes for genome-wide expression studies based on lung tissue have been small due to the invasive nature of sample collection. Increasing evidence for the systemic nature of the disease makes blood a good alternative source to study the disease, but there have also been few large-scale blood genomic studies in COPD. Due to the complexity and heterogeneity of COPD, examining groups of interacting genes may have more relevance than identifying individual genes. Therefore, we used Weighted Gene Co-expression Network Analysis to find groups of genes (modules) that are highly connected. However, module definitions may vary between individual data sets. To alleviate this problem, we used a consensus module definition based on two cohorts, COPDGene and ECLIPSE. We studied the relationship between the consensus modules and COPD phenotypes airflow obstruction and emphysema. We also used these consensus module definitions on an independent cohort (TESRA) and performed a meta analysis involving all data sets. We found several modules that are associated with COPD phenotypes, are enriched in functional categories and are overrepresented for cell-type specific genes. Of the 14 consensus modules, three were strongly associated with airflow obstruction (meta p ≤ 0.0002), and two had some association with emphysema (meta p ≤ 0.06); some associations were stronger in the case-control cohorts, and others in the cases-only subcohorts. Gene Ontology terms that were overrepresented included “immune response” and “defense response.” The cell types whose type-specific genes were overrepresented in modules (p < 0.05) included natural killer cells, dendritic cells, and neutrophils. Together, this is the largest investigation of gene blood expression in COPD with 469 cases in COPDGene, ECLIPSE and TESRA combined, with 6267 genes common to all data sets. Additional, we have 42 and 83 controls in COPDGene and ECLIPSE, respectively.

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.

Tissue expression and source of circulating αKlotho

αKlotho (Klotho), a type I transmembrane protein and a coreceptor for Fibroblast Growth Factor-23, was initially thought to be expressed only in a limited number of tissues, most importantly the kidney, parathyroid gland and choroid plexus. Emerging data may suggest a more ubiquitous Klotho expression pattern which has prompted reevaluation of the restricted Klotho paradigm. Herein we systematically review the evidence for Klotho expression in various tissues and cell types in humans and other mammals, and discuss potential reasons behind existing conflicting data. Based on current literature and tissue expression atlases, we propose a classification of tissues into high, intermediate and low/absent Klotho expression. The functional relevance of Klotho in organs with low expression levels remain uncertain and there is currently limited data on a role for membrane-bound Klotho outside the kidney. Finally, we review the evidence for the tissue source of soluble Klotho, and conclude that the kidney is likely to be the principal source of circulating Klotho in physiology.

Network Analysis of Lung Transcriptomics Reveals a Distinct B-Cell Signature in Emphysema

RATIONALE

Chronic obstructive pulmonary disease (COPD) is characterized by chronic airflow limitation caused by a combination of airways disease (bronchiolitis) and parenchymal destruction (emphysema), whose relative proportion varies from patient to patient.

OBJECTIVES

To explore and contrast the molecular pathogenesis of emphysema and bronchiolitis in COPD.

METHODS

We used network analysis of lung transcriptomics (Affymetrix arrays) in 70 former smokers with COPD to compare differential expression and gene coexpression in bronchiolitis and emphysema.

MEASUREMENTS AND MAIN RESULTS

We observed that in emphysema (but not in bronchiolitis) (1) up-regulated genes were enriched in ontologies related to B-cell homing and activation; (2) the immune coexpression network had a central core of B cell-related genes; (3) B-cell recruitment and immunoglobulin transcription genes (CXCL13, CCL19, and POU2AF1) correlated with emphysema severity; (4) there were lymphoid follicles (CD20(+)IgM(+)) with active B cells (phosphorylated nuclear factor-κB p65(+)), proliferation markers (Ki-67(+)), and class-switched B cells (IgG(+)); and (5) both TNFRSF17 mRNA and B cell-activating factor protein were up-regulated. These findings were by and large reproduced in a group of patients with incipient emphysema and when patients with emphysema were matched for the severity of airflow limitation of those with bronchiolitis.

CONCLUSIONS

Our study identifies enrichment in B cell-related genes in patients with COPD with emphysema that is absent in bronchiolitis. These observations contribute to a better understanding of COPD pathobiology and may open new therapeutic opportunities for patients with COPD.

Functional interactors of three genome-wide association study genes are differentially expressed in severe chronic obstructive pulmonary disease lung tissue

In comparison to genome-wide association studies (GWAS), there has been poor replication of gene expression studies in chronic obstructive pulmonary disease (COPD). We performed microarray gene expression profiling on a large sample of resected lung tissues from subjects with severe COPD. Comparing 111 COPD cases and 40 control smokers, 204 genes were differentially expressed; none were at significant GWAS loci. The top differentially expressed gene was HMGB1, which interacts with AGER, a known COPD GWAS gene. Differentially expressed genes showed enrichment for putative interactors of the first three identified COPD GWAS genes IREB2, HHIP, and FAM13A, based on gene sets derived from protein and RNA binding studies, RNA-interference, a murine smoking model, and expression quantitative trait locus analyses. The gene module most highly associated for COPD in Weighted Gene Co-Expression Network Analysis (WGCNA) was enriched for B cell pathways, and shared seventeen genes with a mouse smoking model and twenty genes with previous emphysema studies. As in other common diseases, genes at COPD GWAS loci were not differentially expressed; however, using a combination of network methods, experimental studies and careful phenotype definition, we found differential expression of putative interactors of these genes, and we replicated previous human and mouse microarray results.

The role of the endothelium in asthma and chronic obstructive pulmonary disease (COPD)

COPD and asthma are important chronic inflammatory disorders with a high associated morbidity. Much research has concentrated on the role of inflammatory cells, such as the neutrophil, in these diseases, but relatively little focus has been given to the endothelial tissue, through which inflammatory cells must transmigrate to reach the lung parenchyma and cause damage. There is evidence that there is an abnormal amount of endothelial tissue in COPD and asthma and that this tissue and its’ progenitor cells behave in a dysfunctional manner. This article reviews the evidence of the involvement of pulmonary endothelium in COPD and asthma and potential treatment options for this.

META-ANALYSIS OF CONTINUOUS PHENOTYPES IDENTIFIES A GENE SIGNATURE THAT CORRELATES WITH COPD DISEASE STATUS

The utility of multi-cohort two-class meta-analysis to identify robust differentially expressed gene signatures has been well established. However, many biomedical applications, such as gene signatures of disease progression, require one-class analysis. Here we describe an R package, MetaCorrelator, that can identify a reproducible transcriptional signature that is correlated with a continuous disease phenotype across multiple datasets. We successfully applied this framework to extract a pattern of gene expression that can predict lung function in patients with chronic obstructive pulmonary disease (COPD) in both peripheral blood mononuclear cells (PBMCs) and tissue. Our results point to a disregulation in the oxidation state of the lungs of patients with COPD, as well as underscore the classically recognized inammatory state that underlies this disease.

RNAi Screen for NRF2 Inducers Identifies Targets That Rescue Primary Lung Epithelial Cells from Cigarette Smoke Induced Radical Stress

Chronic Obstructive Pulmonary Disease (COPD) is a highly prevalent condition characterized by inflammation and progressive obstruction of the airways. At present, there is no treatment that suppresses the chronic inflammation of the disease, and COPD patients often succumb to the condition. Excessive oxidative stress caused by smoke inhalation is a major driving force of the disease. The transcription factor NRF2 is a critical player in the battle against oxidative stress and its function is impaired in COPD. Increasing NRF2 activity may therefore be a viable therapeutic option for COPD treatment. We show that down regulation of KEAP1, a NRF2 inhibitor, protects primary human lung epithelial cells from cigarette-smoke-extract (CSE) induced cell death in an established in vitro model of radical stress. To identify new potential drug targets with a similar effect, we performed a siRNA screen of the ‘druggable’ genome using a NRF2 transcriptional reporter cell line. This screen identified multiple genes that when down regulated increased NRF2 transcriptional activity and provided a survival benefit in the in vitro model. Our results suggest that inhibiting components of the ubiquitin-proteasome system will have the strongest effects on NRF2 transcriptional activity by increasing NRF2 levels. We also find that down regulation of the small GTPase Rab28 or the Estrogen Receptor ESRRA provide a survival benefit. Rab28 knockdown increased NRF2 protein levels, indicating that Rab28 may regulate NRF2 proteolysis. Conversely ESRRA down regulation increased NRF2 transcriptional activity without affecting NRF2 levels, suggesting a proteasome-independent mechanism.

Integrated Genomics Reveals Convergent Transcriptomic Networks Underlying Chronic Obstructive Pulmonary Disease and Idiopathic Pulmonary Fibrosis

RATIONALE

Despite shared environmental exposures, idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary disease are usually studied in isolation, and the presence of shared molecular mechanisms is unknown.

OBJECTIVES

We applied an integrative genomic approach to identify convergent transcriptomic pathways in emphysema and IPF.

METHODS

We defined the transcriptional repertoire of chronic obstructive pulmonary disease, IPF, or normal histology lungs using RNA-seq (n = 87).

MEASUREMENTS AND MAIN RESULTS

Genes increased in both emphysema and IPF relative to control were enriched for the p53/hypoxia pathway, a finding confirmed in an independent cohort using both gene expression arrays and the nCounter Analysis System (n = 193). Immunohistochemistry confirmed overexpression of HIF1A, MDM2, and NFKBIB members of this pathway in tissues from patients with emphysema or IPF. Using reads aligned across splice junctions, we determined that alternative splicing of p53/hypoxia pathway-associated molecules NUMB and PDGFA occurred more frequently in IPF or emphysema compared with control and validated these findings by quantitative polymerase chain reaction and the nCounter Analysis System on an independent sample set (n = 193). Finally, by integrating parallel microRNA and mRNA-Seq data on the same samples, we identified MIR96 as a key novel regulatory hub in the p53/hypoxia gene-expression network and confirmed that modulation of MIR96 in vitro recapitulates the disease-associated gene-expression network.

CONCLUSIONS

Our results suggest convergent transcriptional regulatory hubs in diseases as varied phenotypically as chronic obstructive pulmonary disease and IPF and suggest that these hubs may represent shared key responses of the lung to environmental stresses.

The Relevance of Genomic Signatures at Adhesion GPCR Loci in Humans

Adhesion G protein-coupled receptors (aGPCRs) have a long evolutionary history dating back to very basal unicellular eukaryotes. Almost every vertebrate is equipped with a set of different aGPCRs. Genomic sequence data of several hundred extinct and extant species allows for reconstruction of aGPCR phylogeny in vertebrates and non-vertebrates in general but also provides a detailed view into the recent evolutionary history of human aGPCRs. Mining these sequence sources with bioinformatic tools can unveil many facets of formerly unappreciated aGPCR functions. In this review, we extracted such information from the literature and open public sources and provide insights into the history of aGPCR in humans. This includes comprehensive analyses of signatures of selection, variability of human aGPCR genes, and quantitative traits at human aGPCR loci. As indicated by a large number of genome-wide genotype-phenotype association studies, variations in aGPCR contribute to specific human phenotypes. Our survey demonstrates that aGPCRs are significantly involved in adaptation processes, phenotype variations, and diseases in humans.

Defining the nasal transcriptome in granulomatosis with polyangiitis (Wegener's)

OBJECTIVE

To determine whether disease processes related to granulomatosis with polyangiitis (Wegener's) (GPA) are reflected in gene expression profiles of the nasal mucosa.

METHODS

Nasal brushings of the inferior turbinate were obtained from 32 patients with GPA (10 with active nasal disease, 13 with prior nasal disease, and 9 with no history of nasal disease) and a composite comparator group with and without inflammatory nasal disease (12 healthy people, 15 with sarcoidosis, and 8 with allergic rhinitis). Differential gene expression was assessed between subgroups of GPA and comparators.

RESULTS

A total of 339 genes were differentially expressed between the GPA and comparator groups (absolute fold change >1.5; false discovery rate <0.05). Top canonical pathways up-regulated in nasal brushings from patients with GPA included granulocyte adhesion and diapedesis (P = 8.6(-22) ), agranulocyte adhesion and diapedesis (P = 1.3(-14) ), IL10 signaling (P = 3.0(-11) ), LXR/RXR activation (P = 4.3(-11) ), and TREM1 signaling (P = 9.0(-11) ). A set of genes differentially expressed in GPA independently of nasal disease activity status included genes related to epithelial barrier integrity (fibronectin 1, desmosomal proteins) and several matricellular proteins (e.g., osteonectin, osteopontin). Significant overlap of differentially expressed genes was observed between active and prior nasal disease GPA subgroups. Peripheral blood neutrophil and mononuclear gene expression levels associated with GPA were similarly altered in the nasal gene expression profiles of patients with active or prior nasal disease.

CONCLUSION

Profiling the nasal transcriptome in GPA reveals gene expression signatures related to innate immunity, inflammatory cell chemotaxis, extracellular matrix composition, and epithelial barrier integrity. Thus, airway-based expression profiling is feasible and informative in GPA.

Identification of MMP-9 as a biomarker for detecting progression of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a complex immunological disease with multiple pathological features that is primarily induced by smoking together with additional genetic risk factors. COPD is frequently underdiagnosed; forced expiratory volume in the first second (FEV1) is considered to be the main diagnostic measure for COPD, yet it is insufficiently sensitive to monitor disease progression. Biomarkers capable of monitoring COPD progression and severity are needed. In this report, we evaluated matrix metalloproteinase-9 (MMP-9) as an early marker for the detection and staging of COPD, by assessing the mRNA levels of MMP-9 in peripheral blood samples collected from 22 COPD patients, 6 asymptomatic smokers, and 5 healthy controls. Our results demonstrate that the mRNA levels of MMP-9 increased more than two-fold in severe COPD relative to non-COPD smokers or moderate COPD groups. Moreover, in the very severe COPD group, MMP-9 mRNA levels showed a 4-fold increase relative to the non-COPD smokers or the moderate COPD groups, while there was a mild increase (∼40%) when compared to the severe COPD group. Taken together, our results suggest that MMP-9 serves as a biomarker for the grade and severity of COPD.

A 7-month cigarette smoke inhalation study in C57BL/6 mice demonstrates reduced lung inflammation and emphysema following smoking cessation or aerosol exposure from a prototypic modified risk tobacco product

Modified risk tobacco products (MRTP) are designed to reduce smoking-related health risks. A murine model of chronic obstructive pulmonary disease (COPD) was applied to investigate classical toxicology end points plus systems toxicology (transcriptomics and proteomics). C57BL/6 mice were exposed to conventional cigarette smoke (3R4F), fresh air (sham), or a prototypic MRTP (pMRTP) aerosol for up to 7 months, including a cessation group and a switching-to-pMRTP group (2 months of 3R4F exposure followed by fresh air or pMRTP for up to 5 months respectively). 3R4F smoke induced the typical adaptive changes in the airways, as well as inflammation in the lung, associated with emphysematous changes (impaired pulmonary function and alveolar damage). At nicotine-matched exposure concentrations of pMRTP aerosol, no signs of lung inflammation and emphysema were observed. Both the cessation and switching groups showed a similar reversal of inflammatory responses and no progression of initial emphysematous changes. A significant impact on biological processes, including COPD-related inflammation, apoptosis, and proliferation, was identified in 3R4F-exposed, but not in pMRTP-exposed lungs. Smoking cessation or switching reduced these perturbations to near sham-exposed levels. In conclusion, the mouse model indicated retarded disease progression upon cessation or switching to pMRTP which alone had no adverse effects.

Comprehensive Analysis of Transcriptome Sequencing Data in the Lung Tissues of COPD Subjects

Background and Objectives. Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airflow limitation. Although airway inflammation and oxidative stress are known to be important in the pathogenesis of COPD, the mechanism underlying airflow obstruction is not fully understood. Gene expression profiling of lung tissue was performed to define the molecular pathways that are dysregulated in COPD. Methods. RNA was isolated from lung tissues obtained from 98 subjects with COPD and 91 control subjects with normal spirometry. The RNA samples were processed with RNA-seq using the HiSeq 2000 system. Genes expressed differentially between the two groups were identified using Student's t-test. Results. After filtering for genes with zero counts and noncoding genes, 16,676 genes were evaluated. A total of 2312 genes were differentially expressed between the lung tissues of COPD and control subjects (false discovery rate corrected q < 0.01). The expression of genes related to oxidative phosphorylation and protein catabolism was reduced and genes related to chromatin modification were dysregulated in lung tissues of COPD subjects. Conclusions. Oxidative phosphorylation, protein degradation, and chromatin modification were the most dysregulated pathways in the lung tissues of COPD subjects. These findings may have clinical and mechanistic implications in COPD.

Protocol for a human in vivo model of acute cigarette smoke inhalation challenge in smokers with COPD: monitoring the nasal and systemic immune response using a network biology approach

INTRODUCTION

Cigarette smoke contributes to a diverse range of diseases including chronic obstructive pulmonary disease (COPD), cardiovascular disorders and many cancers. There currently is a need for human challenge models, to assess the acute effects of a controlled cigarette smoke stimulus, followed by serial sampling of blood and respiratory tissue for advanced molecular profiling. We employ precision sampling of nasal mucosal lining fluid by absorption to permit soluble mediators measurement in eluates. Serial nasal curettage was used for transcriptomic analysis of mucosal tissue.

METHODS AND ANALYSIS

Three groups of strictly defined patients will be studied: 12 smokers with COPD (GOLD Stage 2) with emphysema, 12 matched smokers with normal lung function and no evidence of emphysema, and 12 matched never smokers with normal spirometry. Patients in the smoking groups are current smokers, and will be given full support to stop smoking immediately after this study. In giving a controlled cigarette smoke stimulus, all patients will have abstained from smoking for 12 h, and will smoke two cigarettes with expiration through the nose in a ventilated chamber. Before and after inhalation of cigarette smoke, a series of samples will be taken from the blood, nasal mucosal lining fluid and nasal tissue by curettage. Analysis of plasma nicotine and metabolites in relation to levels of soluble inflammatory mediators in nasal lining fluid and blood, as well as assessing nasal transcriptomics, ex vivo blood platelet aggregation and leucocyte responses to toll-like receptor agonists will be undertaken.

IMPLICATIONS

Development of acute cigarette smoke challenge models has promise for the study of molecular effects of smoking in a range of pathological processes.

ETHICS AND DISSEMINATION

This study was approved by the West London National Research Ethics Committee (12/LO/1101). The study findings will be presented at conferences and will be reported in peer-reviewed journals.

Genes related to emphysema are enriched for ubiquitination pathways

BACKGROUND

Increased small airway resistance and decreased lung elasticity contribute to the airflow limitation in chronic obstructive pulmonary disease (COPD). The lesion that corresponds to loss of lung elasticity is emphysema; the small airway obstruction is due to inflammatory narrowing and obliteration. Despite their convergence in altered physiology, different mechanisms contribute to these processes. The relationships between gene expression and these specific phenotypes may be more revealing than comparison with lung function.

METHODS

We measured the ratio of alveolar surface area to lung volume (SA/V) in lung tissue from 43 smokers. Two samples from 21 subjects, in which SA/V differed by >49 cm2/mL were profiled to select genes whose expression correlated with SA/V. Significant genes were tested for replication in the 22 remaining subjects.

RESULTS

The level of expression of 181 transcripts was related to SA/V ( p < 0.05). When these genes were tested in the 22 remaining subjects as a replication, thirty of the 181 genes remained significantly associated with SA/V (P < 0.05) and the direction of association was the same in 164/181. Pathway and network analysis revealed enrichment of genes involved in protein ubiquitination, and western blotting showed altered expression of genes involved in protein ubiquitination in obstructed individuals.

CONCLUSION

This study implicates modified protein ubiquitination and degradation as a potentially important pathway in the pathogenesis of emphysema.

Dynamic gene expressions of peripheral blood mononuclear cells in patients with acute exacerbation of chronic obstructive pulmonary disease: a preliminary study

INTRODUCTION

Acute exacerbation of chronic obstructive pulmonary disease (AECOPD) is a serious event that is responsible for the progress of the disease, increases in medical costs and high mortality.

METHODS

The aim of the present study was to identify AECOPD-specific biomarkers by evaluating the dynamic gene expression profiling of peripheral blood mononuclear cells (PBMCs) from patients with AECOPD on days 1, 3 and 10 after hospital admission and to compare the derived data with data from healthy controls or patients with stable COPD.

RESULTS

We found that 14 genes were co-differentially upregulated and 2 downregulated greater than 10-fold in patients with COPD or AECOPD compared with the healthy individuals. Eight co-differentially upregulated genes and six downregulated genes were identified as a panel of AECOPD-specific genes. Downregulation of TCF7 in PBMCs was found to be associated with the severity of COPD. Dynamic changes of Aminolevulinate-delta-synthase 2 and carbonic anhydrase I had similar patterns of Digital Evaluation Score System scores and may serve as potential genes of interest during the course of AECOPD.

CONCLUSION

Thus, our findings indicate a panel of altered gene expression patterns in PBMCs that can be used as AECOPD-specific dynamic biomarkers to monitor the course of AECOPD.

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.

Overexpression of type VI collagen in neoplastic lung tissues

Type VI collagen (COL6), an extracellular matrix protein, is important in maintaining the integrity of lung tissue. An increase in COL6 mRNA and protein deposition was found in the lungs of patients with pulmonary fibrosis, a chronic inflammatory condition with a strong association with lung cancer. In the present study, we demonstrated overexpression of COL6 in the lungs of non-small cell lung cancers. We hypothesized that excessive COL6 in the lung interstitium may exert stimulatory effects on the adjacent cells. In vitro stimulation of monocytes with COL6 resulted in the production of IL-23, which may promote tumor development in an environment of IL-23-mediated lung inflammation, where tissue modeling occurs concurrently with excessive COL6 production. In addition, COL6 was capable of stimulating signaling pathways that activate focal adhesion kinase and extracellular signal-regulated kinase 1/2 in lung epithelial cells, which may also facilitate the development of lung neoplasms. Taken together, our data suggest the potential role of COL6 in promoting lung neoplasia in diseased lungs where COL6 is overexpressed.

A large lung gene expression study identifying fibulin-5 as a novel player in tissue repair in COPD

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a progressive, incurable lung disease characterised by abnormal tissue repair causing emphysema and small airways fibrosis. Since current therapy cannot modify this abnormal repair, it is crucial to unravel its underlying molecular mechanisms. Unbiased analysis of genome-wide gene expression profiles in lung tissue provides a powerful tool to investigate this.

METHODS

We performed genome-wide gene expression profiling in 581 lung tissue samples from current and ex-smokers with (n=311) and without COPD (n=270). Subsequently, quantitative PCR, western blot and immunohistochemical analyses were performed to validate our main findings.

RESULTS

112 genes were found to be upregulated in patients with COPD compared with controls, whereas 61 genes were downregulated. Among the most upregulated genes were fibulin-5 (FBLN5), elastin (ELN), latent transforming growth factor β binding protein 2 (LTBP2) and microfibrillar associated protein 4 (MFAP4), all implicated in elastogenesis. Our gene expression findings were validated at mRNA and protein level. We demonstrated higher ELN gene expression in COPD lung tissue and similar trends for FBLN5 and MFAP4, and negative correlations with lung function. FBLN5 protein levels were increased in COPD lung tissue and cleaved, possibly non-functional FBLN5 protein was present. Strong coexpression of FBLN5, ELN, LTBP2 and MFAP4 in lung tissue and in silico analysis indicated cofunctionality of these genes. Finally, colocalisation of FBLN5, MFAP4 and LTBP2 with elastic fibres was demonstrated in lung tissue.

CONCLUSIONS

We identified a clear gene signature for elastogenesis in COPD and propose FBLN5 as a novel player in tissue repair in COPD.