Factors influencing withdrawal of life-supporting treatment in cervical spinal cord injury: a large multicenter observational cohort study

Abstract

Background

Traumatic spinal cord injury (SCI) leads to profound neurologic sequelae, and the provision of life-supporting treatment serves great importance among this patient population. The decision for withdrawal of life-supporting treatment (WLST) in complete traumatic SCI is complex with the lack of guidelines and limited understanding of practice patterns. We aimed to evaluate the individual and contextual factors associated with the decision for WLST and assess between-center differences in practice patterns across North American trauma centers for patients with complete cervical SCI.

Methods

This retrospective multicenter observational cohort study utilized data derived from the American College of Surgeons Trauma Quality Improvement Program database between 2017 and 2020. The study included adult patients (> 16 years) with complete cervical SCI. We constructed a multilevel mixed effect logistic regression model to adjust for patient, injury and hospital factors influencing WLST. Factors associated with WLST were estimated through odds ratios with 95% confidence intervals. Hospital variability was characterized using the median odds ratio. Unexplained residual variability was assessed through the proportional change in variation between models.

Results

We identified 5070 patients with complete cervical SCI treated across 477 hospitals, of which 960 (18.9%) had WLST. Patient-level factors associated with significantly increased likelihood of WLST were advanced age, male sex, white race, prior dementia, low presenting Glasgow Coma Scale score, having a pre-hospital cardiac arrest, SCI level of C3 or above, and concurrent severe injury to the head or thorax. Patient-level factors associated with significantly decreased likelihood of WLST included being racially Black or Asian. There was significant variability across hospitals in the likelihood for WLST while accounting for case-mix, hospital size, and teaching status (MOR 1.51 95% CI 1.22–1.75).

Conclusions

A notable proportion of patients with complete cervical SCI undergo WLST during their in-hospital admission. We have highlighted several factors associated with this decision and identified considerable variability between hospitals. Further work to standardize WLST guidelines may improve equity of care provided to this patient population.

Background

Traumatic spinal cord injury (SCI) has been estimated to affect 1.3 million North Americans [1]. Such an injury can be a devastating event with profound consequences for patients, caregivers, and the healthcare system [2, 3]. In the setting of complete cervical SCI, patients suffer from neurologic sequelae including quadriplegia, loss of sensation and potential loss of respiratory control, which can affect a patient’s ability to move, breathe, and perform activities of daily living. The provision of life-supporting therapy is therefore often a central component of medical care in this patient population. Some patients or caregivers choose to explore the option of withdrawal of life-supporting treatment (WLST). The decision for WLST is a complex process that requires careful consideration of patient values, clinical status, prognosis, support networks, and quality of life. The physician’s role is to help guide decision makers through this challenging terrain by providing objective evaluations of different available treatment pathways. Despite this, there are limited guidelines around how to navigate this complex decision for patients with complete cervical SCI.

To date, there are no multicenter studies describing factors associated with WLST in the SCI population. Current literature on this topic is scarce, with most studies being limited to case reports or series [4,5,6]. An understanding of practice patterns with respect to WLST across health care centers treating patients with this condition is equally lacking. These knowledge gaps represent key areas to be addressed to ensure high-quality care for patients suffering from traumatic SCI. The aim of this study was to assess contemporary trends in the practice of WLST and identify individual and contextual factors associated with the decision for WLST among patients with complete cervical SCI. Secondary aims of this study were to assess variation in practice patterns across North American trauma centers with respect to the decision for WLST.

Methods

Study design and data source

We conducted a retrospective observational cohort study. All data was derived from the American College of Surgeons (ACS) Trauma Quality Improvement Program (TQIP) database from 2017 to 2020 [7]. Nearly 900 ACS- and state-verified level I, II, and III trauma centers across North America contribute to the TQIP database for the purpose of quality improvement. Over 300 variables are collected for every trauma admission each year, including patient demographic and injury characteristics, comorbidities, processes of care, and outcomes. It includes all patients from participating centers with at least one severe injury (Abbreviated Injury Scale (AIS) ≥ 3 in at least one body region). Data reliability and quality are maintained through training of data abstractors and inter-rater reliability audits of contributing centers.

Participants

The study cohort was created and reported in accordance with the Reporting of studies Conducted using Observational Routinely collected health Data (RECORD) statement (Fig. 1) [8]. Adult patients (≥ 16 years) with a diagnosis of an acute traumatic American Spinal Injury Association Grade A cervical SCI were included [9]. This definition represents a complete sensory and motor injury and was identified using AIS codes (Additional file 1). Patients identified with complete cervical SCI and additional non-spine injuries with an AIS score of 6 were also excluded, as these are considered non-survivable injuries [10]. Patients without data on WLST were excluded. Moreover, patients with missing data on demographic, injury, and hospital variables used for regression were also excluded. Additionally, patients with pre-existing advanced directive limiting care are excluded from TQIP.

Fig. 1

Study cohort creation diagram in accordance with the RECORD statement. Abbreviations TQIP, Trauma Quality Improvement Program; SCI, spinal cord injury; AIS, Abbreviated Injury Scale; WLST, withdrawal of life-supporting treatment

Variables

Several patient and hospital variables were selected from the TQIP database according to clinical relevance. Patient demographic data included age, sex, race, and medical insurance type. Age was treated as a continuous variable; sex was dichotomized into male and female; and race was categorized as racially Black, White, Asian, and other. These categories were based on coding available in TQIP, and to ensure a minimum of 5 WLST observations per category as recommended by Harrell [11]. Insurance type was categorized as Medicaid, Medicare, private, self-pay, and other. These categories were chosen based on the data coding in TQIP and to ensure a minimum of 5 WLST observations per category to adequately adjust for each variable level in the regression model. We also collected data on comorbidities likely to influence a decision for WLST, such as whether a patient was functionally dependent prior to injury, or had a history of stroke, dementia, disseminated cancer, or chronic renal failure. Each of these data were available as binary variables in TQIP. Data on the characteristics of the hospital presentation and traumatic injury were also collected. This included the presenting Glasgow Coma Scale (GCS); presence of shock (defined as initial Emergency Department blood pressure < 90 mmHg), whether a patient sustained a pre-hospital cardiac arrest; mechanism of injury; and AIS body region scores including head, face, neck, thorax, abdomen, and the upper and lower extremities. We also captured data on the anatomic level of the SCI, specifically if the injury was high cervical, defined to be C3 and above, or lower cervical defined to be C4 and below. The patient’s GCS was categorized as follows: GCS15, GCS13–14, GCS 9–12, GCS 3–8, consistent with categories corresponding to severity of traumatic brain injury [21]. Mechanism of injury was categorized as blunt or penetrating. This type of categorization has been used in prior studies [12,13,14]. AIS body region scores were dichotomized into 1–2 (minor–moderate injuries) and ≥ 3 (serious–critical injuries). This grouping has also been used in prior studies on SCI [13, 14]. Hospital characteristics including the size, and teaching status were also captured. Hospital size was categorized as ≤ 200 beds, 201–400 beds, 400–600, and > 600 beds. Hospital teaching status was categorized as university hospital, community hospital, and non-teaching hospital. We also were able to link each trauma admission to a unique hospital facility key, to cluster patients according to the hospital at which they were treated.

Outcome

The primary outcome was WLST occurring from the time of hospital arrival to discharge, transfer, or death. TQIP defines WLST to be the decision to either limit escalation, remove, or withhold further life-supporting intervention, limited to ventilator support, dialysis, other forms of renal support, medications to support blood pressure or cardiac function, or a specific surgical, interventional, or radiological procedure. Decisions for WLST were required to be documented in the medical record through either a physician order, progress note, case manager note, social services note, nursing note, nursing flowsheet, or discharge summary. A Do Not Resuscitate status was not a requirement nor seen as equivalent to WLST. Furthermore, TQIP excludes the discontinuation of cardiopulmonary resuscitation from the outcome definition. WLST was coded as a binary variable within TQIP.

Statistical analyses

All statistical analyses were performed using R version 4.2.1 with an a priori specified significance level of p = 0.05 for two-tailed tests. Descriptive statistics were reported as mean, standard deviation, minimum, median, and maximum for continuous variables and count with proportions for categorical variables. Covariate balance across the cohort was assessed by comparing the absolute standardized difference between patients with and without WLST, with a pre-specified threshold of 0.1 representing a significant difference between groups [15].

Multilevel logistic regression

A multilevel mixed effects logistic regression model was constructed with WLST as the primary outcome. We followed multilevel logistic regression best practices outlined by Austin and Merlo [16, 17]. Fixed effect model covariates included age, sex, race, insurance type, functionally dependent status, history of stroke, history of dementia, history of disseminated cancer, history of chronic renal failure, presenting GCS, presence of shock, presence of pre-hospital cardiac arrest, mechanism of injury, head, face, neck, thorax, and abdomen AIS scores, spinal cord level of injury, hospital size, hospital teaching status, and year of injury. Age was modeled as a nonlinear function using restricted cubic splines with 4 knots located at the 5th, 35th, 65th, and 95th quantiles (Additional files 2, 3, 4 and 5).

To investigate differences in case-mix and unmeasured variability between hospitals, we modeled hierarchical clustering of data at the hospital level through random intercepts. We treated hospital specific facility identifiers as a random effect in our multilevel logistic regression model. Effect sizes for fixed effects were reported as odds ratios (OR) with 95% confidence intervals (CI). For age, due to the use of non-linear splines, we reported the OR associated with a change in age from 20 to 30 years, and 60–70 years. Point estimates were given using the predict() function, and 95% confidence intervals were computed using a bootstrap with 100 iterations. We computed p-values with Wald’s test [11] for all fixed effects other than age. For age, due to the use of splines, we reported the p-value associated with a likelihood ratio test when comparing a full model, with a nested model not including age as a covariate [11]. Between hospital practice pattern variation with respect to WLST was quantified through the median odds ratio (MOR). A 95% confidence interval was computed for the MOR using a bootstrap with 100 iterations [16, 17]. The MOR represents the odds of WLST for a patient admitted to a hospital relative to the same patient admitted to another randomly selected hospital. It provides a measure of between-cluster variation in the odds of WLST, providing an estimate of the contextual effect for the decision for WLST. Multilevel model analyses have been well described in a series of prior papers by Merlo and Austin [16,17,18,19,20]. In addition, we computed the proportional change in variance (PCV) attributable to measured patient-level and hospital-level covariates by comparing the difference in individual-level variance between a ‘null model’ including only the hospital random effect, with a ‘full model’ including both the full complement of fixed along with the hospital random effects. We also computed the MOR for the null model for comparison.

A sensitivity analysis was done by refitting the full model after excluding all patients with a presenting GCS of 3-8 (Additional file 6). This was done to further address any confounding related to a concurrent severe traumatic brain injury, as this has been shown to influence the decision for WLST [21].

Missing data were analyzed and found to be missing at random, with less than 5% missing data per variable and less than 10% across all variables (Additional files 7 and 8). We conducted a complete cases analysis and did not impute missing data for the model.

Results

Within the 2017–2020 TQIP database, we identified 5070 patients treated across 477 hospitals that that met our study criteria (Fig. 1). Baseline characteristics of the cohort are summarized in Table 1. Mean age for the cohort was 49.6 years. Most patients were male (80.2%) and racially white (63.7%). The highest proportion of patients had private insurance as their payer (39.5%). Within the cohort, 960 (18.9%) patients had WLST across 329 of 477 hospitals (69.0%). Assessment of covariate balance demonstrated significant differences between patients with and without WLST with respect to age, race, insurance type, functional dependence, dementia, presenting GCS, shock, pre-hospital cardiac arrest, mechanism of injury, SCI level, severe concurrent head injuries, and severe concurrent thoracic injuries (Table 1).

Table 1 Baseline characteristics of study cohort with comparisons based on decision for withdrawal of life-supporting treatment

Patient and hospital factors associated with withdrawal of life-supporting treatment

Results for the multilevel logistic regression model for WLST with respect to fixed effects are summarized in Table 2 and Fig. 2. The likelihood for WLST was found to be significantly associated with age with a non-linear relationship. The odds ratio of WLST associated with a 10-year increase of age, from 60 to 70 years was larger than that associated with an increase from 20 to 30 years. Other patient factors found to significantly increase likelihood of WLST, after adjustment, included being male, having Medicare, self-pay, and a prior history of dementia. Conversely, being Black or Asian significantly decreased the likelihood of WLST.

Table 2 Summary of patient- and hospital-level fixed effect predictors associated with withdrawal of life-supporting treatment

Fig. 2

Forest plot of adjusted odds ratios for patient-level and hospital-level fixed effect predictors of WLST. Adjusted odds ratios were estimated using multilevel logistic regression, with point estimates ranked according to effect size. Predictors with point estimates corresponding to a positive association with WLST are shown in blue, while negative associations are shown in red. Abbreviations: AIS, Abbreviated Injury Scale; WLST, withdrawal of life-supporting treatment; SCI, spinal cord injury; GCS, Glasgow Coma Scale; CI, confidence interval

Injury factors associated with significantly increased likelihood for decision for WLST included presenting to hospital in a comatose state with GCS 3–8 (OR 3.29 95% CI 2.60–4.18), with a pre-hospital cardiac arrest (OR 1.92 95% CI 1.48–2.48), or with a penetrating injury (OR 1.79 95% CI 1.25–2.56). An injury level of C3 and above was associated with twice the odds for WLST compared to an injury of C4 and below (OR 2.00 95% CI 1.65–2.44). With respect to polytrauma patients, having a concurrent severe head injury (OR 1.55 95% CI 1.25–1.91) or thoracic injury (OR 1.29 95% CI 1.03–1.62) was associated with significantly increased likelihood for WLST.

A sensitivity analysis of results after excluding patients with GCS 3–8 was consistent with the primary analysis and found that age, male sex, self-pay, lower GCS, concurrent severe thoracic injury, and a higher SCI level were significantly associated with a higher likelihood of WLST; and being black was significantly associated with a lower likelihood of WLST. Complete results from the sensitivity analysis can be found in Additional file 6.

Between hospital variation in practice patterns of withdrawal of life-supporting treatment

Results of the hospital-level clustering analysis of WLST are provided in Table 3. The average intercept and hospital-level variance of the null model on the log-odds scale were − 1.49 and 0.18, respectively. When computing equivalent values for the full model with adjustment for covariates, the average intercept and hospital-level variation were − 4.07 and 0.17, respectively. The MOR of the full model was 1.48 (95% CI 1.22–1.75), highlighting a significant variation in practice patterns between hospitals. The PCV in the null model after inclusion of patient and hospital covariates was 5.8%, suggesting 94.2% of the variation in hospital-level practices for WLST is unmeasured by the covariates included in the model.

Table 3 Random effects at the hospital level, clustering, and model fit for multilevel logistic regression models

Discussion

In this multicenter retrospective cohort study, we identified factors associated with WLST and investigated variation in practice patterns for WLST among patients with complete cervical SCI. In our study, 18.9% of patients had WLST. Patient demographic factors found to be significantly associated with the decision for WLST included age, sex, race, type of medical insurance, and prior history of dementia. Injury factors found to be significantly associated with WLST included low presenting GCS, presence of a pre-hospital cardiac arrest, penetrating injuries, a C3 and above SCI, and a concurrent severe head or thorax injury. We noted significant variation in practice patterns, with a 95% CI MOR of 1.22–1.75 for the adjusted model, suggesting a significant difference in the odds ratio for WLST in a standard patient from any two hospitals across the study.

Despite its relevance, there have been surprisingly limited prior studies investigating WLST in patients with acute SCI. Osterthun et al. conducted a retrospective study on end-of-life decisions among patients treated with acute SCI across hospital in the Netherlands in 2010. Of the 185 SCI patients treated that year, they noted 19 (10.3%) were reported to have WLST [22]. We noted 18.9% of cases across our study population had WLST. However, our study was focused on complete cervical SCI, which represents a severe form of SCI, and a subset of the heterogenous population of patients with acute SCI. Consideration of WLST is likely more relevant in this patient population as they suffer specifically from quadriplegia and often require respiratory support. This difference likely accounts for the higher rate noted in our study. To our knowledge, all other studies reporting cases of WLST in acute SCI are limited to case reports and case series [4,5,6].

Other studies have assessed factors associated with WLST and practice patterns of WLST in the setting of other traumatic injuries. For patients with severe traumatic brain injury, Williamson et al. noted a WLST rate of 20.7% across the cohort of injured patients [21]. Similar to our study, they noted that race and medical insurance type, and prior history of dementia were significantly associated with decision for WLST in the setting of severe traumatic injury. In their study, however, they did note a significant association between hospital teaching status and the decision for WLST. However, in their study they did not account for hospital-level clustering, which may narrow the confidence intervals around their estimates [16,17,18,19,20]. A prospective cohort study conducted by Cooper et al. looked at variation in practice patterns for withdrawal of care orders for patients with any severe traumatic injury (AIS ≥ 3) across various trauma and non-trauma centers distributed among 12 states [23]. They identified 14,190 patients, of which 618 (4.4%) of patients had a withdrawal of care order. This suggests that patients with complete cervical SCI reflect a subset of the trauma patient population where WLST is more common. Similar to our study, Cooper et al. also noted that age, race, comorbidities, and lower presenting GCS were associated with a withdrawal of care order for trauma patients in general.

In our study we found race and medical insurance type were significantly associated with WLST. Black and Asian patients were less likely to have WLST, and patients with Medicaid or Medicare more likely to have WLST. This highlights potential equity issues within decision-making for WLST among patients with complete cervical SCI. Racial disparities in trauma care have previously been identified. Duong et al. investigated mortality among trauma patients in the USA admitted between 2010 and 2016, and found black patients had an OR of mortality as high as 1.30 when compared to white patients [24]. Furthermore, differences in WLST associated with patient medical insurance type may reflect significant equity issues related to cost of treatment and long-term care. Evidence for disparities in trauma care based on insurance status has been found. A single-center study by de Angelis et al. on trauma patients admitted between 2016 and 2018 found uninsured patients more likely to have a shorter hospital length of stay (incidence rate ratio 0.34, 95% CI 0.24–0.51) [25]. Further research addressing these socioeconomic disparities is crucial to ensure all patients receive equitable, high-quality care.

Limitations of our study are largely imparted by the source of data available. TQIP is a rigorously maintained database, with regular audits and training of data abstractors, however the presence of possible misclassification in data elements is always present in secondary data sources. However, there is no strong reason to suspect differential misclassification among patients with and without WLST. We noted less than 5% missing data among variables relevant to our study; as such, we do not suspect this to severely bias the inference from our study. Furthermore, our analysis of predictors was limited to variables available in the database. We did not have any data on patient religious affiliation, attending physician religious affiliation, and family social dynamics, which would be important for better understanding of decisions around WLST among our patient population [26]. Generalizability of findings is another limitation. TQIP derives data from trauma centers with specific focus on quality of care. As such, we have likely underestimated the variation in practice patterns with respect to WLST that would be found across all hospitals in North America.

A strength of this study is the large sample size of 5070 patients, with 960 observations of WLST across 477 North American institutions over 4 years. Therefore, we me had sufficient observations to regress on the 39 model coefficients included in the model. In addition, our use of a multilevel mixed effect model allowed for the estimation of effect sizes, while adjusting for patient- and hospital-level confounders and data clustering.

Conclusions

An understanding of practice patterns for WLST are relevant to providing quality care to patients with complete cervical SCI. In the current study, we have found that 18.9% of patients with complete cervical SCI ultimately had WLST. Both patient demographic and injury-related factors are associated with this decision, with notable practice variation across North American institutions. More studies are needed investigating the source of variability across hospitals, which is important for hospital policy makers aiming to maintain quality of care. Overall, our study points to a need for standardized WLST guidelines to improve quality of care for patients with complete cervical SCI.