Abstract

Background. No epidemiological study has investigated the effect of anion gap (AG) on the prognosis of critically ill patients with acute kidney injury (AKI). Therefore, we aimed to determine the association between serum AG and all-cause mortality in these patients. Methods. From MIMIC III, we extracted demographics, vital signs, laboratory tests, comorbidities, and scoring systems from the first 24 h after patient ICU admission. A generalized additive model was used to identify a nonlinear association between anion gap and 30-day all-cause mortality. We also used the Cox proportional hazards models to measure the association between AG levels and 30-day, 90-day, and 365-day mortality in patients with AKI. Results. A total of 11,573 eligible subjects were extracted from the MIMIC-III. The relationship between AG levels and 30-day all-cause mortality in patients with AKI was nonlinear, with a U-shaped curve. In multivariate analysis, after adjusting for potential confounders, higher AG was a significant predictor of 30-day, 90-day, and 365-day all-cause mortality compared with lower AG (HR, 95% CI: 1.54, 1.33–1.75; 1.55, 1.38–1.73; 1.46, 1.31–1.60). Conclusions. The relationship between AG levels and 30-day all-cause mortality described a U-shaped curve. High-AG levels were associated with increased risk 30-day, 90-day, and 365-day all-cause mortality in critically ill patients with AKI.

1. Introduction

Acute kidney injury (AKI) is a common syndrome characterized by an abrupt, usually reversible decline in glomerular filtration, associated with substantial morbidity and mortality, especially in critically ill patients [1, 2]. In the US, approximately 20% of critically ill patients have AKI in the intensive care unit (ICU) [3]. In the presence of AKI, patient mortality within 1 year after ICU admission is significantly elevated, as high as 60%–70% [4, 5]. Because of the high incidence and poor prognosis of AKI in critically ill patients, researchers are attempting to identify prognostic predictors in AKI [6, 7]. Unfortunately, most of them are not widely used in clinical practice.

The serum anion gap (AG) is a mathematical derivation parameter calculated from the difference in serum cation and anion concentrations. It is the simplest way to assess acid-base status [8, 9] and helps to identify various forms of metabolic acidosis. AG is also one of the most commonly used biomarkers that provides important clues regarding diagnosis or prognosis of various disorders [1012]. There is a positive monotonic relationship between high AG and severity of illness or poor prognosis for sepsis [13], coronary artery disease (CAD) [14], aortic aneurysm [15], and chronic kidney disease (CKD) [16]. To the best of our knowledge, no epidemiological study has investigated the effect of AG on the prognosis of critically ill patients with AKI. Therefore, we aimed to determine the association between serum AG and all-cause mortality in these patients.

2. Methods

2.1. Data Source

We followed the methods of Wang et al. in this study, as we have done previously [1719]. Our study is based on an openly available clinical database called the Multiparameter Intelligent Monitoring in Intensive Care III (MIMIC III) [20]. The database includes more than 40,000 ICU patients admitted to Beth Israel Deaconess Medical Center (Boston, MA, USA) from 2001 to 2012. To apply for access to the database, we completed the National Institutes of Health’s web-based course and passed the Protecting Human Research Participants exam (no. 6182750). We extracted clinical variables, including demographic characteristics, International Classification of Diseases (ICD-9) codes, physiological index, medications, and laboratory tests. The project was approved by the Institutional Review Boards of Beth Israel Deaconess Medical Center and the Massachusetts Institute of Technology (Cambridge, MA). To safeguard patient privacy, data were deidentified; therefore, informed consent was waived.

2.2. Population Selection Criteria

Adult patients (≥18 years) with AKI according to ICD-9 code at first ICU admission for more than two days were included. The exclusion criteria were (1) missing AG at ICU admission and (2) missing >5% individual data.

2.3. Data Extraction

Structured Query Language (SQL) with PostgreSQL (version 9.6) was used to extract data from MIMIC-III. Demographics, vital signs, laboratory tests, comorbidities, scoring systems, and other variables collected within the first 24 h of ICU admission were extracted from MIMIC III. The comorbidities included CAD, congestive heart failure (CHF), atrial fibrillation (AFIB), stroke, renal disease, liver disease, pneumonia, respiratory failure, and acute respiratory distress syndrome (ARDS). Laboratory data were also extracted, including AG, albumin, bicarbonate, bilirubin, creatinine, chloride, glucose, hematocrit, hemoglobin, platelet, sodium, potassium, lactate, blood urea nitrogen (BUN), white blood cell (WBC), prothrombin time (PT), activated partial thromboplastin time (APTT), and international normalized ratio (INR). We also calculated the sequential organ failure assessment (SOFA) score [21] and simplified acute physiology score II (SAPSII) [22] for each patient. Other extracted data included age, gender, systolic blood pressure (SBP), diastolic blood pressure (DBP), mean blood pressure (MBP), heart rate, respiratory rate, temperature, SPO2, AKI stage, renal replacement therapy, and ICU length of stay (LOS). Survival information on vital status was obtained from Social Security Death Index records. The endpoints of our study were 30-day, 90-day, and 365-day all-cause mortality from the date of ICU admission.

2.4. Statistical Analysis

Continuous variables were expressed as or medians and interquartile range (IQR). Categorical data were expressed as frequencies or percentages. Chi-square, 1-way ANOVA, and Kruskal-Wallis tests were used to determine any significant differences between the groups. A generalized additive model (GAM) was used to identify the nonlinear association between AG and 30-day all-cause mortality. We also used the Cox proportional hazards models to determine associations between AG levels and 30-day, 90-day, and 365-day mortality in AKI; these results were expressed as hazard ratios (HRs) with 95% confidence intervals (CIs).

Variables based on epidemiological and biological background were incorporated as potential confounders, and those confounders based on a change in effect estimate of >10% were used to generate an adjusted model [23]. For each endpoint, two multivariate models were constructed on the basis of AG group inclusion according to tertiles. The first tertile was treated as the reference group. In model I, covariates were adjusted for age and gender. In model II, we further adjusted for age, gender, AKI stage, CHF, CAD, liver disease, stroke, respiratory failure, pneumonia, SIRS, potassium, albumin, lactate, platelet, BUN, PT, INR, APTT, WBC, pH, creatinine, bicarbonate, bilirubin, renal replacement therapy, respiration rate, SPO2, heart rate, SBP, DBP, temperature, Elixhauser comorbidity index, SOFA, and SAPSII. Subgroup analysis of the associations between AG and 30-day all-cause mortality was performed using stratified linear regression models. All probability values were 2-sided, and values less than 0.05 were considered statistically significant. R (http://www.R-project.org) and EmpowerStats (http://www.empowerstats.com/en/, X&Y solutions, Inc., Boston, MA) were used for all statistical analysis.

3. Results

3.1. Subject Characteristics

The patients were divided according to AG in tertiles, and baseline characteristics of these patients are summarized in Table 1. A total of 11,573 eligible subjects were extracted from the MIMIC-III. There were 6,626 men and 4,947 women, and patients were generally older. Patients with higher AG () were more likely to report a history of CAD, CHF, AFIB, and renal disease and had higher values of bilirubin, creatinine, potassium, lactate, BUN, WBC, PT, and APTT. SAPSII, SOFA scores, mortality, use of RRT, and ICU LOS were also significantly higher in the high-AG group () than in the lower-AG group (AG< 12).

3.2. AG Levels and All-Cause Mortality

The relationship between AG levels and 30-day all-cause mortality was nonlinear, and a U-shaped curve was observed (Figure 1). We used Cox proportional hazards regression model to determine the associations between AG and 30-day, 90-day, and 365-day all-cause mortality in patients with AKI (Table 2). In model I, high-AG levels () were associated with increased risk of all-cause mortality after adjustment for age and gender. In model II, the lower AG () was used as a reference. After adjustment for confounders (age, gender, acute kidney injury stage, congestive heart failure, coronary artery disease, liver disease, stroke, respiratory failure, pneumonia, SIRS, potassium, albumin, lactate, platelet, BUN, PT, INR, APTT, WBC, pH, creatinine, bicarbonate, sodium, chloride, diabetes, bilirubin, renal replacement therapy, respiration rate, SPO2, heart rate, SBP, DBP, temperature, Elixhauser comorbidity index, SOFA, and SAPSII), higher AG () remained a more significant predictor of 30-day, 90-day, and 365-day all-cause mortality than lower AG (HR, 95% CI: 1.54, 1.33–1.75; 1.55, 1.38–1.73; 1.46, 1.31–1.60). For the purpose of sensitivity analysis, we also handled AG as categorical variable (tertile) and found the same trend ( for trend: <0.0001). In addition, in order to verify that AG was an independent prognostic factor for AKI, we also analyzed the potential associations of bicarbonate, pH, lactate, and urine ketone bodies on all-cause mortality, and the results were included in supplementary material (available here).

3.3. Subgroup Analyses

As shown in Table 3, the test for interactions was statistically significant in several strata ( for interaction <0.05). Among these strata, we observed that patients with higher AGs had significantly higher mortality with hypotension, , , , , , , , , and . Similar trends were observed in patients with CAD and liver disease.

4. Discussion

The relationship between AG and 30-day all-cause mortality among critically ill patients with AKI was nonlinear, and a U-shaped curve was observed. In the fully adjusted model, high-AG levels were associated with increased risk 30-day, 90-day, and 365-day all-cause mortality. To our knowledge, this was the first study to measure the association between serum AG and all-cause mortality in critically ill patients with AKI.

Several studies have explored the relationship between AG and clinical outcomes of various diseases. Yang et al. [14] measured the association between the AG and all-cause mortality in CAD and found that higher AG was associated with worse cardiac function and was a significant predictor of all-cause mortality. Banerjee et al. [24] suggested that greater AG was present among persons with CKD, and AG increased the risk for progression to end-stage renal disease in these patients. AG is a traditional tool for assessing acid-base status, and most previous studies have associated it with acid-base disorders, all of which have significant impacts on morbidity and mortality in critically ill patients [25]. Similarly, our findings showed a positive correlation between serum AG and all-cause mortality in critically ill patients with AKI. Nevertheless, the underlying mechanism requires further research.

AG reflects the unmeasured anion concentration and can be easily calculated from conventional clinical chemical analysis. It is widely used to evaluate the acid-base status and is one of the most commonly used biomarkers that provides important clues for the diagnosis and prognosis of various diseases [13, 15]. Elevated serum AG is usually caused by overproduction of organic acid anions and/or reduction in anion excretion [8]; increased serum lactate and ketoanions are the main reasons for increased AG [26]. AKI is defined as rapid decline in renal function lasting from hours to days, as opposed to chronic kidney disease [27]. Critically ill patients are often exposed to hypoxia and anaerobic tissue conditions, leading to rapid accumulation of pyruvate, which is almost completely converted to lactate [28]. The kidney is known for lactate dehydrogenase dysfunction in the setting of patients with AKI [29], and the renal acid excretion does not fully offset endogenous acid production [30]. Consequently, high AG is common in critically ill patients with AKI.

Similar statistical methods were used to analyze the relationship between lactate, bicarbonate, pH, and urine ketone bodies on the prognosis of AKI patients, as shown in supplementary material. In the fully adjusted model, lactate levels were associated with increased risk 30-day and 90-day all-cause mortality among critically ill patients with AKI, but not with 365-day mortality. As we all know, lactate is closely related to the prognosis of critically ill patients, especially the risk of short-term death [31, 32], which can explain the correlation between lactate and short-term prognosis of critically ill patients with AKI. Moreover, bicarbonate, pH, and urine ketone bodies were not independently associated with 30-day, 90-day, and 365-day all-cause mortality. Interestingly, our results further suggested a positive correlation between AG and all-cause mortality in these patients after adjusting for potential confounders such as lactate, pH, and bicarbonate.

There are several limitations in the present study. First, this was a single-center retrospective observational study, and selection bias was inevitable. Second, we measured serum AG in patients only upon admission to the ICU and did not have laboratory follow-up data. There is the possibility of misclassified measured data that may influence the summary results. Third, although we adjusted for confounding factors, our results may have been influenced by other unknown factors. Finally, we could not determine the underlying mechanism between higher AG and poor prognosis; therefore, further study regarding the mechanism is necessary.

5. Conclusions

We found that the relationship between AG levels and 30-day all-cause mortality was nonlinear, with a U-shaped curve. High-AG levels were associated with increased risk 30-day, 90-day, and 365-day all-cause mortality in critically ill patients with AKI.

Data Availability

The clinical data used to support the findings of this study were supplied by Monitoring in Intensive Care Database III version 1.4 (MIMIC-III v.1.4). Although the database is publicly and freely available, researchers must complete the National Institutes of Health’s web-based course known as Protecting Human Research Participants to apply for permission to access the database.

Disclosure

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

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This research was supported by the Scientific Research Foundation of Wenzhou (grant nos. Y20180515 and Y20150038), the Zhejiang Provincial Natural Science Foundation of China (grant nos. LY19H150002 and LY13H150007), and the Clinical Research Foundation of the 2nd Affiliated Hospital of Wenzhou Medical University (SAHoWMU-CR2019-11-423).

Supplementary Materials

Table 1: relationship between bicarbonate and all-cause mortality in different models. Table 2: relationship between pH and all-cause mortality in different models. Table 3: relationship between lactate and all-cause mortality in different models. Table 4: relationship between urine ketone bodies and all-cause mortality in different models. (Supplementary Materials)