Published ahead of print on August 21, 2008, doi:10.1164/rccm.200803-387OC
© 2008 American Thoracic Society doi: 10.1164/rccm.200803-387OC
Circulating RNA Transcripts Identify Therapeutic Response in Cystic Fibrosis Lung Disease1 Division of Pulmonary Sciences and Critical Care Medicine, University of Colorado School of Medicine, Denver, Colorado; 2 Department of Medicine, National Jewish Health, Denver, Colorado; and 3 Department of Preventive Medicine and Biometrics, and 4 Division of Pediatric Pulmonary Medicine, University of Colorado School of Medicine, Denver, Colorado Correspondence and requests for reprints should be addressed to Milene T. Saavedra, M.D., University of Colorado Health Sciences Center, 4200 E. Ninth Avenue, Box C272, Denver, CO 80262. E-mail: saavedram{at}njc.org
Rationale: Circulating leukocyte RNA transcripts are systemic markers of inflammation, which have not been studied in cystic fibrosis (CF) lung disease. Although the standard assessment of pulmonary treatment response is FEV1, a measure of airflow limitation, the lack of systemic markers to reflect changes in lung inflammation critically limits the testing of proposed therapeutics. Objectives: We sought to prospectively identify and validate peripheral blood leukocyte genes that could mark resolution of pulmonary infection and inflammation using a model by which RNA transcripts could increase the predictive value of spirometry. Methods: Peripheral blood mononuclear cells were isolated from 10 patients with CF and acute pulmonary exacerbations before and after therapy. RNA expression profiling revealed that 10 genes significantly changed with treatment when compared with matched non-CF and control subjects with stable CF to establish baseline transcript abundance. Peripheral blood mononuclear cell RNA transcripts were prospectively validated, using real-time polymerase chain reaction amplification, in an independent cohort of acutely ill patients with CF (n = 14). Patients who responded to therapy were analyzed using general estimating equations and multiple logistic regression, such that changes in FEV1% predicted were regressed with transcript changes. Measurements and Main Results: Three genes, CD64, ADAM9, and CD36, were significant and independent predictors of a therapeutic response beyond that of FEV1 alone (P < 0.05). In both cohorts, receiver operating characteristic analysis revealed greater accuracy when genes were combined with FEV1. Conclusions: Circulating mononuclear cell transcripts characterize a response to the treatment of pulmonary exacerbations. Even in small patient cohorts, changes in gene expression in conjunction with FEV1 may enhance current outcomes measures for treatment response.
Key Words: cystic fibrosis peripheral blood mononuclear cells biomarkers pulmonary exacerbation
Cystic fibrosis (CF) is the most common lethal inherited disease in the Western world, with respiratory failure accounting for more than 80% of deaths from the disease, usually in the third or fourth decade of life. The triad of mucus airway obstruction, chronic Pseudomonas aeruginosa infection, and severe airway inflammation are the major pathogenic factors in CF lung disease (1). There is a critical need for effective antimicrobial and antiinflammatory therapies to mitigate disease in this young population. The design of clinical trials in CF is hampered, in part, by the lack of sensitive measures of treatment response. The most established outcome measure for CF therapies is FEV1, a functional measure of airflow limitation, and a key consideration in the advancement of treatments from phase 2 to phase 3 trials. The process of airway remodeling by inflammatory cells is progressive and may not be detected by FEV1 alone in the timeframe of a typical phase 2 trial. To date, no systemic marker of treatment response has been validated in CF. The gene expression of peripheral leukocytes has never been studied as a predictor of treatment response in CF, yet peripheral blood mononuclear cells (PBMCs) serve as the source of the dense mononuclear infiltrate present at sites of CF airway injury and acquire a unique transcriptional "fingerprint" as a result of repeated passages through the pulmonary vascular bed. Furthermore, circulating leukocytes have demonstrated diagnostic and prognostic utility for multiple pulmonary conditions, including steroid sensitivity in asthma and discrimination between viral and bacterial pulmonary infections (2–6). A condensed model of deterioration in CF lung disease is represented by episodic pulmonary exacerbations. Primarily a clinical diagnosis, exacerbations are treated with antimicrobials, considered the standard of care to decrease pulmonary bacterial burden and host airway inflammation. The aggressive treatment of typical pulmonary exacerbations results in rapid reduction in pulmonary inflammation due to reduced infection (7–10). Response to treatment is assessed by improvements in spirometry and clinical evaluation (11). Because exacerbation frequency is positively correlated with respiratory decline (12–14), the assessment of the effect of therapies on reducing the incidence of pulmonary exacerbations is a key measure in clinical trials for evaluating the most important drugs used in CF care (15–18). In two separate cohorts of patients with CF, we identified three PBMC genes whose changes in expression correlate with the onset and resolution of acute pulmonary exacerbations. The genes consistently added independent predictive power to FEV1, and represent a novel tool to quantify CF therapeutic response noninvasively.
Patient Recruitment Subjects older than 18 years with CF (sweat chloride testing and genotype) were enrolled at the time of intravenous antibiotic initiation for a clinically diagnosed pulmonary exacerbation at a Cystic Fibrosis Foundation (CFF)–accredited adult CF clinic. Patients enrolled met CFF Clinical Practice criteria for an acute pulmonary exacerbation and were treated for a minimum of 2 weeks with intravenous antibiotics (19). This observational trial used within-subject comparisons, such that each subject served as his or her own control, after treatment with antibiotics. Blood was drawn at initiation (±2 d) and completion (±1 wk) of antibiotic therapy. At each time point, the following were sampled: (1) blood for PBMC isolation and whole-blood differentials, (2) sputum microbiology, (3) simple spirometry for FEV1, and (4) plasma cytokines. In addition, transcript abundance in study patients was compared with PBMC array data from eight healthy volunteer subjects and six control subjects with stable CF. Next, a validation study was performed on 14 adult patients with CF suffering from acute pulmonary exacerbations. Enrollment and assessment of blood, microbiology, and spirometry occurred as described above.
PBMC Isolation
Microarray Hybridization and Data Analysis
Quantitative RT-PCR Validation
Plasma Cytokine Measurements
Endpoints
Statistical Methods
Sequential Data Analysis Strategy The stepwise execution of this study is shown in Figure 1. Candidate genes, differentially regulated before and after antibiotic therapy in patients with CF and acute lung infection, were identified in 10 patients on the basis of data generated from high-density oligonucleotide array experiments on patients' PBMC samples. Nonparametric univariate tests ranked genes on their ability to discriminate between the two time points. PBMC transcript abundance from patients with stable CF and normal healthy control subjects, matched in age and sex to the CF population and prepared with identical method and array platforms, was included in the analysis. Samples underwent whole-blood differentials as well as cell counts after each isolation event. Although there was a trend toward decreased total white blood cell count after treatment of pulmonary exacerbations (P = 0.052), there were no statistically significant differences in median percentage values of neutrophils, lymphocytes, monocytes, eosinophils, and basophils, pre- and post-therapy, nor were there significant differences between PBMC absolute numbers pre– and post–antibiotic therapy, after isolation. Second, classifier genes were independently confirmed in the same set of samples using quantitative real-time RT-PCR. Of note, 6 of 10 patients in the development cohort had sufficient remaining RNA, after microarray analysis, for transcript measurements of 10 genes. Third, the resulting "CF gene signature," consisting of genes significantly changed in both microarray and RT-PCR experiments, was then tested for its ability to classify therapeutic response in an independent group of patients with CF undergoing therapy for acute pulmonary infections. Table 1 depicts measurement distributions performed for all 24 samples from patients with CF. Fourth, using FEV1 improvement as a standard measure of clinical response, logistic regression on the change in FEV1% predicted, with gene transcript changes, allowed association between measurements of lung function and genes in the responder groups for both cohorts. Finally, a comparison of the diagnostic value of genes in combination with FEV1 as markers, was depicted using ROC curves.
Baseline Patient Characteristics PBMCs were isolated from 22 adult patients with CF at the beginning and end of a course of antibiotic treatment for an acute pulmonary exacerbation, representing a total of 24 acute pulmonary exacerbations. The development cohort contained 10 patients, and a subsequent validation cohort consisted of 14 patients. Two patients in the development cohort suffered subsequent exacerbations and were included in the validation cohort 2 years later. All patients met CFF guidelines for an acute pulmonary exacerbation (19). The baseline demographics, severity of airflow limitation, genotype, and sputum microbiology for the development and validation cohorts are shown in Table 2. Both groups had moderate to severe airway disease by American Thoracic Society criteria, on the basis of FEV1% predicted measured at the completion of intravenous antibiotic therapy (20). Ninety percent of patients in both groups grew P. aeruginosa on sputum culture at the time of therapy. Staphylococcus aureus was commonly isolated in both groups, although only one strain was methicillin resistant. Whereas no patients were treated with steroids in the development cohort, over half of the validation cohort received steroids (P = 0.002). The median pretreatment FEV1% was lower in the validation group (P = 0.04); however, post-treatment FEV1% predicted did not significantly differ between groups. All patients exhibited FEV1 increases at the conclusion of therapy except for a single patient in the development cohort, whose FEV1% predicted declined by 26% after therapy. This patient concomitantly suffered a severe flare-up of CF-related arthritis at the end of her antibiotic therapy. Given the drop in FEV1, this patient was considered a nonresponder, and her gene copy changes were not used in the regression analyses for genes with FEV1. Clinical nonresponders for the second cohort are described below.
A Transcriptional Signature Characterizes the Effect of Antibiotic Therapy in Patients with CF with Acute Pulmonary Exacerbations Table 3 lists 10 genes significantly changed (P < 0.05) by both microarray and quantitative real-time RT-PCR analysis in the development cohort of patients with CF. After ontology analysis, 19 candidate genes identified by microarray analysis were evaluated by quantitative PCR in the same samples to validate results. Of 19 genes, 10 were significantly changed when measured by both methods. Of note, the real-time RT-PCR validation occurred in six pairs of RNA samples, because the other samples did not have sufficient remaining RNA after microarray to be evaluated by PCR. We designated these 10 PBMC genes, reflecting treatment of CF pulmonary exacerbations, the "CF therapeutic signature." The vast majority of genes (9 of 10) were down-regulated after resolution of the acute pulmonary exacerbation. Only IL-32 increased transcription after treatment, suggesting suppression during the acute exacerbation and a return to normal baseline after treatment, as indicated by transcripts in normal control subjects. In Figure 2, we compare pre- and post-antibiotic therapy median values and distribution for CF genes in the development cohort with array data from eight age- and sex-matched normal control subjects, as well as six control subjects with stable CF. Before antibiotic therapy, expression of all genes, except for TLR2, significantly differed from expression in normal subjects (ANOVA, P < 0.05; pairwise comparison with Fisher's protected least significant difference test [PLSD], P < 0.05) and expression of 5 of 10 genes (PLXND1, ADAM9, CSPG2, CD64, and CD36) differed from that in stable subjects with CF (ANOVA, P < 0.05; pairwise comparison with Fisher's PLSD, P < 0.05). After treatment of the acute pulmonary exacerbation, none of the 10 genes differed significantly between post-therapy patients and control subjects with stable CF. The following genes were significantly different between preantibiotic CF and stable CF and unchanged between postantibiotic CF and stable CF: PLXND1, ADAM9, CSPG2, CD64, and CD36 (using ANOVA and pairwise comparison with Fisher's PLSD). In post-therapy patients with CF, expression of three genes, CD64, ADAM9, and PLXND1, remained significantly different in patients with CF compared with normal control subjects (ANOVA, P < 0.05; pairwise comparison with Fisher's PLSD, P < 0.05).
The CF Therapeutic Signature and Associations with Common Clinical Outcome Measures Measurements of CRP and FEV1 were performed simultaneously with PBMC isolation in the development cohort. As expected, improvement in FEV1% predicted was statistically significant after treatment, (P = 0.02 by paired t test; 95% confidence intervals, –1.8 to –15.78) (Figure 3). The change in plasma CRP levels between pairs did not reach statistical significance (P = 0.09). The correlations between gene expression values and both FEV1 and CRP levels were evaluated with Spearman's rank order correlation coefficient (rs) in Table 4. A significant correlation exists between the development cohort expression values and CRP values for 7 of the 10 genes: TLR2, CD64, CSPG2, HPSE, IL-32, CD163, and ADAM9. No correlation was noted between FEV1 change and the 10-gene CF signature. Furthermore, we conducted correlation analyses between genes and circulating cellular markers of inflammation, namely, neutrophil counts and total white blood cell counts (Table 5). Four genes were highly correlated (P < 0.05) to circulating neutrophil counts before and after therapy: ADAM9, CSPG2, IL32, and CD163.
Plasma Cytokine Changes in the Resolution of the Acute Pulmonary Exacerbations in CF Multiple analytes were measured in aliquots of plasma with the Luminex system (Luminex Corp., Austin, TX), taken from the same blood specimens from which PBMCs were isolated before and after antibiotic therapy. Table 6 depicts mean cytokine measurements at both time points and their associated P values. After the Bonferroni's correction for multiple tests, the change in granulocyte colony–stimulating factor fails to reach significance at the P 0.05 level.
Prospective Evaluation of the CF Therapeutic Signature in a Validation Cohort of Patients with CF Expression of the 10 identified classifier genes was evaluated in an independent population of patients with CF and acute pulmonary exacerbations via real-time RT-PCR (n = 14). PBMC isolation and FEV1 measurements were performed on adult patients with CF at the initiation and termination of antibiotic therapy for acute pulmonary exacerbations using identical methods as for the development cohort. This validation group of patients was heterogeneous in terms of sputum microbiology and variable use of systemic steroids and was ultimately more representative of a realistic clinical setting (Table 2). Although many patients were treated with steroids, as opposed to none in the initial study group, the presence of steroids did not significantly alter gene expression between steroid-treated and steroid-naive patients for 9 of 10 genes (by multi-way ANOVA). Only transcript changes in HPSE significantly differed between steroid-treated and steroid-naive patients (P = 0.03). FEV1 improved significantly between pre- and post-measurement in this cohort (P = 0.003 by paired t test). All patients in the validation cohort had a higher FEV1% predicted at the second measurement. However, two patients in this cohort with severe airway limitation (FEV1 <25% predicted at both measurements) manifested an improvement of 100 ml or less in FEV1% predicted and were considered nonresponders based on clinical parameters. Both patients had poor outcomes: one was readmitted with recurrent pulmonary symptoms in 6 days and the second was never discharged due to a persistently poor clinical response to antibiotics and underwent lung transplant 43 days later. These nonresponders were not included in the logistic regression analysis of the validation cohort, because the model was predicated on identifying genes associated with therapeutic response. In a univariate analysis of individual gene changes within the validation cohort, five genes were significantly changed in responders after antibiotic therapy, based on log-transformed gene expression values measured by real-time RT-PCR: CD36, P = 0.002; CD64, P = 0.007; PLXND1, P = 0.01; CSPG2, P = 0.002; and TLR2, P = 0.05.
Diagnostic Value of the CF Therapeutic Signature in Association with Percentage of Change in FEV1
In CF lung tissue, mononuclear cells that originate from the peripheral blood are abundantly present at sites of airway injury. Histologic studies of CF airways highlight mononuclear cells as the predominant cell population in areas of cartilaginous destruction and identify lymphocytes as the majority cell type infiltrating airway submucosa (21). (22). The consistent change in expression of a small group of PBMC genes among patients with CF with heterogeneous genotypes, lung function, sputum microbiology, and antibiotic regimens creates a strong basis on which to study the participation of these cells in CF lung disease pathogenesis. This study is unique in several aspects. It is the first study of its kind to use circulating leukocyte transcripts to assess response to therapy in CF lung disease. By correlation with known inflammatory markers, CRP and circulating neutrophil counts, the gene predictors are plausible markers of inflammation. Second, the gene transcript changes demonstrate reproducibility across two patient groups and discriminatory capacity to differentiate between acutely ill and subsequently treated patients. Third, our regression model demonstrates that mRNA changes in circulating leukocytes add meaningful diagnostic information to FEV1 in assessing treatment response in acute pulmonary exacerbations.
Genes identified in this study may serve to generate additional hypotheses concerning disease pathogenesis, inflammatory regulation, and treatment response in the CF airway. Seven of the 10 genes (and the proteins they encode) featured in this article (IL32, HPSE, ADAM9, PLXND1, HCA112, CSPG2, and CD163) have not previously been linked to CF lung disease. It is important to note that these genes are not specific to CF and have varying roles in other conditions characterized by pathologic pulmonary inflammation, including asthma and pneumonia (23, 24). The fact that these genes encode for proteins implicated in inflammatory processes lends biologic plausibility to naming them as potential markers of the resolution of CF airway infection and inflammation. As a group, these genes represent functions of immune recognition and response, phagocytosis, and matrix degradation. TLR2 represents a central pattern recognition receptor for the innate response against bacterial infection. IL-32 is a newly described tumor necrosis factor–inducible intracellular cytokine (25). An inducer of proinflammatory cytokines, IL-32 induces blood monocyte differentiation to macrophages, with subsequent phagocytic activity for live bacteria (26). The expression of IL-32 by inflamed luminal epithelia may facilitate differentiation of blood monocytes infiltrating infected lung (26). Three surface-receptor genes participate in phagocytosis. CD64, or Fc This study was designed to measure markers that change with aggressive treatment of a pulmonary exacerbation, the current best therapy for reduction of acute increases in airway infection and inflammation in CF lung disease. This design has several advantages. First, our design closely parallels a clinical trial sequence in which a treatment would be tested for its effect on decreasing inflammation. Second, our gene signature is not pathogen limited. Our study population suffered from infection with a representative variety of bacterial pathogens, typical of the CF population as a whole. When patients with CF were treated for pulmonary exacerbations, expression of most of the candidate genes more closely resembled the normal control subjects (Figure 2), which supports the biologic roles of these genes as markers of decreased infection and inflammation. Furthermore, one-half of the genes were specific for exacerbation among patients with CF, meaning that, when compared with patients with stable CF, values differed significantly before treatment but not after treatment. This includes the three genes (CD64, ADAM9, and CD36) that were most highly diagnostic of therapeutic response. One of the genes that added significant explanatory power to the regression model in both patient groups, ADAM9, is highly representative of immediate postexacerbation transcriptional changes and is the only gene whose postexacerbation expression remained significantly different from both normal PBMC and control subjects with stable CF. Simultaneous measurements of respiratory physiology and plasma markers allowed for statistical comparisons of multiple outcomes measures. The significant correlation of more than one-half of the genes in the CF therapeutic signature with changes in CRP, a well-characterized serum marker of inflammation, in addition to correlations to neutrophil counts, further strengthens the association of these genes with inflammatory processes. The concomitant evaluation of a relatively comprehensive panel of plasma cytokines in the development cohort did not demonstrate statistical significance. Post-treatment cytokine measurements were lower than pretreatment measurements; however, the specific cytokines reduced differed from patient to patient, making none of them broadly applicable as a marker for exacerbation resolution in a cohort of this size. Although this study was not powered to detect a significant difference in plasma cytokine concentrations, our group and others are currently conducting larger-scale longitudinal studies of blood, sputum, and plasma markers in the resolution of acute exacerbations, which will allow simultaneous analysis of markers reflecting compartmental and circulating inflammation. The complexity of CF airway inflammation is such that local events within airways may be heterogeneous and compartmentalized across lung segments. Although PBMCs are easily accessible, they may not completely reflect inflammatory processes occurring in all pulmonary compartments. It would be of great interest to compare the response of PBMC RNA transcripts to direct measures of airway inflammation. Although many markers of airway inflammation can be detected in the sputum, multiple previous studies (particularly of this size) have not been able to detect changes in sputum concentrations of inflammatory markers associated with the treatment of an acute pulmonary exacerbation (33–38). Thus, to directly compare identified changes in PBMC expression with specific inflammatory markers in the sputum, sputum markers of inflammation were not analyzed as a part of this study.. The compartmental quality of local pulmonary events may also reflect, in part, the relative steroid independence of the gene expression change in the steroid-exposed validation cohort. In animal models, glucocorticoids typically suppress local blood flow at sites of inflammation, reduce vascular permeability, and reduce leukocyte migration. These responses would likely have greater effect on local transcription profiles at the site of inflammation rather than on circulating expression (39–41). In addition, whereas corticosteroid therapy did not affect expression of 9 of 10 genes, the study was not powered to detect significant steroid-modulated differences in gene expression. A systemic marker of lung inflammation has many advantages, because blood can be obtained from subjects of any age and disease severity, and may reflect the status of inflammation throughout the lung, rather than one segment, as is assessed by bronchoalveolar lavage. The current study establishes that the use of PBMC expression permits the assessment of leukocyte activity, concomitantly with the functional information represented by FEV1, to reflect treatment of an acute pulmonary exacerbation. This analysis is sensitive, inexpensive, and obtained from tissue that is easily accessible in pediatric and adult populations, and has the potential to be performed in a clinical laboratory. The identification of relevant biomarkers could have more immediate clinical implications for several large subpopulations of patients with CF. In children, airway infection and inflammation can occur as early as 4 weeks of age (42). Computed tomography radiographic studies have demonstrated considerable bronchiectasis and parenchymal abnormalities in children with normal lung function (43). Sensitive markers could allow for a personalized strategy of antiinfectious and antiinflammatory treatment in young children, with the ability to rapidly monitor outcomes from these interventions. Conversely, in patients with severe lung destruction and multiple antibiotic drug–resistant organisms, assessment of response to a particular treatment is often difficult, given day-to-day variability in disease and the degree of irreversibility in airway damage. A sensitive measure of leukocyte activities could be used to gauge response to therapeutics when clinical response lags far behind. It is important to acknowledge the small size of our study. In addition, two patients provided samples in both cohorts, defining two separate pulmonary exacerbations for each individual. A larger trial is underway by our group to define the optimal combination of genes to characterize a successful therapeutic response. Validation will also require determination of intra- and interindividual variability, responders versus nonresponders, and corticosteroid effect. Herein, we have identified potential gene biomarkers, from which various gene combinations with FEV1 serve as accurate predictors of resolution of pulmonary exacerbations, with greater sensitivity, specificity, and discriminatory capacity than FEV1 alone. Future separate studies will be necessary to determine whether genes have applicability as outcomes measures in clinical trials for antiinflammatory drugs. The utilization of leukocyte gene expression to bolster standard physiologic outcome variables has major implications in terms of the conduct of small trials in rare human diseases.
Supported by National Institutes of Health grants K08HL74512 (M.T.S.), U01 HL081335 (F.J.A.), HL090991 (J.A.N.), the Max and Yetta Karasik Foundation (J.A.N.), the Rebecca Runyon Bryan Chair for Cystic Fibrosis (J.A.N.), and the Cystic Fibrosis Foundation (M.T.S. and J.A.N.). This article has an online supplement, which is accessible from this issue's table of contents at www.atsjournals.org Originally Published in Press as DOI: 10.1164/rccm.200803-387OC on August 21, 2008 Conflict of Interest Statement: M.T.S. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. G.J.H. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. L.A.S. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. M.C. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. D.M.R. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. C.D.C. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. M.W.G. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. S.D.S. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. F.J.A. received $2,500 from Medpro, Inc., for teaching a course. J.W. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. J.A.N. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. Received in original form March 7, 2008; accepted in final form August 20, 2008
This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||