|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: TGL FRA JM RP JRH. Performed the experiments: TGL YL JLJ WW KP TC HC JC SLR CL AS. Analyzed the data: TGL FRA RHK WW JM RP DRC JRH. Contributed reagents/materials/analysis tools: TGL FRA RHK WW DRC. Wrote the paper: TGL FRA RHK WW JRH DRC.
Lung function, acute pulmonary exacerbations (APE), and weight are the best clinical predictors of survival in cystic fibrosis (CF); however, underlying mechanisms are incompletely understood. Biomarkers of current disease state predictive of future outcomes might identify mechanisms and provide treatment targets, trial endpoints and objective clinical monitoring tools. Such CF-specific biomarkers have previously been elusive. Using observational and validation cohorts comprising 97 non-transplanted consecutively-recruited adult CF patients at the Intermountain Adult CF Center, University of Utah, we identified biomarkers informative of current disease and predictive of future clinical outcomes. Patients represented the majority of sputum producers. They were recruited March 2004-April 2007 and followed through May 2011. Sputum biomarker concentrations were measured and clinical outcomes meticulously recorded for a median 5.9 (interquartile range 5.0 to 6.6) years to study associations between biomarkers and future APE and time-to-lung transplantation or death. After multivariate modeling, only high mobility group box-1 protein (HMGB-1, mean=5.84 [log ng/ml], standard deviation [SD]=1.75) predicted time-to-first APE (hazard ratio [HR] per log-unit HMGB-1=1.56, p-value=0.005), number of future APE within 5 years (0.338 APE per log-unit HMGB-1, p<0.001 by quasi-Poisson regression) and time-to-lung transplantation or death (HR=1.59, p=0.02). At APE onset, sputum granulocyte macrophage colony stimulating factor (GM-CSF, mean 4.8 [log pg/ml], SD=1.26) was significantly associated with APE-associated declines in lung function (−10.8 FEV1% points per log-unit GM-CSF, p<0.001 by linear regression). Evaluation of validation cohorts produced similar results that passed tests of mutual consistency. In CF sputum, high HMGB-1 predicts incidence and recurrence of APE and survival, plausibly because it mediates long-term airway inflammation. High APE-associated GM-CSF identifies patients with large acute declines in FEV1%, possibly providing a laboratory-based objective decision-support tool for determination of an APE diagnosis. These biomarkers are potential CF reporting tools and treatment targets for slowing long-term progression and reducing short-term severity.
Cystic fibrosis (CF) is the most common lethal genetic disorder among Caucasians, but disease occurs world-wide. Approximately 10 million Americans carry mutations while 25,000 suffer actual disease . Pancreatic enzyme therapy increased survival from infancy to 6 years of age and shifted the commonest cause of death from starvation to lung disease, which remains life-limiting for most patients , . Median predicted survival now approaches 40 years  but remains well short of normal.
Therapeutic successes have greatly increased CF survival, however, this makes survivorship problematic for assessing treatment effects, prompting a search for alternative measurements –. Lung function, acute pulmonary exacerbations (APE) of CF requiring hospitalization and weight best predict survivorship clinically  but may themselves require large trial enrollments and years of observation to discern treatment effects.
In CF, airway inflammation underlies pulmonary deterioration, frequent APE and persistent malnutrition , . Chronic airway infections incite intense, dysfunctional inflammation , , eliciting multiple deleterious signals . However, these inflammatory signals may serve as immediate biomarkers of disease state. A biomarker that additionally predicts clinical outcomes might identify causal mechanisms for airway disease, pinpoint investigational targets, rapidly detect efficacy of new therapies, and distinguish patients most needing urgent interventions , . Such biomarkers have previously been elusive in CF.
We hypothesized that high-mobility group box-1 protein (HMGB-1), a novel inflammatory cytokine, might be informative of future clinical events in CF. HMGB-1 is a highly conserved DNA-shepherding protein. However, when released by necrosis or secretion, HMGB-1 acts as a potent cytokine associated with delayed onset of prolonged inflammation ,  and septic shock in animal models reversible with anti-HMGB-1 antibodies . In this work, we investigated the potential of sputum HMGB-1 measurements to reflect concurrent clinical status and predict future outcomes following a physician-patient encounter .
This study was reviewed and approved by the University of Utah Investigational Review Board prior to enrollment of any patients. All participants were included only after obtaining written informed consent for sputum collection and analysis. Collection of information for the CF Foundation Patient Registry (CFFPR) was reviewed and approved at each participating institution and by the University of Utah Investigational Review Board. Our use of data from the CFFPR was also reviewed and approved by the CF Foundation. All adult patients included in the CFFPR gave written informed consent. All pediatric patients in the CFFPR gave assent if able to understand, and their parents or guardians gave written informed consent and permission to participate. The study was conducted under the supervision of the physicians involved (TGL, WW, HC, RP III, JRH), and patients were advised of risks, benefits and the right to withdraw from further involvement in the study at any point without repercussions. All data, particularly patient identifying data, were physically and electronically secured throughout the study.
From March, 2004 through May, 2011, we recruited sputum-producing outpatients from the Intermountain Adult CF Center at quarterly-regular and unplanned-sick visits. All visits were conducted in compliance with CF Foundation guidelines for outpatient care . We collected sputum, recorded same-day lung function, weight, height, sputum culture results, number of prior-year APE, and recorded whether patients were in stable, mild exacerbation or APE states. We actively followed patients to lung transplantation, death, loss-to-follow-up, or May 31, 2011, the end of study, recording APE numbers and spirometry results . Patients recruited from March, 2004 through September 2006 formed the study groups while remaining patients formed validation groups (Tables 1 and and2).2). Forced expiratory volume in 1 second (FEV1) was measured  and normalized to percent predicted FEV1 (FEV1%) . We assessed generalizability of results by comparing our patients with sputum-producing adults in the 2006 CFFPR.
We diagnosed APE and hospitalized for symptoms and objective evidence of severe acute worsening of CF. Symptoms included increased cough, sputum production or dyspnea, chest pain or tightness, hemoptysis, fever, chills, arthralgias or decreased exercise tolerance. Patients had to meet one objective criterion: 10% decrease in FEV1% or percent predicted forced vital capacity, fever above 38.4°C, documented hemoptysis greater than 100 ml per episode, SaO2 below 90% or PaO2 below 60 mm Hg despite usual oxygen, increased supplemental oxygen requirements or three or more kg unplanned drop in weight within 3 months. Symptomatic patients without objective findings were classified as mild exacerbation and not hospitalized. These criteria were modified from previous APE definitions to allow prospective application at a patient encounter by excluding retrospective criteria such as prior antibiotic treatment , , .
Sputum samples underwent standardized processing (see Text S1). CF Foundation guidelines were followed to collect sputum cultures. Cultures were obtained within 6 months of sputum collections as bacterial infections are typically stable over such a period , a finding that has recently been supported by culture-independent methods . We measured potential biomarkers previously identified as important in CF including granulocyte macrophage colony stimulating factor (GM-CSF) , , – using SearchLight multiplex assay services (Aushon Biosystems, Billerica, MA). We measured HMGB-1 by ELISA using commercially-available antibodies (R&D Systems, Inc, Minneapolis, MN) and previously published protocols  after confirming assay reproducibility on consecutive days (see Text S1). Biomarker measurements reported are log-scale.
Through 5 analyses , we examined (1) associations between biomarkers and concurrent clinical disease measurements, (2) APE effects on biomarkers, (3) biomarker predictions of clinical outcomes using linear, logistic and quasi-Poisson regression , and proportional hazards modeling ,  (See Tables 1 and and22 and Text S1 for details). Throughout this paper, p-values smaller than 0.05 are reported as statistically significant. In Analyses 1–3, we performed multiple univariate tests for associations between biomarker values and various clinical measurements, and we used stringent Bonferroni correction by dividing 0.05 by the number of tests performed to define the corrected p-value required for significance. We repeated Analyses 3–5 using sputum measurements from 4 sets of validation patients chosen so that no patient appeared simultaneously in a study group and comparison validation groups (Table 1). Mutual consistency of results was assessed by a weighted least squares analysis.
During enrollment of the three study groups, 161 patients attended clinic at least once; 26 produced no sputum. Including validation groups, we recruited 97 patients and collected 149 samples: 56 samples were collected during stable, 31 during mild exacerbation and 62 during APE clinical states (Table 3 and Text S1). Every sputum-producing patient approached agreed to participate, but we were occasionally limited by availability of study personnel or by temporary lack of sputum production. All patients were lifetime non-smokers. Study patients were nearly indistinguishable from adult sputum producers in the 2006 CF Foundation Patient Registry (Table 4). Our patients had a weight-for-age z-score 0.24 units less than the national median and used oral steroids less frequently.
Study group patients were followed up to 7.1 years (median 5.9, interquartile range 5.0 to 6.6) through death, lung transplantation, loss to follow-up or end of study. All APE hospitalizations occurred within our center. Nine deaths occurred between 0.6 and 6.5 years, and six lung transplants between 1 and 6.3 years after enrollment. Three patients were lost to follow-up 2.7, 3.4 and 3.6 years after enrollment when they moved.
Biomarkers were measured (Table S1) in three overlapping patient study groups after sample inventory demonstrated feasibility of planned analyses (Table 1). Group 1: the initial sputum sample from each of 80 randomly selected patients underwent SearchLight multiplex assays (Aushon); results for 24 were discarded due to unreliability (as samples thawed in transport to Aushon), leaving results from 56 patients for analysis. Group 2: we measured HMGB-1 by ELISA and obtained a second SearchLight multiplex assay (Aushon) on samples from the 26 patients that during repeated collections resulted in a pair of samples per patient: one from an APE and one from a stable point adjacent in time. Ten gave a clinically-stable-time-point sample first; 16 gave an APE-time-point sample first. The median number of days from stable to APE sample was 98 (interquartile range 91–154) and from APE to stable sample was 119 (interquartile range 84–212). Group 3: 76 patients from Group 1 had sufficient sample volumes to measure HMGB-1.
We looked for associations with concurrent clinical findings to identify plausible biomarkers for CF airway inflammation (Table S2). Backwards model selection using biomarkers with statistically significant or borderline univariate associations (p<0.1) produced three final multivariate models. HMGB-1 and interleukin-(IL)-17A are associated with concurrent FEV1% (Figure 1A and S1A), IL-17A with weight-for-age z-score (Figure S1B), and HMGB-1 and interferon-(IFN)-α with number of prior-year APE (Figure 1B, Figure S1C, Table 5). (Discussion of these results can be found in Text S1.).
We looked for biomarkers that changed in association with a clinical change from stable to APE state. For 19 of 22 biomarkers measured for Study Group 2 patients, the mean change between stable and APE states (Table S1) was positive indicating a non-random inflammatory state change (sign test, p=0.006). C-reactive protein (CRP), IL-1β and IFN-α were statistically significantly higher (paired t-test, p<0.05) during an APE (Figure S1D-F). We investigated a potential effect of sample collection sequences (stable-then-APE and APE-then-stable) in sensitivity analyses. Mean IFN-α change was greater for the sequence stable-then-APE (t-test, p=0.015). No sample-sequence-associated differences were found for CRP or IL-1β. (Discussion of these results can be found in Text S1.).
For the 26 Study Group 2 patients, mean FEV1% decline between stable and APE time points was 7.3 (SD 7.8) percentage points; 5 had no decline. Univariate linear regressions found statistically significant inverse associations between FEV1% decline and GM-CSF, IL-23, CRP and sputum total cell count (TCC) each measured at APE time points and the difference in IL-17 measurements at stable and APE time points (Table S3). Neither biomarker measurements at stable time points nor any other biomarker differences between stable and APE time points were associated with FEV1% decline. After backward selection of multivariate models, GM-CSF and to a lesser extent IL-5, both measured at APE time points, were significantly and independently associated with APE-associated decline in FEV1% (Figure 1C–D, Table 5).
Increased HMGB-1, adjusted for number of prior-year APE and follow-up time, was retained after backwards and forwards selection in a final multivariate quasi-Poisson regression model for the total number of APE occurring within 5 years of clinically-stable-time-point samples (Table 5). All other clinical measurements including P aeruginosa and S aureus infections, FEV1%, 5-year predicted survival  and other potential biomarkers measured at the stable-time point were not statistically significant predictors when forced into the final multivariate model (p=0.23 for FEV1%; p>0.05 for all other covariates).
We have complete follow-up from clinically-stable time points to first APE, enabling proportional hazards modeling of time-to-first APE –. Univariate and multivariate models found HMGB-1 alone predicts time-to-first APE (Table 6, Figure 2A). Because of this result, we explored using median HMGB-1 level as a clinical test for predicting an APE within 6 months and found favorable preliminary results (see Text S1 and Figure S2).
Because of multiple other statistically significant associations, we assessed HMGB-1 ability to predict time-to-lung transplantation or death –. The composite outcome avoids 2 biases: 1) in a study of time-to-death, transplantation of end-stage patients creates an informative censoring bias and 2) in a study of time-to-lung transplantation, use of FEV1% to select transplantation candidates confounds the analysis. However, the composite outcome would still be confounded by modeling factors used to select transplantation candidates , , , thus potential covariates including FEV1% and prior-year APE were necessarily excluded from analysis. Other potential covariates including P aeruginosa and S aureus infection status do not alter lung transplantation candidacy and were included in the analysis. Follow-up time for the composite event was time-to-listing for transplantation, death, loss-to-follow-up or study end. Three patients lost-to-follow-up were treated as alive at the time of censoring. After multivariate proportional hazards modeling with backward and forward selection, we found that higher HMGB-1 values alone are strongly associated with shorter time-to-transplantation or death (Table 6, Figure 2B).
Analyses 3–5 were repeated using data from four validation patient groups (Tables 1 and and2)2) to validate predictive abilities of HMGB-1 and GM-CSF. All validation results were similar to primary results. We evaluated results for mutual consistency using a weighted least-squares analysis and found that differences in coefficients and hazard ratios were insufficiently discrepant to cause concern (Table 7).
Sputum HMGB-1 measurements predicted lung transplantation or death for CF patients during 7.1 years of follow-up (Table 6, Figure 2B). HMGB-1 measured during clinical stability surpassed FEV1% and all other biomarkers alone or in combination in ability to predict number of subsequent APE (Table 5) and time-to-first APE (Table 6, Figure 2A). Simultaneously, HMGB-1 measurements possess immediate clinical relevance by demonstrating strong independent associations with concurrent FEV1% and number of prior-year APE (Table 5) .
Our results confirm published findings from 32 CF patients that HMGB-1 inversely correlates with FEV1% and has an association with patients that suffer an APE. The earlier statistical findings were limited to concurrent clinical status . Our study substantially extends these results finding evidence that HMGB-1 predicts future clinical events including survival (Tables 5 and and66).
HMGB-1 is required for normal cellular function but when released due to necrosis or secretion ,  plays key roles in inflammatory arthritis , sepsis , and acute lung injury . It mediates endotoxin-associated inflammation  through RAGE that is upregulated specifically in CF airways  and other receptors . In rat models, HMGB-1-blocking antibodies reversed arthritis  and produced remarkably improved, durable survival despite delayed intervention following induction of septic shock by cecal ligation and puncture . Increased HMGB-1 is associated with human acute lung injury  and sustained COPD-related inflammation . In CF, HMGB-1 may increase pulmonary inflammation by attracting neutrophils , preventing their efferocytosis  and amplifying the effects of bacterial lipopolysaccharide and cytosine-phosphatidyl-guanosine-DNA constructs . Soluble RAGE, an antagonist to HMGB-1-RAGE mediated inflammation, is absent in CF airways, potentiating HMGB-1 effects . Our successful predictions (Tables 5 and and6,6, Figure 2) strengthen the argument that HMGB-1 may play a causal role in human inflammatory lung disease.
Low lung function shortens CF survival . Sharp FEV1% declines define APE for most patients and may presage permanent lung function loss . In this study, several biomarkers had significant univariate associations with acute APE-associated FEV1% decline, but GM-CSF and to a lesser extent IL-5 measured at APE time points together had the strongest association (Figure 1C–D, Table 5). These associations were independent of FEV1% itself and either cytokine measured during clinical stability. Our findings reflect previously known relationships in structure, expression and function of GM-CSF and IL-5, including the strong correlation between their measurements (Table S4) and the observation that one cytokine, in this case GM-CSF, is often more prominent than the other . For GM-CSF at an APE time point, an approximate 1% increase was associated with a 10 point decline from the clinically-stable FEV1%. A patient’s GM-CSF measurement may quantify APE severity independently of prior clinical state and might provide laboratory-based objective support for hospitalization decisions.
GM-CSF maintains normal alveolar macrophage and innate immune responses to acute P aeruginosa pneumonia in mice . In humans, elevated sputum GM-CSF is associated with increased COPD and asthma severity . However, its role in chronically infected human airways is unclear. P aeruginosa or S aureus in CF airways  elicit increased epithelial GM-CSF secretion leading to prolonged survival and decreased apoptosis of airway neutrophils  and may prolong protease and reactive oxygen species releases .
Our definition of APE, based on experience and bedside utility, played a central role in defining events, sample collections and analysis of our major results with HMGB-1 and GM-CSF. It resembles prior definitions , , , but differs by excluding retrospective criteria such as a prior decision to give antibiotics  and by exchanging quantitative scoring  for bedside facility. Nevertheless our definition includes the key covariates of the existing quantitative model .
Our study focused on adult CF patients with moderately severe disease (Table 3), thus it provides no information easily extrapolated to children or mild or end-stage adult patients. However, we enrolled more than half of our large adult CF center and two-thirds of sputum producers. Our patients were nearly indistinguishable from the 70% of 2006 CFFPR adult patients nationally that produced sputum, a large, identifiable group that would benefit from new insights and therapies (Table 4) . Volume of sputum collection could not be compared to the national population but had no effect in any of our analyses (not shown).
HMGB-1 was the key biomarker identified in multivariate models, proving to be superior to IL-8 and other biomarkers of previously established interest, indicating that HMGB-1 measurement alone is sufficient to report the information encoded in the other biomarkers evaluated (Table S1). However, we could not study all potentially significant biomarkers. To avoid compromising sample collection (see Text S1) , we excluded proteinases such as neutrophil elastase , , matrix metalloproteinase ,  and proteinase-3 ,  which may all cause lung tissue injury. Thus further studies are needed to explore the relative clinical predictive ability of these and still other biomarkers.
Our HMGB-1 and GM-CSF results successfully validated (Table 7) with additional but greatly dissimilar patients (Table 3). In every case, results were unaffected by commonly-prescribed potential modifiers of inflammation such as azithromycin and inhaled steroid therapies or by airway infections with Pseudomonas aeruginosa or Staphylococcus aureus. Nevertheless, this was still a single center study that was exploratory in nature, and collection of additional high-quality data from a larger cohort is needed to confirm our findings and explore directions of causality .
It is not feasible to include FEV1% or several other clinical covariates with HMGB-1 to predict time-to-lung transplantation or death due to biases from confounding and informative censoring. These biases arise because the same clinical covariates are used for selection of lung transplantation candidates. Prediction of either transplant or death alone is equally not feasible due to the same biases. Using a combined endpoint of lung transplant or death and excluding confounding variables avoided these biases but assumed that our clinical decisions  were accurate in identifying patients at highest risk of death . These limitations are not specific to our study size, population or design. Rather, they apply to any population where lung transplantation is utilized. A positive byproduct results from avoiding these biases: this study focuses on causal pathways for airway inflammation leading to end-stage lung disease and avoids confusion from consideration of intermediate clinical findings like FEV1% that may not be causal .
This study identified airway-inflammatory biomarkers reporting on distinct pathophysiologic states that precede, accompany and follow an APE and provided clinically relevant information about permanent lung function losses on short and early mortality on long time scales. GM-CSF and HMGB-1 provide potentially useful new measurements for monitoring short and longer term treatment effects for CF. Both biomarkers provide focal points for additional clinical, epidemiologic and mechanistic investigations of CF airway inflammation. In particular, HMGB-1 provides pathophysiologic information predictive of APE and survival that might guide preemptive treatment and provide a novel target for therapy.
Biomarker Relationships with Concurrent Clinical Outcomes and Biomarker Changes between Stable and APE States. Relationships between A) IL-17A and concurrent FEV1%, B) IL-17A and weight-for-age z-score and C) IFN-α with prior year APE. Three biomarkers, D) CRP, E) IL-1β and F) IFN-α are increased in sputum collected at the time of diagnosis of an APE compared to their measurements during clinical stability using paired t-tests. Red symbols illustrate values for patients that underwent the sputum collection sequence stable-then-APE while blue symbols denote patients that underwent the sequence APE-then-stable. The number of days between stable and APE sputum collections had no effect on differences on measurements (multiple linear regression p>0.5). However, changes in these three biomarker measurements are unsuitable for clinical use to identify an APE due to substantial overlap of APE and stable time point values.
Receiver operator characteristic curves for potential markers of an APE within 6 months. Receiver operator characteristic curves were constructed for 9 biological or clinical markers of CF using APE within 6 months as the clinical discriminator. Log-scale biomarker values were derived from measurements of sputum samples from 26 patients during clinical stability (Study Group 2, Tables 1–2). Each plot has a color bar presenting the range of cut-off values of each biomarker used to create the points on each ROC curve. The area under the curve (AUC or accuracy) and 95% confidence interval for each AUC reported for each curve was the result of boot strapping of the AUC using 1,000 repetitions. HMGB-1 has the highest accuracy of any marker studied.
Sputum Biomarker Measurements.
Univariate Associations of Biomarkers with Clinically Relevant Concurrent Conditions.
Univariate Associations of Biomarkers for APE-associated FEV1% drop.
Correlations (p-values) between biomarkers, Study Group 1.
Detailed Laboratory and Statistical Methods, Additional Results and Discussion.
The authors thank Barbara Cahill, MD, Greg Elliott, MD and Wayne Samuelson, MD, all faculty at the University of Utah, who reviewed early versions of the manuscript and provided valuable critiques without compensation. We thank Hebe Quinton, MS, Dartmouth College, for technical assistance obtaining specific data about the Intermountain Adult CF Center from the US CFFPR for 2006 without compensation.
The national data on CF patients reported were supplied by the U.S. Cystic Fibrosis Foundation. The authors are responsible for interpretation and reporting of the CFFPR data and do not express opinions or recommendations for the Foundation.
This work was supported by grants from the NIH/NHLBI (HL03656 and HL084341), a Harry Shwachman Clinical Investigator Award, a Clinical Research Award from the CF Foundation and by generous support from the Claudia Ruth Goodrich Stevens Research Endowment in Cystic Fibrosis in the Department of Internal Medicine, University of Utah and the Ben B. and Iris M. Margolis Family Foundation of Utah. The funding entities had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.