- Open Access
Effect of short-term oral prednisone therapy on blood gene expression: a randomised controlled clinical trial
Respiratory Research volume 20, Article number: 176 (2019)
Effects of systemic corticosteroids on blood gene expression are largely unknown. This study determined gene expression signature associated with short-term oral prednisone therapy in patients with chronic obstructive pulmonary disease (COPD) and its relationship to 1-year mortality following an acute exacerbation of COPD (AECOPD).
Gene expression in whole blood was profiled using the Affymetrix Human Gene 1.1 ST microarray chips from two cohorts: 1) a prednisone cohort with 37 stable COPD patients randomly assigned to prednisone 30 mg/d + standard therapy for 4 days or standard therapy alone and 2) the Rapid Transition Program (RTP) cohort with 218 COPD patients who experienced AECOPD and were treated with systemic corticosteroids. All gene expression data were adjusted for the total number of white blood cells and their differential cell counts.
In the prednisone cohort, 51 genes were differentially expressed between prednisone and standard therapy group at a false discovery rate of < 0.05. The top 3 genes with the largest fold-changes were KLRF1, GZMH and ADGRG1; and 21 genes were significantly enriched in immune system pathways including the natural killer cell mediated cytotoxicity. In the RTP cohort, 27 patients (12.4%) died within 1 year after hospitalisation of AECOPD; 32 of 51 genes differentially expressed in the prednisone cohort significantly changed from AECOPD to the convalescent state and were enriched in similar cellular immune pathways to that in the prednisone cohort. Of these, 10 genes including CX3CR1, KLRD1, S1PR5 and PRF1 were significantly associated with 1-year mortality.
Short-term daily prednisone therapy produces a distinct blood gene signature that may be used to determine and monitor treatment responses to prednisone in COPD patients during AECOPD.
Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality in the United States (US) and elsewhere. In the US, COPD is now the 3rd leading cause of morbidity and mortality, trailing only ischemic heart disease and lung cancer . Most COPD deaths occur during periods of acute exacerbations (AECOPDs), which are characterised by an abrupt increase in dyspnea and cough related most commonly to an acute respiratory tract infection. AECOPDs are common indications for hospitalisation. In the US each year, there are more than 225,000 hospitalisations related to AECOPD . Most national and international guidelines and strategic documents recommend the use of systemic corticosteroids with or without antibiotics for serious AECOPD. Although on average systemic corticosteroids reduce the length of stay in hospital (by approximately 1 day) and decrease the risk of treatment failure, there is tremendous variation in response across patients with some experiencing benefits while others experiencing only harm . Moreover, even with “successful” therapy of their AECOPDs, a third of patients treated with systemic corticosteroids experience a treatment relapse (i.e. re-exacerbation) within 6 months  and 10–30% die within 1 year of their major AECOPD [5, 6]. The use of systemic corticosteroids during AECOPD is also associated with adverse events including hyperglycemia (which impacts ~ 50% of patients), hypertension (which impacts ~ 15% of patients) and adrenal insufficiency (which impacts 2 to 10% of patients) [3, 4]. Despite these challenges and limitations of systemic corticosteroids, they have been the mainstay of AECOPD management over the past 30 years. In the era of Precision Medicine, it is essential that we identify easily accessible biomarkers that can guide the use of systemic corticosteroids to patients who will benefit the most from these therapies. Here, we report on findings from 2 studies: 1) a clinical trial that determined which genes in peripheral blood were most responsive to short-term oral prednisone therapy in patients with COPD (i.e., discovery of blood biomarkers related to prednisone therapy) and 2) a clinical study which related the most responsive genes from the clinical trial to 1 year mortality in patients who experienced a serious AECOPD and were treated with systemic corticosteroids in hospital (i.e., discovery of predictive biomarkers of prednisone). Through these combined studies, we sought to identify genetic biomarkers that could be used to predict therapeutic benefits of systemic corticosteroids during AECOPD.
Study 1: a clinical trial to determine which genes in peripheral blood of COPD patients are responsive to oral prednisone therapy (NCT02534402)
To determine which genes were responsive to oral prednisone therapy in patients with COPD, we performed a randomised clinical trial (NCT02534402). We enrolled 40 patients from the COPD Clinic at St. Paul’s Hospital (SPH) in Vancouver, Canada. All patients had a clinical diagnosis of COPD from a board-certified pulmonologist and demonstrated a forced expiratory volume in 1 s (FEV1) to forced vital capacity (FVC) ratio of less than 70% (post-bronchodilator). The other enrollment criteria were that participants had to be free of exposure to prednisone or any other systemic corticosteroids and clinically stable (i.e., free of AECOPD) for at least 2 weeks prior to study entry. Following informed consent, patients were randomised in a 3: 1 ratio to either prednisone (30 mg/d for 4 days) + standard therapy or standard therapy alone. Participants attended 5 study visits over 5 consecutive days. The prednisone group received 30 mg/d of prednisone from days 1 through 4. Compliance was directly observed by the study personnel. During the trial period, patients were not permitted to modify any of their inhaled therapies including the dosage of inhaled corticosteroids or be placed on oral antibiotics.
During each of the visits, blood samples were collected in PAXgene (PreAnalytix, Switzerland) and EDTA tubes just before dispensation of the oral prednisone tablets. The PAXgene samples were placed in a -80C freezer. The EDTA blood was spun down within 2 h of collection and plasma and buffy coat were aliquoted and stored at -80C until analysis. Total RNA was extracted from PAXgene tubes using a PAXgene Blood miRNA kit from PreAnalytix, according to the manufacturer’s instructions; mRNA abundance was measured using Human Gene 1.1 ST 96-well array plates (Affymetrix, United States) at the McGill University and Genome Quebec Innovation Centre.
Study 2: rapid transition program cohort (NCT02050022)
To determine which if any of the steroid-responsive genes discovered in study 1 could be replicated and which of these replicated genes related to 1-year mortality among patients who experienced a severe AECOPD, requiring hospitalisation, we used data from the Rapid Transition Program (RTP) cohort. The details of the RTP cohort have been previously reported . Briefly, RTP included patients admitted for AECOPD as the Principal (or Most Responsible) Diagnosis at Vancouver General Hospital (VGH) or St. Paul’s Hospital (SPH) in Vancouver, Canada. The diagnosis of AECOPD was made based on clinical judgement by board-certified pulmonologists or general internists who were the most responsible physicians for these patients, and later confirmed by an independent board-certified pulmonologist, who did not participate in the care of these patients. During the hospitalisation, all patients were treated with at least 4 days of systemic corticosteroids (at the discretion of the attending physicians). Blood samples for this study were collected at baseline (within first 3 days of hospitalization for AECOPD) and between day 30 and 90 post-hospitalisation when the patients were completely stable and their symptoms had returned to baseline levels (i.e., convalescent sample). Similar to study 1, all blood samples were collected in PAXgene and EDTA tubes and were stored at -80C until processing. CBC and differential were performed on all samples collected in EDTA tubes.
Statistical analyses were performed in R. Patient characteristics were assessed for statistical significance using a Student’s t-test and a Kruskal-Wallis test for continuous variables, and a Fisher’s exact test for categorical variables. The oligo package was used to assess the quality of CEL files based on the normalized unscaled standard error (NUSE) and relative log-expression (RLE) metrics, and perform normalization using Robust Multi-array Average. Probes were summarized at the transcript cluster level. The ComBat algorithm was used to correct for plate-to-plate technical variation. The resulting batch-corrected, transcript cluster-level data was annotated with Human Gene 1.1 ST transcript cluster annotations (v36). Transcripts that were not annotated or mapped to more than one gene, were removed. If a gene had multiple transcripts assigned to it, these transcripts were averaged into a single value per sample. The end result of these steps is referred to as the gene expression data.
Diffearential expression was assessed using the limma package, which performs a moderated t-test using empirical Bayes. All analyses were adjusted for white blood cell counts and percentage of lymphocytes, neutrophils, monocytes, basophils, and eosinophils and were blocked by subject (i.e., they were modeled as random effects). Genes with a false discovery rate (FDR) < 0.05 between prednisone-treated versus untreated samples were considered differentially expressed (“prednisone genes”). These were then tested for functional enrichment analysis using the web-based gene set analysis toolkit (WebGestalt) . Gene symbols were used as input for WebGestalt and tested for enrichment in Gene Ontology (GO) processes  and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways . Enrichment was analyzed by a hypergeometric model and adjusted for multiple testing, with a FDR < 0.05 being considered significant.
The prednisone genes were replicated in the RTP cohort by comparing their convalescent gene expression to baseline AECOPD gene expression values. A FDR < 0.05 was considered significant in the replication dataset. The replicated genes were then used in a time-to-event analysis, to determine if any of these were predictive of mortality in the RTP cohort. Baseline AECOPD gene expression was subtracted from convalescent gene expression levels (to obtain “deltas”, or changes in these genes) on a log 2 scale, corresponding to a ratio on the intensity scale. These deltas were tested for their ability to predict mortality in a time-to-event analysis, for each gene individually. This was performed by fitting Cox proportional-hazards models (survival package) to predict number of days between AECOPD and death. All analyses were adjusted for white blood cell counts and percentage of lymphocytes, neutrophils, monocytes, basophils, and eosinophils. Genes with a p-value < 0.05 were considered significant in this time-to-event analysis (“mortality genes”). These mortality genes were combined together in a single Cox proportional-hazards mode.
Study 1: a clinical trial to determine which genes in peripheral blood of COPD patients are responsive to oral prednisone therapy
A flow diagram of patient recruitment is shown in Fig. 1; 30 patients were randomised to the prednisone arm and 10 were allocated to the control group. After excluding 2 patients in the prednisone group and 1 patient in the control group because of incomplete sample data collection or poor microarray data quality, there were 28 in the prednisone and 9 patients in the control group in the final analysis. The demographic and clinical data of the cohort stratified by prednisone treatment status are summarized in Table 1. There were no significant differences between the two groups in terms of age, gender, ethnicity, smoking status, a history of asthma or complete blood count (CBC) including its differential white blood cell count. However, patients in the prednisone group had a lower FEV1, % predicted (p = 0.007), FVC, % predicted (p = 0.020) and were slightly more likely to have used inhaled corticosteroids at study entry (p = 0.041).
There was no significant differences in the (baseline) day 1 transcriptomic expression of any of the measured genes at a FDR < 0.05 between the prednisone and control groups. However, in the prednisone arm, at day 3 and 5, there were 12 and 17 genes, respectively, that were differentially expressed at a FDR < 0.05 when compared with day 1, after adjusting for the total number of white blood cells and its differential cell count (Additional file 1: Table S1 and Additional file 2: Table S2). There were no genes that were differentially expressed between day 3 and 5 among those who received daily prednisone therapy. Therefore, for the remaining analysis, we pooled the gene expression from these two timepoints together to increase statistical power. When we compared days 3 and 5 from the prednisone group to day 1 from both the prednisone and control groups, there were 857 genes were differentially expressed at a FDR < 0.05 without adjusting for total number of white blood cells and their differential cell count, and only 51 genes with adjustments (Fig. 2).
Among these 51 differently expressed genes in the adjusted analysis at a FDR < 0.05, 33 genes were down-regulated and 18 genes were up-regulated. The gene with the greatest magnitude of fold-change was killer cell lectin-like receptor subfamily F, member 1 (KLRF1; fold change [FC], 1.73). The gene with the second greatest fold-change was granzyme H (GZMH; FC, 1.64), followed by adhesion G protein-coupled receptor G1 (ADGRG1; FC, 1.58) and killer cell lectin-like receptor subfamily D, member 1 (KLRD1, FC; 1.54) (Table 2). All 51 genes are shown in Additional file 3: Table S3.
Enrichment analysis was undertaken using prednisone genes at a FDR < 0.05 (Table 3), resulting in 21 genes that were significantly enriched in 11 GO biological processes. All of these gene were related to the immune system except one that was related to the circadian rhythm. The highest regulated GO processes were immune responses (p = 2.57E-07), followed by cytolysis (p = 1.19E-06), and cellular defense response (p = 2.87E-05). KEGG pathway analysis also significantly included 4 pathways; graft-versus-host disease (p = 1.43E-05), natural killer (NK) cell mediated cytotoxicity (p = 1.80E-05), neuroactive ligand-receptor interaction (p = 0.0037) and metabolic pathway (p = 0.0392). All genes in the pathway of NK cell mediated cytotoxicity including KLRD1, PRF1, GZMB and NCR1 were down-regulated by the use of prednisone (Additional file 4: Figure S1).
Study 2: rapid transition program (RTP) cohort
Gene expression data from 218 patients were available at both baseline and convalescence. Baseline characteristics of patients are shown in Table 1. We used a FDR cutoff of 0.05 and determined which of these differentially expressed genes from study 1 (at a FDR of 0.05) were replicable in the Rapid Transition Program cohort. In study 1, there were 51 genes there were differentially expressed between the prednisone and control groups at a FDR < 0.05; of these, 32 genes were replicated at a FDR < 0.05 in study 2 (Additional file 5: Table S4). Enrichment analysis was also undertaken using these 32 genes, resulting in genes significantly enriched in 3 GO process. Similar to study 1, the highest regulated GO process was cytolysis (p = 1.57E-07), followed by cellular defense response (p = 0.0002), and immune response (p = 0.0005). KEGG pathway analysis resulted in 4 pathways including NK cell mediated cytotoxicity (p = 0.0001) (Additional file 6: Table S5).
Twenty-seven patients (12.4%) died within 1 year after hospitalisation related to AECOPD. Among the 32 replicated genes, changes in expression of 10 genes were significantly associated with 1-year mortality in a time-event-analysis (Table 4). Survival curves based on theoretical patients whose gene deltas lie at the 10th and 90th risk percentiles are shown in Fig. 3.
Prednisone is used to treat moderate to severe COPD exacerbations but its effect on blood gene expression in COPD patients is largely unknown. Here, we showed in study 1 that short-term (4 day) prednisone therapy significantly modifies the expression of 51 genes in peripheral blood at a FDR < 0.05 (33 down-regulated and 18 up-regulated) and that these prednisone genes are significantly related to biological processes of immune systems including in the NK cell mediated cytotoxicity. Importantly, in study 2, we replicated 32 of 51 prednisone genes in a cohort of COPD patients hospitalised and treated with systemic corticosteroids for their AECOPD. Of these 32 genes, changes in the expression level of 10 genes were significantly associated with a mortality over 1 year.
The airways of COPD patients are characterized by the activation of inflammatory genes and infiltration of immune cells, which increases during acute exacerbations. The expression of inflammatory genes in the airway of COPD patients is controlled by pro-inflammatory transcription factors such as nuclear factor (NF)-κB . Corticosteroids, on the other hand, repress the expression of inflammatory genes by inducing genes that encode for anti-inflammatory proteins such as inhibitor of NF-κB and by reversing histone acetylation of activated inflammatory genes via histone deacetylases (HDAC) 2 [12, 13].
Interestingly, in the differential analysis (between prednisone and control groups) in study 1, we found that a majority of genes that were significantly down-regulated at a FDR < 0.05 were those that were involved in NK cell function; KLRF1 (NKp80), ADGRG1, KLRD1 (CD94), NKG7, NCR1 (NKp46) and CX3CR1 or were cytotoxic mediators; GZMH, GZMA, GZMB, GNLY and PRF1. Furthermore, in enrichment analysis, prednisone therapy affected the pathway related to cellular immunity and cytolysis including NK cell mediated cytotoxicity. Similarly, in study 2, 11 genes modified by prednisone therapy during AECOPD were enriched in cellular immunity and NK cell mediated cytotoxicity pathways.
NK cells participate in innate immunity as a first line of defense against tumorigenesis or intracellular pathogens including viruses. This cytotoxic activity of NK cells is regulated by diverse combination of activating and inhibitory receptor-ligand interaction [14,15,16,17,18,19,20]. Granzyme is a family of serine proteases of which there are five known members in humans: A, B, H, K and M  and perforin is a pore-forming protein facilitating the entry of granzyme B into the target cell through a polyperforin pore . Cytotoxic T lymphocytes and NK cells release inflammatory cytokines and cytolytic proteins to destroy target cells by inducing them to undergo apoptosis [18, 22, 23]. However, it is well-known that the predominant cells in the airways of COPD patients during AECOPD are neutrophils . Our study did not find any significant effect of prednisone on genes involved in neutrophil-related pathways. This may explain why systemic corticosteroids are only modestly effective in improving patient outcomes in AECOPD.
In the survival analysis, KLRD1 and PRF1 genes, which are involved in the NK cell mediated cytotoxicity pathway was significantly related to 1-year mortality following AECOPD. We also found that CX3CR1 and S1PR5 genes, which encode receptors expressed on a variety of immune cells including NK cell and macrophages and regulate essential cellular processes such as cell proliferation, migration, adhesion and phagocytosis [24,25,26], were also significantly associated with 1-year mortality. CX3CR1 may play an important role in lung inflammation and is involved in emphysema pathogenesis related to cigarette smoke . S1PR5, on the other hand, is expressed on alveolar macrophages of COPD patients and may regulate phagocytosis [24, 27]. Importantly, we showed that prednisone down-regulated the expression levels of all 10 genes, which, in turn, were associated with increased mortality. This raises the possibility that this 9 gene expression signature related to prednisone may be a biomarker of poor responsiveness or even harm related to systemic corticosteroids during AECOPD. However, these latter findings should be considered preliminary owing to the small sample size; additional validation is required in a larger cohort.
There were several limitations of this study. First, in study 1, we observed an up-regulation of NFKB2 which is counter to the known anti-inflammatory effects of prednisone (Additional file 3: Table S3). NFKB2 is a transcription factor which increases expression of genes, which promote inflammation. We postulate that the up-regulation in the level of expression of this gene related to prednisone therapy may be through a negative or indirect feedback mechanism by modulating the mitogen-activated protein kinase pathways or by activating inhibitors of NFKB . Second, gene signature is tissue-specific, and thus our results cannot be generalised to lung tissue or airway epithelial cells. However, the strength of this study included the evaluation of whole blood which is easily accessible, enabling translation of these results to the clinic. Third, our gene signature related to short-term prednisone therapy was not validated at the protein level. Fourth, owing to the small sample size, we could not determine whether smoking modifies the effects of prednisone on gene expression. Previous studies have shown that smoking impairs therapeutic responses of patients to inhaled and systemic corticosteroids [28, 29]. In both study 1 and 2, we did not find any genes that were differentially expressed to patients’ smoking status at a FDR < 0.05.
The current study has shown that short-term systemic corticosteroid therapy produces a distinct gene expression signature in blood of patients with COPD. A subset of these genes significantly relates to total mortality following COPD hospitalisation. Although our findings are preliminary and requires validation in larger studies, they raise the possibility that blood genomic signatures may be useful in identifying subgroups of patients who may experience benefits and others who may experience harm from systemic corticosteroid therapy.
Availability of data and materials
The datasets used for the current study are available from the authors on reasonable request.
Acute exacerbation of chronic obstructive pulmonary disease
Complete blood count
Chronic obstructive pulmonary disease
False discovery rate
Forced expiratory volume in 1 s
Forced vital capacity
The Global Initiative for Chronic Obstructive Lung Disease
Kyoto Encyclopedia of Genes and Genomes
Rapid Transition Program
Mokdad AH, Ballestros K, Echko M, Glenn S, Olsen HE, Mullany E, Lee A, Khan AR, Ahmadi A, Ferrari AJ, et al. The state of US health, 1990-2016: burden of diseases, injuries, and risk factors among US states. Jama. 2018;319:1444–72.
COPD Hospitalizations https://ephtracking.cdc.gov/showIndicatorPages.action?selectedContentAreaAbbreviation=23&selectedIndicatorId=129&selectedMeasureId. Accessed 20 Jan 2019.
Quon BS, Gan WQ, Sin DD. Contemporary management of acute exacerbations of COPD: a systematic review and metaanalysis. Chest. 2008;133:756–66.
Leuppi JD, Schuetz P, Bingisser R, Bodmer M, Briel M, Drescher T, Duerring U, Henzen C, Leibbrandt Y, Maier S, et al. Short-term vs conventional glucocorticoid therapy in acute exacerbations of chronic obstructive pulmonary disease: the REDUCE randomized clinical trial. Jama. 2013;309:2223–31.
Garcia-Sanz MT, Canive-Gomez JC, Senin-Rial L, Aboal-Vinas J, Barreiro-Garcia A, Lopez-Val E, Gonzalez-Barcala FJ. One-year and long-term mortality in patients hospitalized for chronic obstructive pulmonary disease. J Thorac Dis. 2017;9:636–45.
Ho TW, Tsai YJ, Ruan SY, Huang CT, Lai F, Yu CJ. In-hospital and one-year mortality and their predictors in patients hospitalized for first-ever chronic obstructive pulmonary disease exacerbations: a nationwide population-based study. PLoS One. 2014;9:e114866.
Leung JM, Chen V, Hollander Z, Dai D, Tebbutt SJ, Aaron SD, Vandemheen KL, Rennard SI, FitzGerald JM, Woodruff PG, et al. COPD exacerbation biomarkers validated using multiple reaction monitoring mass spectrometry. PLoS One. 2016;11:e0161129.
Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017;45:W130–W37.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:D277–80.
Di Stefano A, Caramori G, Oates T, Capelli A, Lusuardi M, Gnemmi I, Ioli F, Chung KF, Donner CF, Barnes PJ, Adcock IM. Increased expression of nuclear factor-kappaB in bronchial biopsies from smokers and patients with COPD. Eur Respir J. 2002;20:556–63.
Barnes PJ. Corticosteroid effects on cell signalling. Eur Respir J. 2006;27:413–26.
Barnes PJ, Karin M. Nuclear factor-kappaB: a pivotal transcription factor in chronic inflammatory diseases. N Engl J Med. 1997;336:1066–71.
Chang GW, Hsiao CC, Peng YM, Vieira Braga FA, Kragten NA, Remmerswaal EB, van de Garde MD, Straussberg R, Konig GM, Kostenis E, et al. The adhesion G protein-coupled receptor GPR56/ADGRG1 is an inhibitory receptor on human NK cells. Cell Rep. 2016;15:1757–70.
Crooks SW, Bayley DL, Hill SL, Stockley RA. Bronchial inflammation in acute bacterial exacerbations of chronic bronchitis: the role of leukotriene B4. Eur Respir J. 2000;15:274–80.
Guma M, Busch LK, Salazar-Fontana LI, Bellosillo B, Morte C, Garcia P, Lopez-Botet M. The CD94/NKG2C killer lectin-like receptor constitutes an alternative activation pathway for a subset of CD8+ T cells. Eur J Immunol. 2005;35:2071–80.
Kuttruff S, Koch S, Kelp A, Pawelec G, Rammensee HG, Steinle A. NKp80 defines and stimulates a reactive subset of CD8 T cells. Blood. 2009;113:358–69.
Lanier LL. Up on the tightrope: natural killer cell activation and inhibition. Nat Immunol. 2008;9:495–502.
Middleton D, Curran M, Maxwell L. Natural killer cells and their receptors. Transpl Immunol. 2002;10:147–64.
Parham P, Abi-Rached L, Matevosyan L, Moesta AK, Norman PJ, Older Aguilar AM, Guethlein LA. Primate-specific regulation of natural killer cells. J Med Primatol. 2010;39:194–212.
Barry M, Bleackley RC. Cytotoxic T lymphocytes: all roads lead to death. Nat Rev Immunol. 2002;2:401–9.
Ngan DA, Vickerman SV, Granville DJ, Man SF, Sin DD. The possible role of granzyme B in the pathogenesis of chronic obstructive pulmonary disease. Ther Adv Respir Dis. 2009;3:113–29.
Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Innate immunity. In: Molecular biology of the cell. 4th ed. New York: Garland Science; 2002. p. 1616.
Barnawi J, Jersmann H, Haberberger R, Hodge S, Meech R. Reduced DNA methylation of sphingosine-1 phosphate receptor 5 in alveolar macrophages in COPD: a potential link to failed efferocytosis. Respirology. 2017;22:315–21.
McComb JG, Ranganathan M, Liu XH, Pilewski JM, Ray P, Watkins SC, Choi AM, Lee JS. CX3CL1 up-regulation is associated with recruitment of CX3CR1+ mononuclear phagocytes and T lymphocytes in the lungs during cigarette smoke-induced emphysema. Am J Pathol. 2008;173:949–61.
Sanchez T, Hla T. Structural and functional characteristics of S1P receptors. J Cell Biochem. 2004;92:913–22.
Barnawi J, Tran H, Jersmann H, Pitson S, Roscioli E, Hodge G, Meech R, Haberberger R, Hodge S. Potential link between the Sphingosine-1-phosphate (S1P) system and defective alveolar macrophage phagocytic function in chronic obstructive pulmonary disease (COPD). PLoS One. 2015;10:e0122771.
Burge PS, Calverley PM, Jones PW, Spencer S, Anderson JA. Prednisolone response in patients with chronic obstructive pulmonary disease: results from the ISOLDE study. Thorax. 2003;58:654–8.
Bhatt SP, Anderson JA, Brook RD, Calverley PMA, Celli BR, Cowans NJ, Crim C, Martinez FJ, Newby DE, Vestbo J, et al. Cigarette smoking and response to inhaled corticosteroids in COPD. Eur Respir J. 2018;51:1701393.
The study was funded by Genome Canada, Genome British Columbia, Genome Quebec, Canadian Institutes of Health Research, Providence Health Care, St Paul’s Hospital Foundation, and PROOF Centre. Dr. Don Sin is a Tier 1 Canada Research Chair in COPD and holds the DeLazzari Family Chair at HLI.
Ethics approval and consent to participate
Written informed consent was provided by each participant in accordance with the Research Ethics Board. Study 1 was approved by the University of British Columbia (UBC)-Providence Health Care (PHC) Research Ethics Board (H15–01147). Study 2 was approved by the UBC-PHC Research and UBC Clinical Research Ethics Boards (certificate numbers H11–00786 and H13–00790).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Genes differentially expressed by prednisone at a FDR < 0.05 after adjusting for the total number of white blood cells and its differential cell count (day 3 versus day 1 in prednisone group). (DOCX 17 kb)
Table S2. Genes differentially expressed by prednisone at a FDR < 0.05 after adjusting for the total number of white blood cells and its differential cell count (day 5 versus day 1 in prednisone group). (DOCX 16 kb)
Table S3. Fifty-one genes differentially expressed by prednisone at a FDR < 0.05 in study 1. (DOCX 18 kb)
Figure S1 KEGG pathway map of NK cell mediated cytotoxicity. Rectangle and circle represent gene product including RNA and a compound, respectively. The blue boxes are hyperlinked to KEGG orthology entries and the red boxes represent genes which were significantly enriched in this pathway. Among genes which were responsive to short-term prednisone at a FDR < 0.05, KLRD1, PRF1, GZMB and NCR1 were significantly enriched in the pathway of NK cell mediated cytotoxicity at a FDR < 0.05. (PPTX 92 kb)
Table S4. Genes replicated in the Rapid Transition Program (RTP) Cohort at a FDR < 0.05. (DOCX 16 kb)
Table S5. Enrichment analysis: Replicated biological process and pathway in study 2. (DOCX 17 kb)