Abstract

Background. Although serum calcium has been proven to be a predictor of mortality in a wide range of diseases, its prognostic value in critically ill patients with cardiogenic shock (CS) remains unknown. This retrospective observational study is aimed at investigating the association of admission calcium with mortality among CS patients. Methods. Critically ill patients diagnosed with CS in the Medical Information Mart for Intensive Care-III (MIMIC-III) database were included in our study. The study endpoints included 30-day, 90-day, and 365-day all-cause mortalities. First, admission serum ionized calcium (iCa) and total calcium (tCa) levels were analyzed as continuous variables using restricted cubic spline Cox regression models to evaluate the possible nonlinear relationship between serum calcium and mortality. Second, patients with CS were assigned to four groups according to the quartiles (Q1-Q4) of serum iCa and tCa levels, respectively. In addition, multivariable Cox regression analyses were used to assess the independent association of the quartiles of iCa and tCa with clinical outcomes. Results. A total of 921 patients hospitalized with CS were enrolled in this study. A nonlinear relationship between serum calcium levels and 30-day mortality was observed (all values for nonlinear ). Furthermore, multivariable Cox analysis showed that compared with the reference quartile (Q3: ), the lowest serum iCa level quartile (Q1: ) was independently associated with an increased risk of 30-day mortality (Q1 vs. Q3: HR 1.35, 95% CI 1.00-1.83, ), 90-day mortality (Q1 vs. Q3: HR 1.36, 95% CI 1.03-1.80, ), and 365-day mortality (Q1 vs. Q3: HR 1.28, 95% CI 1.01-1.67, ) in patients with CS. Conclusions. Lower serum iCa levels on admission were potential predictors of an increased risk of mortality in critically ill patients with CS.

1. Introduction

Cardiogenic shock (CS) is a severely diminished-cardiac-output state resulting in life-threatening end-organ hypoperfusion and hypoxia [1, 2]. There are numerous causes of CS, including acute myocardial infarction (AMI), severe myocarditis, and end-stage dilated cardiomyopathy [3]. In addition, CS is the most common cause of death for patients hospitalized with AMI [4]. Despite advances in treatment, the in-hospital mortality remains unacceptably high (27%-51%) [57]. As mortality peaks within the first 48 hours after CS onset, it is necessary to find an accurate yet user-friendly predictor for early risk stratification to provide more accurate prognostic information and help implement appropriate treatment [8].

Serum calcium plays an essential role in a range of biological processes related to cardiovascular diseases, including myocardial contraction and relaxation, nerve transmission, vascular smooth muscle contractile activity, platelet adhesion, and blood coagulation [911]. Thus, alterations in serum calcium concentrations might interfere with myocardial function and cause severe cardiovascular complications and organ dysfunctions [12]. Derangement in serum calcium is known to be extremely common in the intensive care unit (ICU) setting, and several previous studies have shown that increased or/and decreased levels of serum calcium were independent risk predictors for mortality in patients with AMI [1318], heart failure [19], acute kidney injury (AKI) [20], and acute stroke [21] or individuals in the general population [2224]; they were also tightly related to cardiovascular risk factors such as hyperlipidemia, hyperglycemia, and hypertension [10, 16].

To the best of our knowledge, there have been no epidemiological studies exploring the prognostic value of serum calcium among critically ill patients with CS. As a common urgent critical illness, patients with CS are at greater risks of kidney injury, impaired gastrointestinal function, or heightened neurohormonal activation, which could affect serum calcium homeostasis [2528]; it remains unclear whether abnormalities in calcium levels could affect the prognosis of CS. Additionally, most previous studies only focused on the serum tCa [1315, 2123, 2931]. Considering the limitations of tCa measurements in the identification of true calcium derangements (i.e., its dependency on serum albumin levels) [3133], the prognostic ability of serum iCa was also explored in this study.

In the present study, we aimed to investigate the possible association of admission serum iCa and tCa levels with the risks of all-cause mortality in patients with CS.

2. Methods

2.1. Study Design

This is a single-center retrospective cohort study, and all the relevant data were collected from the Medical Information Mart for Intensive Care-III (MIMIC-III) database. MIMIC-III is a freely accessible and conveniently sized critical care database covering over 50,000 hospital admissions comprised of 38,645 adults as well as 7,875 neonates admitted to surgical, trauma surgery, coronary, and cardiac surgery recovery ICUs of Beth Israel Deaconess Medical Center (BIDMC) in Boston from 2001 to 2012 [34, 35]. The MIMIC-III database documents contained high-resolution information from hospital monitoring systems (including laboratory data, medication, and hospital administrative data) and bedside monitoring systems (vital signs, caregivers notes, and radiology reports). We passed the “Protecting Human Research Participants” exam and obtained permission to access the dataset (authorization code: 33281932). Furthermore, we conducted this study in accordance with the STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) statement [36].

2.2. Ethical Approval

The establishment of the MIMIC-III database was approved by the Institutional Review Boards (IRB) of the Massachusetts Institute of Technology (No. 0403000206) and BIDMC (2001-P-001699/14). Our study utilized the anonymous data available from this database, and hence, the requirement for informed consent was waived. In summary, the study complied with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.

2.3. Study Population

We included all ICU patients ( years) with the primary diagnosis of CS using International Classification of Diseases, ninth version- (ICD-) 9 diagnosis codes (ICD-9 codes: 785.51) in the MIMIC-III database. Only the data of each patient’s first ICU admission were used in this study. Patients were excluded if they had (1) a secondary diagnosis of hepatic dysfunction, renal failure, acute or chronic pancreatitis, parathyroid diseases, or malignancy on admission; (2) a length of stay in the ICU less than 24 hours; (3) incomplete or unobtainable data of serum iCa and tCa measured during the first 24 hours admission; (4) incomplete follow-up information; or (5) more than 10% of individual data missing.

2.4. Data Extraction, Preparation, and Definitions

Demographics, vital signs, laboratory tests, medications, and others were extracted from the MIMIC-III database using structured query language (SQL) with PostgreSQL (version 9.4.6, http://www.postgresql.org). The code that supports the MIMIC-III documentation and website is publicly available, and contributions from the community of users are encouraged (https://github.com/MIT-LCP/mimic-website).

Baseline demographic variables included age, sex, ethnicity (white or others), and current smoking status (by Natural Language Processing searches in provider notes, categorized as “yes,” or “no/unknown”). We extracted data on the following comorbidities: coronary artery disease (CAD), chronic heart failure (CHF), atrial fibrillation (AF), hypertension, peripheral artery disease (PAD), stroke, diabetes mellitus (DM), and chronic kidney disease (CKD). Vital signs on admission included systolic blood pressure (SBP), diastolic blood pressure (DBP), mean blood pressure (MBP), and heart rate (HR). Laboratory-based data included iCa, tCa, phosphorus, potassium, sodium, chloride, bicarbonate, lactate, anion gap (AG), pH, creatinine, estimated glomerular filtration rate (eGFR), hemoglobin, platelet, and white blood cell count (WBC). The eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-Epi) formula [37]. If patients received a laboratory test more than once during their hospitalization, only the initial test results were included for analysis. Three scoring systems (the Sequential Organ Failure Assessment (SOFA), the Simplified Acute Physiology Score II (SAPS II), and the Glasgow Coma Scale (GCS)) were calculated within the first 24 hours after admission using the values associated with the greatest severity of illness. In addition, treatment information data were also collected, including mechanical ventilation, renal replacement treatment (RRT), and in-hospital medication (inotrope and vasoconstrictor) administration.

2.5. Identification of Cut-Off Values for Serum iCa and tCa Levels

Serum calcium levels were categorized into four groups according to the quartiles (Q1-Q4) of their concentrations.

2.6. Study Outcomes

The primary outcome of our study was 30-day all-cause mortality. Secondary outcomes included 90-day and 365-day all-cause mortality. Patients with missing survival outcome information were excluded from the final cohort.

2.7. Statistical Analysis

The data distribution was examined using the Kolmogorov-Smirnov test. Continuous variables are presented as mean (standardized differences (SD)) or median (interquartile range (IQR)) and categorical variables as total number and percentage. Baseline characteristics of enrolled participants were presented by using either Student -test, Kruskal Wallis rank test, Pearson’s test, or Fisher’s exact test as appropriate.

Restricted cubic spline Cox regression models were used to evaluate the possible nonlinear relationship between serum calcium levels and 30-day all-cause mortality [38]. If the test for nonlinearity was not significant, the test result for overall association and linearity was checked, with significant results indicating linear associations.

The Kaplan-Meier method was used to plot unadjusted survival curves, and the log-rank test was used to compare differences between the quartiles of serum calcium. Moreover, Cox proportional hazards regression analysis was performed to examine the relationship between baseline covariates and each endpoint. We separately included the serum iCa and tCa quartiles in multivariable Cox regression models, adjusting for the potential confounders selected based on in the univariable analysis. The third quartile (Q3: ; ) was used as a reference group, and the results are presented as hazard ratios (HRs) with 95% confidence intervals (CIs). Furthermore, subgroup analyses were performed to investigate the association between serum calcium levels and mortality. Moreover, most commonly, CS is an emergency disease characterized by unacceptably high in-hospital mortality; therefore, we mainly focused on the short-term mortality of CS and performed subgroup analyses only for the 30-day mortality.

As extensive missing data might lead to bias, variables with over 20% missing values were not included in the subsequent analyses. Correspondingly, multivariate imputation (MI) was used for variables with less than 20% missing values [39, 40]. Variables for which MI was adopted included SBP, DBP, MBP, HR, lactate, AG, pH, and GCS.

A two-tailed value of less than 0.050 was considered to be statistically significant. All statistical analyses were performed using SPSS software (version 22.0; IBM Corporation, St. Louis, Missouri, USA) and R software (version.3.6.1; The R Project for Statistical Computing, TX, USA; http://www.r-project.org).

3. Results

3.1. Subject and Variable Characteristics

After application of the inclusion and exclusion criteria, the final study cohort consisted of 921 CS patients (Figure 1). The median age of the study cohort was 72 (62-81) years, and 60.3% (555/921) subjects were male. The median admission serum iCa and tCa were 1.11 (1.04-1.17) mmol/L and 8.3 (7.8-8.9) mg/dL, respectively.

In the current study, serum tCa levels were divided into the Q1 group (), Q2 group (), Q3 group (), and Q4 group (). Similarly, serum iCa levels were divided into the Q1 group (), Q2 group (), Q3 group (), and Q4 group (). A total of 223 patients were in the Q1 group (), 224 patients were in the Q2 group (), 231 patients were in the Q3 group (), and 243 patients were in the Q4 group (). The comparison of baseline characteristics stratified by serum iCa quartiles is summarized in Table 1. Compared to those in the Q2-4 groups, patients in the Q1 group () had lower SBP (), lower tCa concentration (), lower bicarbonate concentration (), and higher lactate concentration () (Table 1); they also were more likely to receive RRT () (Table 1). Characteristics including age, sex, comorbidities, and scoring systems were relatively flat across each group (Table 1).

3.2. Relationship between Serum Calcium Levels and Mortality

Restricted cubic spline analyses showed the nonlinear relationships between serum calcium levels (iCa and tCa) and the risk of 30-day mortality. (all values for nonlinear ; Figure 2). In addition, we also observed that the lowest risk of mortality was associated with approximately 1.10 mmol/L for iCa and 9.0 mg/dL for tCa.

3.3. Survival Analysis

Among the 921 CS patients included, 39.1% (360/921) died during the first 30 days, 47.9% (441/921) died during the first 90 days, and 56.0% (516/921) died during the one-year follow-up period. The 30-day mortality was 48.8% in serum iCa of <1.04 mmol/L, 35.3% in 1.04-1.11 mmol/L, 33.3% in 1.11-1.17 mmol/L, and 39.9% in ≥1.17 mmol/L.

Kaplan–Meier curves for all-cause death according to the quartiles of serum calcium are shown in Figure 3. The curves of the quartiles of calcium differed significantly (log-rank test: for 30-day, 90-day, and 365-day all-cause mortalities), and patients in the lowest serum calcium quartile had the highest cumulative incidence of mortality.

In the Cox regression analysis, we analyzed serum iCa and tCa concentrations stratified by quartiles to determine whether serum calcium was associated with all-cause mortality (Table 2). The univariable Cox regression models showed that the lowest serum iCa level quartile () and the lowest serum tCa level quartile () were significant predictors of 30-day, 90-day, and 365-day mortalities compared with the reference group (iCa: 1.06-1.14 mmol/L; tCa: 7.9-8.7 mg/dL) (Table 2; Table S1-6). Furthermore, after adjusting for more confounding factors including age, SBP, DBP, MBP, phosphorus, potassium, chloride, bicarbonate, lactate, AG, creatinine, eGFR, WBC, SOFA, SAPS II, and vasopressor use, only the lowest serum iCa level () remained an independent predictor of 30-day mortality (HR 1.35, 95% CI 1.00-1.83, ), 90-day mortality (HR 1.36, 95% CI 1.03-1.80, ), and 365-day mortality (HR 1.28, 95% CI 1.01-1.67, ) (Table 2; Table S1-6). Furthermore, the highest serum iCa level quartile () was only associated with 90-day mortality in both the univariable and multivariable Cox regression analyses (Table 2; Table S1-6).

3.4. Sensitivity and Subgroup Analysis

We performed subgroup analyses to assess the association between the serum iCa and tCa concentrations and 30-day all-cause mortality (Table 3). Subgroup analyses showed the lowest serum iCa quartile () was also associated with deteriorative mortality in most strata except in patients with a medical history of CHF (). In addition, the results of subgroup analyses of serum tCa were shown in Table S7. Moreover, we used original data for analysis without using the MI method, and 807 patients remained in the final cohort. After adjustment for more confounding factors including age, SBP, DBP, MBP, phosphorus, potassium, chloride, bicarbonate, lactate, creatinine, eGFR, SOFA, and SAPS II, the lowest serum iCa level () still remained an independent predictor of 30-day mortality (HR 1.36, 95% CI 1.01-1.85, ) (Table S8 and S9).

4. Discussion

In the present study, we evaluated 921 patients to measure the association of admission serum iCa and tCa levels with all-cause mortality in critically ill patients with CS. Our main findings can be summarized as follows. First, a nonlinear relationship between admission serum calcium (iCa and tCa) and 30-day all-cause mortality could be observed. Second, lower iCa levels () and tCa levels () were associated with an increased risk of 30-day, 90-day, and 365-day mortalities. Third, after adjustments for potential confounding factors, the quartile of the lowest iCa level () remained an independent predictor and was associated with an increase in all-cause mortality. To our knowledge, this study is the first to investigate the prognostic value of serum iCa and tCa levels among critically ill patients with CS.

A considerable number of clinical studies have suggested that the reduced serum calcium level was a common electrolyte disturbance among critically ill patients, which was also associated with increased mortality [41]. Our findings were consistent with the results of studies that evaluated the prognostic value of low serum calcium level in other clinical settings including CAD [14, 15, 18, 42], heart failure [43], AKI [20], CKD [44], trauma [45, 46], coronavirus disease 2019 (COVID-19) [47], or unselected emergency department admissions [48]. Lu et al. [15] reported that lower calcium levels were independent predictors for in-hospital mortality in patients with ST-elevation myocardial infarction (STEMI). Similarly, Yan et al. [14] showed that the baseline serum calcium added an incremental predictive value when combined with the Global Registry of Acute Coronary Events (GRACE) score in acute coronary symptom (ACS) patients. This study was the first to demonstrate that the low serum calcium was also associated with mortality in CS patients. In addition, although the most common cardiac cause of CS is ACS, CS can also result from nonischemic cardiac conditions, and few studies have attempted to explore predictors, which could be applicable to non-ACS presentations [8, 49]. In the subgroup analysis, we found that a lower level of iCa concentration () was a significant predictor of poor prognosis in CS caused by nonischemic cardiac conditions. Consequently, we hope the results of this study will supplement the findings of previous studies. Furthermore, decreased serum calcium levels might imply impaired kidney function [50]. In the present study, the adjustment for eGFR, or stratifying for CS patients according to the medical history of CKD, did not change the significant relationship between decreased serum calcium levels and increased risks of mortality. Thus, our findings showed that a lower serum calcium level might be an independent risk factor for the prognosis of CS rather than a surrogate marker of lower eGFR.

Although the exact mechanisms through which serum calcium leads to an elevated mortality rate remain unclear, there might be several possible explanations for this association. First, severe extracellular hypocalcemia could impact cardiac contractility because the sarcoplasmic reticulum is unable to maintain a sufficient amount of calcium content to initiate myocardial contraction [51]. Second, it has been assumed that the low calcium level might indicate an increased calcium consumption, partially reflecting more plaques or thrombi formed and worsening coronary conditions, resulting in poor outcomes through platelet activation [52]. Third, the appearance of low serum iCa was associated with secondary hyperparathyroidism and increased secretion of parathyroid hormone (PTH), which could promote calcium entry via L-type Ca2+ channels with consequent intracellular calcium overloading. Excessive cytosolic Ca2+ would affect the cardiac excitation-contraction coupling function, alter autophagic flux, and induce premature activation of intracellular enzymes, all of which contribute to the pathogenesis of CS [53].

Even in the era of reperfusion therapy, CS remains one of the leading causes of death with in-hospital mortality rates still approaching 50% [6, 54]. Individualized and timely risk assessment for each critically ill patient allows a more precise decision-making for therapeutic strategy and medical resource allocation. The prognostic value of several relatively convenient predictors including neutrophil percentage-to-albumin ratio [55], neutrophil-lymphocyte ratio [56], red blood cell distribution width [57], and low diastolic blood pressure [58] was explored. Similarly, even under conditions without imaging or additional laboratory tests, serum calcium could still serve as an effective marker for quick risk assessments. Our findings might provide additional convenience in some special situations, for example, underdeveloped areas. Moreover, further investigations are needed to explore the therapeutic value of serum calcium and find out whether calcium-supplementation therapy in CS patients with low serum calcium could improve their prognosis.

Several limitations of our study should be noted. First, we used data from a single academic medical center in the USA, with the earliest cases from almost 20 years ago, when care may have been inconsistent with currently accepted standards. The single-center nature of the study may also limit the applicability of our findings to other sites. Therefore, multicenter registry and prospective studies are needed to confirm these findings. Second, we measured serum iCa and tCa levels in patients only upon admission to the ICU and did not assess changes during their ICU stay, which might influence the summary results. Third, accurate calcium state determination depends on blood pH levels, because the binding of calcium to protein is particularly pH-sensitive. As pH decreases, H+ displaces Ca2+ from binding sites, and the amount of iCa increases. Conversely, as the blood pH increases, albumin and the globulins become more negatively charged and bind more calcium, causing the amount of iCa to decrease. Therefore, some sample collection practices (such as prolonged use of a tourniquet or the practice of having the patient clench or pump their fist) can artificially change the pH and cause an inaccurate iCa result, which might influence the results of our study. In addition, although every effort had been made to adjust for confounding factors using multivariate analysis, there remained other unknown factors that confused the prognostic value of serum iCa and tCa.

5. Conclusion

Lower serum iCa concentration was an independent predictor of all-cause mortality in critically ill patients with CS. Further studies, especially large prospective studies, are needed to confirm this relationship and validate its clinical significance.

Abbreviations

CS:Cardiogenic shock
AMI:Acute myocardial infarction
ICU:Intensive care unit
AKI:Acute kidney injury
tCa:Total calcium
iCa:Ionized calcium
MIMIC-III:Medical Information Mart for Intensive Care-III
BIDMC:Beth Israel Deaconess Medical Center
STROBE:STrengthening the Reporting of OBservational studies in Epidemiology
IRB:Institutional Review Boards
ICD:International Classification of Diseases
SQL:Structured query language
CAD:Coronary artery disease
CHF:Chronic heart failure
PAD:Peripheral artery disease
DM:Diabetes mellitus
CKD:Chronic kidney disease
SBP:Systolic blood pressure
DBP:Diastolic blood pressure
MBP:Mean blood pressure
HR:Heart rate
tCa:Total calcium
AG:Anion gap
eGFR:Estimated glomerular filtration rate
WBC:White blood cell
SOFA:Sequential Organ Failure Assessment
SAPS:Simplified Acute Physiology Score
GCS:Glasgow Coma Scale
RRT:Renal replacement treatment
SD:Standardized difference
IQR:Interquartile range
HR:Hazard ratio
CI:Confidence interval
MI:Multivariate imputation
COVID-19:Coronavirus disease 2019
STEMI:ST-elevation myocardial infarction
GRACE:Global Registry of Acute Coronary Events
ACS:Acute coronary symptom
PTH:Parathyroid hormone.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study utilized the anonymous data available in the MIMIC-III database with preexisting institutional review board (IRB) approval.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Authors’ Contributions

Y.Y., J.Y., and R.Y. wrote the original draft; Z.W., Y.Z., and J.X. reviewed and edited the manuscript; P.W. performed the supervision; Z.W. acquired the fund. All authors have read and agreed to the published version of the manuscript. Yue Yu, Jingwen Yu, Renqi Yao, and Pei Wang contributed equally to this work.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (No. 81770244), Medical Science and Technology Youth Cultivation Plan (No. 17QNP013 and No. 20QNPY038), Naval Military University Foreign Student Teaching Research and Reform Project (No. WJYA2018005), Shanghai Municipal Commission of Science and Technology (No. 17ZR1439100), Shanghai Shenkang Medicine Developing Project (No. SHDC12014107), and Shanghai Science and Technology Committee Medicine Leading Project (No. 15411960100).

Supplementary Materials

Table S1: univariable and multivariable Cox regression analysis for serum iCa levels and 30-day mortality. Table S2: univariable and multivariable Cox regression analysis for serum iCa levels and 90-day mortality. Table S3: univariable and multivariable Cox regression analysis for serum iCa levels and 365-day mortality. Table S4: univariable and multivariable Cox regression analysis for serum tCa levels and 30-day mortality. Table S5: univariable and multivariable Cox regression analysis for serum tCa levels and 90-day mortality. Table S6: univariable and multivariable Cox regression analysis for serum tCa levels and 365-day mortality. Table S7: the association between serum tCa levels and 30-day mortality in the subgroup analysis. Table S8: univariable and multivariable Cox regression analysis for serum iCa levels and 30-day mortality using original data. Table S9: univariable and multivariable Cox regression analysis for serum tCa levels and 30-day mortality using original data (Supplementary Materials)