FormalPara Take home message

The early use of dexmedetomidine for sedation of ventilated critically ill patients who are older than 65 years, and in those with an operative diagnosis, across broad range of age categories, has a high probability of reduced mortality. Conversely, younger patients with a non-operative diagnosis have a high probability of increased mortality. Thus, the early use of dexmedetomidine in this group of patients, outside controlled research, is not advised.

 

Introduction

Clinical practice guidelines have recommended dexmedetomidine (DEX), a high-affinity alpha-2 adrenergic agonist, for sedation of critically ill ventilated patients in preference to benzodiazepines [1, 2]. The Sedation Practice in Intensive Care Evaluation (SPICE) III trial, by far the largest study in the field, compared the early use of DEX as primary sedative agent to usual care sedation in mechanically ventilated critically ill patients [3]. While the trial reported no overall difference in mortality, a pre-specified subgroup analysis according to age, dichotomized at the median of 63.7 years, observed a divergent DEX effect on mortality with a higher mortality below and a lower mortality above such median age.

There was no apparent explanation for the observed age effect on mortality. Thus, despite a large sample size, and a wide inclusion criteria in the SPICE III trial, heterogeneity of treatment effect (HTE) on mortality could not be ruled out [4]. It is also possible that HTE for mortality in response to DEX treatment may apply to other groups with different baseline characteristics beyond age.

Accordingly, we conducted a secondary Bayesian analysis of the SPICE III trial population to quantify the HTE of early sedation with DEX compared with usual care. Our aim was to identify patients who have a high probability of lower or higher 90-day mortality, according to age categories, and other specific clinical subgroups (clusters) according to clinically relevant baseline characteristics.

Methods

Study design, population, and intervention

This was a secondary analysis of the SPICE III trial. The SPICE III study design and protocol have been previously published [3, 5]. The study was an open label, multinational, randomized clinical trial conducted in 74 intensive care units (ICU) in Australia, Ireland, Italy, Malaysia, New Zealand, Saudi Arabia, Switzerland, and the United Kingdom between November 2013 and February 2018.

The SPICE III trial enrolled mechanically ventilated patients who were expected remain ventilated after the next calendar day and needed ongoing sedation. Eligible patients were randomly assigned to receive dexmedetomidine infusion started at 1 µg/kg/h, without a loading dose and titrated to effect, or usual care sedatives as determined by the treating clinician. In the present analysis, only patients with available data for the outcomes and for the characteristics considered for group formation were included.

Definitions and subgroups

Outcomes were defined as per the primary manuscript [3]. HTE was assessed according to the following subgroups:

  • Age groups: patients were divided in two groups according to World Health Organization conventional definitions of older age > 65 years.

  • Clusters according to pre-defined baseline characteristics described below.

Clustering process

Clustering is an approach, which aims to separate patients into groups with similar traits (called “clusters”) [6, 7]. This is a non-supervised approach, meaning that no information of predefined classes is used. The number of clusters and its characteristics are completely determined by the computational algorithm. For cluster detection, based on clinical relevance, the following variables were selected a priori on the basis of clinical judgment: (1) Acute Physiology and Chronic Health Evaluation (APACHE) II score, without the age component [8]. (2) Baseline partial pressure of arterial oxygen divided by percentage of inspired oxygen (PaO2/FiO2 ratio). (3) Admission diagnosis (defined in mutually exclusive categories as sepsis, respiratory, gastrointestinal, cardiovascular, trauma, neurological or other). (4) Type of admission (defined as non-operative, elective surgery or emergency surgery). (5) ICU source of admission (defined in mutually exclusive categories as emergency department, ward, transferred from other services or operating room). (6) Gender (male or female). (7) Baseline use of opioids (morphine or fentanyl, yes/no). (8) Baseline use of propofol (yes/no). (9) Baseline use of midazolam (yes/no). (10) Baseline use of ketamine (yes/no), and (11) Baseline use of DEX. Sensitivity analysis excluding missing baseline sedatives and imputation of PaO2/FiO2 by the median showed consistent clusters, Supplementary Information [e-Table 2 and eFigure 2B].

As age is a contributing component to the APACHE II score, to create a marker of patient severity that is independent of age, the age component was removed from the APACHE II score.

The K-means for mixed large data (kamila) method was used to detect the clusters [6]. The best number of clusters was defined by inspecting the prediction strength of clusters after 1000 cross-validations. More information on the clustering process is provided in the Supplementary Information eMethods section.

Outcomes

The study primary outcome was 90-day all-cause mortality. Additional secondary outcomes were the number of days alive, free from coma or delirium and the number of ventilator-free days at day 28 after randomization. We included secondary outcomes that showed a difference between DEX and usual care in the primary intention to treat analysis.

Statistical analysis

Baseline and outcome data are presented according to the subgroups assessed. Continuous data are presented as median (quartile 25%–quartile 75%) and compared with the Wilcoxon rank-sum test; while, categorical binary data are presented as number and percentage and compared with Fisher’s exact test.

HTE was assessed independently in each of the two subgroups proposed following a Bayesian hierarchical logistic model for the primary outcome and a Bayesian hierarchical linear model for the secondary outcomes [9]. All hierarchical models were modeled as a simple regression and shrinkage model. The hierarchical models partially pool the data and shrink the estimates in each subgroup towards the overall estimate, with shrinkage proportional to the size of the subgroup limiting the risk of type 1 error. For all analyses, weakly informative priors were used, aiming to encompass all plausible effect sizes. Since the sample size of the original study is large, it is expected that the likelihood will dominate the posteriors. The models used and priors’ definitions are described in detail in the Supplementary Information statistical analysis section. The effect estimate was reported as odds ratio (OR) where an OR > 1 represents harm and an OR < 1 represents benefit with early DEX use. As the recorded mortality according to age categories was different, for clarity, a 10% risk difference represents a 4% and 2% change in mortality in older and in younger patients, respectively.

For the secondary outcome, the mean difference (MD) and 95% credible interval (CrI) of the distribution are presented, as the probability of a MD < 0.00 or > 0.00 (null) and < 2 or > 2 (2-day reduction or increase).

A mixed-effect Bayesian logistic and linear regression model with centers as random effect and considering weakly informative priors was used to assess the interaction between the treatment groups and age on the continuous scale. The results are shown in marginal effect plots. Additional analyses considering pessimistic and optimistic priors were conducted as sensitivity analyses (full description in Supplementary Information).

To further understand the interaction according to age and severity of disease on HTE for DEX, we assessed the within-age association between use of DEX and 90-day mortality in a mixed-effect Bayesian logistic regression model according to APACHE II (without the age component). In this model, interactions between APACHE II groups (stratified into six groups) and allocation groups, age subgroups (≤ 65 vs. > 65 years) and allocation groups, and age subgroups and APACHE II quintiles were included. Also, to assess the within-cluster variation for the effect of DEX according to age, a new model including interactions between age groups (stratified into six groups) and allocation groups, clusters and allocation groups, and cluster and age quintiles were included. The models considered a Bernoulli distribution, with centers as random effect and with starting values randomly generated. All priors were drawn from normal distributions and were weakly informative (full definition of priors available in the statistical section in the statistical section in Supplementary Information).

All effect estimates were drawn from the median of the posterior distribution and the 95% CrI from the 95% percentiles of the distribution. Number of missing data is shown in the Supplementary Information [eTable 1]. All analyses were performed using the R (R, version 4.0.2, Core Team, Vienna, Austria, 2016) software with the beanz package and Stan through brms [10, 11].

Results

Patients

Out of a total of 3918 patients consented and randomized in the SPICE III trial, 3904 (99.6%) with known primary outcome were included. From this population, 1825 (46.7%) were older than 65 years old [Supplementary Information Study Flow Diagram, eFigure 1]. Older patients had a similar severity of illness as measured by APACHE II scores (after removing the age component) and lower PaO2/FiO2 ratios than younger patients. In addition, a higher proportion of older patients were admitted with a cardiovascular diagnosis or from the general ward, and fewer were admitted with trauma or from an emergency department. Detailed baseline characteristics are shown in Table 1.

Table 1 Baseline characteristics of the patients according to the categories of age and the clusters identified

Clusters identified

From patients 3322/3904 (85.1%) with a complete dataset, cluster analysis identified two distinct clusters, one with 976 (29.4%) patients (cluster 1) and another with 2346 (70.6%) patients (cluster 2) [Supplementary Information Study Flow Diagram, eFigure 1]. The optimal number of clusters was confirmed to be two by the prediction strength and the plot of Gower’s distance as shown in the Supplementary Information [eFigure 2 A–B]. While patients in cluster 1 were admitted mainly from the operating room due to cardiovascular or gastrointestinal diseases, patients in cluster 2 were predominantly admitted due to non-operative reasons and from the emergency department due to sepsis or respiratory disease (Table 1). Patients in cluster 2 also had a higher median APACHE II scores and lower PaO2/FiO2 ratios and were more often sedated at baseline with midazolam and/or ketamine, but less often sedated with propofol compared to cluster 1. The number of patients allocated to DEX or usual care sedation were similar among the groups [Supplementary Information Study Flow Diagram, eFigure 1].

Age related heterogeneity of dexmedetomidine effect on mortality

The overall 90-day mortality was higher in older patients, 700/1825 (38.4%) vs 435/2079 (20.9%) in younger patients (Supplementary Information eFigure 3). However, the probability that assignment to DEX group results in lower OR for mortality in older patients (more beneficial), compared to young patients, was 99.8% (Table 2). There was an interaction between age and the use of DEX and a 99.3% probability of lower mortality with the use of DEX with increasing age (Fig. 1a).

Table 2 Primary and secondary outcomes—effect estimates according to age and clusters
Fig. 1
figure 1

Age-related heterogeneity of treatment effect—dexmedetomidine and mortality. OR odds ratio. Values less than 1 indicate lower mortality. a Marginal effect plot for the interaction between the allocation group and age, as a continuous variable, for 90-day mortality. b The posterior distribution of mortality, depicted as odds ratios. The probability of benefit (OR < 1) is 97.7% in patients > 65 years old with 98.5% probability of harm in patients ≤ 65 years old (OR > 1)

The use of DEX resulted in lower mortality compared to usual care in patients older than 65 years (OR 0.83 [95% CrI 0.68–1]), with a probability of benefit of 97.7% and of a 10% mortality benefit of 81.7% (Fig. 1b, Table 2). On the other hand, DEX led to higher mortality in patients younger than 65 years (OR 1.26 [95% CrI 1.02–1.56), with a 98.5% probability of harm and an 89.9% probability of > 10% increased mortality (Fig. 1b). The use of different priors did not materially change these findings (Table 2 and Supplementary Information eFigures 4 and 5).

Cluster related heterogeneity of dexmedetomidine effect on mortality

The overall 90-day mortality rate was higher in cluster 2 at 731/2346 (31.2%) vs 240/976 (24.6%) in cluster 1. There was a higher probability of a mortality benefit with DEX in cluster 1 patients (OR 0.86 [95% CrI 0.65–1.14], 84.8% probability of benefit, compared to cluster 2 (OR 1.04 [95% CrI 0.87–1.24], (Table 2 and Supplementary Information eFigure 6). The probability that DEX would result in lower OR of death in patients in cluster 1 compared to patients in cluster 2 was 86.4%. The use of different priors did not materially change these findings (Table 2 and Supplementary Information eFigure 7).

Interaction between age and severity of illness

The results of the model assessing interactions between age, severity of illness and use of DEX are shown in Fig. 2a–d. In older patients, there was a high probability of lower mortality (OR < 1) with the use of DEX which declined with increasing APACHE II, panels a, b. A high probability of > 10% benefit with the use of DEX was seen in older patients with APACHE II < 20, panel c. In young patients, > 10% harm associated with the use of DEX was seen with high probability in those with an APACHE II ≥ 25, panel d.

Fig. 2
figure 2

Risk of death and interaction between age, severity of illness and dexmedetomidine treatment. APACHE II Acute Physiology and Chronic Health Evaluation II, CrI credible interval. The effect estimates, odds ratio (OR) for the interaction between dexmedetomidine allocation, age category and six different cut-offs of APACHE II are presented. OR < 1.0 represents a favorable outcome and > 1.0 represents unfavorable outcome with the use of dexmedetomidine. a Odd ratios according to age category, age group > 65 years depicts a high probability of OR < 1.0 with increased APACHEII. b Probability of benefit with the allocation to dexmedetomidine was higher in older age group. c Probability of > 10% benefit with the allocation to dexmedetomidine was higher in patients > 65 years but declined with increasing APACHEII. d Probability of > 10% harm with the allocation to dexmedetomidine was higher in patients ≤ 65 years and increased with rising APACHEII

Interaction between clusters and age

The results of the model assessing interactions between cluster, age and use of DEX are shown in Fig. 3a–d. The probability of reduced mortality (OR < 1) with DEX was high when age was greater than 50 years in patients in cluster 1 (panels a, b). In cluster 2, however, a high probability of reduced mortality with DEX was found only when age was greater than 80 years (panel b). A high probability of > 10% benefit in cluster 1 was seen with DEX with age above 60 years, panel c. In cluster 1, the probability of > 10%, harm, was low with age greater than 40 years, but in cluster 2, this probability remained high up to age of 70 years (panels c, d).

Fig. 3
figure 3

Risk of death and interaction between clusters, age and dexmedetomidine treatment. CrI credible interval. The effect estimates, odds ratio (OR) for the interaction between dexmedetomidine allocation, cluster assignment and six different cut-offs of age categories are presented. OR < 1.0 represents a favorable outcome and > 1.0 represents unfavorable outcome with the use of dexmedetomidine. a Odd ratios according to cluster, operative cluster 1 depicts a high probability of benefit with increased age. b Probability of benefit with the allocation to dexmedetomidine was higher in cluster 1 but mainly in those > 50 years old. c Probability of > 10% benefit with the allocation to dexmedetomidine was higher in cluster 1 but mainly in those > 60 years. d Probability of > 10% harm with the allocation to dexmedetomidine was higher in cluster 2, non-operative, mainly in patients younger than 50 years

Age and cluster heterogeneity of treatment effect and secondary outcomes

The probability of more days alive, and free of delirium and coma was 97.4% in patients older than 65 years, allocated to DEX, with a neutral effect in younger patients (Table 2 and Supplementary Information eFigure 8). Conversely, in clusters 1 and 2, the probability was lower, at 88.5% and 77.3%, respectively (Table 2).

A high probability of greater number of ventilator-free days with DEX treatment was seen in older patients and in those assigned to cluster 1 (Table 2 and Supplementary Information eFigure 9).

Adverse events

Adverse and serious adverse events, reported by site investigators, were distributed at comparable rates across age and cluster categories. [Supplementary Information. eTable 3].

Discussion

Key findings

In this Bayesian secondary analysis of the SPICE III trial, we assessed the HTE of early sedation with DEX according to age and two distinct clusters, based on 12 baseline characteristics. Using this approach, we found evidence of HTE on 90-day mortality and key secondary outcomes. In particular, early sedation with DEX reduced 90-day mortality in patients older than 65 years with high (97.7%) probability. In addition, early DEX was more likely to improve 90-day mortality in the mainly operative cluster. Conversely, we observed a high probability of increased 90-day mortality in younger patients and in the cluster of non-operative patients, which increased with increasing APACHE II scores. Importantly, median age in both clusters was similar, suggesting that other factors such as sedative drug interaction or neuroendocrine response, independent of age, may have accounted for HTE of DEX in these patients.

Relationship to previous studies

In the SPICE III trial, a pre-specified subgroup analysis of age, dichotomized at the cohort median, suggested a significant interaction with DEX treatment [3]. This interaction, however, was not observed within the subgroup of operative or non-operative patients [3]. In addition, there was no detailed analysis of treatment effects in any of the subgroups for the secondary outcomes or occurrence of adverse events. In this analysis, we adopted a Bayesian approach to evaluate in detail the HTE of DEX on the 90-day mortality, relevant secondary outcomes and adverse events. In addition to the age category, we identified two distinct clusters, operative and non-operative, and assessed the interaction according to different six age and six APACHEII categories.

Despite an expected age effect on the pharmacokinetics and efficacy of DEX [12, 13], no previous sedation trials have assessed age related HTE for mortality or other sedation related outcomes. This is possibly due to a significantly smaller sample size in previous randomized sedation trials compared with SPICE III [14,15,16,17]. In addition, Bayesian and/or cluster analysis was not a feature of previous sedation trials.

In the SLEAP sedation interruption trial, HTE was reported in subgroup analysis, suggesting a beneficial effect of sedative interruption for surgical/trauma patients when compared with medical non-operative patients [18]. HTE according to age, however, was not assessed. More recently, the differential effect of vasopressors in septic patients older than 65 years was investigated [19], suggesting more awareness of potential and important age related HTE in critically ill patients.

Implications of study findings

Our findings imply that the early use of DEX in ventilated critically ill patients is likely beneficial in patients older than 65 years across broad diagnostic categories and illness severity. Conversely, early DEX-based sedation in younger patients, appears likely to increase mortality, especially in non-operative critically ill patients with high severity of illness. These findings also imply that the positive effects of early DEX on delirium-, coma- and ventilator-free time are preserved in older patients irrespective of cluster assignment. Our findings also suggest a comparable distribution of adverse or serious adverse events according to age or cluster categories.

Due to the post hoc secondary nature of this analysis, no inference of causality or mediation of effect can be claimed. Thus, caution must be taken in extrapolating these findings into practice. Nonetheless, the use of early DEX as sedative in the older and operative population, appears safe and may be efficacious. Its use in younger non-operative patients, however, appears unwarranted. The mechanism for the underlying effect remains uncertain; further studies and analysis are needed to address the cause/s of the observed mortality.

Strengths and limitations

The SPICE III trial is the largest in the field of sedation, conducted in 74 ICUs in 8 countries with diverse health systems and different models of care, and virtually complete 90-day follow-up. The plan of the present analysis (including prior choices) was based on available evidence and pre-defined priors. Also, the clusters were not arbitrarily defined but rather derived from an unsupervised approach. Conventional subgroup analysis, as used in the primary trial, is at increased risk of type I and II errors. In the present paper, the use of hierarchical models according to a simple regression and shrinkage model limits the degree of type I error. Moreover, the posterior probabilities derived from the Bayesian analysis provide more power than the traditional interaction tests used for subgroup analyses [20, 21]. In addition, we assessed factors other than age as potential sources of HTE, including the interaction of age with illness severity, adding to the findings of the original trial. Finally, due to the large sample size in the original trial, the clusters assessed still had a considerable number of patients, lending statistical robustness to our observations.

The present study has some limitations. First, this is an unplanned secondary post hoc analysis. Second, the limitations of the original trial influence the present analysis. They include the lack of blinding, the inclusion of patients needing deep sedation, and the absence of a strictly protocolized strategy for managing sedation or delirium. Third, in defining the clusters, missingness in some baseline characteristics was relatively high. Nonetheless, robust sensitivity analysis revealed, exactly, the same clusters. Finally, our results are not generalizable outside the study population and should not be extrapolated to all critically ill patients. In addition, with possible selection bias in the clustering process, HTE could be present in different groups considering characteristics not assessed in this study. Due to these limitations, and because the clustering procedure is mainly descriptive, our findings can only provide the statistical underpinning and rationale for further investigations and stratification in clinical trials.

Conclusion

This secondary analysis of the SPICE III trial identified HTE of early use of DEX, vs. usual care, as primary sedative in critically ill patients needing mechanical ventilation. In particular, a high probability of reduced 90-day mortality was observed in older patients regardless of operative or non-operative cluster status. Conversely, a high probability of increased 90-day mortality was observed in younger patients of non-operative cluster status. These findings warrant confirmation in future randomized trials.