Data source and study design
A retrospective cohort study was performed using data sourced from version 2.2 of the Medical Information Mart for Intensive Care IV database. This database comprises two in-house systems: a comprehensive hospital-wide electronic health record (EHR) and an ICU-specific clinical information system, covering data from 2008 to 202418. This study was performed in compliance with institutional guidelines. The MIMIC-IV database received approval from the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center. This study was reviewed and approved by the Human Research Ethics Committee of Ningbo No. 2 Hospital (approval number: YJ-NBEY-KY-2022-107-01, dated September 7, 2022). Due to the retrospective design and the anonymized nature of the patient health information, the same ethics committee waived the requirement for informed consent under the same approval (YJ-NBEY-KY-2022-107-01). Access to the database was secured by one of the authors (JQ L), who received the necessary authentication and completed the Collaborative Institutional Training Initiative examination (authentication number 60691748).
Given the retrospective nature of this study, traditional sample size calculations or power analyses were not performed. However, the volume of available records in the database was deemed adequate to meet the study’s objectives. This assessment is based on investigating the relationship between the albumin-corrected anion gap (ACAG) and mortality outcomes within a large and diverse patient population. The convenience samples drawn from the database effectively represent the intensive care population, thus supporting the statistical analyses employed in this study.
Participants
The study encompassed all sepsis patients from the MIMIC-IV v2.2 database. Sepsis was defined according to the Sepsis-3.0 criteria, implemented using structured query language (SQL) queries to identify patients with suspected infection and a sequential organ failure assessment (SOFA) score change of ≥2 within 24 h19. All patients in this cohort were admitted to the ICU as the first point of hospital entry; patients admitted to general wards and later transferred to the ICU were excluded to ensure a homogeneous population with early ICU admission. Data extraction was performed using SQL via a database management tool20, with all utilized software, databases, and analytical tools listed in an alphabetically sorted Table of Materials. Inclusion criteria specified sepsis patients aged over 18 years who were admitted to the ICU for the first time. The exclusion criteria included: (1) patients under 18 years old; (2) ICU stays of less than 48 h; (3) multiple ICU admissions due to sepsis; and (4) insufficient data, such as missing records for albumin or anion gap (Figure 1).
Research procedures and definitions
Data extraction from MIMIC-IV was performed using SQL queries. The extracted data encompassed patient demographics (age, height, weight, gender, insurance status, race, and marital status), medical history (including hypertension, type 2 diabetes, heart failure, myocardial infarction, malignant tumors, chronic kidney disease, cirrhosis, hepatitis, pneumonia, chronic obstructive pulmonary disease, and hyperlipidemia), and initial laboratory test results. These results included white and red blood cell counts, sodium, potassium, calcium, chloride, glucose, anion gap, blood pH, thrombin time, fibrinogen, partial thromboplastin time, international normalized ratio, total cholesterol, aspartate aminotransferase, and alanine aminotransferase. Additionally, information on special treatments (such as continuous renal replacement therapy and mechanical ventilation), clinical scores (SOFA, acute physiology and chronic health evaluation III [APACHE III], simplified acute physiology score II [SAPS II], Oxford acute severity of illness score [Oasis], Charlson comorbidity index [Charlson], and Glasgow coma scale [GCS] scores), and clinical outcomes (length of hospital stay, in-hospital mortality, ICU stay, and ICU mortality) was collected. During data cleaning, predictors with >30% missing data were excluded.
The albumin-corrected anion gap (ACAG) was calculated using the first available measurements of serum anion gap and serum albumin obtained within 24 h of ICU admission. If multiple measurements were available within the first 24 h, the values closest to the admission time were used. The formula for ACAG is:
Albumin - corrected anion gap = Anion gap + [40 - Serum albumin (g/dL)] × 0.2521 (1)
Outcomes and Measures
The primary outcomes of this study were 30-day in-hospital mortality and 30-day ICU mortality. The index date for survival analysis was the time of ICU admission. Patients were followed for 30 days from this index date, regardless of hospital discharge status. Mortality events occurring after discharge were captured using the social security death index (SSDI) records integrated within the MIMIC-IV database, ensuring accurate ascertainment of out-of-hospital deaths.
Statistical analysis
Continuous variables were expressed as mean (standard deviation) or median (interquartile range), while categorical variables were presented as frequency and percentage. For data adhering to a normal distribution, t-tests or analysis of variance (ANOVA) were utilized. Data not following a normal distribution were analyzed using the Mann-Whitney U test or Kruskal-Wallis test. Kaplan-Meier survival analysis was employed to assess the incidence of endpoint events across different levels of ACAG, with differences evaluated using the log-rank test (Figure 2).
The Cox proportional hazards model was used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs) for the association between ACAG and mortality. ACAG was analyzed both as a continuous variable and by quartiles. Model adjustments were made as follows: model 1 involved univariate analysis; model 2 adjusted for age, gender, height, weight, and race; and model 3 further adjusted for insurance status, marital status, WBC and RBC counts, RDW, albumin, chloride, ALT, AST, SOFA score, APACHE III score, SAPS II score, Oasis score, Charlson score, and a range of comorbidities including hypertension, type 2 diabetes mellitus, heart failure, myocardial infarction, malignant tumors, chronic kidney disease, acute renal failure, cirrhosis, hepatitis, pneumonia, stroke, hyperlipidemia, acute kidney injury, and chronic obstructive pulmonary disease. Collinearity among covariates in model 3 was assessed using the variance inflation factor (VIF), with VIF > 5 indicating significant collinearity. No significant collinearity was detected (all VIF < 3).
A restricted cubic spline (RCS) regression model with four knots was applied to examine the nonlinear relationship between baseline ACAG and mortality in both hospital and ICU settings. The RCS analysis specifically tested for an inflection point where the slope of the mortality risk curve changes. Subgroup analyses were conducted to explore potential differences across strata based on age (≤70 years vs. >70 years), sex, body mass index (BMI; <21.4, 21.4–31.2, ≥31.2 kg/m2), presence of hypertension, type 2 diabetes, myocardial infarction, chronic kidney disease, stroke, acute kidney injury, as well as use of continuous renal replacement therapy (CRRT) and mechanical ventilation. These analyses evaluated the consistency of ACAG’s prognostic value for the primary outcomes. Cox models were also used in subgroup analyses to adjust for all baseline variables. Interaction p-values were calculated to assess effect modification. Data processing and analysis were performed using R version 4.3.0 (RRID: SCR_001905), with statistical significance set at two-tailed p < 0.05. For handling missing values, multiple imputation by chained equations (MICE) was applied using the R package “mice” under the missing at random assumption. Variables with missing values exceeding 30% were excluded prior to imputation.