Abstract

Background. Several studies have suggested that serum ionized calcium (iCa) is associated with mortality in critical illness. However, evidence regarding the predictive significance of serum total calcium (tCa) in critical illness remains scarce. The aim of this study was to assess the association of tCa levels with mortality in critical illness. Methods. We employed the MIMIC-III v1.3 database. tCa was measured upon ICU admission and its relationship with mortality was determined using smooth curve fitting. The association between admission tCa levels and hospital mortality was determined using logistic regression. Results. Inclusion criteria were met by 44,886 critically ill patients. A U-shaped pattern was observed between tCa and hospital mortality. Similar trends were observed for hospital mortality when quintiles were used to group patients according to tCa. In multivariate analysis, adjusted for age and sex, the model indicated that admission tCa levels ⩽7.6mg/dl, 7.7-8.1mg/dl, and 9.0mg/dl were associated with an increase in mortality when compared to the reference level (8.6-9.0mg/dl). However, adjusted for more clinical characteristics, tCa was not associated with hospital mortality. Conclusions. The relationship between tCa and hospital mortality followed a ‘‘U’’ shaped curve. tCa had certain prognostic value in critically ill patients, but it had no independent association with hospital mortality.

1. Introduction

Electrolyte disturbances, including hypokalemia and hypocalcemia, affect many fundamental physiologic regulatory mechanisms [1, 2]. Particular aberrations in calcium metabolism can lead to serious cardiovascular complications and organ dysfunctions [35]. Derangement of ionized calcium (iCa) is widespread in patients with critical illness [68]. Several previous studies have studied the prognostic value of iCa in critical illness [7, 9, 10]. Most of these studies demonstrated that hypocalcemia was associated with increased mortality in critical illness [1114], whereas there were few studies of hypercalcemia [15, 16].

Of total body calcium distribution, 98% is found in the skeleton, and the remaining 2% is found in the circulation. There are two main physicochemical states of serum total calcium (tCa). Approximately 50% of tCa is iCa, and the remainder is bound to albumin [17]. Thus, there is a close relationship between the tCa and iCa. Roland’s study found that critically ill patients with lower tCa concentration have a high rate of lower iCa [18]. However, we do not know whether derangement in tCa was associated with mortality in critical illness.

Therefore, we extracted and analyzed the Multiparameter Intelligent Monitoring in Intensive Care Database III version 1.3 (MIMIC-III v1.3) to assess the impact of tCa on mortality in critical illness.

2. Methods

2.1. Data Source

MIMIC-III v1.3 comprises more than 40,000 intensive care unit (ICU) patients admitted to Beth Israel Deaconess Medical Center (Boston, MA, USA) from 2001 to 2012 [19]. To apply for access to the database, we completed the National Institutes of Health’s web-based course and passed the Protecting Human Research Participants exam (NO. 6182750). We extracted clinical data, including patient demographics and laboratory test results. The establishment of this database was subject to approval by the institutional review boards of Massachusetts Institute of Technology and Beth Israel Deaconess Medical Center. To safeguard patient privacy, they were deidentified.

2.2. Population Selection Criteria

The MIMIC-III database recorded 58,976 ICU admissions. Adult patients (18 years) were eligible. Patients were excluded if there was no tCa measured during ICU stay or if there were more than 10% data missing.

2.3. Data Extraction

Structured Query Language (SQL) with PostgreSQL (version 9.6) was used to extract data from MIMIC-III. The data extracted were clinical parameters, age, sex, temperature, respiratory rate, heart rate, blood pressure, blood oxygen saturation, comorbidities, laboratory parameters, and scoring systems. We extracted the following comorbidities: chronic heart failure, coronary artery disease, stroke, atrial fibrillation, chronic kidney diseases (CKD), acute kidney injury (AKI), pneumonia, liver disease, chronic obstructive pulmonary disease (COPD), malignancy, chronic pancreatitis, acute pancreatitis, and acute respiratory distress syndrome (ARDS). The following laboratory parameters were extracted: bicarbonate, creatinine, chloride, blood urea nitrogen (BUN), glucose, hematocrit, hemoglobin, platelet, sodium, and white blood cell (WBC). We also extracted sequential organ failure assessment (SOFA) and simplified acute physiology score II (SAPS II). Hospital mortality was the primary endpoint. tCa was measured on ICU admission. The baseline characteristic data were extracted within the first 24 h after patient ICU admission.

2.4. Statistical Analysis

Continuous variables are presented as mean ± standard deviation (SD) and compared using the analysis of variance or the Kruskal-Wallis test. Categorical variables were expressed as percentage and were compared using the chi-squared test.

The association between the tCa and hospital mortality was determined using a multivariate logistic regression and was expressed as the adjusted odds ratio with pertinent 95% confidence interval (CI). Two multivariate models were constructed on the basis of tCa according to quintiles derived with curve-fitting methods based on hospital mortality. The fourth quintile was considered as the control group. In model I, covariates were adjusted only for age and sex. In model II, covariates were adjusted for age, sex, oxygen saturation (SPO2), SOFA, systemic inflammatory response syndrome (SIRS), overall anxiety severity and impairment scale (OASIS), simplified acute physiology score II (SAPSII), systolic blood pressure (SBP), diastolic blood pressure (DBP), heart rate, respiratory rate, chronic heart failure, AKI, liver disease, acute pancreatitis, renal replacement therapy (RRT), anion gap, creatinine, chloride, glucose, and hemoglobin.

We performed stratification analyses to determine whether the effect of tCa differed across various subgroups. AKI, CKD, chronic heart failure, pneumonia, liver disease, acute pancreatitis, acute respiratory distress syndrome (ARDS), and RRT were also covered. All statistical analyses were performed using the EmpowerStats version 2.17.8 (http://www.empowerstats.com/cn/) and R software version 3.42. Two-tailed P value <.05 was considered to be statistically significant.

3. Results

3.1. Subject Characteristics

A total of 44,886 patients met inclusion criteria. Quintiles were used to group patients according to the tCa: 7,701 patients were in the ⩽7.6 mg/dl group, 8,382 patients were in the 7.7- 8.1 mg/dl group, 8,824 patients were in the 8.2- 8.5 mg/dl group, 10,219 patients were in the 8.6-9.0 mg/dl group, and 9,760 patients were in the 9.0 mg/dl group. Characteristics and hematologic laboratory data of the study patients according to admission tCa levels are displayed in Table 1.

Characteristics including age, sex, temperature, and SPO2 were relatively flat across each group. Blood pressure was lowest in severe hypocalcemia (⩽7.6 mg/dl), but heart rate and respiratory rate were highest. Patients with severe hypocalcemia were more likely to report a history of AKI, pneumonia, malignancy, chronic pancreatitis, and acute pancreatitis. Moreover, the proportion of comorbidities, such as stroke, CKD, and chronic obstructive pulmonary disease (COPD), were higher in those with severe hypercalcemia (9.0 mg/dl). Participants with severe hypercalcemia also had higher levels of bicarbonate, creatinine, glucose, hematocrit, hemoglobin, platelets, sodium, and BUN, than did other groups. In addition, Elixhauser Comorbidity Index, SOFA, and SAPSII scores were significantly higher in severe hypocalcemia than in other groups.

3.2. tCa Levels and Hospital Mortality

We divided initial tCa into variable categories. Figure 1 shows the relationship between tCa and logit-transformed hospital mortality. The relationships were nonlinear. We observed a U-shaped relationship, suggesting that both hyper- and hypocalcemia were associated with equally increased risk of hospital mortality.

In multivariate analysis, adjusted for clinical characteristics, including age and sex, the adjusted ORs (95% CIs) for admission tCa levels ⩽7.6mg/dl, 7.7-8.1mg/dl, and 9.0mg/dl compared to the reference level (8.6-9.0mg/dl) were 1.84 (1.69, 2.00), 1.23 (1.13, 1.35), and 1.12 (1.02, 1.22), respectively. These were linked to an increase in mortality. However, after adjustment for additional clinical characteristics including age, sex, SPO2, SOFA, SIRS, OASIS, SAPSII, SBP, DBP, heart rate, respiratory rate, chronic heart failure, AKI, liver disease, acute pancreatitis, renal replacement therapy, anion gap, creatinine, chloride, glucose, and hemoglobin, admission tCa levels were not independently associated with hospital mortality (Figure 2, Table 2).

3.3. Subgroup Analyses

The association between tCa and the risk of hospital mortality was similar for most strata (Table 3). Patients with AKI, CKD, chronic heart failure, pneumonia, liver disease, acute pancreatitis, and RRT had a significantly higher risk of hospital mortality for hypocalcemia, especially in severe hypocalcemia (⩽7.6 mg/dl). Patients with these diseases showed no difference in the risk of mortality for hypercalcemia. Neither hyper- nor hypocalcemia was linked with mortality in patients with ARDS. No matter what kind of the acid-base state, severe hypocalcemia increased the risk of hospital mortality. Hypercalcemia was associated with mortality only in acidosis. Although several tCa levels were related to the mortality in hypoalbuminemia, albumin had little effect on tCa in prognostic value of mortality.

4. Discussion

We demonstrated a U-shaped relationship between initial tCa and hospital mortality. Adjusted for clinical characteristics including age and sex, both hyper- and hypocalcemia were linked with increased hospital mortality in critically ill patients. When adjusted for several clinical characteristic, tCa was not associated with hospital mortality. Although many previous studies showed that iCa was associated with mortality in critical illness [2023], the relationship between tCa and mortality had not been studied. To our knowledge, this was the first study to establish an association between initial tCa and mortality in critically ill patients.

Numerous risk factors appear to contribute to poor outcomes in critical illness. Although tCa was not an independent risk factor for mortality after adjustment for various clinical characteristics, it has a certain prognostic implication in critically ill patients. As observed in present study, a U-shaped relationship and a model that has been calibrated for age and gender indicated that tCa <8.2mg/l or >9.0mg/l were associated with hospital mortality. The greater the deviation from the normal value, the higher the death rate. In subgroup analysis, we showed that there was a correlation between tCa and hospital mortality for most critical illness disease categories. Thus, on ICU admission, the initial tCa level assessment could serve as a preliminary prognostic marker for mortality. However, we can not yet establish a predictive model of serum calcium for mortality. The precise prognostic value of tCa for mortality requires further investigation.

Calcium homeostasis is regulated by the vitamin D-parathyroid-calcium axis. Dysfunction of parathyroid hormone is highly prevalent in critical illness [24], and hypocalcemia is common in critically ill patients [25]. Potential mechanisms include proinflammatory cytokines impairing parathyroid hormone, hypersecretion of catecholamine, and redistribution of serum calcium [26, 27]. Nair’s study [24] showed that vitamin D insufficiency or deficiency was prevalent among ICU patients. This prevalent vitamin D insufficiency may explain the high incidence of hypocalcemia. Severe hypercalcemia is associated with increase in mortality risk, attributable to either primary hyperparathyroidism or malignancy under most circumstances [28]. In addition, severe hypercalcemia is a risk factor for AKI in critically ill patients, and the occurrence of AKI would result in increased mortality [29, 30].

The major strength of this study is that we determine the association of tCa and hospital mortality for the first time. Mortality as primary outcome is clearly important in critically ill patients. In addition, the extracted data are from a large clinical database (MIMIC-III) and the number of patients enrolled is very large.

This study has a few limitations. First, the study has a single-center retrospective design. It is therefore subject to selection bias. Second, tCa was not an independent risk factor for mortality, but critically ill patients are heterogeneous. Whether increasing the specificity of the study will improve the prognostic value of tCa requires further investigation. Third, individual data are missing and outliers are present, possibly influencing the summary results. Finally, the follow-up length of mortality varied. Generally, 30-day mortality is the leading choice. We use hospital mortality, and this may result in too many instances of censored data.

5. Conclusions

In a large, heterogeneous group of critically ill patients, we found that the relationship between tCa and hospital mortality in critically ill patients followed a ‘‘U’’ shaped curve. tCa had a certain prognostic implication in critically ill patients but had no independent association with hospital mortality.

Data Availability

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

Disclosure

The funders of the project were not involved in study design, collection, data analysis, writing of the report, and publication.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Benji Wang and Bihuan Cheng designed the experiment. Yuqiang Gong and Binyu Ying extracted data. Benji Wang and Yuqiang Gong performed all statistical analyses. Benji Wang and Bihuan Cheng drafted the paper. All authors revised and approved the final manuscript.

Acknowledgments

This research was supported by Zhejiang Provincial Natural Science Foundation of China under Grant no. LY14H150011.