Patients with Chronic Obstructive Pulmonary Disease (COPD) suffer from periodical exacerbations, characterized by recurrent episodes of worsening respiratory symptoms. Exacerbations result in further decreases in lung function, impairing the patient’s quality of life, and increasing the use of healthcare resources. In addition, development of Community Acquired Pneumonia (CAP) is also common in COPD exacerbations. Gene expression profiling (GEP) in cells obtained from peripheral blood has proved to be a very useful non-invasive approach for biomarker discovery and for studying molecular pathogenesis of prevalent diseases. While there is previous work identifying gene expression markers associated with COPD using Peripheral Blood Mononuclear Cells (PBMCs), the transcriptomic picture associated with critical respiratory illness in this disease is not known at the present moment. The objective was to describe the gene expression signatures associated with critical COPD respiratory illness compared to non-critical COPD exacerbations, as a preliminary study aimed to identify the genes associated with the severity of this disease.
Patients and samples
12 patients with pre-existing diagnosis of COPD in need of admission to the Intensive Care Unit due to acute respiratory failures were compared to a group of 16 COPD patients with non critical disease exacerbation admitted to the Respiratory Medicine Service. 4 healthy voluntary donors of similar ages to the patients were recruited for gene expression data normalization. A sample of 2.5 ml of blood was collected in the first 24 hours following admission to the ICU or to the Respiratory Medicine Service. Acute exacerbation of COPD was defined as a patient showing 2 symptoms (at least 1 major) for 2 consecutive days. Major symptoms were defined as increased dyspnoea, sputum volume, or sputum purulence and minor symptoms were increased cough, wheeze, sore throat, or coryzal symptoms. Definition of CAP was based on current American Thoracic Society and Infectious Disease Society of America guidelines . Informed consent was obtained directly from each patient before enrolment. Patient’s identification remained anonymous. The protocol was approved by the Ethics Committee on Clinical Research of each one of the participating centers.
1. Investigación Biomédica del Clínico (ibC), Hospital Clínico Universitario de Valladolid, Avda Ramón y Cajal 3, 47005 Valladolid, Spain
2. Infection & Immunity Medical Investigation Unit (IMI), Microbiology and Immunology Service, Hospital Clínico Universitario, Valladolid, Spain - IECSCYL, 3. Critical Care Department, Hospital Son Llatzer, Palma de Mallorca, Spain - SEMICYUC 4. Respiratory Medicine Service, Hospital Clínico Universitario-SACYL / SEPAR, Valladolid, Spain 5. Critical Care Department. Parc Tauli, Hospital Universitari de Sabadell, Spain. 6. Critical Care Department, Hospital Clínico Universitario-SACYL/SEMICYUC, Valladolid, Spain. 7. Critical Care Department, Hospital Universitari Joan XXIII. Tarragona, Spain. 8. Internal Medicine Service. Hospital Clínico Universitario - SACYL, Valladolid, Spain. 9. Microbiology Service, Hospital Son Llatzer, Palma de Mallorca, Spain.
Sputum samples were routinely Gram stained and cultured on general purpose media (blood agar, chocolate agar, and the differentials mediums McConkey agar and Chapman agar). Fungal infections were screened by culturing sputum samples on Sabouraud agar containing chloramphenicol. Viral diagnosis was performed on RNA from pharyngeal swabs collected in the first 24 hours following admission to the hospital by reverse transcription-polymerase chain reaction–based methods using reagents purchased from Roche™ (Swine Inf A/H1N1 detection set). These samples were also assessed by multiplex polymerase chain reaction (Luminex) with xTAG RVP kit from Luminex-Abbott for infection with Respiratory Syncytial Virus, Influenza B virus, Parainfluenza viruses 1–4, Human metapneumovirus, Enteroviruses, Rhinovirus, Adenovirus, Bocavirus and Coronaviruses NL63, HKU1, 229E and OC43, in accordance with manufacturer’s instructions.
Total RNA was extracted from blood samples using the PAXgene Blood RNA System (PreAnalytix, Hombrechtikon, Switzerland). RNA was quantified by spectrometry (NanoDrop ND1000, NanoDrop Technologies, Wilminton, Delaware USA) and quality confirmed by RNA Experion Bioanalyzer (BioRad, California USA) assay. Up to 1750 ng of each RNA sample was concentrated with the RNeasy MinElute Cleanup kit (QIAGEN, Hilden, Germany). RNA was eluted with 10 microliter of RNase-free H2O. 300 ng of purified total RNA were used to produce Cyanine 3-CTP-labeled cRNA using the Quick Amp Labeling kit (Agilent p/n 5190–0442) according to the manufacturer’s instructions. Following ‘One-Color Microarray-Based Gene Expression Analysis’ protocol Version 5.7 (Agilent p/n 4140–90040), 3 μg of labeled cRNA was hybridized with Whole Human Genome Oligo Microarray Kit (Agilent p/n G2519F-014850) containing 41,000+ unique human genes and transcripts. Arrays were scanned in an Agilent Microarray Scanner (Agilent G2565BA) according to the manufacturer’s protocol and data extracted using Agilent Feature Extraction Software 9.5.3 following the Agilent protocol GE1-v5_95_Feb07 and the QC Metric Set GE1_QCMT_Jan08. Resulting microarray data sets have been uploaded to the ArrayExpress microarray data repository [accession number: E-MEXP-3589]. Changes were verified in microarray gene expression for representative genes of our analysis by qPCR using Real-time Ready plates purchased to Roche, using β-actin and β2-microglobulin as reference (housekeeping) genes.
Data analysis was carried out by using GeneSpring GX 11.0 software. The original data was cleansed and normalized using the robust multichip average (RMA) algorithm consists of three steps: background correction, p75 normalization and expression calculation. Subsequent to logarithms transformation, baseline transformation of the data was performed using the median of control samples. Before statistical analyses all microarrays were subjected to quality and filtering criteria. Quality of the microarray data was assessed on Principal Component analysis (PCA) plots. All the 32 arrays passed these criteria and were included in the analyses. Student T tests (GeneSpring GX11.0) were used to identify genes differentially expressed between ICU and non-ICU groups at the level of significance p < 0.01 with Benjamini-Hochberg multiple testing corrections. GeneSpring GX 11.0 was used also for performing gene hierarchical clustering. Ingenuity pathway analysis 8.5 (IPA) (Ingenuity Systems, Redwood City, CA) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to select, annotate and visualize genes by function and pathway (gene ontology).
Gene expression analysis
T-test evidenced 1501 genes differentially expressed between ICU and non ICU patients. ICU patients showed higher expression levels of 558 genes and lower expression levels of 943 than non critically ill patients.
Genes showing higher levels of RNAm transcripts in the ICU group
IPA analysis of the most representative biological functions revealed that those genes showing higher levels of their corresponding transcripts in the ICU group participate principally in four biological functions: [immune response] (Table 2), [severe acute respiratory syndrome], [respiratory infection] (Table 3) and [inflammatory response] (Table 4). Patients in the ICU group not only showed increased neutrophil counts in blood, but also higher levels of RNAm corresponding to a group of proteins known to participate in neutrophil-mediated antimicrobial defense and tissue injury (Figure 1).
Genes showing lower levels of RNAm transcripts in the ICU group
IPA analysis revealed that those genes showing lower levels of their corresponding transcripts in the ICU group participate in the immune response against pathogens: [Natural Killer Cell Signaling], [T Cell Receptor Signaling], [Regulation of IL-2 Expression in Activated and Anergic T Lymphocytes], [Transforming Growth FActor-β Signaling], [CCR5 Signaling in Macrophages], [CD28 Signaling in T Helper Cells] (Additional file 1: Table S1 and Additional file 2: Table S2).
Real time qPCR results
Differential expression was confirmed in samples from all patients by qPCR. Differences between gene expression levels in the qPCR assays were assessed by performing the Mann–Whitney U test and by using SPSS 20.0 software analysis (significance was fixed at the level p < 0.05), as shown in Additional file 3: Figure S1.
Neutrophil related genes
ELANE (Elastase, neutrophil expressed), CTSG (cathepsin G), and PRTN3 (proteinase 3) are neutrophil serine proteases (NSPs) produced during neutrophil development in bone marrow and stored in the azurophilic granules of mature neutrophils. NSPs contribute to the neutrophil oxygen-independent system-mediated protection of the host against invading pathogens. On the other hand, NSPs play a critical role in neutrophil-associated lung inflammatory and tissue-destructive diseases, including COPD. NSPs have a broad substrate specificity and degrade a variety of extra-cellular matrix proteins including elastin, collagen (type I–IV), fibronectin, laminin and proteoglycans. High concentrations of NSPs are found in purulent secretions of COPD patients. Intra-tracheal instillation of mice with human neutrophil elastase or proteinase 3 leads to tissue destruction and airspace enlargement. Neutrophil elastase may also impair host defence interfering with muco-ciliary clearance of bacteria and phagocytosis of pathogens. In turn, both elastase and cathepsin G impair T-cell function through cleavage of CD2, CD4 and CD8 on the surface of T-cells. Patients lacking α1-Pi, the main physiological inhibitor of neutrophil elastase, are at greater risk of developing emphysema. CTSD (Cathepsin D) is also present in azurophil granules. By using a mice model, Bracke et al.. have recently demonstrated that cigarette smoke induces the expression of CTSD in pulmonary macrophages, supporting the role of cathepsins in the respiratory compromise associated to COPD . MPO (Myeloperoxidase) constitutes the major component of neutrophil azurophilic granules. MPO is important in bacterial killing, but also drives inflammatory reactions and tissue oxidation. By using a mice model of influenza infection, Sugamata R et al.. found that the absence of MPO reduced inflammatory damage with suppression of leakage of total proteins in bronchoalveolar lavage fluid associated with alteration of claudins in the lung. Activated neutrophils (MPO + cells) are found in severe COPD. In this sense, inhibition of MPO may be a novel and useful therapeutic treatment for COPD. In turn, AZU1 (azurocidin 1), BPI (bactericidal/permeability-increasing protein) and alpha-defensins (such as DEFA3, defensin, alpha 3, neutrophil-specific), are as NSPs, major constituents of neutrophil azurophilic granules. The protein encoded by AZU1 is an antibiotic protein, with monocyte chemotactic and antibacterial activity. It is also an important multifunctional inflammatory mediator. Bactericidal/permeability-increasing protein (BPI) encodes a lipopolysaccharide binding protein with bactericidal activity on gram-negative organisms. Defensins show activities against Gram-positive and Gram-negative bacteria, fungi, yeast, and enveloped viruses, but, as NSPs, also play important roles in promoting inflammation in the lungs, potentially contributing to lung injury . Along with NSPs, defensins may also be released upon neutrophil stimulation. Both have been described to affect the integrity of the epithelial layer, decrease the frequency of ciliary beat, increase the secretion of mucus, and induce the synthesis of epithelium-derived mediators that may influence the amplification and resolution of neutrophil-dominated inflammation.
A number of other gene transcripts related to neutrophils were identified in our analysis as more represented in the ICU group: LTF (lactotransferrin), TCN1 (transcobalamin I (vitamin B12 binding protein, R binder family)) as proteins which form part of secondary granules in neutrophils. S100A12 (S100 calcium binding protein A12) which codifies a protein proposed to be involved in specific calcium-dependent signal transduction pathways and which regulatory effect on cytoskeletal components may modulate various neutrophil activities. STXBP2 (Syntaxin binding protein 2) is involved in neutrophil degranulation. MPP1 (Membrane protein, palmitoylated 1, 55 kDa), which participates in the regulation of neutrophil chemotaxis. FCAR (Fc fragment of IgA, receptor for): this protein interacts with IgA-opsonized targets and triggers several immunologic defense processes, including phagocytosis. ITGAX (Homo sapiens cDNA, FLJ99683) which mediates adherence of neutrophils and monocytes to stimulated endothelium cells, and in the phagocytosis of complement coated particles. FCAR: this gene encodes a receptor for the Fc region of IgA, present on the surface of myeloid lineage cells such as neutrophils, monocytes, macrophages, and eosinophils, and triggers several immunologic defence processes, including phagocytosis, antibody-dependent cell-mediated cytotoxicity, and stimulation of the release of inflammatory mediators.
This “neutrophil signature” supports the notion of the existence at the systemic level of neutrophils ready to fight against infection but also ready to produce mediators which are able to induce tissue injury. Participation of these molecules in COPD pathogenesis had already been described at the respiratory level, but this work is the first suggesting an important activity of these enzymes at the systemic level. Migration of these cells from the blood vessels could certainly help to clear the causative microbe, but also to damage the pulmonary parenchyma, contributing to the explanation of the severity of respiratory conditions of these patients. Finding this “neutrophil signature” was possible since we used whole blood instead of PBMCs for the transcriptomic study. This signature is probably related to the fact that the vast majority of patients in the ICU group showed bacterial infection, and neutrophils are major actors in antibacterial defence. In fact, critically ill patients in our cohort showed significantly higher total counts of neutrophils in blood than non-critically ill ones. On the other hand, a limitation of our work is that we did not isolated neutrophils for the gene expression profiling assays. Further works targeting specific cell types would contribute to the definition of the exact contribution of neutrophils and other leukocytes to the gene expression signatures linked to severe respiratory illness in patients with COPD.
Other genes showing increased expression in the ICU group
C8G (Complement component 8, gamma polypeptide) forms part of C8, protein which plays a central role in assembly of the "membrane attack complex" (MAC) of complement. The MAC is a macromolecular pore that targets and lyses pathogens that challenge the host. Deposition of these pores on human cells contributes to tissue damage. Interestingly, MPO products seem to be able to activate complement, indicating a potential link between MPO over-expression and induction of MAC formation. CD24: This gene encodes a sialoglycoprotein that is expressed on mature granulocytes and in many B cells. Other interesting gene relatively up-regulated in critical patients was VEGFA (Vascular endothelial growth factor A). This protein is a glycosylated mitogen that specifically acts on endothelial cells and has various effects, including mediating increased vascular permeability, inducing angiogenesis, vasculogenesis and endothelial cell growth, promoting cell migration, and inhibiting apoptosis. Circulating levels of VEGF are up-regulated in patients with acutely exacerbated COPD and decrease after recovery from exacerbation. Respiratory levels of VEGF have been shown to negatively correlate with pulmonary function in stable COPD, which suggests its important role in COPD airway remodelling. ELN (Elastin): This gene encodes a protein that is one of the two components of elastic fibbers, and participates in extracellular matrix organization, such as those taking place after severe tissue injury. MUC1 (Mucin 1, cell surface associated) play an essential role in forming protective mucous barriers on epithelial surfaces and in the response to hypoxia. Finally, we identified increased transcript levels in the ICU group of a group of genes involved in antibacterial defence: PGLYRP1 (peptidoglycan recognition protein 1), CEBPE (CCAAT/enhancer binding protein (C/EBP), epsilon), CHIT1 (chitinase 1 (chitotriosidase)) and CYBA (cytochrome b-245, alpha polypeptide).
Genes showing decreased expression in the ICU group
On the other hand, the significantly lower counts of lymphocytes and monocytes in the blood of critically ill patients compared to non-critically ill ones probably explains the relatively depressed expression of immune related genes found in the former group. Migration of these cells to the site of infection (lung) or increased apoptosis in this group of patients could contribute to the observed lymphopenia and monocytopenia. Alternatively, the presence of an increased nitro-oxidative stress environment and the critical condition of the most severe patients could down-modulate the expression of genes involved in the response against pathogens.
Availability of supporting data
The data sets supporting the results of this article are available in the ArrayExpress repository; reference number E-MEXP-3589.
COPD: Chronic obstructive pulmonary disease; IPA: Ingenuity pathway analysis; KEGG: Kyoto encyclopedia of genes and genomes; CTSG: Cathepsin G; ELANE: Elastase, neutrophil expressed; PRTN3: Proteinase 3; MPO: Myeloperoxidase; CTSD: Cathepsin D; DEFA3: Defensin, alpha 3, neutrophil-specific; AZU1: Azurocidin 1; BPI: Bactericidal/permeability-increasing protein; NSPs:Neutrophil serine proteases.
The authors declare that they have no competing interests
LS, MSG, MDO, DAO, FB, AH, SR, VR, JB participated in patient recruitment and assisted in the analysis, interpretation of data, and writing the report. RA developed the bioinformatics analysis. LR performed the gene expression profiling works. IML and JFBM assisted in the design of the study, coordinated patient recruitment, analyzed and interpreted the data, and assisted in writing the paper. CASAJ introduced the clinical data and helped with data interpretation. MCG, VF performed the microbiology diagnosis and assisted in writing the paper. ROL and CD assisted in writing the paper. All authors read and approved the final manuscript.