Time series analysis Between Serum Calcium Levels and 60-Day Mortality in Critically ill Patients with Multiple Myeloma: A Retrospective Cohort Study | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Time series analysis Between Serum Calcium Levels and 60-Day Mortality in Critically ill Patients with Multiple Myeloma: A Retrospective Cohort Study Yan Zeng, Xiaomin Cheng, Yongran Si, Jun Wang, Jingwei Zhang This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8812562/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background Hypercalcemia is a well-recognized adverse prognostic factor in patients with multiple myeloma (MM), particularly in critically ill subgroups where its incidence is substantially higher than in the general MM population. However, the prognostic value of dynamic serum calcium fluctuations during intensive care unit (ICU) hospitalization remains underexplored, limiting the ability to refine risk stratification and personalized management. Objective This study aims to identify heterogeneous serum calcium trajectories in critically ill MM patients and investigate their association with 60-day mortality. Methods A retrospective cohort study was conducted using data from 254 critically ill MM patients extracted from the MIMIC-IV database. Latent Class Growth Modeling (LCGM) was applied to classify serum calcium trajectories over the first 28 days of hospitalization, with measurements collected at predefined time points-. Multivariable Cox proportional hazards regression models (adjusted for demographic, clinical, and disease severity covariates) and Kaplan-Meier (K-M) curves were used to assess the association between trajectory patterns and 60-day mortality. Sensitivity analyses included trajectory reclassification and E-value assessment to validate result robustness. Results Four distinct serum calcium trajectories were identified: low-level (n = 90, 35.43%), medium-level (n = 138, 57.48%), U-shaped (n = 9, 3.54%), and high-level (n = 17, 6.68%). Compared with the high-level trajectory (reference group), the low-level (HR = 0.17; 95% CI 0.07–0.41; P < 0.001) and medium-level (HR = 0.26; 95% CI 0.12–0.56; P = 0.001) trajectories were independently associated with significantly lower 60-day mortality. The U-shaped trajectory showed no significant mortality difference from the high-level trajectory (HR = 0.56; 95% CI 0.14–2.25; P = 0.416). K-M analysis confirmed the prognostic gradient across trajectories, and sensitivity analysis (E-value = 2.49) validated the robustness of these associations to unmeasured confounding. Conclusion Persistently high or rebound (U-shaped) serum calcium trajectories are independent predictors of increased 60-day mortality in critically ill MM patients, reflecting underlying aggressive disease biology or treatment resistance. Dynamic serum calcium monitoring via trajectory analysis refines early risk stratification and provides actionable insights for guiding personalized ICU management, addressing an unmet need in current clinical practice. Multiple myeloma (MM) serum calcium mortality Latent Class Growth Modeling (LCGM) Trajectory Figures Figure 1 Figure 2 Introduction Multiple myeloma (MM) is a hematologic malignancy characterized by the clonal proliferation of plasma cells within the bone marrow. Globally, MM incidence is on the rise, making it the second most common hematologic cancer after Hodgkin's lymphoma. In Europe, the annual incidence rate ranges from 4.5 to 6.0 cases per 100,000 individuals, with a corresponding mortality rate of 4.1 per 100,000 [1]. Complications associated with MM include anemia, hypercalcemia, renal dysfunction, and bone disease. Hypercalcemia is observed in 19.5% of MM patients at diagnosis, with severity ranging from mild to life-threatening, and is frequently linked to bone lesions, kidney dysfunction, and high-risk chromosomal abnormalities in MM [2, 3]. Notably, emerging evidence suggests that critically ill MM patients admitted to the intensive care unit (ICU) exhibit a significantly higher incidence of hypercalcemia (up to 30%-40%) compared to the general MM population. In such patients, hypercalcemia is not only a consequence of bone resorption driven by myeloma cell-derived cytokines, but also intertwined with ICU-specific insults such as acute kidney injury (AKI) and sepsis, forming a vicious cycle that accelerates organ dysfunction and poor outcomes. This unique pathophysiological feature highlights the need to explore dynamic calcium changes rather than static baseline values in this population[4]. Hypercalcemia at initial diagnosis of MM is widely recognized as an adverse prognostic factor [2, 3]. Given the complexity of critical illness and the variability of medical interventions (e.g., anti-myeloma therapy that may alter bone metabolism, renal replacement therapy that regulates calcium excretion), serum calcium levels change dynamically over time. The persistence or exacerbation of hypercalcemia in critically ill MM patients may further increase mortality risk. Guidelines, including the EHA-ESMO clinical practice guidelines for MM, recommend monthly or at least quarterly monitoring of core laboratory parameters (complete blood count, creatinine, and serum calcium) during follow-up, but lack guidance on interpreting dynamic fluctuations, leaving clinicians unable to leverage serial calcium data for real-time risk assessment in ICU settings[5, 6]. A recent study further demonstrated a nonlinear relationship between serum calcium levels and in-hospital mortality in critically ill MM patients, with a serum calcium level of approximately 8.40 mg/dL associated with the lowest mortality risk and risk increasing with rising levels, an important consideration for ICU clinicians[6]. Despite this, most current studies focus on baseline serum calcium data from newly diagnosed MM patients at admission [2, 6], failing to account for the dynamic fluctuations of serum calcium levels during ICU stay. Given the complexity of critical illness and the variability of medical interventions, serum calcium levels change dynamically over time, and investigating these time-series changes and their prognostic implications is crucial for optimizing risk stratification and treatment adjustment [7]. Medical time series analysis (TSA) has been widely applied in prognostic studies [8]. Among TSA methods, the Latent Class Growth Modeling (LCGM) is particularly valuable for identifying heterogeneous trajectory subgroups in heterogeneous populations such as critically ill MM patients who may exhibit diverse calcium change patterns, by partitioning the population into latent classes based on individual characteristics (e.g., intercepts and slopes) and modeling class-specific change patterns over time[9]. Trajectory analysis using LCGM has been successfully applied to characterize dynamic biomarkers in various hematologic and critical diseases[10], but to date, no studies have utilized this approach to explore the relationship between dynamic serum calcium trajectories and mortality in MM patients, especially in the high-risk critically ill subgroup. To address this research gap, the present study leverages data from the MIMIC-IV database to construct LCGM of serum calcium trajectories during the first 28 days of hospitalization, integrating both baseline measurements and dynamic time-series data. The primary objectives are to identify distinct serum calcium trajectories in critically ill MM patients, evaluate their association with 60-day mortality, and provide evidence for refining early risk stratification and guiding personalized clinical decisions. Materials and Methods Population and Data Collection This retrospective cohort study utilized the Medical Information Mart for Intensive Care-IV (MIMIC-IV v3.0) database, comprising de-identified clinical data from >70,000 intensive care unit (ICU) admissions at Beth Israel Deaconess Medical Center (Boston, MA) between 2008 and 2019[6]. Access to the MIMIC-IV database was obtained in accordance with the official access protocol specified by PhysioNet (https://physionet.org/content/mimiciv/3.0/). Access to the database required completion of the "CITI Data or Specimens Only Research" training course on the National Institutes of Health (NIH) website, with approval granted for data extraction for research purposes (Certificate No for Zhang: ID: 13813672). Ethical oversight was provided by institutional review boards of Massachusetts Institute Technology and Beth Israel Deaconess Medical Center. Informed consent was waived given the retrospective, de-identified nature of the data[6]. Structured Query Language (SQL), PostgreSQL (version 13.0), and Navicat software (version 16.0) were used to identify the cohort and extract relevant clinical data [11]. Inclusion criteria were defined as: (1) adult patients aged ≥18 years; (2) a confirmed diagnosis of MM based on International Classification of Diseases (ICD) coding: ICD-9 codes (203.00, 203.01, 203.02) or ICD-10 codes (C90.00, C90.01, C90.02)[6]; (3) first-time ICU admission (to avoid duplicate patient data); and (4) availability of at least two serum calcium measurements during hospitalization. Exclusion criteria were: (1) absence of baseline serum calcium measurement on ICU admission; (2) fewer than two serum calcium measurements throughout the ICU hospitalization period; (3) incomplete data on the primary outcome (60-day mortality). Demographics, comorbidities, and laboratory indicators were extracted from the MIMIC-IV database. Serum calcium levels were systematically collected at predefined time points: ICU days 1, 3, 5, and 7, with weekly averages calculated for weeks 2, 3, and 4 of hospitalization. Additionally, the minimum, maximum, and mean serum calcium values during the entire ICU stay were computed. For missing serum calcium data on days 3, 5, or 7, a last-observation-carried-forward approach was applied, wherein the nearest adjacent available measurement was imputed. For other laboratory parameters, maximum, minimum, or mean values were selected based on clinical relevance and prior literature[10, 12]. The primary outcome of interest was 60-day mortality. Variables with missing data exceeding 10% were excluded from the analysis, consistent with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline[13]. For remaining variables with minor missingness (<10%), multiple imputation was performed using chained equations with 5 imputed datasets via the mice package in R [14]. Identification of Serum Calcium Trajectory Subgroups and Model Selection Heterogeneous serum calcium trajectories were identified using the R package "lcmm" (version 2.0.2) via LCGM [7]. A latent class mixed model was applied for longitudinal serum calcium data analysis, and a parametric survival model with a class-specific Weibull distribution was incorporated to characterize the baseline hazard function, enabling simultaneous modeling of longitudinal trajectories and time-to-60-day-mortality outcomes [7].Trajectory modeling was based on serum calcium data systematically collected within the initial 28 days after ICU admission to be consistent with the study’s predefined observation window, including measurements on ICU days 1, 3, 5, and 7, as well as weekly averages for weeks 2-4. Both linear and quadratic functional forms for the time variable were evaluated to optimize model fit. Quadratic terms were included to capture potential acceleration or deceleration in serum calcium changes commonly observed in critically ill patients[15]. Model selection was guided by a comprehensive set of criteria to ensure statistical robustness and clinical interpretability:Lower values of the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), and sample-adjusted BIC (SABIC) indicated superior model fit; A significant P-value (< 0.05) from the Lo-Mendell–Rubin (LMR) adjusted likelihood ratio test supported adding an additional trajectory class; Visual inspection of trajectory plots was conducted to confirm the clinical plausibility of identified patterns; Higher maximum log-likelihood and entropy (≥ 0.80, indicating strong classification accuracy) were prioritized; Robust subgroup classification required an average posterior probability ≥ 70% for each class; Only models with each subgroup containing > 3% of the total study population were retained to avoid unstable estimates from small subgroups[10, 16, 17].The optimal number of trajectory subgroups was determined by balancing all aforementioned criteria, with the final model selected based on the best overall fit and clinical interpretability. Statistical analysis Continuous variables were described as mean ± standard deviation (SD) for normally distributed data (assessed via the Shapiro-Wilk test) and median (interquartile range, IQR) for skewed data. Group comparisons across latent trajectory classes were performed using the Kruskal-Wallis test for continuous variables. Categorical variables were presented as counts (percentages, %), with group comparisons conducted using Pearson’s chi-squared test; Fisher’s exact test was used when the expected frequency of any cell was < 5. Univariate and multivariate Cox proportional hazards regression analyses were performed to evaluate the association between candidate variables (including serum calcium trajectories) and 60-day mortality. The selection of covariates for adjustment in multivariate models was guided by three principles: (1) clinical relevance based on prior literature[18], (2) expert consensus on risk factors for mortality in critically ill MM patients[6], and (3) statistical significance (p < 0.05) in univariate Cox regression. Four hierarchical multivariate Cox models were constructed to assess the robustness of the association between serum calcium trajectories and 60-day mortality. Model 1: Unadjusted (only serum calcium trajectory as the predictor);Model 2: Adjusted for age and gender (basic demographic confounders);Model 3: Adjusted for age, gender, mean systolic blood pressure (SBP_mean), maximum blood urea nitrogen (Bun_max), minimum platelet count (Plt_min), and congestive heart failure (clinical and laboratory confounders);Model 4: Fully adjusted for age, gender, SBP_mean, Bun_max, Plt_min, congestive heart failure, Acute Physiology Score III (APSIII), Simplified Acute Physiology Score II (SAPSII), Oxford Acute Severity of Illness Score (OASIS), and Sequential Organ Failure Assessment (SOFA) score (adding disease severity indices). Kaplan–Meier (KM) survival curves were generated to visualize survival probabilities across serum calcium trajectory classes, with between-group differences evaluated using the log-rank test. The E-value calculation Quantifies the potential impact of unmeasured confounding on the observed associations. The E-value was defined as the minimum strength of association (on the risk ratio scale) that an unmeasured confounder would need to have with both the exposure (serum calcium trajectory) and the outcome (60-day mortality), conditional on measured covariates, to fully explain the observed association [19]. All statistical analyses were conducted using R Statistical Software (version 4.2.2; http://www.R-project.org, The R Foundation) and the Free Statistics Analysis Platform (version 2.1.1; Beijing, China). Specific R packages included survival for Cox regression and KM analysis, and car for testing model assumptions. A two-tailed p-value < 0.05 was considered statistically significant[20]. Results Baseline characteristics A total of 254 patients were included in this study after screening according to predefined inclusion and exclusion criteria ( Supplementary Figure 1 ). The overall 60-day mortality rate was 24.02% (61/254). The baseline characteristics of the study population stratified by 60-day mortality status are summarized in Table 1 . Compared with the surviving group (n = 193), the non-surviving group (n = 61) had a significantly older age and a higher prevalence of congestive heart failure. Regarding vital signs, the non-surviving group exhibited a lower mean diastolic blood pressure and a higher mean respiratory rate. For laboratory parameters, the non-surviving group had a significantly lower minimum platelet count and a higher maximum blood urea nitrogen level. In terms of illness severity scores, the non-surviving group had significantly higher scores on the APSIII, SAPSII, OASIS, and SOFA score. Identification of subpopulations using the LCGM The optimal number of serum calcium trajectory subgroups was determined based on multiple model fit metrics, including log-likelihood, AIC, BIC, SABIC, entropy, and class composition ratio ( Table 2 ). As the number of classes increased from 1 to 5, the log-likelihood, AIC, and SABIC values initially decreased and then increased, with the 4-class model yielding the lowest AIC (2262.018) and SABIC (2271.196) values. Additionally, the 4-class model exhibited a moderate log-likelihood (-1106.008962), a relatively high entropy (0.683), and each class accounted for >3% of the total population (ranging from 3.544% to 57.480%), meeting the criteria for robust subgroup classification.Thus, the 4-class model was selected as the optimal solution, and the trajectories were labeled as follows ( Figure 1 ). Class 1(Low-level trajectory, 35.433%, n = 90) of patients, characterized by a slight increase from an initially low serum calcium level; Class 2 (Medium-level trajectory, 57.480%, n = 138) of patients, characterized by a slight decrease from a moderate baseline serum calcium level; Class 3 (U-shaped trajectory, 3.544%, n = 9) of patients, characterized by serum calcium levels decreasing from an initial high level to a low point and then rising back to a medium level; Class 4 (High-level trajectory, 6.683%, n = 17) of patients, characterized by persistently high serum calcium levels maintained throughout the observation period. The baseline characteristics of the study population stratified by serum calcium trajectories are summarized in Supplementary Table 1 . Univariable Cox Regression Analysis As shown in Table 3 , univariable Cox regression analysis was performed with the high-level trajectory (Class 4) as the reference group. The results revealed that patients in the low-level trajectory (Class 1) and medium-level trajectory (Class 2) had significantly lower 60-day mortality risk (Class 1: HR = 0.28, 95% CI = 0.13-0.63, p = 0.002; Class 2: HR = 0.34, 95% CI = 0.16-0.71, p = 0.004). In contrast, no significant difference in 60-day mortality was observed between the U-shaped trajectory (Class 3) and the reference group (HR = 0.51, 95% CI = 0.14-1.89, p = 0.313).Additionally, univariable analysis identified several independent risk factors for 60-day mortality, including older age, lower DBP-mean, higher RR-mean, lower PLT-min, higher BUN-max, presence of congestive heart failure, and higher APSIII, SAPSII, OASIS, and SOFA scores (all p < 0.001). Multivariable Cox Regression Analysis To adjust for potential confounding factors, four hierarchical multivariable Cox regression models were constructed ( Table 4 ). In the fully adjusted model (Model 4), which included age, gender, SBP-mean, BUN-max, PLT-min, congestive heart failure, APSIII, SAPSII, OASIS, and SOFA score, the associations between serum calcium trajectories and 60-day mortality remained robust. Compared with the high-level trajectory (Class 4), the low-level trajectory (Class 1) and medium-level trajectory (Class 2) were still significantly associated with lower 60-day mortality risk (Class 1: HR = 0.17, 95% CI = 0.07-0.41, p < 0.001; Class 2: HR = 0.26, 95% CI = 0.12-0.56, p = 0.001). The U-shaped trajectory (Class 3) still showed no significant difference in 60-day mortality compared with the reference group (HR = 0.56, 95% CI = 0.14-2.25, p = 0.416). Kaplan-Meier curves analysis Kaplan-Meier (KM) survival curves stratified by serum calcium trajectory subgroups are presented in Figure 2 . The results showed that the high-level trajectory (Class 4) had the lowest 60-day survival probability, while the low-level trajectory (Class 1) had the highest, followed by the medium-level trajectory (Class 2). Although the survival curve of the U-shaped trajectory (Class 3) intersected with those of Class 1 and Class 2 around weeks 2–3, the log-rank test indicated significant difference in survival probabilities among the four subgroups (p = 0.0082). Sensitivity Analysis Given the lack of significant difference in 60-day survival rates between Class 3 and Class 4, we compared serum calcium levels across the trajectory classes at various time points. As shown in Table 5 , Class 1 consistently exhibited lower serum calcium levels on days 1, 3, 5, and 7, as well as lower mean, minimum, and maximum values during hospitalization compared to the other three classes. The initial values for Class 2 were lower than those for Class 3 during the first three days, became comparable by day 5, and subsequently dropped below Class 3 by week three. Class 3 initially had higher values than Class 4 over the first three days, followed by a gradual decline and subsequent rise, while remaining consistently lower than Class 4. When stratifying the population by 60-day mortality and comparing serum calcium levels between survivors and non-survivors at multiple time points, no significant differences were observed (all p > 0.05) ( Supplementary Table 2 ). These findings suggest that Class 3 and Class 4 share similar prognostic characteristics attributable to their consistently elevated calcium levels. Based on the similar prognostic characteristics between the U-shaped trajectory (Class 3) and the high-level trajectory (Class 4) identified in the primary analysis, we merged these two classes into "Group 3" and reclassified the low-level trajectory (Class 1) and medium-level trajectory (Class 2) as "Group 1" and "Group 2", respectively. Multivariable Cox regression analysis with Group 3 as the reference group ( Supplementary Table 3 ) showed that after adjusting for the same covariates, Group 1 (HR = 0.21, 95% CI = 0.10–0.45, p < 0.001) and Group 2 (HR = 0.30, 95% CI = 0.15–0.61, p = 0.001) still had significantly lower 60-day mortality risk. KM survival curves for the three reclassified groups confirmed that Group 3 had the lowest 60-day survival probability ( Supplementary Figure 2 ), with a statistically significant difference among the three groups (overall p = 0.0091). To evaluate the robustness of the observed associations to unmeasured confounding, an E-value analysis was performed. The E-value for the point estimate of the association between serum calcium trajectories and 60-day mortality was 2.494 (95% CI: 1.492–4.051). This relatively large E-value indicates that an unmeasured confounder would need to have a strong association with both serum calcium trajectory and 60-day mortality (with a risk ratio of at least 2.494) to fully explain the observed association, suggesting that the study results are robust to potential unmeasured confounding. Discussion To the best of our knowledge, this is the first study to identify heterogeneous serum calcium trajectories in critically ill MM patients using Latent Class Growth Modeling (LCGM) and to systematically explore the association between these dynamic trajectories and 60-day mortality. This work addresses a critical gap in the existing literature, which has predominantly focused on static baseline serum calcium levels in newly diagnosed or stable MM patients[ 6 ] [ 21 ], rather than capturing the dynamic fluctuations inherent to critical illness, characterized by complex interactions between myeloma-related bone resorption, acute kidney injury, sepsis, and intensive care interventions[ 4 , 22 ]. Our findings indicate that the overall 60-day mortality rate among critically ill MM patients was 24.02%, which is lower than the 30-day mortality rate (35.2%) reported by Hampshire et al. [ 23 ] in a cohort of hematologic malignancy patients admitted to UK ICUs, but slightly higher than the in-hospital mortality rate (19.4%) reported by Mao et al. in critically ill MM patients[ 6 ]. This discrepancy is likely multifactorial: first, the choice of endpoint (60-day mortality vs. 30-day/in-hospital mortality) captures a longer follow-up period, including post-discharge deaths that may be associated with unresolved organ dysfunction or disease progression; second, variations in cohort severity; third, differences in healthcare settings and treatment patterns (e.g., utilization of novel anti-myeloma agents or renal replacement therapy) between study populations may further explain the variance[ 21 ].A key novel finding of this study is the identification of four distinct serum calcium trajectories over the first 28 days of ICU hospitalization, with clear prognostic stratification. Specifically, patients with low-level (Class 1) and medium-level (Class 2) trajectories exhibited significantly lower 60-day mortality (HR = 0.17 and 0.26, respectively, in the fully adjusted model) compared to those with high-level (Class 4) trajectories. Notably, the U-shaped trajectory (Class 3), characterized by an initial high calcium level, transient decline, and subsequent rise, showed no significant mortality difference from the high-level trajectory (HR = 0.56, p = 0.416). This observation is crucial. It suggests that even a transient normalization of serum calcium does not mitigate the adverse prognosis portended by the occurrence of hypercalcemia in the critical care setting[ 21 ] [ 24 ], but extends these findings by demonstrating that dynamic trajectory patterns, rather than single-point measurements, provide a more robust prognostic assessment. Elevated serum calcium levels are inversely correlated with survival in MM patients [ 25 ], a relationship deeply rooted in the disease’s biological aggressiveness and organ-damaging sequelae. This correlation is primarily driven by the intimate link between hypercalcemia and skeletal-related events (SREs)—including pathological fractures, extensive osteolytic lesions, renal impairment, and increased infection risk-all of which independently worsen prognosis [ 25 , 26 ]. Mechanistically, MM-related hypercalcemia is predominantly mediated by osteoclast hyperactivation, a process amplified by myeloma cell-derived cytokines such as IL-6 and TNF-α that disrupt the bone remodeling balance, leading to unregulated bone resorption and calcium release into the circulation [ 24 , 27 ]. Notably, this cytokine-driven bone destruction is not only a cause of hypercalcemia but also a hallmark of high-risk MM biology, as these cytokines simultaneously promote myeloma cell proliferation and immune suppression[ 27 ]. Although the full spectrum of mechanisms by which hypercalcemia impairs survival remains incompletely elucidated, accumulating evidence highlights its association with adverse biological features. The presence of hypercalcemia at diagnosis is an established independent poor prognostic factor, with nearly 20% of newly diagnosed symptomatic MM patients presenting with this complication[ 28 ]. Critically, hypercalcemia in MM is strongly linked to high-risk cytogenetic abnormalities—including del 13q (detected by fluorescence in situ hybridization) and amp1q21, which are themselves associated with rapid disease progression, treatment resistance, and shortened overall survival [ 3 ]. The prognostic significance of hypercalcemia persists across age groups and treatment modalities, including conventional chemotherapy and novel agents such as thalidomide, lenalidomide, and bortezomib, suggesting that hypercalcemia is not merely a clinical manifestation but a surrogate marker of inherently aggressive disease biology and potential drug tolerance[ 3 ]. This biological context provides critical insights into the four serum calcium trajectories identified in our study. The high-level trajectory (Class 4), characterized by persistently elevated calcium levels—likely reflects uncontrolled myeloma activity, driven by high-risk cytogenetic profiles and sustained osteoclast-mediated bone destruction[ 3 ] [ 24 ]. These patients may harbor clones with inherent resistance to anti-myeloma therapies, as evidenced by the failure to normalize calcium levels despite intensive care interventions. Similarly, the U-shaped trajectory (Class 3), with initial hypercalcemia, transient reduction, and subsequent rebound, suggests a temporary response to acute interventions (e.g., hydration, bisphosphonates, or short-course anti-myeloma therapy) followed by disease progression or acquired drug resistance [30]. The recurrence of hypercalcemia in this trajectory aligns with the notion that cytogenetically driven bone destruction is not fully reversible with standard therapies, leading to persistent organ damage and poor outcomes comparable to those with sustained hypercalcemia[ 3 , 22 ].In contrast, the low-level (Class 1) and medium-level (Class 2) trajectories, characterized by stable calcium levels within or below the normal range, are indicative of better disease control and treatment sensitivity. These trajectories likely represent patients with non-aggressive disease biology (e.g., absence of high-risk cytogenetics), effective suppression of myeloma-driven cytokine production, and responsive bone remodeling following treatment. The stability of calcium levels in these groups suggests durable control of both myeloma cell proliferation and osteoclast activity. Notably, the segmented linear distribution of average calcium levels across trajectories (Class 1 < Class 2 < Class 3 < Class 4) mirrors the gradient of prognostic risk, consistent with Mao et al.’s observation of a nonlinear relationship between calcium levels and mortality [ 6 ]. This continuity further validates that serum calcium trajectories are not arbitrary groupings but reflect a biological spectrum of disease aggressiveness, from well-controlled to refractory disease. Collectively, our findings highlight the direct clinical applicability of dynamic serum calcium trajectory monitoring. In the ICU setting, these trajectories can function as an early-warning biomarker, prompting clinicians to intensify management of myeloma-specific complications, such as escalating bone-targeted therapy or reassessing anti-myeloma regimens, when a persistently high or rebounding pattern is observed. Incorporation of trajectory analysis with validated severity scores may enhance risk stratification and support timely clinical decision-making. This study has several notable strengths. First, it innovatively applies LCGM to characterize serum calcium trajectories in critically ill MM patients, shifting the focus from static to dynamic biomarker assessment, a paradigm increasingly recognized in critical care research. Second, the use of hierarchical multivariable Cox regression and E-value analysis (E = 2.494) enhances the robustness of the observed associations against confounding factors. Third, the identification of prognostic trajectories provides a practical tool for risk stratification, addressing the unmet need for dynamic monitoring guidance in critically ill MM management. Despite these insights, our study has several limitations that warrant consideration. First, the reliance on a single-center public database (MIMIC-IV) precluded access to detailed cytogenetic data (e.g., del 13q, amp1q21), bone marrow biopsy results, and long-term treatment response data—critical variables needed to directly validate the link between trajectories, high-risk biology, and treatment resistance. Future prospective studies should integrate these variables to explicitly test whether trajectory patterns correlate with specific genetic signatures and treatment outcomes. Second, the small sample size of the U-shaped trajectory group limits the statistical power to explore subgroup-specific mechanisms, such as differences in drug resistance patterns or organ damage trajectories. Larger multi-center cohorts are needed to confirm the clinical significance of the U-shaped pattern and identify potential interventions to prevent calcium rebound. Third, while E-value analysis (E = 2.494) supports the robustness of our findings to unmeasured confounding, residual factors such as bisphosphonate use, nutritional status, and fluid balance—all of which can influence serum calcium levels—were not fully accounted for. Future analyses should incorporate these variables to refine trajectory modeling and strengthen causal inferences. Finally, the retrospective design limits our ability to assess the impact of targeted interventions (e.g., early intensification of anti-myeloma therapy) on trajectory modification and survival. Prospective interventional studies are needed to determine whether altering unfavorable trajectories (e.g., preventing calcium rebound in Class 3) improves outcomes in high-risk critically ill MM patients. Conclusion In conclusion, this study demonstrates that dynamic serum calcium trajectories, modeled via LCGM, are robust, independent prognostic markers for 60-day mortality in critically ill MM patients. A key insight is the heterogeneous prognostic significance of hypercalcemia: sustained high-level or U-shaped (rebound) trajectories identify patients with the poorest outcomes, likely reflecting aggressive disease biology or treatment resistance. In contrast, stable low-to-medium trajectories are associated with significantly improved survival. These findings support a paradigm shift from static calcium assessment to dynamic trajectory monitoring, positioning serial calcium measurement as a real-time biomarker for disease activity in the ICU. Integrating this approach with conventional severity scores may refine early risk stratification and guide personalized management. Abbreviations Multiple myeloma (MM);Latent Class Growth Modeling (LCGM);Latent Class Growth Modeling (LCGM);Intensive Care Unit (ICU);International Classification of Diseases (ICD);Structured Query Language (SQL);Strengthening the Reporting of Observational Studies in Epidemiology (STROBE);Akaike Information Criterion (AIC);Bayesian Information Criterion (BIC);Sample-adjusted BIC (SABIC);Lo-Mendell–Rubin (LMR);Acute Physiology Score III (APSIII);Simplified Acute Physiology Score II (SAPSII);Oxford Acute Severity of Illness Score (OASIS);Sequential Organ Failure Assessment (SOFA);Hazard Ratio (HR);Confidence Interval (CI);Kaplan-Meier (KM);Skeletal-related events (SREs);White blood cell (WBC);Blood urea nitrogen (BUN);Platelet (PLT);International Normalized Ratio (INR);Prothrombin Time (PT);Activated Partial Thromboplastin Time (APTT);Heart rate (HR);Beats per minute (bpm);Systolic blood pressure (SBP);Diastolic blood pressure (DBP);Respiratory Rate (RR);Oxygen saturation (SpO₂);Creatinine (Cr);Sodium (Na);Potassium (K);Anion gap (AG);Bicarbonate (HCO₃⁻);Chloride (Cl⁻);Glucose (Glu);Hematocrit (HCT);Hemoglobin (Hb);Interleukin-6 (IL-6);Macrophage inflammatory protein-1α (MIP-1α);Tumor necrosis factor-α (TNF-α);Fluorescence in situ hybridization (FISH). Declarations Author contributions Yan Zeng: Writing-original draft, Writing-review & editing, Data curation, Methodology. Xiaomin Cheng: Writing-original draft, Data curation, Methodology. Yongran Si: Data curation, Validation. Jun Wang: Formal analysis, Software, Visualization. Jingwei Zhang: Writing-review & editing, Conceptualization, Investigation, Supervision. All authors read and approved the final manuscript. Funding This work was supported by the Technological Innovation Project of Chengdu Municipal Science and Technology Bureau (No. 2024-YF05-00936-SN). Data accessibility The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author. Data in the article can be obtained from mimic-IV data- base (https://mimic.mit.edu/) Ethical considerations The study utilized data from the MIMIC (Medical Information Mart for Intensive Care) database, which is a publicly available dataset that does not require individual ethical approval. The data are de-identified and anonymized, ensuring that no personally identifiable information can be linked to individual participants. Therefore, ethical approval was not required for this study. There was no requirement of individual informed consent to extract data from the mimic-IV database because mimic-IV database information was publicly available and all patient data were deidentified Declaration of generative AI and AI-assisted technologies in the writing process During the preparation of this work, the authors used ChatGPT in order to polish. After using this tool or service, the authors reviewed and edited the content as needed and took full responsibility for the content of the publication. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Acknowledgements We express our profound gratitude to the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center for their invaluable contribution to the MIMIC project. Publisher’s note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher. References Palumbo, A., et al., Personalized therapy in multiple myeloma according to patient age and vulnerability: a report of the European Myeloma Network (EMN). Blood, 2011. 118 (17): p. 4519-29. Cowan, A.J., et al., Diagnosis and Management of Multiple Myeloma: A Review. Jama, 2022. 327 (5): p. 464-477. Zagouri, F., et al., Hypercalcemia remains an adverse prognostic factor for newly diagnosed multiple myeloma patients in the era of novel antimyeloma therapies. Eur J Haematol, 2017. 99 (5): p. 409-414. Mousseaux, C., et al., Epidemiology, clinical features, and management of severe hypercalcemia in critically ill patients. Ann Intensive Care, 2019. 9 (1): p. 133. Dimopoulos, M.A., et al., Multiple myeloma: EHA-ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up(†). Ann Oncol, 2021. 32 (3): p. 309-322. Mao, Y., S. Zhu, and Y. Geng, Association between serum calcium and in-hospital mortality in critical patients with multiple myeloma: a cohort study. Hematology, 2022. 27 (1): p. 795-801. Ye, Q., et al., Serial platelet count as a dynamic prediction marker of hospital mortality among septic patients. Burns & Trauma, 2024. 12 . Song, M., et al., Trajectory of body shape in early and middle life and all cause and cause specific mortality: results from two prospective US cohort studies. Bmj, 2016. Kyheng, M., et al., Joint latent class model: Simulation study of model properties and application to amyotrophic lateral sclerosis disease. BMC Med Res Methodol, 2021. 21 (1): p. 198. Si, Y., et al., Time series analysis between platelet counts and 60-day mortality in sepsis patients with thrombocytopenia: a retrospective cohort study. BMC Infectious Diseases, 2025. 25 (1). Tang, S., et al., Association between red blood cell distribution width-platelet ratio (RPR) and mortality in patients with heart failure from the MIMIC IV database: A retrospective cohort study. Heliyon, 2024. 10 (16). Zeng, Y., et al., Association between Platelet Count and In-Hospital Mortality in Critical Patients with Multiple Myeloma: A Cohort Study. PLOS One, 2025. 20 (6). von Elm, E., et al., The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet, 2007. 370 (9596): p. 1453-7. Li, J., et al., Predicting Mortality in Intensive Care Unit Patients With Heart Failure Using an Interpretable Machine Learning Model: Retrospective Cohort Study. Journal of Medical Internet Research, 2022. 24 (8). Nagin, D.S., B.L. Jones, and J. Elmer, Recent Advances in Group-Based Trajectory Modeling for Clinical Research. Annu Rev Clin Psychol, 2024. 20 (1): p. 285-305. Kim, S.Y., Determining the Number of Latent Classes in Single- and Multi-Phase Growth Mixture Models. Struct Equ Modeling, 2014. 21 (2): p. 263-279. Ding, X., et al., Associations between sepsis occurrence, hemoglobin level and mortality in patients with non-trauma hemorrhagic brain injuries: trajectory-based analysis. European Journal of Medical Research, 2025. 30 (1). Zhang, Q., et al., The prognostic value of the platelet-to-lymphocyte ratio in multiple myeloma patients treated with a bortezomib-based regimen. Scientific Reports, 2025. 15 (1). Chung, W.T. and K.C. Chung, The use of the E-value for sensitivity analysis. J Clin Epidemiol, 2023. 163 : p. 92-94. Cui, K., et al., The association between blood urea nitrogen to albumin ratio and the 28 day mortality in tuberculosis patients complicated by sepsis. Scientific Reports, 2024. 14 (1). Tothong, W., et al., Prevalence, Outcomes and Impact of Disease-Related Complications in the Survival of Multiple Myeloma Patients. Hematol Rep, 2024. 16 (1): p. 89-97. Rees, M.J. and S. Kumar, High-risk multiple myeloma: Redefining genetic, clinical, and functional high-risk disease in the era of molecular medicine and immunotherapy. Am J Hematol, 2024. 99 (8): p. 1560-1575. Hampshire, P.A., et al., Admission factors associated with hospital mortality in patients with haematological malignancy admitted to UK adult, general critical care units: a secondary analysis of the ICNARC Case Mix Programme Database. Crit Care, 2009. 13 (4): p. R137. Jeon, J.E., et al., A Pathogen-Responsive Gene Cluster for Highly Modified Fatty Acids in Tomato. Cell, 2020. 180 (1): p. 176-187.e19. Weidle, U.H., et al., Molecular Mechanisms of Bone Metastasis. Cancer Genomics Proteomics, 2016. 13 (1): p. 1-12. Kanellias, N., et al., Newly Diagnosed Multiple Myeloma Patients with Skeletal-Related Events and Abnormal MRI Pattern Have Poor Survival Outcomes: A Prospective Study on 370 Patients. J Clin Med, 2022. 11 (11). Freire-de-Lima, L., et al., Multiple Myeloma Cells Express Key Immunoregulatory Cytokines and Modulate the Monocyte Migratory Response. Front Med (Lausanne), 2017. 4 : p. 92. Bao, L., et al., Hypercalcemia caused by humoral effects and bone damage indicate poor outcomes in newly diagnosed multiple myeloma patients. Cancer Med, 2020. 9 (23): p. 8962-8969. Tables Tables 1 to 5 are available in the Supplementary Files section. Additional Declarations No competing interests reported. Supplementary Files Table1to5.docx supp.docx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8812562","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":599216096,"identity":"723403d9-923f-4006-a180-0d6201f74438","order_by":0,"name":"Yan Zeng","email":"","orcid":"","institution":"Sichuan University, Sichuan University affiliated Chengdu Second People's Hospital, Chengdu Second People's Hospital","correspondingAuthor":false,"prefix":"","firstName":"Yan","middleName":"","lastName":"Zeng","suffix":""},{"id":599216100,"identity":"2f2b713e-086e-4436-a2ac-85bf271dcb3e","order_by":1,"name":"Xiaomin Cheng","email":"","orcid":"","institution":"Sichuan University, Sichuan University affiliated Chengdu Second People's Hospital, Chengdu Second People's Hospital","correspondingAuthor":false,"prefix":"","firstName":"Xiaomin","middleName":"","lastName":"Cheng","suffix":""},{"id":599216101,"identity":"4ab854a6-24ba-49fb-9552-d93594d35f49","order_by":2,"name":"Yongran Si","email":"","orcid":"","institution":"Sichuan University, Sichuan University affiliated Chengdu Second People's Hospital, Chengdu Second People's Hospital","correspondingAuthor":false,"prefix":"","firstName":"Yongran","middleName":"","lastName":"Si","suffix":""},{"id":599216103,"identity":"cde0d96b-fc76-4d4a-a23f-ef5a64277400","order_by":3,"name":"Jun Wang","email":"","orcid":"","institution":"Sichuan University, Sichuan University affiliated Chengdu Second People's Hospital, Chengdu Second People's Hospital","correspondingAuthor":false,"prefix":"","firstName":"Jun","middleName":"","lastName":"Wang","suffix":""},{"id":599216105,"identity":"d3b4f4e6-e9d3-4806-b52a-17b88a360e7c","order_by":4,"name":"Jingwei Zhang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAwUlEQVRIiWNgGAWjYBACfvbG5sc/DP7X8xOtRbLncJsxQwFzgmQDsVoMbqQ3SDN8YE4wOEC0NQcSG4wLDNjyjI8nb2D4UbGNsA7GhoMNj2cY8BSbnXlWwNhz5jZhLcyMjQ0GPAYSjNtu5BgwM7YRoYWNmbFBgsfAgHHzDGK18LAxNkjzGCQkbpAgVosED2Ob4QyDA8YSQL8cJMov9vefP37w4c8BOf725I0PflQQoQUJkBI1cC2k6hgFo2AUjIIRAgDnwz8UQmAXQwAAAABJRU5ErkJggg==","orcid":"","institution":"Sichuan University, Sichuan University affiliated Chengdu Second People's Hospital, Chengdu Second People's Hospital","correspondingAuthor":true,"prefix":"","firstName":"Jingwei","middleName":"","lastName":"Zhang","suffix":""}],"badges":[],"createdAt":"2026-02-07 05:53:16","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8812562/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8812562/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":104404020,"identity":"f1eb76ae-e733-4993-94ea-47c18e7e6c2b","added_by":"auto","created_at":"2026-03-11 12:19:37","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":262770,"visible":true,"origin":"","legend":"\u003cp\u003e\u0026nbsp;See image above for figure legend.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-8812562/v1/3aa46e4de7c0b4485d83804c.png"},{"id":104177949,"identity":"a608beca-6f14-4ab7-89b0-3f23808f6118","added_by":"auto","created_at":"2026-03-08 16:50:56","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":309481,"visible":true,"origin":"","legend":"\u003cp\u003e\u0026nbsp;See image above for figure legend.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-8812562/v1/4fc8994f2bc9d49081ec7178.png"},{"id":105216134,"identity":"d1a2c9c2-3fd2-4e0f-bcbc-298c59c35e19","added_by":"auto","created_at":"2026-03-23 14:42:22","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1187643,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8812562/v1/9743073e-854a-4fa0-9eda-48bc585a84bf.pdf"},{"id":104404606,"identity":"f8401128-55f8-4d48-b854-90e932708a1b","added_by":"auto","created_at":"2026-03-11 12:20:37","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":30147,"visible":true,"origin":"","legend":"","description":"","filename":"Table1to5.docx","url":"https://assets-eu.researchsquare.com/files/rs-8812562/v1/0ac1cc45c7aa7708fb1b9f93.docx"},{"id":104177952,"identity":"2ad65b06-810e-4ab5-8ae6-65ab09dca1b1","added_by":"auto","created_at":"2026-03-08 16:50:57","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":340551,"visible":true,"origin":"","legend":"","description":"","filename":"supp.docx","url":"https://assets-eu.researchsquare.com/files/rs-8812562/v1/6e3d697098a50141efda17ee.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Time series analysis Between Serum Calcium Levels and 60-Day Mortality in Critically ill Patients with Multiple Myeloma: A Retrospective Cohort Study","fulltext":[{"header":"Introduction","content":"\u003cp\u003eMultiple myeloma (MM) is a hematologic malignancy characterized by the clonal proliferation of plasma cells within the bone marrow. Globally, MM incidence is on the rise, making it the second most common hematologic cancer after Hodgkin\u0026apos;s lymphoma. In Europe, the annual incidence rate ranges from 4.5 to 6.0 cases per 100,000 individuals, with a corresponding mortality rate of 4.1 per 100,000 [1]. Complications associated with MM include anemia, hypercalcemia, renal dysfunction, and bone disease. Hypercalcemia is observed in 19.5% of MM patients at diagnosis, with severity ranging from mild to life-threatening, and is frequently linked to bone lesions, kidney dysfunction, and high-risk chromosomal abnormalities in MM [2, 3]. Notably, emerging evidence suggests that critically ill MM patients admitted to the intensive care unit (ICU) exhibit a significantly higher incidence of hypercalcemia (up to 30%-40%) compared to the general MM population. In such patients, hypercalcemia is not only a consequence of bone resorption driven by myeloma cell-derived cytokines, but also intertwined with ICU-specific insults such as acute kidney injury (AKI) and sepsis, forming a vicious cycle that accelerates organ dysfunction and poor outcomes. This unique pathophysiological feature highlights the need to explore dynamic calcium changes rather than static baseline values in this population[4]. Hypercalcemia at initial diagnosis of MM is widely recognized as an adverse prognostic factor [2, 3].\u0026nbsp;Given the complexity of critical illness and the variability of medical interventions (e.g., anti-myeloma therapy that may alter bone metabolism, renal replacement therapy that regulates calcium excretion), serum calcium levels change dynamically over time.\u0026nbsp;The persistence or exacerbation of hypercalcemia in critically ill MM patients may further increase mortality risk.\u003c/p\u003e\n\u003cp\u003eGuidelines, including the EHA-ESMO clinical practice guidelines for MM, recommend monthly or at least quarterly monitoring of core laboratory parameters (complete blood count, creatinine, and serum calcium) during follow-up, but lack guidance on interpreting dynamic fluctuations, leaving clinicians unable to leverage serial calcium data for real-time risk assessment in ICU settings[5, 6]. A recent study further demonstrated a nonlinear relationship between serum calcium levels and in-hospital mortality in critically ill MM patients, with a serum calcium level of approximately 8.40 mg/dL associated with the lowest mortality risk and risk increasing with rising levels, an important consideration for ICU clinicians[6]. Despite this, most current studies focus on baseline serum calcium data from newly diagnosed MM patients at admission [2, 6], failing to account for the dynamic fluctuations of serum calcium levels during ICU stay. Given the complexity of critical illness and the variability of medical interventions, serum calcium levels change dynamically over time, and investigating these time-series changes and their prognostic implications is crucial for optimizing risk stratification and treatment adjustment [7].\u003c/p\u003e\n\u003cp\u003eMedical time series analysis (TSA) has been widely applied in prognostic studies [8]. Among TSA methods, the Latent Class Growth Modeling (LCGM) is particularly valuable for identifying heterogeneous trajectory subgroups in heterogeneous populations such as critically ill MM patients who may exhibit diverse calcium change patterns, by partitioning the population into latent classes based on individual characteristics (e.g., intercepts and slopes) and modeling class-specific change patterns over time[9]. Trajectory analysis using LCGM has been successfully applied to characterize dynamic biomarkers in various hematologic and critical diseases[10], but to date, no studies have utilized this approach to explore the relationship between dynamic serum calcium trajectories and mortality in MM patients, especially in the high-risk critically ill subgroup.\u003c/p\u003e\n\u003cp\u003eTo address this research gap, the present study leverages data from the MIMIC-IV database to construct LCGM of serum calcium trajectories during the first 28 days of hospitalization, integrating both baseline measurements and dynamic time-series data. The primary objectives are to identify distinct serum calcium trajectories in critically ill MM patients, evaluate their association with 60-day mortality, and provide evidence for refining early risk stratification and guiding personalized clinical decisions.\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cp\u003e\u003cstrong\u003ePopulation and Data Collection\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis retrospective cohort study utilized the Medical Information Mart for Intensive Care-IV (MIMIC-IV v3.0) database, comprising de-identified clinical data from \u0026gt;70,000 intensive care unit (ICU) admissions at Beth Israel Deaconess Medical Center (Boston, MA) between 2008 and 2019[6]. Access to the MIMIC-IV database was obtained in accordance with the official access protocol specified by PhysioNet (https://physionet.org/content/mimiciv/3.0/). Access to the database required completion of the \u0026quot;CITI Data or Specimens Only Research\u0026quot; training course on the National Institutes of Health (NIH) website, with approval granted for data extraction for research purposes (Certificate No for Zhang: ID: 13813672). \u0026nbsp;Ethical oversight was provided by institutional review boards of Massachusetts Institute Technology and Beth Israel Deaconess Medical Center. Informed consent was waived given the retrospective, de-identified nature of the data[6]. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003eStructured Query Language (SQL), PostgreSQL (version 13.0), and Navicat software (version 16.0) were used to identify the cohort and extract relevant clinical data [11]. Inclusion criteria were defined as: (1) adult patients aged \u0026ge;18 years; (2) a confirmed diagnosis of MM based on International Classification of Diseases (ICD) coding: ICD-9 codes (203.00, 203.01, 203.02) or ICD-10 codes (C90.00, C90.01, C90.02)[6]; (3) first-time ICU admission (to avoid duplicate patient data); and (4) availability of at least two serum calcium measurements during hospitalization. Exclusion criteria were: (1) absence of baseline serum calcium measurement on ICU admission; (2) fewer than two serum calcium measurements throughout the ICU hospitalization period; (3) incomplete data on the primary outcome (60-day mortality).\u003c/p\u003e\n\u003cp\u003eDemographics, comorbidities, and laboratory indicators were extracted from the MIMIC-IV database. Serum calcium levels were systematically collected at predefined time points: ICU days 1, 3, 5, and 7, with weekly averages calculated for weeks 2, 3, and 4 of hospitalization. Additionally, the minimum, maximum, and mean serum calcium values during the entire ICU stay were computed. For missing serum calcium data on days 3, 5, or 7, a last-observation-carried-forward approach was applied, wherein the nearest adjacent available measurement was imputed. For other laboratory parameters, maximum, minimum, or mean values were selected based on clinical relevance and prior literature[10, 12]. The primary outcome of interest was 60-day mortality. Variables with missing data exceeding 10% were excluded from the analysis, consistent with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline[13]. For remaining variables with minor missingness (\u0026lt;10%), multiple imputation was performed using chained equations with 5 imputed datasets via the\u0026nbsp;mice\u0026nbsp;package in R [14].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of Serum Calcium Trajectory Subgroups and Model Selection\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eHeterogeneous serum calcium trajectories were identified using the R package \u0026quot;lcmm\u0026quot; (version 2.0.2) via LCGM [7]. A latent class mixed model was applied for longitudinal serum calcium data analysis, and a parametric survival model with a class-specific Weibull distribution was incorporated to characterize the baseline hazard function, enabling simultaneous modeling of longitudinal trajectories and time-to-60-day-mortality outcomes [7].Trajectory modeling was based on serum calcium data systematically collected within the initial 28 days after ICU admission to be consistent with the study\u0026rsquo;s predefined observation window, including measurements on ICU days 1, 3, 5, and 7, as well as weekly averages for weeks 2-4. Both linear and quadratic functional forms for the time variable were evaluated to optimize model fit. Quadratic terms were included to capture potential acceleration or deceleration in serum calcium changes commonly observed in critically ill patients[15].\u003c/p\u003e\n\u003cp\u003eModel selection was guided by a comprehensive set of criteria to ensure statistical robustness and clinical interpretability:Lower values of the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), and sample-adjusted BIC (SABIC) indicated superior model fit; A significant P-value (\u0026lt; 0.05) from the Lo-Mendell\u0026ndash;Rubin (LMR) adjusted likelihood ratio test supported adding an additional trajectory class; Visual inspection of trajectory plots was conducted to confirm the clinical plausibility of identified patterns; Higher maximum log-likelihood and entropy (\u0026ge; 0.80, indicating strong classification accuracy) were prioritized; Robust subgroup classification required an average posterior probability \u0026ge; 70% for each class; Only models with each subgroup containing \u0026gt; 3% of the total study population were retained to avoid unstable estimates from small subgroups[10, 16, 17].The optimal number of trajectory subgroups was determined by balancing all aforementioned criteria, with the final model selected based on the best overall fit and clinical interpretability.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eStatistical analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eContinuous variables were described as mean \u0026plusmn; standard deviation (SD) for normally distributed data (assessed via the Shapiro-Wilk test) and median (interquartile range, IQR) for skewed data. Group comparisons across latent trajectory classes were performed using the Kruskal-Wallis test for continuous variables. Categorical variables were presented as counts (percentages, %), with group comparisons conducted using Pearson\u0026rsquo;s chi-squared test; Fisher\u0026rsquo;s exact test was used when the expected frequency of any cell was \u0026lt; 5. Univariate and multivariate Cox proportional hazards regression analyses were performed to evaluate the association between candidate variables (including serum calcium trajectories) and 60-day mortality. The selection of covariates for adjustment in multivariate models was guided by three principles: (1) clinical relevance based on prior literature[18], (2) expert consensus on risk factors for mortality in critically ill MM patients[6], and (3) statistical significance (p \u0026lt; 0.05) in univariate Cox regression. Four hierarchical multivariate Cox models were constructed to assess the robustness of the association between serum calcium trajectories and 60-day mortality. Model 1: Unadjusted (only serum calcium trajectory as the predictor);Model 2: Adjusted for age and gender (basic demographic confounders);Model 3: Adjusted for age, gender, mean systolic blood pressure (SBP_mean), maximum blood urea nitrogen (Bun_max), minimum platelet count (Plt_min), and congestive heart failure (clinical and laboratory confounders);Model 4: Fully adjusted for age, gender, SBP_mean, Bun_max, Plt_min, congestive heart failure, Acute Physiology Score III (APSIII), Simplified Acute Physiology Score II (SAPSII), Oxford Acute Severity of Illness Score (OASIS), and Sequential Organ Failure Assessment (SOFA) score (adding disease severity indices). Kaplan\u0026ndash;Meier (KM) survival curves were generated to visualize survival probabilities across serum calcium trajectory classes, with between-group differences evaluated using the log-rank test. The E-value calculation Quantifies the potential impact of unmeasured confounding on the observed associations. The E-value was defined as the minimum strength of association (on the risk ratio scale) that an unmeasured confounder would need to have with both the exposure (serum calcium trajectory) and the outcome (60-day mortality), conditional on measured covariates, to fully explain the observed association [19]. All statistical analyses were conducted using R Statistical Software (version 4.2.2; http://www.R-project.org, The R Foundation) and the Free Statistics Analysis Platform (version 2.1.1; Beijing, China). Specific R packages included survival for Cox regression and KM analysis, and car for testing model assumptions. A two-tailed p-value \u0026lt; 0.05 was considered statistically significant[20].\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eBaseline characteristics\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA total of 254 patients were included in this study after screening according to predefined inclusion and exclusion criteria (\u003cu\u003eSupplementary Figure 1\u003c/u\u003e). The overall 60-day mortality rate was 24.02% (61/254). The baseline characteristics of the study population stratified by 60-day mortality status are summarized in\u0026nbsp;\u003cu\u003eTable 1\u003c/u\u003e.\u003c/p\u003e\n\u003cp\u003eCompared with the surviving group (n = 193), the non-surviving group (n = 61) had a significantly older age and a higher prevalence of congestive heart failure. Regarding vital signs, the non-surviving group exhibited a lower mean diastolic blood pressure \u0026nbsp; and a higher mean respiratory rate. For laboratory parameters, the non-surviving group had a significantly lower minimum platelet count and a higher maximum blood urea nitrogen level. In terms of illness severity scores, the non-surviving group had significantly higher scores on the APSIII, SAPSII, OASIS, and SOFA score.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of subpopulations using the LCGM\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe optimal number of serum calcium trajectory subgroups was determined based on multiple model fit metrics, including log-likelihood, AIC, BIC, SABIC, entropy, and class composition ratio (\u003cu\u003eTable 2\u003c/u\u003e). As the number of classes increased from 1 to 5, the log-likelihood, AIC, and SABIC values initially decreased and then increased, with the 4-class model yielding the lowest AIC (2262.018) and SABIC (2271.196) values. Additionally, the 4-class model exhibited a moderate log-likelihood (-1106.008962), a relatively high entropy (0.683), and each class accounted for \u0026gt;3% of the total population (ranging from 3.544% to 57.480%), meeting the criteria for robust subgroup classification.Thus, the 4-class model was selected as the optimal solution, and the trajectories were labeled as follows (\u003cu\u003eFigure 1\u003c/u\u003e). Class 1(Low-level trajectory, 35.433%, n = 90) of patients, characterized by a slight increase from an initially low serum calcium level; Class 2 (Medium-level trajectory, 57.480%, n = 138) of patients, characterized by a slight decrease from a moderate baseline serum calcium level;\u003c/p\u003e\n\u003cp\u003eClass 3 (U-shaped trajectory, 3.544%, n = 9) of patients, characterized by serum calcium levels decreasing from an initial high level to a low point and then rising back to a medium level; Class 4 (High-level trajectory, 6.683%, n = 17) of patients, characterized by persistently high serum calcium levels maintained throughout the observation period. The baseline characteristics of the study population stratified by serum calcium trajectories are summarized in\u0026nbsp;\u003cu\u003eSupplementary Table\u003c/u\u003e\u003cu\u003e\u0026nbsp;1\u003c/u\u003e.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eUnivariable\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eCox Regression Analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAs shown in\u0026nbsp;\u003cu\u003eTable 3\u003c/u\u003e, univariable Cox regression analysis was performed with the high-level trajectory (Class 4) as the reference group. The results revealed that patients in the low-level trajectory (Class 1) and medium-level trajectory (Class 2) had significantly lower 60-day mortality risk (Class 1: HR = 0.28, 95% CI = 0.13-0.63, p = 0.002; Class 2: HR = 0.34, 95% CI = 0.16-0.71, p = 0.004). In contrast, no significant difference in 60-day mortality was observed between the U-shaped trajectory (Class 3) and the reference group (HR = 0.51, 95% CI = 0.14-1.89, p = 0.313).Additionally, univariable analysis identified several independent risk factors for 60-day mortality, including older age, lower DBP-mean, higher RR-mean, lower PLT-min, higher BUN-max, presence of congestive heart failure, and higher APSIII, SAPSII, OASIS, and SOFA scores (all p \u0026lt; 0.001).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMultivariable Cox Regression Analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo adjust for potential confounding factors, four hierarchical multivariable Cox regression models were constructed (\u003cu\u003eTable 4\u003c/u\u003e). In the fully adjusted model (Model 4), which included age, gender, SBP-mean, BUN-max, PLT-min, congestive heart failure, APSIII, SAPSII, OASIS, and SOFA score, the associations between serum calcium trajectories and 60-day mortality remained robust. Compared with the high-level trajectory (Class 4), the low-level trajectory (Class 1) and medium-level trajectory (Class 2) were still significantly associated with lower 60-day mortality risk (Class 1: HR = 0.17, 95% CI = 0.07-0.41, p \u0026lt; 0.001; Class 2: HR = 0.26, 95% CI = 0.12-0.56, p = 0.001). The U-shaped trajectory (Class 3) still showed no significant difference in 60-day mortality compared with the reference group (HR = 0.56, 95% CI = 0.14-2.25, p = 0.416).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eKaplan-Meier curves analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eKaplan-Meier (KM) survival curves stratified by serum calcium trajectory subgroups are presented in\u0026nbsp;\u003cu\u003eFigure 2\u003c/u\u003e. The results showed that the high-level trajectory (Class 4) had the lowest 60-day survival probability, while the low-level trajectory (Class 1) had the highest, followed by the medium-level trajectory (Class 2). Although the survival curve of the U-shaped trajectory (Class 3) intersected with those of Class 1 and Class 2 around weeks 2\u0026ndash;3, the log-rank test indicated significant difference in survival probabilities among the four subgroups (p = 0.0082).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSensitivity Analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eGiven the lack of significant difference in 60-day survival rates between Class 3 and Class 4, we compared serum calcium levels across the trajectory classes at various time points. As shown in \u003cu\u003eTable 5\u003c/u\u003e, Class 1 consistently exhibited lower serum calcium levels on days 1, 3, 5, and 7, as well as lower mean, minimum, and maximum values during hospitalization compared to the other three classes. The initial values for Class 2 were lower than those for Class 3 during the first three days, became comparable by day 5, and subsequently dropped below Class 3 by week three. Class 3 initially had higher values than Class 4 over the first three days, followed by a gradual decline and subsequent rise, while remaining consistently lower than Class 4. When stratifying the population by 60-day mortality and comparing serum calcium levels between survivors and non-survivors at multiple time points, no significant differences were observed (all p \u0026gt; 0.05) (\u003cu\u003eSupplementary Table 2\u003c/u\u003e). These findings suggest that Class 3 and Class 4 share similar prognostic characteristics attributable to their consistently elevated calcium levels.\u003c/p\u003e\n\u003cp\u003eBased on the similar prognostic characteristics between the U-shaped trajectory (Class 3) and the high-level trajectory (Class 4) identified in the primary analysis, we merged these two classes into \u0026quot;Group 3\u0026quot; and reclassified the low-level trajectory (Class 1) and medium-level trajectory (Class 2) as \u0026quot;Group 1\u0026quot; and \u0026quot;Group 2\u0026quot;, respectively. Multivariable Cox regression analysis with Group 3 as the reference group (\u003cu\u003eSupplementary Table 3\u003c/u\u003e) showed that after adjusting for the same covariates, Group 1 (HR = 0.21, 95% CI = 0.10\u0026ndash;0.45, p \u0026lt; 0.001) and Group 2 (HR = 0.30, 95% CI = 0.15\u0026ndash;0.61, p = 0.001) still had significantly lower 60-day mortality risk. KM survival curves for the three reclassified groups confirmed that Group 3 had the lowest 60-day survival probability (\u003cu\u003eSupplementary Figure 2\u003c/u\u003e), with a statistically significant difference among the three groups (overall p = 0.0091).\u003c/p\u003e\n\u003cp\u003eTo evaluate the robustness of the observed associations to unmeasured confounding, an E-value analysis was performed. The E-value for the point estimate of the association between serum calcium trajectories and 60-day mortality was 2.494 (95% CI: 1.492\u0026ndash;4.051). This relatively large E-value indicates that an unmeasured confounder would need to have a strong association with both serum calcium trajectory and 60-day mortality (with a risk ratio of at least 2.494) to fully explain the observed association, suggesting that the study results are robust to potential unmeasured confounding.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eTo the best of our knowledge, this is the first study to identify heterogeneous serum calcium trajectories in \u003cb\u003ecritically ill MM patients\u003c/b\u003e using Latent Class Growth Modeling (LCGM) and to systematically explore the association between these dynamic trajectories and 60-day mortality. This work addresses a critical gap in the existing literature, which has predominantly focused on static baseline serum calcium levels in newly diagnosed or stable MM patients[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e] [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e], rather than capturing the dynamic fluctuations inherent to critical illness, characterized by complex interactions between myeloma-related bone resorption, acute kidney injury, sepsis, and intensive care interventions[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eOur findings indicate that the overall 60-day mortality rate among critically ill MM patients was 24.02%, which is lower than the 30-day mortality rate (35.2%) reported by Hampshire et al. [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e] in a cohort of hematologic malignancy patients admitted to UK ICUs, but slightly higher than the in-hospital mortality rate (19.4%) reported by Mao et al. in critically ill MM patients[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. This discrepancy is likely multifactorial: first, the choice of endpoint (60-day mortality vs. 30-day/in-hospital mortality) captures a longer follow-up period, including post-discharge deaths that may be associated with unresolved organ dysfunction or disease progression; second, variations in cohort severity; third, differences in healthcare settings and treatment patterns (e.g., utilization of novel anti-myeloma agents or renal replacement therapy) between study populations may further explain the variance[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e].A key novel finding of this study is the identification of four distinct serum calcium trajectories over the first 28 days of ICU hospitalization, with clear prognostic stratification. Specifically, patients with low-level (Class 1) and medium-level (Class 2) trajectories exhibited significantly lower 60-day mortality (HR\u0026thinsp;=\u0026thinsp;0.17 and 0.26, respectively, in the fully adjusted model) compared to those with high-level (Class 4) trajectories. Notably, the U-shaped trajectory (Class 3), characterized by an initial high calcium level, transient decline, and subsequent rise, showed no significant mortality difference from the high-level trajectory (HR\u0026thinsp;=\u0026thinsp;0.56, p\u0026thinsp;=\u0026thinsp;0.416). This observation is crucial. It suggests that even a transient normalization of serum calcium does not mitigate the adverse prognosis portended by the occurrence of hypercalcemia in the critical care setting[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e] [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e], but extends these findings by demonstrating that dynamic trajectory patterns, rather than single-point measurements, provide a more robust prognostic assessment.\u003c/p\u003e \u003cp\u003eElevated serum calcium levels are inversely correlated with survival in MM patients [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e], a relationship deeply rooted in the disease\u0026rsquo;s biological aggressiveness and organ-damaging sequelae. This correlation is primarily driven by the intimate link between hypercalcemia and skeletal-related events (SREs)\u0026mdash;including pathological fractures, extensive osteolytic lesions, renal impairment, and increased infection risk-all of which independently worsen prognosis [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e, \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Mechanistically, MM-related hypercalcemia is predominantly mediated by osteoclast hyperactivation, a process amplified by myeloma cell-derived cytokines such as IL-6 and TNF-α that disrupt the bone remodeling balance, leading to unregulated bone resorption and calcium release into the circulation [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. Notably, this cytokine-driven bone destruction is not only a cause of hypercalcemia but also a hallmark of high-risk MM biology, as these cytokines simultaneously promote myeloma cell proliferation and immune suppression[\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eAlthough the full spectrum of mechanisms by which hypercalcemia impairs survival remains incompletely elucidated, accumulating evidence highlights its association with adverse biological features. The presence of hypercalcemia at diagnosis is an established independent poor prognostic factor, with nearly 20% of newly diagnosed symptomatic MM patients presenting with this complication[\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. Critically, hypercalcemia in MM is strongly linked to high-risk cytogenetic abnormalities\u0026mdash;including del 13q (detected by fluorescence in situ hybridization) and amp1q21, which are themselves associated with rapid disease progression, treatment resistance, and shortened overall survival [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. The prognostic significance of hypercalcemia persists across age groups and treatment modalities, including conventional chemotherapy and novel agents such as thalidomide, lenalidomide, and bortezomib, suggesting that hypercalcemia is not merely a clinical manifestation but a surrogate marker of inherently aggressive disease biology and potential drug tolerance[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThis biological context provides critical insights into the four serum calcium trajectories identified in our study. The high-level trajectory (Class 4), characterized by persistently elevated calcium levels\u0026mdash;likely reflects uncontrolled myeloma activity, driven by high-risk cytogenetic profiles and sustained osteoclast-mediated bone destruction[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e] [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. These patients may harbor clones with inherent resistance to anti-myeloma therapies, as evidenced by the failure to normalize calcium levels despite intensive care interventions. Similarly, the U-shaped trajectory (Class 3), with initial hypercalcemia, transient reduction, and subsequent rebound, suggests a temporary response to acute interventions (e.g., hydration, bisphosphonates, or short-course anti-myeloma therapy) followed by disease progression or acquired drug resistance [30]. The recurrence of hypercalcemia in this trajectory aligns with the notion that cytogenetically driven bone destruction is not fully reversible with standard therapies, leading to persistent organ damage and poor outcomes comparable to those with sustained hypercalcemia[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e].In contrast, the low-level (Class 1) and medium-level (Class 2) trajectories, characterized by stable calcium levels within or below the normal range, are indicative of better disease control and treatment sensitivity.\u003c/p\u003e \u003cp\u003eThese trajectories likely represent patients with non-aggressive disease biology (e.g., absence of high-risk cytogenetics), effective suppression of myeloma-driven cytokine production, and responsive bone remodeling following treatment. The stability of calcium levels in these groups suggests durable control of both myeloma cell proliferation and osteoclast activity. Notably, the segmented linear distribution of average calcium levels across trajectories (Class 1\u0026thinsp;\u0026lt;\u0026thinsp;Class 2\u0026thinsp;\u0026lt;\u0026thinsp;Class 3\u0026thinsp;\u0026lt;\u0026thinsp;Class 4) mirrors the gradient of prognostic risk, consistent with Mao et al.\u0026rsquo;s observation of a nonlinear relationship between calcium levels and mortality [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. This continuity further validates that serum calcium trajectories are not arbitrary groupings but reflect a biological spectrum of disease aggressiveness, from well-controlled to refractory disease. Collectively, our findings highlight the direct clinical applicability of dynamic serum calcium trajectory monitoring. In the ICU setting, these trajectories can function as an early-warning biomarker, prompting clinicians to intensify management of myeloma-specific complications, such as escalating bone-targeted therapy or reassessing anti-myeloma regimens, when a persistently high or rebounding pattern is observed. Incorporation of trajectory analysis with validated severity scores may enhance risk stratification and support timely clinical decision-making.\u003c/p\u003e \u003cp\u003eThis study has several notable strengths. First, it innovatively applies LCGM to characterize serum calcium trajectories in critically ill MM patients, shifting the focus from static to dynamic biomarker assessment, a paradigm increasingly recognized in critical care research. Second, the use of hierarchical multivariable Cox regression and E-value analysis (E\u0026thinsp;=\u0026thinsp;2.494) enhances the robustness of the observed associations against confounding factors. Third, the identification of prognostic trajectories provides a practical tool for risk stratification, addressing the unmet need for dynamic monitoring guidance in critically ill MM management.\u003c/p\u003e \u003cp\u003eDespite these insights, our study has several limitations that warrant consideration. First, the reliance on a single-center public database (MIMIC-IV) precluded access to detailed cytogenetic data (e.g., del 13q, amp1q21), bone marrow biopsy results, and long-term treatment response data\u0026mdash;critical variables needed to directly validate the link between trajectories, high-risk biology, and treatment resistance. Future prospective studies should integrate these variables to explicitly test whether trajectory patterns correlate with specific genetic signatures and treatment outcomes. Second, the small sample size of the U-shaped trajectory group limits the statistical power to explore subgroup-specific mechanisms, such as differences in drug resistance patterns or organ damage trajectories. Larger multi-center cohorts are needed to confirm the clinical significance of the U-shaped pattern and identify potential interventions to prevent calcium rebound. Third, while E-value analysis (E\u0026thinsp;=\u0026thinsp;2.494) supports the robustness of our findings to unmeasured confounding, residual factors such as bisphosphonate use, nutritional status, and fluid balance\u0026mdash;all of which can influence serum calcium levels\u0026mdash;were not fully accounted for. Future analyses should incorporate these variables to refine trajectory modeling and strengthen causal inferences. Finally, the retrospective design limits our ability to assess the impact of targeted interventions (e.g., early intensification of anti-myeloma therapy) on trajectory modification and survival. Prospective interventional studies are needed to determine whether altering unfavorable trajectories (e.g., preventing calcium rebound in Class 3) improves outcomes in high-risk critically ill MM patients.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eIn conclusion, this study demonstrates that dynamic serum calcium trajectories, modeled via LCGM, are robust, independent prognostic markers for 60-day mortality in critically ill MM patients. A key insight is the heterogeneous prognostic significance of hypercalcemia: sustained high-level or U-shaped (rebound) trajectories identify patients with the poorest outcomes, likely reflecting aggressive disease biology or treatment resistance. In contrast, stable low-to-medium trajectories are associated with significantly improved survival. These findings support a paradigm shift from static calcium assessment to dynamic trajectory monitoring, positioning serial calcium measurement as a real-time biomarker for disease activity in the ICU. Integrating this approach with conventional severity scores may refine early risk stratification and guide personalized management.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eMultiple myeloma (MM);Latent Class Growth Modeling (LCGM);Latent Class Growth Modeling (LCGM);Intensive Care Unit (ICU);International Classification of Diseases (ICD);Structured Query Language (SQL);Strengthening the Reporting of Observational Studies in Epidemiology (STROBE);Akaike Information Criterion (AIC);Bayesian Information Criterion (BIC);Sample-adjusted BIC (SABIC);Lo-Mendell\u0026ndash;Rubin (LMR);Acute Physiology Score III (APSIII);Simplified Acute Physiology Score II (SAPSII);Oxford Acute Severity of Illness Score (OASIS);Sequential Organ Failure Assessment (SOFA);Hazard Ratio (HR);Confidence Interval (CI);Kaplan-Meier (KM);Skeletal-related events (SREs);White blood cell (WBC);Blood urea nitrogen (BUN);Platelet (PLT);International Normalized Ratio (INR);Prothrombin Time (PT);Activated Partial Thromboplastin Time (APTT);Heart rate (HR);Beats per minute (bpm);Systolic blood pressure (SBP);Diastolic blood pressure (DBP);Respiratory Rate (RR);Oxygen saturation (SpO₂);Creatinine (Cr);Sodium (Na);Potassium (K);Anion gap (AG);Bicarbonate (HCO₃⁻);Chloride (Cl⁻);Glucose (Glu);Hematocrit (HCT);Hemoglobin (Hb);Interleukin-6 (IL-6);Macrophage inflammatory protein-1\u0026alpha; (MIP-1\u0026alpha;);Tumor necrosis factor-\u0026alpha; (TNF-\u0026alpha;);Fluorescence in situ hybridization (FISH).\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eYan Zeng: Writing-original draft, Writing-review \u0026amp; editing, Data curation, Methodology. Xiaomin Cheng: Writing-original draft, Data curation, Methodology.\u003c/p\u003e\n\u003cp\u003eYongran Si: Data curation, Validation. Jun Wang: Formal analysis, Software, Visualization. Jingwei Zhang: Writing-review \u0026amp; editing, Conceptualization, Investigation, Supervision. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the Technological Innovation Project of Chengdu Municipal Science and Technology Bureau (No. 2024-YF05-00936-SN).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData accessibility\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author. Data in the article can be obtained from mimic-IV data-\u003c/p\u003e\n\u003cp\u003ebase (https://mimic.mit.edu/)\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthical considerations\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe study utilized data from the MIMIC (Medical Information Mart for Intensive Care) database, which is a publicly available dataset that does not require individual ethical approval. The data are de-identified and anonymized, ensuring that no personally identifiable information can be linked to individual participants. Therefore, ethical approval was not required for this study. There was no requirement of individual informed consent to extract data from the mimic-IV database because mimic-IV database information was publicly available and all patient data were deidentified\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDeclaration of generative AI and AI-assisted technologies in the writing process\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDuring the preparation of this work, the authors used ChatGPT in order to polish. After using this tool or service, the authors reviewed and edited the content as needed and took full responsibility for the content of the publication.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDeclaration of competing interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe express our profound gratitude to the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center for their invaluable contribution to the MIMIC project.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePublisher\u0026rsquo;s note\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eclaim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003ePalumbo, A., et al., \u003cem\u003ePersonalized therapy in multiple myeloma according to patient age and vulnerability: a report of the European Myeloma Network (EMN).\u003c/em\u003e Blood, 2011. \u003cstrong\u003e118\u003c/strong\u003e(17): p. 4519-29.\u003c/li\u003e\n\u003cli\u003eCowan, A.J., et al., \u003cem\u003eDiagnosis and Management of Multiple Myeloma: A Review.\u003c/em\u003e Jama, 2022. \u003cstrong\u003e327\u003c/strong\u003e(5): p. 464-477.\u003c/li\u003e\n\u003cli\u003eZagouri, F., et al., \u003cem\u003eHypercalcemia remains an adverse prognostic factor for newly diagnosed multiple myeloma patients in the era of novel antimyeloma therapies.\u003c/em\u003e Eur J Haematol, 2017. \u003cstrong\u003e99\u003c/strong\u003e(5): p. 409-414.\u003c/li\u003e\n\u003cli\u003eMousseaux, C., et al., \u003cem\u003eEpidemiology, clinical features, and management of severe hypercalcemia in critically ill patients.\u003c/em\u003e Ann Intensive Care, 2019. \u003cstrong\u003e9\u003c/strong\u003e(1): p. 133.\u003c/li\u003e\n\u003cli\u003eDimopoulos, M.A., et al., \u003cem\u003eMultiple myeloma: EHA-ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up(\u0026dagger;).\u003c/em\u003e Ann Oncol, 2021. \u003cstrong\u003e32\u003c/strong\u003e(3): p. 309-322.\u003c/li\u003e\n\u003cli\u003eMao, Y., S. Zhu, and Y. Geng, \u003cem\u003eAssociation between serum calcium and in-hospital mortality in critical patients with multiple myeloma: a cohort study.\u003c/em\u003e Hematology, 2022. \u003cstrong\u003e27\u003c/strong\u003e(1): p. 795-801.\u003c/li\u003e\n\u003cli\u003eYe, Q., et al., \u003cem\u003eSerial platelet count as a dynamic prediction marker of hospital mortality among septic patients.\u003c/em\u003e Burns \u0026amp; Trauma, 2024. \u003cstrong\u003e12\u003c/strong\u003e.\u003c/li\u003e\n\u003cli\u003eSong, M., et al., \u003cem\u003eTrajectory of body shape in early and middle life and all cause and cause specific mortality: results from two prospective US cohort studies.\u003c/em\u003e Bmj, 2016.\u003c/li\u003e\n\u003cli\u003eKyheng, M., et al., \u003cem\u003eJoint latent class model: Simulation study of model properties and application to amyotrophic lateral sclerosis disease.\u003c/em\u003e BMC Med Res Methodol, 2021. \u003cstrong\u003e21\u003c/strong\u003e(1): p. 198.\u003c/li\u003e\n\u003cli\u003eSi, Y., et al., \u003cem\u003eTime series analysis between platelet counts and 60-day mortality in sepsis patients with thrombocytopenia: a retrospective cohort study.\u003c/em\u003e BMC Infectious Diseases, 2025. \u003cstrong\u003e25\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eTang, S., et al., \u003cem\u003eAssociation between red blood cell distribution width-platelet ratio (RPR) and mortality in patients with heart failure from the MIMIC IV database: A retrospective cohort study.\u003c/em\u003e Heliyon, 2024. \u003cstrong\u003e10\u003c/strong\u003e(16).\u003c/li\u003e\n\u003cli\u003eZeng, Y., et al., \u003cem\u003eAssociation between Platelet Count and In-Hospital Mortality in Critical Patients with Multiple Myeloma: A Cohort Study.\u003c/em\u003e PLOS One, 2025. \u003cstrong\u003e20\u003c/strong\u003e(6).\u003c/li\u003e\n\u003cli\u003evon Elm, E., et al., \u003cem\u003eThe Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies.\u003c/em\u003e Lancet, 2007. \u003cstrong\u003e370\u003c/strong\u003e(9596): p. 1453-7.\u003c/li\u003e\n\u003cli\u003eLi, J., et al., \u003cem\u003ePredicting Mortality in Intensive Care Unit Patients With Heart Failure Using an Interpretable Machine Learning Model: Retrospective Cohort Study.\u003c/em\u003e Journal of Medical Internet Research, 2022. \u003cstrong\u003e24\u003c/strong\u003e(8).\u003c/li\u003e\n\u003cli\u003eNagin, D.S., B.L. Jones, and J. Elmer, \u003cem\u003eRecent Advances in Group-Based Trajectory Modeling for Clinical Research.\u003c/em\u003e Annu Rev Clin Psychol, 2024. \u003cstrong\u003e20\u003c/strong\u003e(1): p. 285-305.\u003c/li\u003e\n\u003cli\u003eKim, S.Y., \u003cem\u003eDetermining the Number of Latent Classes in Single- and Multi-Phase Growth Mixture Models.\u003c/em\u003e Struct Equ Modeling, 2014. \u003cstrong\u003e21\u003c/strong\u003e(2): p. 263-279.\u003c/li\u003e\n\u003cli\u003eDing, X., et al., \u003cem\u003eAssociations between sepsis occurrence, hemoglobin level and mortality in patients with non-trauma hemorrhagic brain injuries: trajectory-based analysis.\u003c/em\u003e European Journal of Medical Research, 2025. \u003cstrong\u003e30\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eZhang, Q., et al., \u003cem\u003eThe prognostic value of the platelet-to-lymphocyte ratio in multiple myeloma patients treated with a bortezomib-based regimen.\u003c/em\u003e Scientific Reports, 2025. \u003cstrong\u003e15\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eChung, W.T. and K.C. Chung, \u003cem\u003eThe use of the E-value for sensitivity analysis.\u003c/em\u003e J Clin Epidemiol, 2023. \u003cstrong\u003e163\u003c/strong\u003e: p. 92-94.\u003c/li\u003e\n\u003cli\u003eCui, K., et al., \u003cem\u003eThe association between blood urea nitrogen to albumin ratio and the 28 day mortality in tuberculosis patients complicated by sepsis.\u003c/em\u003e Scientific Reports, 2024. \u003cstrong\u003e14\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eTothong, W., et al., \u003cem\u003ePrevalence, Outcomes and Impact of Disease-Related Complications in the Survival of Multiple Myeloma Patients.\u003c/em\u003e Hematol Rep, 2024. \u003cstrong\u003e16\u003c/strong\u003e(1): p. 89-97.\u003c/li\u003e\n\u003cli\u003eRees, M.J. and S. Kumar, \u003cem\u003eHigh-risk multiple myeloma: Redefining genetic, clinical, and functional high-risk disease in the era of molecular medicine and immunotherapy.\u003c/em\u003e Am J Hematol, 2024. \u003cstrong\u003e99\u003c/strong\u003e(8): p. 1560-1575.\u003c/li\u003e\n\u003cli\u003eHampshire, P.A., et al., \u003cem\u003eAdmission factors associated with hospital mortality in patients with haematological malignancy admitted to UK adult, general critical care units: a secondary analysis of the ICNARC Case Mix Programme Database.\u003c/em\u003e Crit Care, 2009. \u003cstrong\u003e13\u003c/strong\u003e(4): p. R137.\u003c/li\u003e\n\u003cli\u003eJeon, J.E., et al., \u003cem\u003eA Pathogen-Responsive Gene Cluster for Highly Modified Fatty Acids in Tomato.\u003c/em\u003e Cell, 2020. \u003cstrong\u003e180\u003c/strong\u003e(1): p. 176-187.e19.\u003c/li\u003e\n\u003cli\u003eWeidle, U.H., et al., \u003cem\u003eMolecular Mechanisms of Bone Metastasis.\u003c/em\u003e Cancer Genomics Proteomics, 2016. \u003cstrong\u003e13\u003c/strong\u003e(1): p. 1-12.\u003c/li\u003e\n\u003cli\u003eKanellias, N., et al., \u003cem\u003eNewly Diagnosed Multiple Myeloma Patients with Skeletal-Related Events and Abnormal MRI Pattern Have Poor Survival Outcomes: A Prospective Study on 370 Patients.\u003c/em\u003e J Clin Med, 2022. \u003cstrong\u003e11\u003c/strong\u003e(11).\u003c/li\u003e\n\u003cli\u003eFreire-de-Lima, L., et al., \u003cem\u003eMultiple Myeloma Cells Express Key Immunoregulatory Cytokines and Modulate the Monocyte Migratory Response.\u003c/em\u003e Front Med (Lausanne), 2017. \u003cstrong\u003e4\u003c/strong\u003e: p. 92.\u003c/li\u003e\n\u003cli\u003eBao, L., et al., \u003cem\u003eHypercalcemia caused by humoral effects and bone damage indicate poor outcomes in newly diagnosed multiple myeloma patients.\u003c/em\u003e Cancer Med, 2020. \u003cstrong\u003e9\u003c/strong\u003e(23): p. 8962-8969.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003eTables 1 to 5 are available in the Supplementary Files section.\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"
[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Multiple myeloma (MM), serum calcium, mortality, Latent Class Growth Modeling (LCGM), Trajectory","lastPublishedDoi":"10.21203/rs.3.rs-8812562/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8812562/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eHypercalcemia is a well-recognized adverse prognostic factor in patients with multiple myeloma (MM), particularly in critically ill subgroups where its incidence is substantially higher than in the general MM population. However, the prognostic value of dynamic serum calcium fluctuations during intensive care unit (ICU) hospitalization remains underexplored, limiting the ability to refine risk stratification and personalized management.\u003c/p\u003e\u003ch2\u003eObjective\u003c/h2\u003e \u003cp\u003eThis study aims to identify heterogeneous serum calcium trajectories in critically ill MM patients and investigate their association with 60-day mortality.\u003c/p\u003e\u003ch2\u003eMethods\u003c/h2\u003e \u003cp\u003eA retrospective cohort study was conducted using data from 254 critically ill MM patients extracted from the MIMIC-IV database. Latent Class Growth Modeling (LCGM) was applied to classify serum calcium trajectories over the first 28 days of hospitalization, with measurements collected at predefined time points-. Multivariable Cox proportional hazards regression models (adjusted for demographic, clinical, and disease severity covariates) and Kaplan-Meier (K-M) curves were used to assess the association between trajectory patterns and 60-day mortality. Sensitivity analyses included trajectory reclassification and E-value assessment to validate result robustness.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eFour distinct serum calcium trajectories were identified: low-level (n\u0026thinsp;=\u0026thinsp;90, 35.43%), medium-level (n\u0026thinsp;=\u0026thinsp;138, 57.48%), U-shaped (n\u0026thinsp;=\u0026thinsp;9, 3.54%), and high-level (n\u0026thinsp;=\u0026thinsp;17, 6.68%). Compared with the high-level trajectory (reference group), the low-level (HR\u0026thinsp;=\u0026thinsp;0.17; 95% CI 0.07\u0026ndash;0.41; P\u0026thinsp;\u0026lt;\u0026thinsp;0.001) and medium-level (HR\u0026thinsp;=\u0026thinsp;0.26; 95% CI 0.12\u0026ndash;0.56; P\u0026thinsp;=\u0026thinsp;0.001) trajectories were independently associated with significantly lower 60-day mortality. The U-shaped trajectory showed no significant mortality difference from the high-level trajectory (HR\u0026thinsp;=\u0026thinsp;0.56; 95% CI 0.14\u0026ndash;2.25; P\u0026thinsp;=\u0026thinsp;0.416). K-M analysis confirmed the prognostic gradient across trajectories, and sensitivity analysis (E-value\u0026thinsp;=\u0026thinsp;2.49) validated the robustness of these associations to unmeasured confounding.\u003c/p\u003e\u003ch2\u003eConclusion\u003c/h2\u003e \u003cp\u003ePersistently high or rebound (U-shaped) serum calcium trajectories are independent predictors of increased 60-day mortality in critically ill MM patients, reflecting underlying aggressive disease biology or treatment resistance. Dynamic serum calcium monitoring via trajectory analysis refines early risk stratification and provides actionable insights for guiding personalized ICU management, addressing an unmet need in current clinical practice.\u003c/p\u003e","manuscriptTitle":"Time series analysis Between Serum Calcium Levels and 60-Day Mortality in Critically ill Patients with Multiple Myeloma: A Retrospective Cohort Study","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-03-08 16:50:52","doi":"10.21203/rs.3.rs-8812562/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"
[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"52bd5a09-9823-4a33-bbeb-76fd1c134206","owner":[],"postedDate":"March 8th, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2026-03-23T14:41:31+00:00","versionOfRecord":[],"versionCreatedAt":"2026-03-08 16:50:52","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-8812562","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8812562","identity":"rs-8812562","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.