- Open Access
Dynamic gene expressions of peripheral blood mononuclear cells in patients with acute exacerbation of chronic obstructive pulmonary disease: a preliminary study
Critical Care volume 18, Article number: 508 (2014)
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.
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.
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.
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.
Chronic obstructive pulmonary disease (COPD) is an inflammation-based syndrome characterized by progressive deterioration of pulmonary function and increasing airway obstruction . COPD is a major and growing public health burden, ranking as the fourth leading cause of death in the world . In China, it is the fourth leading cause of mortality in urban areas and the third leading cause in rural areas . Patients with COPD often experience a sudden deterioration, termed acute exacerbations of chronic obstructive pulmonary disease (AECOPD), along with a progressive decline in lung function; AECOPD becomes more frequent and severe when the severity of disease increases ,. There is a great need for early and sensitive diagnosis and novel therapeutic targets for the disease, especially for patients with AECOPD in whom COPD is diagnosed in the late phase of disease, when they have significant or irreversible impairment .
The progress of COPD is accelerated by the occurrence of the exacerbation induced by multiple factors, including infection. AECOPD is a serious event that is related to decreased health status, increased medical and social costs and increased mortality . Inflammatory cells (for example, lymphocytes, monocytes or macrophages, and their products) could interact with each other or with structural cells in the airways and the lung parenchymal and pulmonary vasculature, leading to the worsening of COPD . Increased numbers of CD8+ lymphocytes were suggested as one of COPD’s characteristics, being present only in smokers who develop the disease . Increased pulmonary inflammatory mediators in patients with COPD could attract inflammatory cells from the circulation, amplify the inflammatory process and induce structural changes .
Peripheral blood mononuclear cells (PBMCs) act as a critical component in the immune system to fight infection and adapt to intruders and play an important role in the development of AECOPD. Gene expression profiles of PBMCs were found to be disease-specific and associated with severity . PBMC samples were suggested as easy to gather and important to the discovery of biomarkers for diagnosis and therapeutic management of COPD ,, although gene expression changes in lung tissues were noted to be associated with COPD -. The aim of the present study was to determine AECOPD-specific biomarkers of PBMCs using the concept of clinical bioinformatics and integrating genomics, bioinformatics, clinical informatics and systems biology -. We translated all clinical measures, including patient complaints, history, therapies, clinical symptoms and signs, physician’s examinations, biochemical analyses, imaging profiles, pathologies and other measurements, into digital format using a digital evaluation scoring system. PBMCs were isolated from healthy volunteers and patients with stable COPD or AECOPD, and we investigated the disease specificity that we inferred from clinical informatics analysis to search for COPD- or AECOPD-specific genes and dynamic biomarkers for AECOPD.
Material and methods
The present study was approved by the Ethical Evaluation Committee of Zhongshan Hospital and designed using a case–control approach. From among 220 candidates comprising blood donors (60 healthy controls), inpatients (80 patients with AECOPD) and outpatients (80 patients with stable COPD) in Zhongshan Hospital, patients with AECOPD, patients with stable COPD and healthy controls matched for age and sex were recruited into the study between October 2011 and March 2012. The inclusion criteria for patients with COPD were as follows: (1) forced expiratory volume in 1 second (FEV1) <80% of predicted value adjusted for age, weight and height, and (2) an improvement in FEV1 following bronchodilator inhalation <12% of baseline FEV1. Patients with asthma who had a persistent airflow obstruction were excluded. Stable COPD was defined according to American Thoracic Society/European Respiratory Society consensus criteria as no requirement for increased treatment above maintenance therapy, other than bronchodilators, for 30 days . AECOPD was the reason for hospital admission and was characterized as a worsening of the patient’s respiratory symptoms that was beyond normal day-to-day variations and led to a change in medication ,. Healthy controls enrolled were blood donors at Zhongshan Hospital. Subjects with respiratory diseases, or any family history of lung disease, were excluded. PBMCs were harvested once from healthy controls and patients with stable COPD, as well as from patients with AECOPD, on the admission day and 3 and 10 days after the admission. Informed consent was given by the subjects themselves before they underwent lung function tests, high-resolution computed tomography and blood collection. The time points used in the present study were selected on the basis of our previous study for collecting plasma samples from healthy controls and from patients with stable COPD or AECOPD. The details of the study design are explained in Figure 1.
Digital evaluation score system
The Digital Evaluation Score System (DESS) is a score index used to translate clinical descriptions and information into clinical informatics, as described previously . Using this instrument, we took into account patient symptoms and signs, biochemical analyses and clinical imaging for patients with stable COPD or AECOPD. Briefly, for the assessment of severity, each component was assigned a score of 0, 1, 2 or 4. The score of 4 as the maximum value indicates far above normal range or much severer condition, and 0 as the minimum value indicates within normal physiological range. After compiling patient data, we added the points for each variable. The DESS scores ranged from 0 to 256 points, with a higher score indicating a severer condition. Patients were scored on the day when their blood samples were collected.
Isolation of PBMC RNA
PBMCs were isolated by using BD Vacutainer CPT cell preparation tubes (Becton Dickinson, Franklin Lakes, NJ, USA) according to the manufacturer’s instructions. Approximately 4 ml of whole blood was collected from each subject. Following centrifugation, cells were lysed for RNA isolation. DNase-free total RNA preparation was performed using TRIzol reagent (Life Technologies, Carlsbad, CA, USA) and the RNeasy kit (QIAGEN, Valencia, CA, USA) according to the manufacturers’ recommendations. RNA concentrations were determined by using a NanoDrop ND-1000 spectrophotometer (NanoDrop, Wilmington, DE, USA). RNA quality was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), and samples with an RNA integrity number >6.0 were used.
The Human 12×135K Gene Expression Array (Roche NimbleGen Systems, Madison, WI, USA), with about 45,000+ human genes and transcripts represented with public domain annotations, was applied for this study. Sample labeling and array hybridization were performed according to the one-color microarray-based gene expression analysis protocol (Roche NimbleGen Systems). Double-stranded cDNA (ds-cDNA) was synthesized from 5 μg of total RNA using an Invitrogen SuperScript reverse transcriptase ds-cDNA synthesis kit (Life Technologies) in the presence of 100 pmol oligo(dT) primers. ds-cDNA was cleaned and labeled in accordance with the NimbleGen gene expression analysis protocol. Briefly, ds-cDNA was incubated with 4 μg of RNase A at 37°C for 10 minutes and cleaned using phenol:chloroform:isoamyl alcohol, followed by ice-cold absolute ethanol precipitation. The purified cDNA was quantified using the NanoDrop ND-1000 spectrophotometer. For Cy3 labeling of cDNA, the NimbleGen one-color DNA labeling kit was used according to the manufacturer’s guidelines as detailed in its gene expression analysis protocol. One microgram of ds-cDNA was incubated for 10 minutes at 98°C with 1 optical density of Cy3-9mer primer. Next, 100 pmol of deoxynucleoside triphosphates and 100 U of the Klenow fragment (New England Biolabs, Ipswich, MA, USA) were added, and the mix was incubated at 37°C for 2 hours. The reaction was stopped by adding 0.1 vol of 0.5 M ethylenediaminetetraacetic acid, and the labeled ds-cDNA was purified by isopropanol/ethanol precipitation. Microarrays were hybridized at 42°C for 16 to 20 hours with 4 μg of Cy3-labeled ds-cDNA in NimbleGen hybridization buffer/hybridization component A in a hybridization chamber. Following hybridization, washing was performed using the NimbleGen wash buffer kit. After being washed in an ozone-free environment, the slides were scanned using an Axon GenePix 4000B microarray scanner (Molecular Devices, Sunnyvale, CA, USA).
For clinical data, all values were expressed as mean ± SE. Analyses were performed using SPSS software (SPSS 18.0; SPSS, Chicago, IL, USA). For microarray analysis, slides were scanned at 5 μm/pixel resolution using the Axon GenePix 4000B microarray scanner piloted by GenePix Pro 6.0 software (Molecular Devices). Scanned images (in TIFF file format) were then imported into NimbleScan software (version 2.5) files for grid alignment and expression data analysis. Expression data were normalized through quantile normalization and the Robust Multi-array Average (RMA) algorithm included in the NimbleScan software. The probe-level (*_norm_RMA.pair) files and gene-level (*_RMA.calls) files were generated after normalization. All gene-level files were imported into GeneSpring GX software (version 11.5.1; Agilent Technologies) for further analysis. Differentially expressed genes between two samples were identified by fold change filtering. Hierarchical clustering was performed using the GeneSpring GX software. Gene Ontology (GO) database analysis and pathway analysis were performed using the standard enrichment computation method. The GO database covers three domains: biological process, cellular component and molecular function. Fisher’s exact test was used to find more overlaps between the descriptive list and the GO annotation list than would be expected by chance. The P-value denoted the significance of GO term enrichment in the descriptive genes. The gene expression data are publicly available in the Gene Expression Omnibus database [GEO:GSE60399] .
Clinical informatics analysis
Clinical phenotypes are described in Table 1, including age, sex, smoking status, lung function test results and emphysema scores of the subjects. Control subjects were nonsmokers, and patients with stable COPD or AECOPD were ex-smokers. Because of the severity of disease, lung function tests were not performed at the onset of AECOPD; however, the baseline FEV1/forced vital capacity (FVC%) and FEV1/predicted percentage of patients with AECOPD were similar to those of patients with stable COPD. In addition, there was no significant difference in the extent of emphysema between patients with stable COPD and those with AECOPD (P = 0.47). DESS scores of subjects from each group are shown in Additional file 1. DESS values of patients with stable COPD or AECOPD were significantly higher than those of control subjects (P < 0.01), as shown in Table 2. DESS scores represented the severity of COPD and declined as the patient’s condition improved. DESS values of patients with AECOPD on day 1 of hospital admission (AE-1) were significantly higher than those on day 3 (AE-3) and day 10 (AE-10) (P < 0.05 and P < 0.01, respectively) (Table 2).
Gene expression profiles
The quality of the genetic data obtained after filtering and the distribution of data sets were assessed and visualized by creating box plots, which showed that there were no significant differences in the distributions of log2 ratios among the groups (see Additional file 2: Figure S1). The variation or reproducibility of gene expression between arrays of different groups was visualized and assessed by creating scatterplots, which are shown in Figure 2. There was a significant variation in gene arrays between healthy controls and patients with stable COPD or AECOPD (Figures 2A to 2D) and between patients with stable COPD and AECOPD (Figures 2E to 2G). The variation in gene array data at AE-1 and AE-3 was significantly different from that at AE-10 (Figures 2I and 2J), whereas there was no difference between AE-1 and AE-3 (Figure 2H). The results of hierarchical clustering showed gene expression profiles similar to those revealed by the scatterplots shown in Figure S2 of Additional file 2.
To identify differentially expressed genes, a fold change filtering between each group pair was performed with a threshold fold change ≥2.0. There were ten comparison pairs with information for fold changes and regulation (that is, SEQ-ID, log fold change, log or absolute fold change, or regulation), normalized intensities or annotations (that is, GENE_NAME, synonyms, description, NCBI_GENE_ID, chromosome, GO, UniGene ID, The Institute of Genomic Research Database-TDB (TIGRID) or Ensembl ID), as shown in Additional file 3. Table 3 shows the number of genes overexpressed more than twofold, (for example, 4,508, 3,899, 4,167 and 3,488 genes of stable, AE-1, AE-3 and AE-10, respectively, above controls; 4,067, 5,063 or 5,451 genes of AE-1, AE-3 and AE-10, respectively, above stable COPD; 586 genes of AE-3 above AE-1; and 1,735 and 1,706 genes of AE-10, respectively, above AE-1 and AE-3). Tables 4, 5 and 6, respectively, list the genes overexpressed (above controls) in PBMCs from patients with stable COPD, AE-1, AE-3 or AE-10 by more than 30-fold (Table 4), between 20- and 30-fold (Table 5) and between 15- and 20-fold (Table 6). Tables 7, 8 and 9 list the genes overexpressed (above patients with stable COPD) in PBMCs from patients with AE-1, AE-3 or AE-10 by more than 30-fold (Table 7), between 20- and 30-fold (Table 8) and between 15- and 20-fold. Table 10 presents upregulated genes in PBMCs of patients at AE-1, AE-3 or AE-10.
Table 11 lists the number of genes downregulated more than twofold, including 4,516, 2,975, 3,426 and 2,798 genes of PBMCs from patients with stable COPD on AE-1, AE-3 and AE-10, respectively, below controls; 3,207, 4,510 and 5288 genes on AE-1, AE-3 and AE-10, respectively, below stable COPD; 598 genes from AE-3 below AE-1; and 2,162 and 1,918 genes from AE-10 below those from AE-1 and AE-3, respectively. Downregulated genes of PBMCs from patients with stable COPD, AE-1, AE-3 or AE-10 greater than tenfold, between 10- and 8-fold or between 8- and 6-fold below healthy control subjects are listed in Tables 12, 13 and 14, respectively. Downregulated genes of PBMCs from patients at AE-1, AE-3 or AE-10 compared to stable COPD, or among patients with AECOPD, are shown in Tables 15 and 16.
To search for COPD-specific genes, co–differentially expressed genes of PBMCs from patients with stable COPD or AECOPD were compared with those from control subjects (listed in Additional file 4). There were five groups and four comparison pairs with information regarding fold changes and regulation (that is, SEQ-ID, fold change, log or absolute fold change, or regulation), normalized intensities or annotations (that is, GENE_NAME, synonyms, description, NCBI_GENE_ID, chromosome, GO, UniGene ID, TIGRID or Ensembl ID). Seventy-nine genes were upregulated and 23 genes downregulated in PBMCs from patients with COPD, including both stable COPD and AECOPD, as compared to the healthy control subjects, as shown in Table 17. Of them, 14 genes were upregulated and 2 were downregulated more than tenfold, as compared to control subjects, including carcinoembryonic antigen–related cell adhesion molecule 1, collagen type VIα3(VI), collagen type I(α)2(I), nucleolar protein 3 (apoptosis repressor with CARD domain), melanophilin, cell surface–associated mucin 1, nuclear protein 1, chemokine (C-X-C motif) ligand 17, claudin 4, ribonuclease 1, imprinted maternally expressed transcript, defensin α1, transcription factor CP2-like 1 and sterol carrier protein 2 (SCP2).
To search for AECOPD-specific genes, co–differentially expressed genes of PBMCs from patients with AECOPD on days 1, 3 and 10 were compared to those from either patients with stable COPD or healthy control subjects (listed in Additional file 4). There were five groups and six comparison pairs with information regarding fold changes and regulation (that is, SEQ-ID, fold change, log or absolute fold change, or regulation), normalized intensities or annotations (that is, GENE_NAME, synonyms, description, NCBI_GENE_ID, chromosome, GO, UniGene ID, TIGRID or Ensembl ID). As compared with both patients with stable COPD and healthy control subjects, 58 genes were upregulated more than fivefold and 238 downregulated more than twofold in patients with AECOPD. Of them, eight upregulated (more than tenfold) and eight downregulated (more than threefold) genes are listed in Table 18. These genes include FBJ murine osteosarcoma viral oncogene homologue (FOS); interferon α-inducible protein 27 (IFI27); cysteine-rich angiogenic inducer 61 (CYR61), connective tissue growth factor (CTGF); G protein–coupled receptor family C group 5 member A (GPRC5A); FBJ murine osteosarcoma viral oncogene homologue B (FOSB); decorin (DCN); hypothetical LOC387763 (LOC387763); killer cell immunoglobulin-like receptor, two domains, short cytoplasmic tail, 2 (KIR2DS2); SH2 domain containing 1B (SH2D1B); CD8b molecule (CD8B); olfactory receptor family 2, subfamily W, member 5 (OR2W5); fibroblast growth factor binding protein 2 (FGF2); and transcription factor 7 (TCF7).
Dynamic change in gene expression in patients with AECOPD
Dynamic changes (down–down, down–up, up–down and up–up) of co–differentially expressed genes of PBMCs from patients with AECOPD are listed in Additional file 4, including fold changes and regulation (that is, SEQ-ID, fold change, log or absolute fold change, or regulation), normalized intensities or annotations (that is, GENE_NAME, synonyms, description, NCBI_GENE_ID, chromosome, GO, UniGene ID, TIGRID or Ensembl ID). Table 19 shows the dynamic changes in the patterns of down–down (52 genes), down–up (131 genes), up–down (238 genes) and up–up (8 genes) more than twofold, as compared with the gene expression on the previous day. The major genes of PBMCs from patients with AECOPD were aminolevulinate, delta-, synthase 2 (ALAS2), erythrocyte membrane protein band 4.2 (EPB42) and carbonic anhydrase I (CA1) in a down–down pattern; selenium-binding protein 1 (SELENBP1) and myosin heavy chain 9, non-muscle (MYH9), in a down–up pattern; HLA complex group 27 (HCG27), BCL2-related protein A1 (BCL2A1), G protein–coupled receptors 109A and 109B (GPR109A and GPR109B) in an up–down pattern; and zeta protein kinase C (PRKCZ), ATP-binding cassette, subfamily A, member 8 (ABCA8), and folate receptor 1 (adult) (FOLR1) in an up–up pattern (Table 19). Levels of genes from patients with AECOPD were also compared with those from patients with stable COPD, as shown in Figure 3, where positive or negative values indicate up- or downregulation as compared with those from patients with stable COPD. When correlated with DESS, ALAS2 and CA1 had similar patterns of change with DESS.
Gene ontology analysis and pathway analysis
Within ten comparison pairs, up- or downregulated genes mainly involved in the biological process are shown in Figures S3 and S4 of Additional file 2, those in cellular components are shown in Figures S5 and S6 of Additional file 2 and those in molecular functions are shown in Figures S7 and S8 of Additional file 2. Additional file 5 lists gene numbers for ten comparison pairs with certain GO terms and different ranges of enrichment scores.
In the biological process, COPD-specific upregulated genes were involved mainly in peptide cross-linking, blood vessel development, biological adhesion or cell adhesion (Figure 4A). COPD-specific downregulated genes were involved mainly in T cell receptor signaling pathways, antigen receptor–mediated signaling pathways, immune response–activating cell surface receptor signaling pathways or steroid biosynthetic process (Figure 4B). AECOPD-specific genes upregulated in response to organic substance, response to wounding, multicellular organismal process or response to chemical stimulus are shown in Figure 4C. AECOPD-specific downregulated genes were involved mainly in the regulation of immune response and the immune system process or in the immune response and immune system process themselves (Figure 4D). In the cellular component, COPD-specific upregulated genes were involved mainly in the extracellular region, the extracellular matrix part, the proteinaceous extracellular matrix or the extracellular matrix (Figure 5A). COPD-specific downregulated genes were involved mainly in the major histocompatibility complex class II (MHC II) protein complex, microbody lumen, peroxisomal matrix or MHC II protein complex (Figure 5B). AECOPD-specific upregulated genes were involved mainly in the extracellular region part, the extracellular matrix, the extracellular space or the extracellular region (Figure 5C). AECOPD-specific downregulated genes were involved mainly in the cell periphery and the plasma membrane and were integral to the plasma membrane and intrinsic to the plasma membrane (Figure 5D). In molecular function, COPD-specific upregulated genes participated mainly in extracellular matrix structural constituent, platelet-derived growth factor binding, serine-type endopeptidase activity and protein binding (Figure 6A). COPD-specific downregulated genes were involved mainly in nucleoside kinase activity, MHC class II receptor activity, C-acyltransferase activity and ephrin receptor activity (Figure 6B). AECOPD-specific upregulated genes were involved mainly in protein binding, growth factor binding, calcium ion binding and polysaccharide binding (Figure 6C). AECOPD-specific downregulated genes were involved mainly in receptor activity, signaling receptor activity, molecular transducer activity and signal transducer activity (Figure 6D).
COPD-specific upregulated genes also participated in extracellular matrix receptor interaction, protein digestion and absorption, focal adhesion and the phosphatidylinositol 3-kinase-Akt signaling pathway (Figure 7A). AECOPD-specific upregulated genes participated in Chagas disease, complement and coagulation cascades, pertussis and Staphylococcus aureus infection (Figure 7B). AECOPD-specific downregulated genes participated in antigen processing and presentation, natural killer cell–mediated cytotoxicity, graft-versus-host disease and thyroid cancer (Figure 7C).
PBMCs play a critical and important role in the occurrence of AECOPD, owing to less capacity for balancing the proinflammatory immune response caused by infection and for secreting adequate amounts of anti-inflammatory cytokines . The fact that patients with COPD are more susceptible to acute exacerbation has been suggested to be associated with PBMC dysfunction and failure of adaptation to infection, stimuli or hypoxia, although there have been not yet studies on the phenotypes of PBMCs in AECOPD. For example, PBMCs from patients with COPD could not induce hypoxia-inducible factor 1 and vascular endothelial growth factor, owing to a reduction in histone deacetylase 7 under hypoxic condition . It was suggested that overproduction of proinflammatory cytokines (CXCL6 and interleukin 6 (IL-6)) from human PBMCs could be stimulated by the infection through activation of Toll-like receptor 4, nicotinamide adenine dinucleotide phosphate oxidase phosphatidylinositol 3-kinase and nuclear factor κB , at least as partial mechanisms by which PBMCs may be involved in the occurrence of AECOPD. The present study provides initial evidence that dynamic alterations of PBMC genetic phenotypes occurred in patients with AECOPD after their hospital admission and during their hospital stay.
Gene expression profiles of PBMCs were investigated in patients with COPD, compared with healthy controls and correlated with lung function measurement . Differential expression of 45 known genes was identified, of which 16 markers had significant correlation with quantitative traits and differential expression between cases and controls and 2 genes, RP9 and NAPE-PLD, were identified as decreased in patients with COPD, as compared to controls, in both lung tissue and blood. Gene expression profiles of PBMCs were recently identified and validated in smokers with and without COPD and corrected with clinical phenotypes such as sex, age, body mass index, family history, smoking status and pack-years of smoking . Of them, 16 candidate genes were found to be associated with airflow obstruction and secondary clinical phenotypes, 12 with emphysema, 13 with gas trapping and 8 with distance walked. Both previous studies demonstrated the gene expression profiles of PBMCs from patients with stable COPD and addressed the potential significance of smoking. In the present study, we selected healthy control subjects and patients who were not current smokers and demonstrated gene expression profiles of PBMCs from patients with COPD, including stable COPD and AECOPD. We addressed COPD-specific gene expression profiles that should appear in both stable COPD and COPD exacerbation conditions and found COPD-specific 79 genes were upregulated and 23 genes down-regulated more than fivefold as compared with gene expression in controls. In the present study, we selected consistent up- or downregulated gene expression on days 1, 3 and 10 of AECOPD-specific as compared with gene expression in both healthy controls and patients with stable COPD, as AECOPD-specific gene expression profiles. We found that 58 AECOPD-specific genes were upregulated more than fivefold and 238 genes were downregulated more than twofold, as compared to both control subjects and patients with stable COPD.
Variation of gene expression profiles is dependent upon multiple uncontrollable factors, such as study population, age, history, genetic background and treatment. In addition, gene expression profiles vary between harvested sample types, such as sputum, bronchoalveolar lavage fluid, blood or lung tissues. For example, 102 genes were identified to distinguish between non- or mild emphysema and severe emphysema in lung tissue  and to distinguish 70 microRNAs and 2,667 mRNAs between smoking patients with or without COPD . In the present study, we investigated gene expression profiles of PBMCs from control subjects, patients with stable COPD, and patients with AECOPD on day 1, day 3 and day 10 of hospital admission, and we found about 3,000 overexpressed genes and 2,000 downregulated genes in patients with stable COPD or AECOPD, as compared with control subjects. These findings indicate that those COPD-specific genes exist in the stable COPD condition and during acute exacerbations of COPD.
Of the COPD-specific genes we studied, CEACAM1, COL6A3, NOL3, COL1A2, MLPH, MUC1, P8, UNQ473, CLDN4, RNASE1, H19, DEFA1 and LOC653600 were upregulated more than tenfold, mainly related to nuclear proteins, collagens or molecular structure. We noted that transcription factor CP2 (TFCP2L1) and SCP2 were downregulated more than tenfold. In previous studies, these genes, including CEACAM1, TFCP2L1 and SCP2, were not found to be associated with COPD. The SCP2 gene is located within chromosome 1 and encodes the nonspecific lipid transfer protein SCP2, which is involved in organellar fatty acid metabolism , and the translocation of cytoplasmic free cholesterol to the mitochondria . Our results indicate that PBMCs from patients with stable COPD or AECOPD had downregulated SCP2, which might point to severe metabolic disorder and thus that SCP2 downregulation might contribute to one of the common comorbidities of COPD . TFCP2 is a member of a family of transcription factors that regulate genes involved in events from early development to terminal differentiation . PBMCs with downregulated TFCP2 of patients with COPD might have less capacity of the transcriptional switch of globin gene promoters, many other cellular and viral gene promoters, or interaction with certain inflammatory response factors, although the exact mechanism and pathological role remain unclear.
AECOPD-specific gene expression profiles were selected by comparing them with both healthy control subjects and patients with stable COPD, including 647 upregulated genes and 238 downregulated genes (greater than twofold upregulation). Of them, FOS, IFI27, CYR61, CTGF, GPRC5A, FOSB, DCN and LOC387763 were upregulated more than tenfold and KIR2DS2, SH2D1B, CD8B, OR2W5, KSP37 and TCF7 were downregulated more than threefold.
We noticed that some genes, such as FOS, CYR61 and CTGF, were upregulated in PBMCs from patients with either stable COPD or AECOPD, consistent with the lung tissue gene expression profiles of patients with COPD or smokers, in whom the genes were expressed mainly in alveolar epithelial cells, airway epithelial cells and stromal and inflammatory cells . Other genes, including GPRC5A, LOC387763 and KIR2DS2, were not found to be associated with AECOPD in previous publications. CTGF is a cysteine-rich peptide implicated in several biological processes, such as cell proliferation, survival and migration, and involved in pulmonary vascular remodeling and hypertension in COPD. It was evidenced by the experimental finding that CTGF short-hairpin RNA could significantly prevent CTGF and cyclin D1 expression, arrest cell cycle at the G0/G1 phase, suppress cell proliferation in smoking-exposed pulmonary smooth muscle cells and ameliorate pulmonary vascular remodeling . Another study demonstrated that some inflammatory genes (IL-1β, IL-6, IL-8, CCL2 and CCL8) were upregulated, whereas some growth factor receptor genes (BMPR2, CTGF, FGF1, KDR and TEK) were downregulated in lung tissue samples from patients who were current smokers or had moderate COPD .
Downregulation of TCF7 was found in PBMCs of patients with COPD and current smoking and was correlated with some clinical phenotypes, such as emphysema, gas trapping and distance walked . In the present study, we also found that TCF7 was downregulated in ex-smokers with COPD by about an absolute threefold compared with control subjects, and, in patients with AECOPD, TCF7 was downregulated by about an absolute tenfold compared with both control subjects and patients with stable COPD. These findings indicate that TCF7 not only is a COPD-specific gene but also is associated with the severity of the disease. TCF7 is a member of a family of HMG box containing factors associated with β-catenin to mediate Wnt signaling, controls the switch between cell self-renewal and differentiation and plays a role in B cell and T cell development. TCF7 was found to be the most downregulated transcription factor when CD34+ cells switched into CD34− cells through a coordinated regulation of the binding between TCF7 and the short isoforms of RUNX1 . It is possible the downregulation of TCF7 and associated regulation may be one part of molecular mechanism of PBMC incapacity during AECOPD.
Dynamic alterations of gene expression profiles in patients with AECOPD were evaluated with dynamic DESS scores. ALAS2, EPB42 and CA1 were co–differentially expressed with a down–down type in patients with AECOPD. Among these three genes, the CA1 gene encodes a protein which is important in respiratory function, fluid secretion and maintenance of cellular acid–base homeostasis . The genes with a down–up type included SELENBP1, MYH9 and an unnamed gene in chromosome 19, both of which are associated with psychotic disorders ,. One limitation of the present study is the small sample size, which detracts from the generalizability of the results presented.
Dynamic alterations of PBMC gene expression profiles were initially investigated in patients with AECOPD, as compared with healthy control subjects or patients with stable COPD. A panel of genes, including eight that were upregulated and eight that were downregulated, were recommended as AECOPD-specific dynamic biomarkers. AECOPD-specific up- or downregulated genes in the biological process, cellular components or molecular function were defined and participated in complement and coagulation cascades, infection, antigen processing and presentation, natural killer cell–mediated cytotoxicity, and/or cancer-causing potential. The integration of dynamic bioinformatics with clinical phenotypes helped us to identify and validate AECOPD-specific biomarkers to help define the severity, duration and response of the disease to therapies.
Circulating dynamic biomarkers were identified for the specificity and severity of AECOPD.
A panel of 16 genes were selected as AECOPD-specific biomarkers.
This is an initial study designed to examine gene expression profiles of peripheral blood mononuclear cells and identify dynamic changes of AECOPD-specific biomarkers.
Acute exacerbations of chronic obstructive pulmonary disease on day 1
Acute exacerbations of chronic obstructive pulmonary disease on day 3
Acute exacerbations of chronic obstructive pulmonary disease on day 10
Acute exacerbation of chronic obstructive pulmonary disease
Aminolevulinate, delta-, synthase 2
Carbonic anhydrase I
Chronic obstructive pulmonary disease
Chemokine (C-X-C motif) ligand 8
Digital evaluation score system
Erythrocyte membrane protein band 4.2
Forced expiratory volume in 1 second
Forced vital capacity
Major histocompatibility complex
Myosin, heavy polypeptide 9, non-muscle
Peripheral blood mononuclear cell
Sterol carrier protein 2
Selenium-binding protein 1
Transcription factor 7
Transcription factor CP2-like 1
Celli BR, MacNee W, Agusti A, Anzueto A, Berg B, Buist AS, Calverley PMA, Chavannes N, Dillard T, Fahy B, Fein A, Heffner J, Lareau S, Meek P, Martinez F, McNicholas W, Muris J, Austegard E, Pauwels R, Rennard S, Rossi A, Siafakas N, Tiep B, Vestbo J, Wouters E, ZuWallack R: Standards for the diagnosis and treatment of patients with COPD: a summary of the ATS/ERS position paper. Eur Respir J. 2004, 23: 932-946. 10.1183/09031936.04.00014304.
Pauwels RA, Rabe KF: Burden and clinical features of chronic obstructive pulmonary disease (COPD). Lancet. 2004, 364: 613-620. 10.1016/S0140-6736(04)16855-4.
Fang X, Wang X, Bai C: COPD in China: the burden and importance of proper management. Chest. 2011, 139: 920-929. 10.1378/chest.10-1393.
Celli BR, Barnes PJ: Exacerbations of chronic obstructive pulmonary disease. Eur Respir J. 2007, 29: 1224-1238.
Donaldson GC, Seemungal TA, Patel IS, Lloyd-Owen SJ, Wilkinson TM, Wedzicha JA: Longitudinal changes in the nature, severity and frequency of COPD exacerbations. Eur Respir J. 2003, 22: 931-936. 10.1183/09031936.03.00038303.
Murtagh E, Heaney L, Gingles J, Shepherd R, Kee F, Patterson C, MacMahon J: Prevalence of obstructive lung disease in a general population sample: the NICECOPD study. Eur J Epidemiol. 2005, 20: 443-453. 10.1007/s10654-005-1248-8.
Seemungal TA, Donaldson GC, Bhowmik A, Jeffries DJ, Wedzicha JA: Time course and recovery of exacerbations in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2000, 161: 1608-1613. 10.1164/ajrccm.161.5.9908022.
Hogg JC, Chu F, Utokaparch S, Woods R, Elliott WM, Buzatu L, Cherniack RM, Rogers RM, Sciurba FC, Coxson HO, Paré PD: The nature of small-airway obstruction in chronic obstructive pulmonary disease. N Engl J Med. 2004, 350: 2645-2653. 10.1056/NEJMoa032158.
Barnes PJ, Shapiro SD, Pauwels RA: Chronic obstructive pulmonary disease: molecular and cellular mechanisms. Eur Respir J. 2003, 22: 672-688. 10.1183/09031936.03.00040703.
Melguizo C, Prados J, Luque R, Ortiz R, Rama AR, Caba O, Rodríguez-Serrano F, Álvarez PJ, Aránega A: Modulation of multidrug resistance gene expression in peripheral blood mononuclear cells of lung cancer patients and evaluation of their clinical significance. Cancer Chemother Pharmacol. 2013, 71: 537-541. 10.1007/s00280-012-2036-2.
Pinto-Plata V, Toso J, Lee K, Park D, Bilello J, Mullerova H, De Souza MM, Vessey R, Celli B: Profiling serum biomarkers in patients with COPD: associations with clinical parameters. Thorax. 2007, 62: 595-601. 10.1136/thx.2006.064428.
Bhattacharya S, Tyagi S, Srisuma S, Demeo DL, Shapiro SD, Bueno R, Silverman EK, Reilly JJ, Mariani TJ: Peripheral blood gene expression profiles in COPD subjects. J Clin Bioinforma. 2011, 1: 12-10.1186/2043-9113-1-12.
Golpon HA, Coldren CD, Zamora MR, Cosgrove GP, Moore MD, Tuder RM, Geraci MW, Voelkel NF: Emphysema lung tissue gene expression profiling. Am J Respir Cell Mol Biol. 2004, 31: 595-600. 10.1165/rcmb.2004-0008OC.
Ning W, Li CJ, Kaminski N, Feghali-Bostwick CA, Alber SM, Di YP, Otterbein SL, Song R, Hayashi S, Zhou Z, Pinsky DJ, Watkins SC, Pilewski JM, Sciurba FC, Peters DG, Hogg JC, Choi AM: Comprehensive gene expression profiles reveal pathways related to the pathogenesis of chronic obstructive pulmonary disease. Proc Natl Acad Sci U S A. 2004, 101: 14895-14900. 10.1073/pnas.0401168101.
Spira A, Beane J, Pinto-Plata V, Kadar A, Liu G, Shah V, Celli B, Brody JS: Gene expression profiling of human lung tissue from smokers with severe emphysema. Am J Respir Cell Mol Biol. 2004, 31: 601-610. 10.1165/rcmb.2004-0273OC.
Wang X, Liotta L: Clinical bioinformatics: a new emerging science. J Clin Bioinforma. 2011, 1: 1-10.1186/2043-9113-1-1.
Wang X: Role of clinical bioinformatics in the development of network-based Biomarkers. J Clin Bioinforma. 2011, 1: 28-10.1186/2043-9113-1-28.
Chen H, Wang X: Significance of bioinformatics in research of chronic obstructive pulmonary disease. J Clin Bioinforma. 2011, 1: 35-10.1186/2043-9113-1-35.
Rodriguez-Roisin R: Toward a consensus definition for COPD exacerbations. Chest. 2000, 117 (5 Suppl 2): 398S-401S. 10.1378/chest.117.5_suppl_2.398S.
Chen H, Song Z, Qian M, Bai C, Wang X: Selection of disease-specific biomarkers by integrating inflammatory mediators with clinical informatics in AECOPD patients: a preliminary study. J Cell Mol Med. 2012, 16: 1286-1297. 10.1111/j.1582-4934.2011.01416.x.
Wang X, Wu X: Expression analysis of stable chronic obstructive pulmonary disease and acute exacerbation of chronic obstructive pulmonary disease. [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE60399] (published 14 August 2014).
Rupp J, Kothe H, Mueller A, Maass M, Dalhoff K: Imbalanced secretion of IL-1β and IL-1RA in Chlamydia pneumoniae-infected mononuclear cells from COPD patients. Eur Respir J. 2003, 22: 274-279. 10.1183/09031936.03.00007303.
To M, Yamamura S, Akashi K, Charron CE, Barnes PJ, Ito K: Defect of adaptation to hypoxia in patients with COPD due to reduction of histone deacetylase 7. Chest. 2012, 141: 1233-1242. 10.1378/chest.11-1536.
Ngkelo A, Meja K, Yeadon M, Adcock I, Kirkham PA: LPS induced inflammatory responses in human peripheral blood mononuclear cells is mediated through NOX4 and G iα dependent PI-3kinase signalling. J Inflamm (Lond) 2012, 9:1.,
Bahr TM, Hughes GJ, Armstrong M, Reisdorph R, Coldren CD, Edwards MG, Schnell C, Kedl R, LaFlamme DJ, Reisdorph N, Kechris KJ, Bowler RP: Peripheral blood mononuclear cell gene expression in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2013, 49: 316-323. 10.1165/rcmb.2012-0230OC.
Ezzie ME, Crawford M, Cho JH, Orellana R, Zhang S, Gelinas R, Batte K, Yu L, Nuovo G, Galas D, Diaz P, Wang K, Nana-Sinkam SP: Gene expression networks in COPD: microRNA and mRNA regulation. Thorax. 2012, 67: 122-131. 10.1136/thoraxjnl-2011-200089.
Szulakowski P, Crowther AJ, Jiménez LA, Donaldson K, Mayer R, Leonard TB, MacNee W, Drost EM: The effect of smoking on the transcriptional regulation of lung inflammation in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2006, 174: 41-50. 10.1164/rccm.200505-725OC.
Morris HR, Larsen BS, Billheimer JT: A mass spectrometric study of the structure of sterol carrier protein SCP 2from rat liver. Biochem Biophys Res Commun 1988, 154:476–482.,
Yanase T, Hara T, Sakai Y, Takayanagi R, Nawata H: Expression of sterol carrier protein 2 (SCP2) in human adrenocortical tissue. Eur J Endocrinol. 1996, 134: 501-507. 10.1530/eje.0.1340501.
Mannino DM, Thorn D, Swensen A, Holguin F: Prevalence and outcomes of diabetes, hypertension and cardiovascular disease in COPD. Eur Respir J. 2008, 32: 962-969. 10.1183/09031936.00012408.
Swendeman SL, Spielholz C, Jenkins NA, Gilbert DJ, Copeland NG, Sheffery M: Characterization of the genomic structure, chromosomal location, promoter, and development expression of the α-globin transcription factor CP2. J Biol Chem. 1994, 269: 11663-11671.
Wang R, Xu YJ, Liu XS, Zeng DX, Xiang M: Knockdown of connective tissue growth factor by plasmid-based short hairpin RNA prevented pulmonary vascular remodeling in cigarette smoke-exposed rats. Arch Biochem Biophys. 2011, 508: 93-100. 10.1016/j.abb.2011.01.019.
Llinàs L, Peinado VI, Goñi JR, Rabinovich R, Pizarro S, Rodriguez-Roisin R, Barberà JA, Bastos R: Similar gene expression profiles in smokers and patients with moderate COPD. Pulm Pharmacol Ther. 2011, 24: 32-41. 10.1016/j.pupt.2010.10.010.
Wu JQ, Seay M, Schulz VP, Hariharan M, Tuck D, Lian J, Du J, Shi M, Ye Z, Gerstein M, Snyder MP, Weissman S: Tcf7 is an important regulator of the switch of self-renewal and differentiation in a multipotential hematopoietic cell line. PLoS Genet. 2012, 8: e1002565-10.1371/journal.pgen.1002565.
Davis MB, West LF, Barlow JH, Butterworth PH, Lloyd JC, Edwards YH: Regional localization of carbonic anhydrase genes CA1 and CA3 on human chromosome 8. Somat Cell Mol Genet. 1987, 13: 173-178. 10.1007/BF01534697.
Amar S, Ovadia O, Maier W, Ebstein R, Belmaker RH, Mishmar D, Agam G: Copy number variation of the SELENBP1 gene in schizophrenia. Behav Brain Funct. 2010, 6: 40-10.1186/1744-9081-6-40.
Kanazawa T, Chana G, Glatt SJ, Mizuno H, Masliah E, Yoneda H, Tsuang MT, Everall IP: The utility of SELENBP1 gene expression as a biomarker for major psychotic disorders: replication in schizophrenia and extension to bipolar disorder with psychosis. Am J Med Genet B Neuropsychiatr Genet. 2008, 147B: 686-689. 10.1002/ajmg.b.30664.
The work was supported by Shanghai Leading Academic Discipline Project (project B115), Zhongshan Distinguished Professor Grant (to XDW), the National Nature Science Foundation of China (91230204, 81270099, 81320108001, 81270131, 81300010), the Shanghai Committee of Science and Technology (12JC1402200, 12431900207, 11410708600), the Zhejiang Provincial Natural Science Foundation (Z2080988), the Zhejiang Provincial Science Technology Department Foundation (2010C14011) and the Ministry of Education, Academic Special Science and Research Foundation for PhD Education (20130071110043).
The authors declare that they have no competing interests.
XW carried out the study, participated in the data analysis and drafted the manuscript. XRS participated in the data mining and analysis. CSC and CXB participated in the study design and data analysis and helped to revise the manuscript. XDW conceived of the study, participated in its design and coordination and finalized the manuscript. All authors read and approved the final manuscript.
Xiaodan Wu, Xiaoru Sun contributed equally to this work.
Electronic supplementary material
Additional file 1: DESS scores. This file lists Digital Evaluation Score System (DESS) scores of subjects from each group. (DOC 44 KB)
Additional file 2: A box plot showing distributions of log2 ratios among groups. They reflect our assessment of the quality of genetic data after the filtering and distribution of data sets. Figure S2. Hierarchical clustering shows distinguishable gene expression profiles and relationships between different groups. Figure S3. Co–differentially upregulated genes within 10 comparison pairs mainly involved in the biological process. Stable vs Con (A); AE-1 vs Con (B); AE-3 vs Con (C); AE-10 vs Con (D); AE-1 vs Stable (E); AE-3 vs Stable (F); AE-10 vs Stable (G); AE-3 vs AE-1 (H); AE-10 vs AE-1 (I); AE-10 vs AE-3 (J). Figure S4. Co–differentially downregulated genes within 10 comparison pairs mainly involved in the biological process. Figure S5. Co–differentially upregulated genes within 10 comparison pairs mainly involved in the cellular component. Figure S6. Co–differentially downregulated genes within 10 comparison pairs mainly involved in the cellular component. Figure S7. Co–differentially upregulated genes within 10 comparison pairs mainly involved in the molecular function. Figure S8. Co–differentially downregulated genes within 10 comparison pairs mainly involved in the molecular function. (PDF 603 KB)
Additional file 3: Differentially expressed genes. This file lists 10 comparison pairs with information of fold changes and regulation, normalized intensities or annotations. (PDF 104 KB)
Additional file 4: Co–differentially expressed genes. This file lists COPD-specific and AECOPD-specific genes, as well as dynamically changed genes, in patients with AECOPD. (PDF 264 KB)
Additional file 5: Gene Ontology database. This file lists gene numbers for 10 comparison pairs with certain GO (Gene Ontology) terms and different enrichment score ranges. (XLS 422 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Wu, X., Sun, X., Chen, C. et al. Dynamic gene expressions of peripheral blood mononuclear cells in patients with acute exacerbation of chronic obstructive pulmonary disease: a preliminary study. Crit Care 18, 508 (2014). https://doi.org/10.1186/s13054-014-0508-y
- Chronic Obstructive Pulmonary Disease
- Gene Ontology
- Connective Tissue Growth Factor
- Downregulated Gene
- Stable Chronic Obstructive Pulmonary Disease