Abstract
Shorter peripheral blood leukocyte (PBL) telomere length (TL) has been associated with poor outcomes in various chronic lung diseases. Whether PBL-TL is associated with survival from critical illness was tested in this study.
We analysed data from a prospective observational cohort study of 937 critically ill patients at Vanderbilt University Medical Center (VUMC). PBL-TL was measured using quantitative PCR of DNA isolated from PBLs. Findings were validated in an independent cohort of 394 critically ill patients with sepsis admitted to the University of California San Francisco (UCSF).
In the VUMC cohort, shorter PBL-TL was associated with worse 90-day survival (adjusted hazard ratio (aHR) 1.3, 95% CI 1.1–1.6 per 1 kb TL decrease; p=0.004); in subgroup analyses, shorter PBL-TL was associated with worse 90-day survival for patients with sepsis (aHR 1.5, 95% CI 1.2–2.0 per 1 kb TL decrease; p=0.001), but not trauma. Although not associated with development of acute respiratory distress syndrome (ARDS), among ARDS subjects, shorter PBL-TL was associated with more severe ARDS (OR 1.7, 95% CI 1.2–2.5 per 1 kb TL decrease; p=0.006). The associations of PBL-TL with survival (adjusted HR 1.6, 95% CI 1.2–2.1 per 1 kb TL decrease; p=0.003) and risk for developing severe ARDS (OR 2.5, 95% CI 1.1–6.3 per 1 kb TL decrease; p=0.044) were validated in the UCSF cohort.
Short PBL-TL is strongly associated with worse survival and more severe ARDS in critically ill patients, especially patients with sepsis. These findings suggest that telomere dysfunction may contribute to outcomes from critical illness.
Abstract
Shorter peripheral blood leukocyte telomere length is strongly associated with worse survival and more severe ARDS in critically ill patients with sepsis http://bit.ly/2mWHKxf
Introduction
Critical illness is commonly associated with a state of catabolic stress, during which the patient exhibits signs of a systemic inflammatory response. This response is associated with complications, such as the acute respiratory distress syndrome (ARDS), systemic infection, multiorgan failure and death [1, 2]. Sepsis and multisystem organ failure remain leading causes of death in critically ill patients, occurring in up to 50% of patients [3–5].
Telomeres are nucleoprotein structures that protect the ends of chromosomes from degradation [6, 7]. Telomere length progressively shortens with cell division, ageing and exposure to various genotoxic stresses such as ionising radiation or reactive oxygen species [8]. When critically shortened, telomeres lose their shelterin complex, a protective protein-capping structure [7, 9]. Uncapping exposes a telomere's blunt end, activating a DNA damage response, which leads to either cellular senescence or cell death via apoptosis [10]. Although each cell has 92 telomeres, uncapping a single telomere may be sufficient to activate a DNA damage response [11].
Telomere length may be measured using a variety of methods including Southern blots, flow-FISH (fluorescent in situ hybridisation) or quantitative (q)PCR [12]. Telomere lengths reported by these methods represent the average telomere length within the population of analysed cells. The shorter the average telomere length, the greater the probability that cells within the population have critically short telomeres. Furthermore, average telomere length within a tissue can be a surrogate marker of the probability of short telomeres in other tissues within the same individual [13]. Considering these relationships, peripheral blood leukocyte (PBL) telomere lengths (TLs) have been analysed in large epidemiological studies in which shorter telomere lengths have been shown to be associated with mental illness, atherosclerotic cardiovascular disease and pulmonary fibrosis [12].
Short PBL-TL is associated with worse survival in patients with idiopathic pulmonary fibrosis (IPF) [14], where telomere dysfunction is increasingly recognised as a molecular driver of the disease [15]. Patients with IPF can suffer from acute exacerbations of their underlying disease. These exacerbations are characterised pathologically by diffuse alveolar damage [16], which is also the pathologic finding in ARDS [17, 18]. In addition, telomere dysfunction isolated to alveolar epithelial cells increases susceptibility of mice to influenza infection, and telomere dysfunction isolated to mesenchymal cells causes pulmonary oedema in mice [19]. Based on the murine studies of lung telomere dysfunction and the histopathological relationship between acute exacerbations of IPF and ARDS, we pursued the hypothesis that PBL-TL may also be associated with survival of critically ill patients, including those with ARDS.
Material and methods
Study design and patient populations
We analysed data from two prospective observational cohort studies of critically ill patients, one from the Vanderbilt University Medical Center (VUMC, Nashville, TN, USA) and one from the University of California San Francisco (UCSF, San Francisco, CA, USA).
The VUMC cohort comprised patients aged ≥18 years admitted to the medical, surgical, trauma or cardiovascular intensive care units (ICUs) from January 2006 to September 2013. All patients who were not discharged from the ICU by ICU day 2 were eligible for enrolment on the morning of ICU day 2, unless they had chronic lung disease, uncomplicated overdose, imminent death or had undergone cardiothoracic surgery. DNA was isolated from PBLs collected the morning after ICU admission. In the UCSF cohort, patients with sepsis admitted to the ICU from the emergency department from September 2009 to January 2017 were eligible, excluding patients with trauma or isolated neurological conditions, with the whole-blood samples collected within 24 h of admission to the ICU.
In both cohorts, the diagnosis of ARDS was made using clinical data including two-physician consensus review of chest radiographs for bilateral opacities, in accordance with the Berlin definition of ARDS [17]. The severity of ARDS was classified based on degree of hypoxaemia: mild (arterial oxygen tension (PaO2)/inspiratory oxygen fraction (FIO2) ratio >200–300 mmHg), moderate (PaO2/FIO2 >100–200 mmHg) and severe (PaO2/FIO2 ≤100 mmHg) [17]. Patients were followed for ARDS for the first four ICU days at Vanderbilt and the first five ICU days at UCSF. For all participants, the presence of sepsis was determined prospectively by a trained study clinician based on the presence of suspected or documented infection plus two or more systemic inflammatory response syndrome criteria. Demographic and clinical information, including APACHE (Acute Physiology and Chronic Health Evaluation) II score and ARDS risk factors were collected. Outcome data included 90-day mortality and 1-year mortality (60-day in-hospital mortality for UCSF sepsis cohort), need for dialysis, duration of mechanical ventilation, duration of ICU stay and duration of hospital stay. Patients with DNA of insufficient quality were excluded from the study.
Ethics approval was obtained from institutional committees at VUMC and at UCSF for use of the source databases and biospecimens. Written informed consent was obtained for all patients with the following exceptions. At VUMC, the institutional review board granted a waiver of consent in cases where the patient was unable to consent and no surrogate could be identified. At UCSF, the institutional review board granted a waiver of consent for patients who died before consent could be obtained [20].
PBL-TL measurements
In both cohorts, PBL genomic DNA was isolated for patients using Gentra Puregene cell kit (Qiagen, Valencia, CA, USA). DNA was visualised on agarose gels to determine quality, and degraded DNA samples were excluded from analysis. Average individual PBL-TLs were measured in triplicate by personnel blinded to all clinical data using quantitative uniplex qPCR and the acidic ribosomal phosphoprotein 36B4 gene as a reference housekeeping gene, as described previously [14, 21–23]. Average telomere lengths were determined by subtracting telomere and reference median cycle threshold (CT) values. Telomere lengths were calculated by comparison to reference samples [24]. Using these methods, the intra-assay coefficient of variation (CV) was <1%, the interassay CV was <3% and the intraclass correlation coefficient for triplicate measurements was 0.981 (95% CI 0.973–0.986).
Statistical analysis
Clinical characteristics of the study cohorts were compared using Kruskal–Wallis rank sum test for continuous variables and Fisher's exact tests for categorical variables. Logistic regression models were utilised to examine the relationship between PBL-TL and hospital mortality controlling for prespecified cofounders reported to influence telomere lengths including age, sex, ethnicity, smoking history and APACHE II. We used ordinal regression models to evaluate the association between PBL-TL and severity of ARDS. Kaplan–Meier plots were constructed to compare patient survival by quartile of PBL-TL, and survival rates in different groups were compared by a log-rank test. The association between PBL-TL and survival time was estimated using Cox proportional hazards regression adjusting for age, sex, ethnicity, smoking history and APACHE II. The proportionality of hazards assumption was checked by Schoenfeld residuals plots against time for each covariate. Cubic splines were applied to permit nonlinear associations and the best regression models were selected based upon the Akaike information criterion to avoid overfitting. Missing values were imputed using transcan function in the R Hmisc package. Variance inflation factor was used to check multicollinearity among variables in the models. All analyses were performed using the R software version 3 (www.r-project.org).
Results
Discovery cohort characteristics
The VUMC cohort consisted of 937 critically ill patients. The mean±sd age was 53±17 years, and 60% were male (table 1). 72% of patients were enrolled within 24 h of hospital admission. Comparison of clinical characteristics between ARDS patients and non-ARDS patients is shown in supplementary table S1. In subgroup analysis, several subcohorts were analysed including a trauma cohort consisting of 281 patients, an ARDS cohort consisting of 249 patients, a sepsis cohort consisting of 368 patients, and a sepsis with ARDS cohort consisting of 151 patients. Some patients were included in more than one subgroup (supplementary figure S1).
Clinical characteristics of critically ill patients in the Vanderbilt University Medical Center cohort
Relationship between PBL-TL and survival of critically ill patients
In the VUMC cohort, the median PBL-TL was 6.96 kb pairs (interquartile range (IQR) 6.47–7.41 kb pairs) (supplementary figure S2); shorter PBL-TLs were independently associated with worse 90-day survival and 1-year survival for all critically ill patients (for 90-day survival: unadjusted hazard ratio (HR) 1.6, 95% CI 1.3–1.8 per 1 kb TL decrease; p<0.0001; adjusted HR 1.3, 95% CI 1.1–1.5 per 1 kb TL decrease; p=0.004; for 1-year survival: unadjusted HR 1.5, 95% CI 1.3–1.8 per 1 kb TL decrease; p<0.0001; adjusted HR 1.3, 95% CI 1.1–1.5 per 1 kb TL decrease; p=0.004) (table 2). Among patients with PBL-TL below the 25th percentile, 34% died within 90 days of their critical illness compared to 15% of patients with PBL-TL above the 75th percentile (figure 1). Shorter PBL-TL was significantly associated with increased 90-day mortality, and 1-year mortality in adjusted models both in the overall VUMC cohort and sepsis subgroup (supplementary table S2). The 90-day mortality rate increased 36% for every 1 kb decrease in PBL-TL in VUMC overall cohort, after adjusting for other covariates (supplementary table S2 and figure 2).
Relationship of peripheral blood leukocytes telomere length to survival in critically ill patients across both Vanderbilt University Medical Center (VUMC) cohort and University of California San Francisco (UCSF) sepsis cohort using Cox proportional hazards models
Shorter telomere lengths are associated with worse survival for critically ill patients in the Vanderbilt University Medical Center (VUMC) cohort, VUMC sepsis subgroup and University of California San Francisco (UCSF) sepsis cohort. a) VUMC cohort (n=937); b) VUMC sepsis subgroup (n=368); c) UCSF sepsis cohort (n=394). Estimated survival for critically ill patients stratified by telomere length (TL) quartiles for survival. Q4: longest quartile; Q3: 50–75%; Q2: 25–50%; Q1: shortest quartile.
Relationship between peripheral blood leukocyte telomere lengths and mortality for critically ill patients in a logistic regression model. a) Vanderbilt University Medical Center (VUMC) cohort; b) VUMC sepsis subgroup; c) University of California San Francisco sepsis cohort. Model controlled for age, male sex, ethnicity, smoking status and Acute Physiology and Chronic Health Evaluation (APACHE) II score.
Relationship between PBL-TL and survival of trauma patients
In the VUMC cohort, critically ill trauma patients had significantly longer PBL-TLs compared to non-trauma patients after adjustment for age (p<0.0001). There was no relationship between PBL-TL and the survival of trauma patients (adjusted HR 0.9, 95% CI 0.8–1.7 per 1 kb TL decrease; p=0.59) (supplementary table S3).
Relationship between PBL-TL and the development of ARDS or organ failure
In the overall VUMC cohort, PBL-TL was not significantly different between ARDS and non-ARDS patients (supplementary table S1), and PBL-TL was not associated with the development of ARDS (p=0.87). In multivariate ordinal logistic regression models, among ARDS patients in the overall VUMC cohort, shorter PBL-TLs were associated with increased risk for developing severe ARDS (table 3). In exploratory univariate analysis, PBL-TL was associated with the number of failed organs (p=0.02) and vasopressor use (p=0.01), but not the need for mechanical ventilation (p=0.78) (supplementary figure S3). These associations did not persist after adjusting for covariates.
Association of acute respiratory distress syndrome (ARDS) severity with peripheral blood leukocyte telomere length in multivariate ordinal logistic regression models
Relationship between PBL-TL and outcomes of sepsis patients
In the VUMC cohort, after adjustment for relevant individual covariates (age, male sex, ethnicity, smoking status and APACHE II score), sepsis patients had significantly shorter PBL-TLs than nonsepsis patients (p=0.0007); shorter PBL-TL was independently associated with worse 90-day survival and 1-year survival of sepsis patients (for 90-day survival: unadjusted HR 1.6, 95% CI 1.3–2.1 per 1 kb TL decrease; p<0.0001; adjusted HR 1.5, 95% CI 1.2–2.0 per 1 kb TL decrease; p=0.001; for 1-year survival: unadjusted HR 1.5, 95% CI 1.2–1.9 per 1 kb TL decrease; p=0.0001; adjusted HR 1.4, 95% CI 1.1–1.8 per 1 kb TL decrease; p=0.003) (table 2 and figure 2). In the VUMC sepsis subgroup, 42% of patients with PBL-TL below the 25th percentile died within 90 days of enrolment compared to 24% of patients with PBL-TL above the 75th percentile (figure 1). Using multivariate ordinal logistic regression models, among ARDS patients in the VUMC sepsis cohort, shorter PBL-TLs were associated with an increased risk of developing severe ARDS (table 3). There was no significant association between PBL-TL and the number of failed organs, vasopressor use or need for mechanical ventilation in the sepsis subgroup (supplementary figure S4).
Validation cohort
Because PBL-TL was most prominently associated with survival in the VUMC sepsis subgroup, and was not associated with outcomes in the trauma cohort, we next sought to validate these findings in a separate sepsis cohort from UCSF. The UCSF sepsis validation cohort consisted of 394 critically ill patients with sepsis; comparison of the UCSF sepsis patients to the VUMC sepsis patients is shown in supplementary table S4. In the UCSF sepsis validation cohort, the median PBL-TL was 6.25 (IQR 5.87–6.58) kb pairs (supplementary figure S2). Shorter PBL-TL was independently associated with worse 60-day survival (unadjusted HR 1.5, 95% CI 1.2–2.0 per 1 kb TL decrease; p=0.001; adjusted HR 1.6, 95% CI 1.2–2.1 per 1 kb TL decrease; p=0.003) (table 2 and figure 2). 39% of patients with PBL-TL below the 25th percentile died within 60 days of their critical illness compared to 20% of patients with PBL-TL above the 75th percentile (figure 1). In addition, shorter PBL-TLs were associated with an increased risk of developing severe ARDS in the UCSF validation cohort (adjusted OR 2.5, 95% CI 1.1–6.3 per 1 kb TL decrease; p=0.044).
Discussion
This study reports that shorter PBL-TL is associated with higher mortality among critically ill patients and that this association is strongest among patients with sepsis. Additionally, among patients with ARDS, shorter PBL-TL was associated with more severe lung injury. These findings suggest that telomere dysfunction may contribute to outcomes from critical illness.
Older adults (age >65 years) are at increased risk of death due to critical illness including sepsis [25, 26]. Molecular explanations for this increased risk remain poorly defined. This study is novel, because it examined the relationship between PBL-TL and outcomes in critically ill patients. PBL-TL is a marker of the molecular age of the patient, with the length being influenced by an individual's chronological age [8], genetic background [27] and lifetime exposure to environmental insults that contribute to telomere attrition. Therefore, these findings suggest that telomere dysfunction in individuals with an advanced molecular age may contribute to critical illness outcomes. The relevant cell types and molecular mediators of this telomere dysfunction that contribute to the poor outcomes will require future study.
The association between short telomeres and increased mortality of critically ill patients was strongest in patients with sepsis. One explanation for the association is that it may be due in part to an impaired capacity of patients with short telomeres to control infections. This hypothesis would be consistent with prior reports that shorter PBL-TL is associated with a higher risk of acquiring infections such as pneumonia [28]. This increased risk of infection may be due to immune cell senescence in leukocytes [29] or T-cell defects such as depleted naive T-cells, increased apoptosis and restricted T-cell repertoire, which have been reported in patients with short telomeres [30]. Alternatively, oxidative stress generated during sepsis could lead to telomere shortening in a progressive positive feedback loop in sepsis, culminating in immune defects and multiple organ failure, which may contribute to poor sepsis survival [31]. Another potential explanation is that the findings are due to differences in the proportion of leukocyte subsets in septic patients. However, when neutrophil count was included as a covariate in the Cox model for survival in a subset of patients from the VUMC cohort with available neutrophil count data, neutrophil count was not found to be a confounder (the p-value for neutrophil count is 0.93).
In the subset of critically ill patients with ARDS, PBL-TL was not associated with the onset of ARDS, but was significantly associated with ARDS severity, with severe ARDS patients having shorter PBL-TLs compared to mild ARDS patients. PBL-TL may be a surrogate marker of an individual patient's risk for short telomeres within cells residing in lung tissue [13]. Telomere dysfunction in alveolar type II cells of mice has been shown to predispose the animals to death from viral pneumonia and mesenchymal cell telomere dysfunction in the lung causes pulmonary oedema [19]. Similarly, patients with IPF, a disease associated with short telomeres in alveolar type II cells, develop acute exacerbations, which are characterised pathologically by diffuse alveolar damage [16]. Putman et al. [32] found that critically ill patients with interstitial lung abnormalities were significantly more likely to develop ARDS. Although the pathogenesis of acute lung injury remains incompletely understood, these findings suggest that telomere dysfunction in alveolar type II cells could increase susceptibility of those cells to external stressors and therefore to the development of lung oedema.
There are limitations to this study. Average PBL-TL was measured on genomic DNA samples isolated from blood leukocytes using quantitative PCR. This method has been reported to have a higher interassay coefficient of variation than other methods of telomere length measurements [33]. Nevertheless, the test characteristics of telomere length measurements in our lab are consistent with test characteristics of other methods for measuring telomeres, such as flow-FISH [34]. Additional factors that may influence PBL-TL measurements include DNA integrity or the population of leukocytes from which DNA is extracted [35]. There was a difference in the average PBL-TL between the VUMC and UCSF sepsis cohorts, probably due to the 11-year mean age difference between the two cohorts. We did not find a significant difference in leukocyte subpopulations between groups of critically ill patients of these two cohorts, suggesting that differences in leukocyte populations do not explain the findings. Furthermore, these differences in age range would decrease the probability that the findings would be replicated, and consequently should increase the generalisability of the findings. The majority of patients were enrolled prior to the current Sepsis-3 guidelines [36]. However, the majority of patients had severe sepsis, which is consistent with Sepsis-3 guidelines. Finally, we were unable to test the correlation between PBL-TL and telomere lengths of cellular subtypes in other organs. It may be useful in future studies to examine patient survival in relationship to telomere lengths within cellular subtypes of other organs.
In conclusion, shorter PBL-TL is associated with mortality and ARDS severity in two independent cohorts of critically ill patients, suggesting that telomere dysfunction is associated with the pathobiology of critical illness. Measures of PBL-TL may provide a unique opportunity to advance understanding of the interaction between ageing and the pathobiology of ARDS and sepsis.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material ERJ-01044-2019.SUPPLEMENT
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-01044-2019.Shareable
Acknowledgements
We are grateful to all the patients and coordinators who participated in this study at all sites.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Author contributions: P.J. Wolters, C.S. Calfee, L.B. Ware and M.A. Matthay designed the study. S. Liu and G. Green performed telomere length measurements. C. Wang and S. Liu analysed the data. H. Zhuo, L.B. Ware, C.S. Calfee, M.A. Matthay, K.D. Liu, K.N. Kangelaris, A. Gomez, A. Jauregui, K. Vessel, S. Ke and C. Hendrickson collected data for the study. S. Liu drafted the manuscript. All authors contributed to study conduct, data interpretation and manuscript preparation.
Conflict of interest: S. Liu reports grants from National Natural Science Foundation of China, during the conduct of the study.
Conflict of interest: C. Wang has nothing to disclose.
Conflict of interest: G. Green has nothing to disclose.
Conflict of interest: H. Zhuo has nothing to disclose.
Conflict of interest: K.D. Liu has nothing to disclose.
Conflict of interest: K.N. Kangelaris has nothing to disclose.
Conflict of interest: A. Gomez has nothing to disclose.
Conflict of interest: A. Jauregui has nothing to disclose.
Conflict of interest: K. Vessel has nothing to disclose.
Conflict of interest: S. Ke has nothing to disclose.
Conflict of interest: C. Hendrickson has nothing to disclose.
Conflict of interest: M.A. Matthay reports grants from NIH/NHLBI (for research and clinical trials of ARDS and sepsis), Department of Defense (to carry out a clinical trial of ARDS) and Bayer Pharmaceuticals (for an observational research study of ARDS), and personal fees from Cerus Therapeutics (ARDS consultation) and CSL Behring (ARDS consultation), outside the submitted work.
Conflict of interest: C.S. Calfee reports grants from NIH, during the conduct of the study; grants from GlaxoSmithKline (for an observational study of sepsis and ARDS), grants and personal fees from Bayer (for an observational study of ARDS and consultancy about ARDS), and personal fees from Prometic (consultancy), CSL Behring (medical advisory board), Roche/Genentech (consultancy) and Quark (consultancy), outside the submitted work.
Conflict of interest: L.B. Ware reports grants from National Institutes of Health, during the conduct of the study; personal fees from Quark Pharmaceuticals (advisory board), CSL Behring (advisory board) and Bayer (advisory board), outside the submitted work.
Conflict of interest: P.J. Wolters reports grants from MedImmune, Genentech and Boehringer Ingelheim, and personal fees from Roche, Boehringer Ingelheim, Blade Therapeutics and Pliant, outside the submitted work.
Support statement: The study was funded in part by NIH grants HL103836 and HL135849 (L.B. Ware), HL131621 and HL140026 (C.S. Calfee), HL51856 (M.A. Matthay), National Natural Science Foundation of China 81700063 (S. Liu), and the Nina Ireland Program for Lung Health. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received May 25, 2019.
- Accepted September 21, 2019.
- Copyright ©ERS 2020