Do sputum or circulating blood samples reflect the pulmonary transcriptomic differences of COPD patients? A multi-tissue transcriptomic network META-analysis

Partial correlation network analysis identifies coordinated gene expression within a regional cluster of COPD genome-wide association signals

Chronic obstructive pulmonary disease (COPD) is a complex disease influenced by well-established environmental exposures (most notably, cigarette smoking) and incompletely defined genetic factors. The chromosome 4q region harbors multiple genetic risk loci for COPD, including signals near HHIP, FAM13A, GSTCD, TET2, and BTC. Leveraging RNA-Seq data from lung tissue in COPD cases and controls, we estimated the co-expression network for genes in the 4q region bounded by HHIP and BTC (~70MB), through partial correlations informed by protein-protein interactions. We identified several co-expressed gene pairs based on partial correlations, including NPNT-HHIP, BTC-NPNT and FAM13A-TET2, which were replicated in independent lung tissue cohorts. Upon clustering the co-expression network, we observed that four genes previously associated to COPD: BTC, HHIP, NPNT and PPM1K appeared in the same network community. Finally, we discovered a sub-network of genes differentially co-expressed between COPD vs controls (including FAM13A, PPA2, PPM1K and TET2). Many of these genes were previously implicated in cell-based knock-out experiments, including the knocking out of SPP1 which belongs to the same genomic region and could be a potential local key regulatory gene. These analyses identify chromosome 4q as a region enriched for COPD genetic susceptibility and differential co-expression.

Transcriptomic analysis reveals distinct effects of cigarette smoke on murine airspace and bone-marrow derived macrophages

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is an inflammatory airway disease characterized by emphysema and chronic bronchitis and a leading cause of mortality worldwide. COPD is commonly associated with several comorbid diseases which contribute to exacerbated patient outcomes. Cigarette smoke (CS) is the most prominent risk factor for COPD development and progression and is known to be detrimental to numerous effector functions of lung resident immune cells, including phagocytosis and cytokine production. However, how CS mediates the various pathologies distant from the lung in COPD, and whether CS has a similar biological effect on systemic immune cells remains unknown.

METHODS

C57BL/6 mice were exposed to 8 weeks of CS as an experimental model of COPD. Bone marrow cells were isolated from both CS-exposed and room air (RA) control mice and differentiated to bone marrow-derived macrophages (BMDMs). Airspace macrophages (AMs) were isolated from the same CS-exposed and RA mice and bulk RNA-Seq performed. The functional role of differentially expressed genes was assessed through gene ontology analyses. Ingenuity Pathway Analysis was used to determine the activation states of canonical pathways and upstream regulators enriched in differentially expressed genes in both cell types, and to compare the differences between the two cell types.

RESULTS

CS induced transcriptomic changes in BMDMs, including an upregulation of genes in sirtuin signalling and oxidative phosphorylation pathways and a downregulation of genes involved in histone and lysine methylation. In contrast, CS induced decreased expression of genes involved in pathogen response, phagosome formation, and immune cell trafficking in AMs. Little overlap was observed in differentially expressed protein-coding genes in BMDMs compared to AMs and their associated pathways, highlighting the distinct effects of CS on immune cells in different compartments.

CONCLUSIONS

CS exposure can induce transcriptomic remodelling in BMDMs which is distinct to that of AMs. Our study highlights the ability of CS exposure to affect immune cell populations distal to the lung and warrants further investigation into the functional effects of these changes and the ensuing role in driving multimorbid disease.

Transcriptomic characterization of the human segmental endotoxin challenge model

Segmental instillation of lipopolysaccharide (LPS) by bronchoscopy safely induces transient airway inflammation in human lungs. This model enables investigation of pulmonary inflammatory mechanisms as well as pharmacodynamic analysis of investigational drugs. The aim of this work was to describe the transcriptomic profile of human segmental LPS challenge with contextualization to major respiratory diseases. Pre-challenge bronchoalveolar lavage (BAL) fluid and biopsies were sampled from 28 smoking, healthy participants, followed by segmental instillation of LPS and saline as control. Twenty-four hours post instillation, BAL and biopsies were collected from challenged lung segments. Total RNA of cells from BAL and biopsy samples were sequenced and analysed for differentially expressed genes (DEGs). After challenge with LPS compared with saline, 6316 DEGs were upregulated and 241 were downregulated in BAL, but only one DEG was downregulated in biopsy samples. Upregulated DEGs in BAL were related to molecular functions such as "Inflammatory response" or "chemokine receptor activity", and upregulated pro-inflammatory pathways such as "Wnt-"/"Ras-"/"JAK-STAT" "-signaling pathway". Furthermore, the segmental LPS challenge model resembled aspects of the five most prevalent respiratory diseases chronic obstructive pulmonary disease (COPD), asthma, pneumonia, tuberculosis and lung cancer and featured similarities with acute exacerbations in COPD (AECOPD) and community-acquired pneumonia. Overall, our study provides extensive information about the transcriptomic profile from BAL cells and mucosal biopsies following LPS challenge in healthy smokers. It expands the knowledge about the LPS challenge model providing potential overlap with respiratory diseases in general and infection-triggered respiratory insults such as AECOPD in particular.

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.

Recent evidence from omic analysis for redox signalling and mitochondrial oxidative stress in COPD

COPD is driven by exogenous and endogenous oxidative stress derived from inhaled cigarette smoke, air pollution and reactive oxygen species from dysregulated mitochondria in activated inflammatory cells within the airway and lung. This is compounded by the loss in antioxidant defences including FOXO and NRF2 and other antioxidant transcription factors together with various key enzymes that attenuate oxidant effects. Oxidative stress enhances inflammation; airway remodelling including fibrosis and emphysema; post-translational protein modifications leading to autoantibody generation; DNA damage and cellular senescence. Recent studies using various omics technologies in the airways, lungs and blood of COPD patients has emphasised the importance of oxidative stress, particularly that derived from dysfunctional mitochondria in COPD and its role in immunity, inflammation, mucosal barrier function and infection. Therapeutic interventions targeting oxidative stress should overcome the deleterious pathologic effects of COPD if targeted to the lung. We require novel, more efficacious antioxidant COPD treatments among which mitochondria-targeted antioxidants and Nrf2 activators are promising.

Clinical characteristics and gene mutation profiles of chronic obstructive pulmonary disease in non-small cell lung cancer

Purpose

The coexistence of chronic obstructive pulmonary disease (COPD) often leads to a worse prognosis in patients with non-small cell lung cancer (NSCLC). Meanwhile, approaches targeting specific genetic alterations have been shown to significantly improve the diagnosis and treatment outcomes of patients with NSCLC. Herein, we sought to evaluate the impact of COPD on the clinical manifestations and gene mutation profiles of NSCLC patients with both circulating tumor (ctDNA) and tumor DNA (tDNA).

Materials and methods

The influence of COPD on clinical features was observed in 285 NSCLC cohorts suffering from NSCLC alone, NSCLC coexisting with COPD, or NSCLC coexisting with prodromal changes in COPD (with emphysema, bullae, or chronic bronchitis). The gene mutation profiles of specific 168 NSCLC-related genes were further analyzed in the NSCLC sub-cohorts with formalin-fixed and paraffin-embedded tumor DNA (FFPE tDNA) samples and plasma circulating tumor DNA (PLA ctDNA) samples. Moreover, mutation concordance was assessed in tDNA and paired ctDNA of 110 NSCLC patients.

Results

Relative to patients with NSCLC alone, patients with NSCLC coexisting with COPD and prodromal changes presented with worse lung functions, more clinical symptoms, signs and comorbidities, and inconsistent gene mutation profiles. In addition, patients in the latter two groups exhibited a higher average frequency of gene mutation. Lastly, mutation concordance between tDNA and ctDNA samples was significantly reduced in NSCLC patients coexisting with COPD.

Conclusions

Collectively, our findings revealed that coexistence of COPD leads to worse clinical manifestations and altered gene mutation profiles in patients with NSCLC. Additionally, for NSCLC patients with COPD, the use of ctDNA instead of tDNA may not be the most efficient approach to identifying gene mutations.

COPD is Associated with Epigenome-wide Differential Methylation in BAL Lung Cells

DNA methylation patterns in chronic pulmonary obstructive disease (COPD) might offer new insights into disease pathogenesis. To assess methylation profiles in the main COPD target organ, we performed an epigenome-wide association study on bronchoalveolar lavage (BAL) cells. Bronchoscopies were performed in 18 COPD subjects and 15 controls (ex- and current smokers). DNA methylation was measured with Illumina MethylationEPIC BeadChip covering >850,000 CpGs. Differentially methylated positions (DMPs) were examined for 1) enrichment in pathways and functional gene relationships using Kyoto Encyclopedia of Genes and Genomes and Gene Ontology; 2) accelerated aging using Horvath's epigenetic clock; 3) correlation with gene expression; and 4) co-localization with genetic variation. We found 1,155 Bonferroni significant (P < 6.74 × 10^-8) DMPs associated with COPD, many with large effect sizes. Functional analysis identified biologically plausible pathways and gene relationships, including enrichment for transcription factor activity. Strong correlation was found between COPD and chronological age, but not with accelerated epigenetic aging. For 79 unique DMPs, DNA methylation correlated significantly with gene expression in BAL cells. Thirty-nine percent of DMPs were co-localized with COPD-associated SNPs. To the best of our knowledge, this is the first EWAS of COPD on BAL cells, and our analyses revealed many differential methylation sites. Integration with mRNA data showed a strong functional readout for relevant genes, identifying sites where DNA methylation might directly impact expression. Almost half of DMPs were co-located with SNPs identified in previous GWAS of COPD, suggesting joint genetic and epigenetic pathways related to disease.

Expression Profiling Suggests Loss of Surface Integrity and Failure of Regenerative Repair as Major Driving Forces for Chronic Obstructive Pulmonary Disease Progression

Chronic obstructive pulmonary disease (COPD) poses a major risk for public health, yet remarkably little is known about its detailed pathophysiology. Definition of COPD as nonreversible pulmonary obstruction revealing more about spatial orientation than about mechanisms of pathology may be a major reason for this. We conducted a controlled observational study allowing for simultaneous assessment of clinical and biological development in COPD. Sixteen healthy control subjects and 104 subjects with chronic bronchitis, with or without pulmonary obstruction at baseline, were investigated. Using both the extent of and change in bronchial obstruction as main scoring criteria for the analysis of gene expression in lung tissue, we identified 410 genes significantly associated with progression of COPD. One hundred ten of these genes demonstrated a distinctive expression pattern, with their functional annotations indicating participation in the regulation of cellular coherence, membrane integrity, growth, and differentiation, as well as inflammation and fibroproliferative repair. The regulatory pattern indicates a sequentially unfolding pathology that centers on a two-step failure of surface integrity commencing with a loss of epithelial coherence as early as chronic bronchitis. Decline of regenerative repair starting in Global Initiative for Chronic Obstructive Lung Disease stage I then activates degradation of extracellular-matrix hyaluronan, causing structural failure of the bronchial wall that is only resolved by scar formation. Although they require independent confirmation, our findings provide the first tangible pathophysiological concept of COPD to be further explored.Clinical trial registered with www.clinicaltrials.gov (NCT00618137).

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.

Inhaled Steroids, Circulating Eosinophils, Chronic Airway Infection, and Pneumonia Risk in Chronic Obstructive Pulmonary Disease. A Network Analysis

Rationale: Treatment of chronic obstructive pulmonary disease (COPD) with inhaled corticosteroids (ICS) is controversial, because it can reduce the risk of future exacerbations of the disease at the expense of increasing the risk of pneumonia.Objectives: To assess the relationship between the presence of chronic bronchial infection (CBI), reduced number of circulating eosinophils, ICS treatment, and the risk of pneumonia in patients with COPD.Methods: This was a post hoc long-term observational study of an historical cohort of 201 patients with COPD (Global Initiative for Chronic Obstructive Lung Disease II-IV) who were carefully characterized (including airway microbiology) and followed for a median of 84 months. Results were analyzed by multivariate Cox regression and network analysis.Measurements and Main Results: Mean age was 70.3 years, 90.5% of patients were male, mean FEV₁ was 49%, 71.6% of patients were treated with ICS, 57.2% of them had bronchiectasis, and 20.9% had <100 blood eosinophils/μl. Pathogenic microorganisms were isolated in 42.3% of patients, and 22.4% of patients fulfilled the definition of CBI. During follow-up, 38.8% of patients suffered one or more episodes of pneumonia, with CBI (hazard ratio [HR], 1.635) and <100 eosinophils/μl (HR, 1.975) being independently associated with the risk of pneumonia, particularly when both coexist (HR, 3.126). ICS treatment increased the risk of pneumonia in those patients with <100 eosinophils/μl and CBI (HR, 2.925).Conclusions: Less than 100 circulating eosinophils/μl combined with the presence of CBI increase the risk of pneumonia in patients with COPD treated with ICS.

Transcriptional analysis identifies potential biomarkers and molecular regulators in pneumonia and COPD exacerbation

Lower respiratory infections, such as community-acquired pneumonia (CAP), and chronic obstructive pulmonary disease (COPD) rank among the most frequent causes of death worldwide. Improved diagnostics and profound pathophysiological insights are urgent clinical needs. In our cohort, we analysed transcriptional networks of peripheral blood mononuclear cells (PBMCs) to identify central regulators and potential biomarkers. We investigated the mRNA- and miRNA-transcriptome of PBMCs of healthy subjects and patients suffering from CAP or AECOPD by microarray and Taqman Low Density Array. Genes that correlated with PBMC composition were eliminated, and remaining differentially expressed genes were grouped into modules. One selected module (120 genes) was particularly suitable to discriminate AECOPD and CAP and most notably contained a subset of five biologically relevant mRNAs that differentiated between CAP and AECOPD with an AUC of 86.1%. Likewise, we identified several microRNAs, e.g. miR-545-3p and miR-519c-3p, which separated AECOPD and CAP. We furthermore retrieved an integrated network of differentially regulated mRNAs and microRNAs and identified HNF4A, MCC and MUC1 as central network regulators or most important discriminatory markers. In summary, transcriptional analysis retrieved potential biomarkers and central molecular features of CAP and AECOPD.