Rationale: Acute exacerbations of chronic obstructive pulmonary disease (AE-COPDs) are associated with a significant disease burden. Blood immune phenotyping may improve our understanding of a COPD endotype at increased risk of exacerbations. Objective: To determine the relationship between the transcriptome of circulating leukocytes and COPD exacerbations. Methods: Blood RNA sequencing data (n = 3,618) from the COPDGene (Genetic Epidemiology of COPD) study were analyzed. Blood microarray data (n = 646) from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study were used for validation. We tested the association between blood gene expression and AE-COPDs. We imputed the abundance of leukocyte subtypes and tested their association with prospective AE-COPDs. Flow cytometry was performed on blood in SPIROMICS (Subpopulations and Intermediate Outcomes in COPD Study) (n = 127), and activation markers for T cells were tested for association with prospective AE-COPDs. Measurements and Main Results: Exacerbations were reported 4,030 and 2,368 times during follow-up in COPDGene (5.3 ± 1.7 yr) and ECLIPSE (3 yr), respectively. We identified 890, 675, and 3,217 genes associated with a history of AE-COPDs, persistent exacerbations (at least one exacerbation per year), and prospective exacerbation rate, respectively. In COPDGene, the number of prospective exacerbations in patients with COPD (Global Initiative for Chronic Obstructive Lung Disease stage ⩾2) was negatively associated with circulating CD8+ T cells, CD4+ T cells, and resting natural killer cells. The negative association with naive CD4+ T cells was replicated in ECLIPSE. In the flow-cytometry study, an increase in CTLA4 on CD4+ T cells was positively associated with AE-COPDs. Conclusions: Individuals with COPD with lower circulating lymphocyte counts, particularly decreased CD4+ T cells, are more susceptible to AE-COPDs, including persistent exacerbations.

Background

Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF.

Methods

We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets.

Results

We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts.

Conclusions

We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.