Trajectories of mean arterial pressure/norepinephrine equivalent dose index in patients with septic shock receiving low-dose hydrocortisone: a retrospective cohort study in Thailand
Article information
Abstract
Background
We aimed to analyze the trajectories of the mean arterial pressure/norepinephrine equivalent dose (MAP/NEQ) index in patients with septic shock treated with low-dose hydrocortisone and to determine the association of these trajectories with mortality.
Methods
A retrospective cohort of 203 patients with septic shock receiving low-dose hydrocortisone was examined. MAP and NEQ data were collected from electronic health records, and group-based trajectory modeling was employed to identify distinct patterns in the MAP/NEQ index over the initial 72 hours of treatment. Univariable and multivariable logistic regression analyses were conducted to assess the associations of MAP/NEQ index trajectories with clinical variables and in-hospital mortality.
Results
The overall in-hospital mortality rate was 56.2%. Three MAP/NEQ index trajectory patterns were identified: unchanged (76.8%), gradually increased (14.3%), and rapidly increased (8.9%). The shock reversal rates were 50% for the unchanged group, 89.7% for the gradually increased group, and 100% for the rapidly increased group. Compared with the unchanged group, both gradually and rapidly increased groups were associated with significantly lower mortality, with adjusted odds ratios of 0.15 (95% CI, 0.05–0.40; P<0.001) and 0.29 (95% CI, 0.09–0.92; P=0.035), respectively.
Conclusions
In patients with septic shock treated with low-dose hydrocortisone, gradually and rapidly increased MAP/NEQ index trajectories were associated with significantly lower mortality risks and higher rates of shock reversal compared to those with unchanged trajectories. These findings highlight the importance of monitoring the MAP/NEQ index to guide treatment and improve septic shock outcomes.
INTRODUCTION
Septic shock, a medical emergency characterized by severe infection, leads to reduced blood flow and oxygen delivery to tissues. This condition is frequently encountered in intensive care units (ICUs) and is recognized by the World Health Organization as a significant global public health issue. Approximately 85% of sepsis cases and sepsis-related deaths occur in low-to-middle-income countries [1]. A systematic review reported a sepsis-related hospital mortality rate of 27% globally [2]. In Thailand, the prevalence of shock among patients with sepsis is 30.5%, with a 28-day mortality rate of 38% [3].
The use of corticosteroids in septic shock has been debated and researched for more than four decades. The benefits of corticosteroids remain controversial owing to varying patient selection criteria, types and doses of corticosteroids, and differences in administration duration, leading to inconsistent primary outcomes regarding mortality. However, secondary outcomes have consistently shown that corticosteroids reduce the duration of vasopressor use [4,5]. Current treatment guidelines for septic shock, based on the 2021 guidelines of the Surviving Sepsis Campaign, recommend the use of intravenous corticosteroids in patients who require escalating doses of vasopressors [6].
Recently, the ratio of mean arterial pressure (MAP) to norepinephrine equivalent dose (NEQ), known as the MAP/NEQ index, has been introduced as a novel marker for assessing septic shock severity. The MAP/NEQ demonstrated improved prognostic validity compared with commonly used sepsis severity scores and vasopressor dose alone [7]. The MAP/NEQ has been used to predict early mortality in patients with shock treated with vasopressors [8] and to determine the appropriate timing for initiating enteral nutrition in patients with shock [9].
As corticosteroids reduce the duration of vasopressor use, this can be expected to influence the trajectory of the MAP/NEQ index over time after the initiation of hydrocortisone administration. In this study, we aimed to characterize hemodynamic response trajectories in patients with septic shock receiving low-dose hydrocortisone in real-world clinical practice rather than to test treatment efficacy in a trial setting. We hypothesized that tracking the MAP/NEQ index trajectories after starting low-dose intravenous hydrocortisone treatment would predict septic shock recovery and in-hospital mortality.
MATERIALS AND METHODS
The study was approved by the Human Ethics Committee of the Faculty of Medicine at Prince of Songkla University (No. 66-281-14-4). All methods were performed in accordance with the relevant guidelines and regulations, and the study adhered to the Declaration of Helsinki. The requirement for informed consent was waived because of the retrospective study design.
Study Design
This was a single-center, retrospective cohort study of patients with septic shock admitted to the medical intensive care unit (MICU) of Songklanagarind Hospital from January 1, 2017, to December 31, 2022. The hospital is an 850-bed academic tertiary care facility with 12 MICU beds. Sepsis management at our hospital followed a sepsis protocol adapted from the 2016 Surviving Sepsis Campaign guidelines [10]. Vasoactive agents available at our hospital include norepinephrine, epinephrine, dopamine, and dobutamine. The criteria for initiating low-dose hydrocortisone treatment in patients with septic shock receiving high-dose vasopressors were left to the physician's discretion.
This study included adult patients aged ≥18 years diagnosed with septic shock according to the Sepsis-3 definition [11] and treated with hydrocortisone at a dose of 50 mg every 6 hours. Another inclusion criterion was the availability of the vital signs and vasopressor use data at 0 hours and at two or more of five time points at 6, 12, 24, 48, and 72 hours after receiving hydrocortisone, with an allowance of ±1 hour. Patients with preexisting adrenal insufficiency and those who received steroids before developing septic shock were excluded. Patients who received hydrocortisone treatment exceeding 200 mg/day for septic shock or for other indications were also excluded. Additionally, pregnant patients were excluded from the study.
Data Collection
Patient data collected from electronic medical records were sex, age, height, weight, underlying diseases, date of hospital admission, initial lactate level, and severity scores, including Sequential Organ Failure Assessment (SOFA) and non-cardiovascular SOFA scores. Treatment information, including the administration of antibiotics, vasopressors, inotropes (dosage and duration), and hydrocortisone (dosage and time of initiation), was also obtained. Vital signs were recorded at the time of shock diagnosis and at 0, 6, 12, 24, 48, and 72 hours after receiving hydrocortisone. Hyperglycemia or hypernatremia was also noted, and information on the number of days spent in the hospital and the discharge status was collected.
NEQ was defined as follows: norepinephrine dose (µg/kg/min) + epinephrine dose (µg/kg/min) + 1/10 × phenylephrine dose (µg/kg/min) + 1/100 × dopamine dose (µg/kg/min) + 1/8 × metaraminol dose (µg/kg/min) + 2.5 × vasopressin dose (U/min) + 10 × angiotensin II dose (µg/kg/min) [12].
Shock reversal was defined as recovery to MAP ≥65 mm Hg with no vasopressor administration for at least 24 hours [13]. Hyperglycemia was defined as blood sugar ≥150 mg/dl on two consecutive occasions within 72 hours of starting hydrocortisone administration [13]. Hypernatremia was defined as a blood sodium level ≥150 mEq/L at any time within 72 hours of starting hydrocortisone administration [14].
Outcome Assessment
The primary outcome was the in-hospital mortality rate. Secondary outcomes included the shock reversal rate, time from septic shock to shock reversal, hyperglycemia or hypernatremia, and length of hospital stay.
Statistical Analyses
Group-based trajectory models (GBTMs) are useful for identifying clusters of individuals with distinct developmental trajectories within a cohort. In this study, individuals within the same cluster followed similar progressions in the MAP/NEQ index measured repeatedly at six time points. For the GBTMs, we used the censored normal model because the MAP/NEQ index is a continuous variable. Initially, we assumed the trajectory patterns to be cubic, with the MAP/NEQ index as the dependent variable and time (hour) as the independent variable.
Model selection involved determining both the optimal number of groups and the shape of each trajectory. Several models with varying group numbers and trajectory forms were compared to identify the best fit. The model fit was evaluated using several statistical criteria, including the Bayesian Information Criterion (BIC), the log Bayes factor (2 logeB10), and the average posterior probability for group membership. A smaller BIC and a value of 2 logeB10 >10 were considered indicators of superior model fit. Each trajectory group was required to have an average posterior probability >0.7 and a minimum group membership of 5% to ensure adequate classification and model stability [15,16].
For descriptive analyses, continuous variables are presented as the median with interquartile range, whereas categorical variables are expressed as percentage. Baseline variables were compared between groups using the Kruskal-Wallis test or the Mann-Whitney U-test for continuous variables and the chi-square test or Fisher’s exact test for categorical variables. When the overall trend was significant, post hoc pairwise comparisons were conducted using the Bonferroni correction.
Patients were classified into three MAP/NEQ index trajectory groups based on individual linear slope parameters as unchanged, gradually increased, and rapidly increased. The relationship between the MAP/NEQ index trajectory group and the primary outcome was examined using multivariable logistic regression models. The covariates included in the model for adjustment were age, nosocomial infection, intravenous fluid resuscitation, appropriate antibiotics, non-cardiovascular SOFA score, and initial lactate level. Odds ratios (ORs) and 95% CIs were estimated using logistic regression. The software used for GBTM was Stata version 16 (StataCorp.), with an added user-written Stata command called "Traj," developed by Jones and Nagin for trajectory modeling [17].
RESULTS
In this study, 265 patients with septic shock admitted to the MICU were eligible for inclusion. After excluding 62 patients (4 with adrenal insufficiency, 32 with preexisting steroid treatment, 4 receiving a different steroid dose, and 22 with incomplete data), 203 patients were included in the study cohort. A flowchart of the patient selection process is shown in Figure 1. The baseline characteristics according to primary outcome are presented in Supplementary Table 1. Type of infection, positive hemoculture, amount of intravenous fluid, appropriate antibiotic use, and severity score significantly affected mortality.
Model comparison using the BIC, the log Bayes factor, and minimum group membership favored the model with three trajectories as having the optimal fit. By contrast, the models with two- or four-groups yielded poorer BIC values and included implausibly small subgroups (Supplementary Table 2). The assessment of model adequacy for the three-trajectory model is presented in Supplementary Table 3. The average posterior probabilities for all trajectory groups exceeded the recommended cutoff of 0.7, indicating acceptable classification quality.
Figure 2 illustrates the MAP/NEQ trajectories measured at 0, 6, 12, 24, 48, and 72 hours after hydrocortisone administration. The study population was categorized into three distinct patterns of unchanged (n=156, 76.8%), gradually increased (n=29, 14.3%), and rapidly increased (n=18, 8.9%). Distinct trajectories emerged within the first 6 hours and were clearly separated by 24 hours. Patients whose MAP/NEQ index increased above 500 mm Hg/µg/kg/min during this early period generally followed favorable (gradually or rapidly increasing) trajectories, whereas those remaining below 300 mm Hg/µg/kg/min typically belonged to the unchanged, non-responder group.
Trajectories of the MAP/NEQ index measured at 0, 6, 12, 24, 48, and 72 hours after hydrocortisone initiation. Three distinct patterns were identified: group 1 (blue), unchanged (non-responders, 76.8%); group 2 (red), gradually increased (delayed responders, 14.3%); and group 3 (green), rapidly increased (early responders, 8.9%). Divergence between trajectories became evident at 6 hours and was clearly pronounced by 24 hours. MAP: mean arterial pressure; NEQ: norepinephrine equivalent dose.
The baseline characteristics of the MAP/NEQ trajectory groups are presented in Table 1. The median patient age was 69 years, and 53.2% of the patients were male. Significant differences among the three groups were observed in terms of catheter-associated urinary tract infection, baseline SOFA score, non-cardiovascular SOFA score, initial lactate level, MAP/NEQ index at the time of starting hydrocortisone treatment, shock reversal rate, time to shock reversal, occurrence of hypernatremia, and in-hospital mortality.
The hemodynamic variables and vasoactive agents used at the time of shock diagnosis and initiation of hydrocortisone treatment are detailed in Table 2. The initial median MAP/NEQ index was 1,103.1 mm Hg/µg/kg/min, with no significant differences among the three groups. However, after initiating hydrocortisone treatment, the median MAP/NEQ indices for the unchanged and gradually increased groups were 245.4 mm Hg/µg/kg/min and 432.2 mm Hg/µg/kg/min, respectively, which were significantly lower than the median MAP/NEQ index of 1,318.5 mm Hg/µg/kg/min for the rapidly increased group.
Hemodynamic variables at the time of shock diagnosis and of starting hydrocortisone according to the MAP/NEQ index trajectories
Table 3 shows the associations between MAP/NEQ index trajectory groups and mortality. Compared with the unchanged group, the gradually increased and rapidly increased groups were significantly associated with a lower risk of mortality after adjusting for confounders (OR, 0.15; 95% CI, 0.05–0.40; P<0.001 and OR, 0.29; 95% CI, 0.09-0.92; P=0.035, respectively).
DISCUSSION
This study focused on analyzing MAP/NEQ index trajectories in patients with septic shock receiving low-dose hydrocortisone as part of clinical care. Hydrocortisone was considered adjunctive therapy rather than an intervention under investigation. We identified three MAP/NEQ index trajectories after initiating hydrocortisone treatment. The groups with an increased MAP/NEQ index, whether gradual or rapid, were significantly associated with a lower risk of in-hospital mortality and higher rates of shock reversal than the unchanged group.
To the best of our knowledge, this is the first study to evaluate the trajectories of hemodynamic responses to low-dose bolus hydrocortisone using the MAP/NEQ index. Previous studies have examined hemodynamic responses to various hydrocortisone preparations. For instance, Bargate et al. found that, in patients with septic shock, a bolus of hydrocortisone with fludrocortisone increased the systolic blood pressure, diastolic blood pressure, and MAP within 2–4 hours of steroid initiation [18]. Similarly, Keh et al. [19] observed an increase in MAP on day 1 and a decrease in norepinephrine requirements with continuous hydrocortisone administration. Our study used the MAP/NEQ index, which represents the ratio of blood pressure to vasopressor dose and is analogous to the PaO2/FiO2 ratio used in respiratory physiology [7]. Our results showed that the MAP/NEQ index trajectories were notably different 24 hours after starting hydrocortisone treatment. The increase in MAP and decrease in NEQ can be attributed to the effects of hydrocortisone, which include direct effects on vascular tone, left ventricular afterload, immunologic effects, and levels of endothelial glucocorticoid receptors that reduce vasoplegia [20].
Our study reflects the real-life practice of using low-dose hydrocortisone for septic shock, with a median of 7 hours from septic shock onset to steroid administration. This aligns with guidelines recommending the initiation of steroids in patients treated with vasopressors for a minimum of 4 hours. The median NEQ dose was 0.25 µg/kg/min, which is consistent with the guideline-recommended dose for corticosteroids [6]. The median time from septic shock onset to shock reversal was 3 days, similar to that in the hydrocortisone group in the Adjunctive Corticosteroid Treatment in Critically Ill Patients with Septic Shock trial [21]. We used intermittent bolus doses of hydrocortisone for ease of administration, which have been associated with higher shock reversal rates on day 7 than continuous infusion [22].
We identified three patterns of MAP/NEQ index response to hydrocortisone. Unfortunately, most patients (76%) exhibited an unchanged pattern, potentially because the patients were steroid non-responders or were severely ill. This lack of response might reflect tissue resistance to glucocorticoids at the receptor level [23], genetic variability in glucocorticoid receptor sensitivity [24], or irreversible vasoplegia occurring in the late phase of septic shock. Previous studies have provided information on corticosteroid responders. For example, the corticosteroid therapy of septic shock study reported that 46.7% of patients with septic shock were steroid non-responders, which was defined as the absence of a response to a corticotropin test [14]. Similarly, a Thai study reported a non-responsiveness rate of 55.2% among patients with septic shock [25]. No benefit was found in using steroids in patients without a sufficient response in the adrenocorticotropic hormone (ACTH) stimulation test, defined as an increase ≤ 9 μg/dl in cortisol 30–60 minutes after ACTH administration. Therefore, steroid responsiveness might not be an appropriate index for hemodynamic response to hydrocortisone. Patients in our study appeared to be more severely affected than those in previous studies, as indicated by a mortality rate of 56.2% in patients with septic shock, compared to ≤50% in two recent major randomized trials of steroid use in septic shock [13,21].
An original study by Bosch et al. [7] reported the prognostic validity of the MAP/NEQ index using MIMIC-III and an electronic ICU database. The MAP/NEQ index within 24 hours of ICU admission ranged from 39 to 8,396 mm Hg/µg/kg/min in the MIMIC-III database. Tertiles were used to create MAP/NEQ cutoffs of < 136, 136–324, and ≥325 mm Hg/µg/kg/min, classified as severe, moderate, and mild severity, respectively. A previous study by Yang et al. [8] identified the MAP/NEQ index as an early mortality predictor in patients with shock using the MIMIC-IV database. Higher MAP/NEQ indices were associated with shorter hospital stays and lower mortality rates. A MAP/NEQ index ≤183 mm Hg/µg/kg/min was linked to a significantly increased mortality risk.
The MAP/NEQ index in our study was higher than those reported in previous studies by Bosch et al. [7] and Yang et al. [8] because we used a different NEQ formula. Previous studies used a formula developed in 2008 that included only norepinephrine, dopamine, epinephrine, and phenylephrine. The units for most vasoactive agents were µg/min, except for dopamine, which was measured in µg/kg/min with a conversion ratio of 1:2. Recently, Kotani et al. [12] updated the NEQ formula to include new vasopressors such as vasopressin and angiotensin II, suggesting its extension for use with the MAP/NEQ index. The updated formula uses units of µg/kg/min, with a conversion ratio of 1:100 for dopamine. This approach yields smaller denominator values and consequently higher MAP/NEQ ratios than those reported in earlier studies. Subsequent studies have shown that the updated NEQ formula demonstrated favorable predictive performance for mortality risk in patients with sepsis in ICUs [26] and supports its use in settings involving alternative vasoactive agents, such as angiotensin II [27] and methylene blue [28]. Our results suggest that using a single MAP/NEQ index as a prognostic marker in patients with septic shock receiving hydrocortisone can be more effective when measured at the time of starting hydrocortisone treatment than at the time of septic shock diagnosis.
The gradually and rapidly increased trajectory groups demonstrated better outcomes in terms of shock reversal, supporting the MAP/NEQ index as a reliable indicator of hemodynamic recovery in patients receiving hydrocortisone. In our study, MAP/NEQ trajectories were associated with both in-hospital mortality and shock reversal outcomes. Hydrocortisone use can lead to metabolic complications such as hyperglycemia and hypernatremia, both of which can worsen patient outcomes. In our study, hyperglycemia (≥150 mg/dl) occurred in approximately 60% of patients after initiation of hydrocortisone administration. Previous evidence suggests that a stable mean blood glucose level between 120 and 140 mg/dl is associated with lower ICU mortality and reduced hypoglycemia risk in patients with sepsis but without diabetes [29]. Hypernatremia (≥150 mEq/L) was observed mainly in the unchanged MAP/NEQ group and has been linked to poor prognosis in septic shock [30]. These findings suggest that, in patients who remain hemodynamically unchanged after 24–48 hours, the development of hypernatremia or persistent hyperglycemia should prompt reconsideration or discontinuation of hydrocortisone to avoid harm. Serial MAP/NEQ monitoring every 6–12 hours can help identify early steroid responders, whereas a persistently flat trajectory, characterized by MAP/NEQ values less than 300 mm Hg/µg/kg/min after 24–48 hours, should prompt reassessment of therapy. This trajectory-based approach could be incorporated into sepsis protocols to guide the timely escalation or de-escalation of corticosteroid treatment.
Our study has several limitations. First, its nature as a single-center study might introduce selection bias. Second, the limited sample size might have reduced the statistical power of the association analysis. Third, rapid hemodynamic fluctuations might have been missed because of the 6-hour data collection intervals. Fourth, data on potential confounders influencing mortality outcomes might not have been fully adjusted. Although fluid balance is an important confounder, only the total volume of fluid resuscitation was recorded because net balance was not preplanned in the study protocol. Moreover, other inflammatory markers, such as C-reactive protein and procalcitonin, were not routinely measured in our hospital. We acknowledge that the observed associations might reflect residual confounding or reverse causality; therefore, the trajectory–mortality relationship should be interpreted as prognostic rather than causal, given the retrospective design and lack of randomization. Fifth, only patients receiving bolus hydrocortisone were included; hence, the findings should not be directly extrapolated to continuous infusion regimens. Finally, the rapidly increased trajectory group, representing less than 10% of the cohort, might have reduced statistical stability and widened CIs, although model adequacy indices support the robustness of the trajectory classification.
In conclusion, three MAP/NEQ index trajectories in patients receiving hydrocortisone were identified using GBTMs. The rapidly and gradually increased trajectory patterns were associated with lower mortality and higher shock reversal than the unchanged trajectory pattern. The MAP/NEQ index trajectory patterns demonstrated the utility of MAP/NEQ index measurement in predicting outcomes in patients with septic shock receiving hydrocortisone. However, our results require further investigation with larger sample sizes and various populations.
KEY MESSAGES
▪ Three distinct mean arterial pressure/norepinephrine equivalent dose (MAP/NEQ) index trajectory patterns were identified in septic shock patients treated with low-dose hydrocortisone: unchanged, gradually increased, and rapidly increased.
▪ Gradually and rapidly increased MAP/NEQ trajectories were associated with significantly lower in-hospital mortality and higher shock reversal rates compared to unchanged trajectories.
Notes
CONFLICT OF INTEREST
No potential conflict of interest relevant to this article was reported.
FUNDING
None.
ACKNOWLEDGMENTS
We appreciate the prompt responses of the Division of Digital Innovation and Data Analytics, Faculty of Medicine, Prince of Songkhla University; and the Medical Record Section, Songklanagarind Hospital.
AUTHOR CONTRIBUTIONS
Conceptualization: AC, SS, VV. Data curation: AC, VV. Formal analysis: AC, SS, VV. Methodology: AC, SS, VV. Project administration: VV. Visualization: VV. Writing original draft: VV. Writing - review & editing: AC, SS, VV. All authors read and agreed to the published version of the manuscript.
SUPPLEMENTARY MATERIALS
Supplementary materials can be found via https://doi.org/10.4266/acc.000125.
Characteristics of patients grouped by survival or death during the study period
acc-000125-Supplementary-Table-1.pdfModel fit statistics and group adequacy for trajectory models with 1–4 groups
acc-000125-Supplementary-Table-2.pdfModel adequacy results for the final three-trajectory model
acc-000125-Supplementary-Table-3.pdf