- 1Department of Cardiology, Qinhuangdao First Hospital, Qinhuangdao, Hebei, China
- 2Medical Department, Qinhuangdao First Hospital, Qinhuangdao, Hebei, China
- 3Department of Cardiology, Tianjin Medical University General Hospital, Tianjin, China
Background: Acute kidney injury (AKI) following acute myocardial infarction (AMI) notably affects patient outcomes. The impact of KDIGO AKI staging on post-discharge short- and long-term outcomes, particularly early-stage AKI, is not well understood. This study evaluates the prognostic implications of various KDIGO stages in AMI patients.
Methods: Utilizing the Medical Information Mart for Intensive Care IV (version 3.0) database, this retrospective cohort study included adult patients primarily diagnosed with AMI. Statistical analyses, including doubly robust estimation, propensity score matching, logistic regression, and Cox regression, were performed. The study compared Non-AKI (KDIGO stage 0) with Mild-AKI (maximum KDIGO stage 1 during hospitalization), and Normal-or-mild AKI (KDIGO stages 0–1) with Moderate-to-severe AKI (KDIGO stages 2–3).
Results: Among 5,715 patients analyzed, 4,306 (75.36%) developed AKI. Doubly robust analysis revealed no significant differences in outcomes between Non-AKI and Mild-AKI groups (28-day mortality: OR 0.97, 95% CI 0.68–1.38; 180-day mortality: HR 0.94, 95% CI 0.76–1.18; 1-year mortality: HR 0.98, 95% CI 0.81–1.20). However, Moderate-to-severe AKI was significantly associated with worse outcomes compared to Normal-or-mild AKI (28-day mortality: OR 1.67, 95% CI 1.36–2.05; 180-day mortality: HR 1.06, 95% CI 1.02–1.10; 1-year mortality: HR 1.22, 95% CI 1.07–1.38; all p < 0.001). Subgroup analyses revealed that patients under 65 years with Mild-AKI showed higher risks of 180-day and 1-year mortality compared to Non-AKI, while Moderate-to-severe AKI consistently demonstrated worse outcomes across all subgroups (age, SOFA score, heart failure status, and renal disease status). These findings were robust across multiple sensitivity analyses.
Conclusions: Patients with Mild-AKI can be considered as having “subclinical AKI,” with prognoses similar to Non-AKI patients. In contrast, Moderate-to-severe AKI significantly worsens prognosis compared to Normal-or-mild AKI.
Introduction
In recent years, growing attention has been directed toward the interplay between cardiovascular conditions, including acute heart failure, acute myocardial infarction (AMI), and cardiovascular surgery, and the onset of acute kidney injury (AKI). This is due to the profound impact AKI has on patient outcomes and prognosis (1). Research has suggested a strong correlation between cardiac and renal function, giving rise to the term “cardiorenal syndrome (CRS)” in the context of heart failure and acute coronary syndrome (ACS) (2).
The pathophysiology of CRS is characterized by a complex interplay of hemodynamic and non-hemodynamic factors that cause mutual cardiac and renal damage. Key contributors include common risk factors such as hypertension, diabetes, atherosclerosis, and chronic inflammation, which drive disease progression (3). Hemodynamic disturbances like venous congestion and increased intra-abdominal pressure reduce renal blood flow, impair glomerular filtration, and activate the renin-angiotensin system (RAAS), worsening renal function (4, 5). Non-hemodynamic mechanisms involve neurohormonal dysregulation, oxidative stress, and inflammation, which contribute to chronic renal hypoxia, tissue injury, and fibrosis (6–8). Inflammatory mediators such as TNF-α, IL-1, and IL-6 play crucial roles, leading to both cardiac and renal remodeling (9). Endothelial dysfunction further exacerbates the cycle of damage by impairing vasodilation, increasing vascular permeability, and promoting thrombosis and atherosclerosis (10). Together, these mechanisms create a self-perpetuating cycle of organ dysfunction, contributing to the progression of CRS (2).
The occurrence of AKI following AMI significantly prolongs hospital stay, increases medical costs, and elevates both short- and long-term mortality. Reported incidence rates of AKI after AMI range from 5.2% to 59% across studies, primarily due to variations in the criteria used to define AKI and differences in study populations (11).
Currently, the most widely used criteria for diagnosing AKI include the RIFLE criteria (12), AKIN criteria (a later version of the RIFLE classification) (13), and KDIGO guidelines (14). Researchers applied both the KDIGO and RIFLE criteria to AMI patients, revealing that KDIGO detects AKI more effectively than RIFLE, with detection rates of 36.6% vs. 14.8% (15). KDIGO integrates elements from both RIFLE and AKIN, combining their strengths to offer a more standardized and comprehensive definition of AKI, thereby minimizing discrepancies between diagnostic frameworks (11). KDIGO is notably more sensitive in detecting AKI, particularly in the early stages (16). Its three-stage classification simplifies clinical application while preserving the diagnostic accuracy of more complex systems like RIFLE. Although AKI diagnosis relies on acute increases in serum creatinine (SC) and reduced urine output (UO), UO measurement is underutilized in clinical practice, despite evidence suggesting its significant diagnostic and prognostic value (17). The use of the more sensitive KDIGO criteria, which incorporate UO, may thus provide an advantage in detecting AKI in patients with AMI. However, few studies have examined the impact of KDIGO staging on short- and long-term outcomes in AMI patients after discharge, and it remains unclear whether even early-stage AKI affects prognosis. Therefore, this study aims to comprehensively analyze outcomes in AMI patients across different KDIGO AKI stages to reveal the prognostic implications of each stage.
Methods
Study design
This study conforms to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines, as outlined in the Supplementary materials. It aims to investigate the short-term and long-term impacts of mild and moderate-to-severe AKI on the prognosis of ICU patients with acute myocardial infarction, utilizing real-time monitoring of AKI based on KDIGO criteria. The KDIGO criteria model in MIMIC-IV dynamically evaluated AKI stages through serum creatinine changes over the past seven days and 48 h, alongside hourly urine output monitored over 6, 12, and 24-h intervals. This approach improved the sensitivity of AKI assessment, facilitating earlier detection and more precise classification of kidney injury. The project received approval from the institutional review boards at both the Massachusetts Institute of Technology (MIT) and Beth Israel Deaconess Medical Center (BIDMC), with informed consent being waived.
This retrospective observational study utilized data from the Medical Information Mart for Intensive Care IV (MIMIC-IV, version 3.0) database. This updated version of MIMIC-III includes critical care information for ICU patients at BIDMC from 2008 to 2022. The database contains comprehensive records from patient hospitalizations, such as laboratory tests, medications given, vital signs, and other details. Author PG gained access to the database after fulfilling the data usage agreement and obtaining Collaborative Institutional Training Initiative (CITI) certification. Since all patient information is de-identified, informed consent was not necessary (17).
Study population
Inclusion criteria: (1) patients aged 18 years or older; (2) AMI listed among the top three discharge diagnoses. Exclusion criteria: (1) not a first hospitalization; (2) absence of ICU records; (3) ICU stay time < 1 day.
The study is divided into two parts. In the first part, patients classified under KDIGO AKI stages 0 and 1 were grouped as Non-AKI and Mild-AKI. In the second part, all patients were included, with those in KDIGO AKI stages 0 and 1 categorized as Normal-or-mild-AKI, while those in stages 2 and 3 were classified as Moderate-to-severe-AKI.
Data extraction and preprocessing
Data extraction was performed using PostgreSQL 14 and SQL queries (Berkeley, California, USA). The dataset extracted included demographics, ICU length of stay, complications, laboratory test results, treatments, and other pertinent clinical information. Laboratory results were taken from the first tests conducted upon ICU admission, as these initial values are available quicker and support timely patient assessment with clinical prediction models. The estimated glomerular filtration rate (eGFR) was calculated using the Cockcroft-Gault (CG) equation (18): eGFR = 175 × (standardized serum creatinine)−1.154 × (age)−0.203 × 1.212 (if Black) × 0.742 (if female).
Endpoint
The endpoints were 28-day mortality, 180-day mortality, and 1-year all-cause mortality.
Statistical methods
The normality of continuous variables was evaluated using the Kolmogorov-Smirnov test. Variables with a normal distribution were presented as mean ± standard deviation, whereas variables that did not follow a normal distribution were expressed as median and interquartile range (IQR) [M (P25, P75)]. The homogeneity of variances for continuous variables across groups was evaluated using Levene's test. For comparisons between two cohorts, continuous variables that followed an independent normal distribution and demonstrated homogeneity of variances were analyzed using Student's t-test. If these assumptions were not met, the Mann-Whitney U test was used to assess differences between groups. For categorical variables, Fisher's exact test was applied when the sample size was < 40. Otherwise, the Chi-square test was used to assess differences between groups. Categorical data were presented as frequencies and percentages. Multiple imputation was performed using the ‘mice' package in R for variables with missing data. Variables with more than 20% missing values were excluded from imputation and not included in model construction. To ensure robust imputation results, the number of imputations was set to 100.
The doubly robust estimation approach was utilized to determine the independent associations between the occurrence of AKI in patients with myocardial infarction and their prognosis. This method combines outcome modeling and propensity score weighting to provide reliable estimates, even if one of the models (outcome or propensity score) is misspecified. This method, also known as survey-weighted generalized linear models, amalgamates a multivariate regression model with a propensity score model to evaluate both the correlation and the causal influence of an exposure on an outcome (19, 20). Typically, unbiased estimation of causal effects using either a regression model or a propensity score model individually is possible only when the respective statistical model is correctly specified. In contrast, the doubly robust estimator combines these two approaches, ensuring that an unbiased effect estimate can still be obtained if at least one of the models is correctly specified.
The gradient boosted model (GBM) was applied to estimate propensity scores for AKI, with the aim of minimizing covariate imbalance between the Non-AKI and Mild-AKI groups, as well as the Normal-or-mild-AKI and Moderate-to-severe-AKI groups. GBM, a machine learning algorithm, iteratively builds and combines models into an ensemble to enhance the accuracy of response variable estimates. Its main principle involves constructing new models that are highly correlated with the negative gradient of a predefined loss function. In this study, regression trees were used as the base learners for the GBM, incorporating 39 covariates in total (21).
An inverse probability of treatment weighting (IPTW) approach was applied to construct a weighted cohort, utilizing the estimated propensity scores as weights. To evaluate the performance of the propensity score model in achieving balance between the groups, covariate imbalances were analyzed for both the unadjusted and weighted cohorts. Standardized mean differences (SMDs) were computed to measure discrepancies between the groups. Subsequently, logistic regression or Cox regression was conducted on this weighted cohort, adjusting for variables that remained unbalanced between groups in the propensity score model. This approach is referred to as a doubly robust analysis using ‘survey' package. Logistic regression analyses utilized the ‘stats' package. In our study, the survival package was employed to fit Cox proportional hazards models and assess the proportional hazards (PH) assumption. For time-dependent covariates that violated the PH assumption, appropriate transformations, such as time-dependent covariate effects or stratification, were applied. These adjustments allowed for more accurate estimation of hazard ratios and improved overall model fit. The survival package provided essential functions for testing the PH assumption (e.g., cox.zph) and incorporating time-dependent effects (e.g., coxph with time-dependent covariates).
Statistical analyses were conducted using R software (version 4.4.1; R Foundation for Statistical Computing, Vienna, Austria). All tests were two-tailed, with a significance level set at P < 0.05.
Sensitivity analysis
We conducted a series of sensitivity analyses to evaluate the robustness of the study's findings and to determine how our conclusions might be influenced by using different association inference models. In these analyses, we applied additional models. For the outcome of 28-day mortality, we used a Log-rank test model, a Multivariate Several models were utilized in the analysis, including a logistic regression model adjusted for all covariates, a multivariate logistic regression model adjusted for unbalanced covariates, a survey-weighted Generalized Linear Model (GLM) incorporating IPTW and adjusted for all covariates, and a survey-weighted GLM with IPTW adjusted for unbalanced covariates. For the outcomes of 180-day and 1-year mortality, the analysis employed a log-rank test, a multivariate Cox proportional hazards model adjusted for all covariates, a multivariate Cox model adjusted for unbalanced covariates, a survey-weighted Cox model with IPTW adjusted for all covariates, and a survey-weighted Cox model with IPTW adjusted for unbalanced covariates. The effect sizes and corresponding p-values derived from these models were reported and compared.
Results
Baseline characteristics
A total of 5,715 patients were included in this study, as illustrated in Figure 1. In the first phase, 1,409 patients were categorized into the Non-AKI group and 1,175 into the Mild-AKI group. After PSM, both groups comprised 656 patients, as detailed in Table 1 and Supplementary Tables S1–S4; In the second phase, 2,584 patients were classified into the Normal-or-mild-AKI group, while 3,131 were assigned to the Moderate-to-severe-AKI group. Following PSM, both groups contained 1,507 patients, as shown in Table 2. In total, 4,306 patients, accounting for 75.36% of the cohort, developed AKI.

Figure 1. Flow chart of inclusion and exclusion criteria for the target study population of patients with acute myocardial infarction who developed acute kidney injury.

Table 1. Baseline characteristics before and after propensity score matching of the Non-AKI and Mild-AKI cohorts.

Table 2. Baseline characteristics before and after propensity score matching of the Normal-or-mild-AKI and Moderate-to-severe-AKI cohorts.
Doubly robust analysis
A propensity score model was initially developed using 39 covariates through GBM. Figure 2 illustrates the relative contributions of each covariate to the calculated propensity scores. Figure 2A highlights that the most significant covariates distinguishing the Non-AKI and Mild-AKI groups include the use of loop diuretics, eGFR, SOFA score, BUN, and vasopressor use, all closely associated with the onset of AKI; Figure 2B shows that the key covariates differentiating the Normal-or-mild-AKI group from the Moderate-to-severe-AKI group are vasopressor use, SOFA score, loop diuretics, eGFR, and WBC, all of which are strongly linked to the progression to Moderate-to-severe AKI.

Figure 2. The relative influence factors evaluate the discriminative power of the 39 covariates within the propensity score model in predicting acute kidney injury (AKI) outcomes (the higher the value in the bar graph, the more important the variable). (A) For Non-AKI and Mild-AKI; (B) For Normal-or-mild-AKI and Moderate-to-severe-AKI; Change in standardized mean difference (SMD) of cohorts before and after propensity score matching: the red curve represents pre-matching, the green curve represents post-matching, and the blue curve represents inverse probability of treatment weighting (IPTW) adjustment. The smaller the curve fluctuation, the better the data quality. (C) For Non-AKI and Mild-AKI; (D) For Normal-or-mild-AKI and Moderate-to-severe-AKI.
Using the estimated propensity scores, IPTW was applied to standardize differences between the Non-AKI and Mild-AKI groups, as well as between the Normal-or-mild-AKI and Moderate-to-severe-AKI groups. Details are presented in Table 1 and Figures 2C, D. In the first analysis, most covariates in the weighted cohorts were comparable or balanced between the Non-AKI and Mild-AKI groups, with some exceptions: SOFA score, loop diuretics, vasopressor use, renal disease, hemoglobin, sodium, potassium, BUN, creatinine, and eGFR; In the second analysis, most covariates were similarly balanced between the Normal-or-mild-AKI and Moderate-to-severe-AKI groups, with exceptions for SOFA score, loop diuretics, vasopressor use, heart failure, renal disease, WBC, bicarbonate, BUN, creatinine, and eGFR.
To address the residual imbalance in covariates within the weighted cohorts, several regression models were constructed using doubly robust estimation.
Outcomes and sensitivity studies
The doubly robust analysis revealed no significant differences in short- or long-term outcomes between patients with Non-AKI and Mild-AKI (28-day mortality: OR 0.97, 95% CI 0.68–1.38, p = 0.854; 180-day mortality: HR 0.94, 95% CI 0.76–1.18, p = 0.618; 1-year mortality: HR 0.98, 95% CI 0.81–1.20, p = 0.857); However, when comparing the Normal-or-mild-AKI group with the Moderate-to-severe-AKI group, patients with Moderate-to-severe-AKI had a significantly worse prognosis (28-day mortality: OR 1.67, 95% CI 1.36–2.05, p < 0.001; 180-day mortality: HR 1.06, 95% CI 1.02–1.10, p < 0.001; 1-year mortality: HR 1.22, 95% CI 1.07–1.38, p < 0.001). As shown in Table 3, Supplementary Tables S5–S28 and Figure 3, sensitivity analyses consistently confirmed these findings across all estimation models.

Figure 3. Unadjusted Kaplan-Meier survival curves for 28-day mortality (A, D), 180-day mortality (B, E), and 1-year mortality (C, F). (A–C) represent comparisons between Non-AKI and Mild-AKI groups, while (D–F) represent comparisons between Normal-or-mild-AKI and Moderate-to-severe-AKI groups. The shaded areas in the graphs indicate the 95% confidence interval.
Subgroup analysis
We performed a subgroup analysis stratified by age (< 65 or ≥65 years), SOFA score (< 5 or ≥5), heart failure status, and renal disease status. Within the Non-AKI and Mild-AKI group, patients under 65 years with Mild-AKI exhibited a significantly higher risk of 180-day and 1-year mortality compared to those without AKI. However, no significant differences in outcomes were observed between Non-AKI and Mild-AKI patients within the other subgroups; In the comparison between the Normal-or-mild-AKI and Moderate-to-severe-AKI groups, patients with Moderate-to-severe-AKI consistently showed significantly worse outcomes compared to those with Normal-or-mild-AKI across all subgroups. These findings are illustrated in Figure 4.

Figure 4. Forest plot of subgroup analysis for 28-day mortality (A, D), 180-day mortality (B, E), and 1-year mortality (C, F). (A–C) represent comparisons between Non-AKI and Mild-AKI groups, while (D–F) represent comparisons between Normal-or-mild-AKI and Moderate-to-severe-AKI groups. The dashed lines in the plots indicate the null effect line, and any result intersecting the null effect line suggests no significant difference.
Discussion
AMI is a leading cause of global mortality and morbidity. Patients with AMI are highly susceptible to complications such as AKI, which is common in critically ill patients and linked to worse clinical outcomes, including higher morbidity and mortality (22). While the association between AKI and adverse outcomes is well-documented, the impact of varying severities of AKI on the prognosis of AMI patients remains not fully understood. Addressing this crucial knowledge gap can inform targeted interventions and improve clinical decision-making. Our study uses advanced statistical methods, such as doubly robust estimation and rigorous sensitivity analyses, to explore the short- and long-term outcomes of ICU patients with different severities of AKI following AMI. This comprehensive approach aims to enhance clinical practices and guide future research by providing a detailed understanding of AKI severity's prognostic implications.
Cardiorenal syndrome (CRS) is a complex disorder characterized by bidirectional interactions between cardiac and renal dysfunction, mediated by multiple molecular mechanisms. It is a major cause of AKI in patients with AMI (20). Oxidative stress and inflammation represent central pathways in CRS pathogenesis. Activation of the NF-κB signaling pathway promotes the production of pro-inflammatory cytokines (e.g., IL-6, TNF-α) and oxidative stress markers (e.g., NOX2, iNOS), contributing to tissue damage in both organs (23, 24). Concurrently, impairment of the Nrf2 antioxidant pathway reduces the expression of cytoprotective enzymes (e.g., HO-1, MnSOD), exacerbating oxidative injury (20). These processes are evident across all CRS subtypes, as demonstrated in animal models of myocardial infarction (CRS1) and chronic kidney disease (CRS2) (25, 26). The renin-angiotensin-aldosterone system (RAAS) is hyperactivated in CRS, leading to vasoconstriction, sodium retention, and fibrosis. Upregulation of ACE and AT1R, coupled with downregulation of protective AT2R and MasR, has been observed in experimental models (27). RAAS inhibitors (e.g., ACEIs/ARBs) reduce urinary podocin loss in CRS2 patients, indicating glomerular protection (28). Moreover, the TGF-β1/Smad pathway mediates fibrosis through collagen deposition and epithelial-mesenchymal transition. Studies in CRS rats show elevated TGF-β1 and phosphorylated Smad3 in cardiac and renal tissues, which are attenuated by empagliflozin and dapagliflozin (29). Aberrant Wnt/β-catenin signaling contributes to cardiac hypertrophy and renal fibrosis. In CRS2 models, activation of β-catenin promotes pro-fibrotic gene expression (e.g., Twist, Snail1), while its inhibition with ICG-001 ameliorates organ damage (30). Gut microbiota dysbiosis further exacerbates CRS by producing uremic toxins (e.g., TMAO), which enhance inflammation and fibrosis via NF-κB and TGF-β1 pathways (31, 32). Noncoding RNAs, such as miR-21 and lncRNA ANRIL, also play roles by modulating fibrosis (e.g., targeting PPARα) and inflammasome activation (33, 34).
Research indicates that the incidence of AKI among ICU patients ranges from 12.1% to 60.93% (1). However, our study found a significantly higher AKI incidence of 75.35% in patients with AMI. We believe this discrepancy may stem from our more sensitive method of detecting AKI. In our study, each laboratory test and every fluctuation in fluid input/output were dynamically monitored throughout the hospitalization period, allowing for earlier detection of AKI. This increased sensitivity, as described in our study design, likely contributed to the higher incidence we observed. Supporting this, a study involving 1,050 AMI patients demonstrated that using the KDIGO criteria identified significantly more cases of AKI compared to the RIFLE criteria, suggesting that KDIGO is more sensitive for detecting AKI in AMI patients (15). Additionally, Kanic et al. (35) found that even minor rises in serum creatinine and progressive increases in AKI severity, as evaluated by the KDIGO criteria, were associated with poorer long-term outcomes in AMI patients. This underscores the importance of employing more sensitive methods, like KDIGO, to detect AKI in this population.
Compared to patients with Mild-AKI, those Non-AKI used loop diuretics and vasopressors more frequently, and they exhibited lower eGFR, higher SOFA scores, and elevated BUN levels. These five variables were identified as the most significant during the PSM process. The use of loop diuretics, lower eGFR, higher BUN levels, and elevated SOFA scores all indicate poorer renal function. Additionally, the frequent use of vasopressors suggests a higher incidence of hypotensive states, which can lead to renal ischemia and further kidney function deterioration (36).
In the PSM process comparing the Normal-or-mild-AKI and Moderate-to-severe-AKI groups, the five most significant variables were vasopressor use, SOFA score, loop diuretics, eGFR, and WBC count. The first four variables were consistent with the findings in the Non-AKI group. However, the WBC count was notably higher in the Moderate-to-severe-AKI group. Elevated leukocyte levels, particularly WBC, play a critical role in the pathophysiology of AKI, involving complex immunopathological interactions. These include mechanisms such as damage-associated molecular patterns (DAMPs), pathogen-associated molecular patterns (PAMPs), oxidative stress, hypoxia-inducible factors, the complement system, and various immune cells like dendritic cells, neutrophils, lymphocytes, and macrophages (37). A study by Chen et al. found that in AMI patients, the percentage of neutrophils in peripheral blood (NEUT%) was positively correlated with both the incidence of AKI and short-term all-cause mortality (38). Elevated serum calcium can induce renal vasoconstriction, reducing renal blood flow and causing tubular injury, while hypocalcemia may indicate the severity of cardiac dysfunction and renal impairment. Previous research (39) has demonstrated that acute kidney injury occurs more frequently in patients with ST-elevation myocardial infarction complicated by cardiogenic shock, leading to poor short-term clinical outcomes. A recent study (40) utilizing the MIMIC database developed a predictive model for AKI risk in AMI patients. Their model identified estimated glomerular filtration rate, creatinine, blood urea nitrogen, cardiogenic shock, and creatine-kinase myocardial band as the five most significant predictors. These findings are largely consistent with our results.
It is important to note that certain variables remained unbalanced between groups despite PSM. To minimize the impact of these variables on study outcomes, we performed corrections in subsequent modeling, though it remains necessary to discuss why statistical methods struggled to eliminate imbalances in these factors. Seven variables showed persistent intergroup imbalance in both analyses (Non-AKI vs. Mild-AKI and Normal/mild-AKI vs. Moderate-severe-AKI): SOFA score, loop diuretic use, vasopressor administration, preexisting renal disease, BUN, creatinine, and eGFR. These metrics are inherently linked to renal function, explaining why baseline disparities for such variables persisted even after PSM—patients with divergent renal profiles inherently exhibit unequal baselines for these parameters. The renal subcomponent of the SOFA score directly assesses renal function, while a history of preexisting renal disease characterizes chronic renal status. Loop diuretic use often reflects fluid overload, a hallmark of renal dysfunction, whereas vasopressor therapy typically indicates hypotension that may compromise renal perfusion. BUN, creatinine, and eGFR serve as direct renal function biomarkers: elevations in BUN, creatinine and declines in eGFR are sensitive indicators of deteriorating renal health.
Additionally, imbalance persisted in several indices across groups, presumably because statistical methods like PSM—being baseline characteristic-matched—cannot eliminate differences in indices strongly associated with the disease itself. In the Non-AKI vs. Mild-AKI group, renal concentrating dysfunction in AKI causes hyponatremia, while tubular potassium excretion impairment leads to hyperkalemia, contributing to data imbalance (41). Anemia (low hemoglobin) acts as both a risk factor for AKI (e.g., reduced oxygen-carrying capacity exacerbates renal injury during ischemia) and a consequence (e.g., decreased renal erythropoietin secretion), forming a bidirectional relationship (42).
In the Normal/mild-AKI vs. Moderate-severe-AKI group, heart failure impairs cardiac pumping, reducing renal perfusion, triggering renal vasoconstriction, and activating the renin-angiotensin-aldosterone system (RAAS) to induce/worsen AKI (43), making these patients more prone to moderate-severe AKI. WBC counts increase with infection, inflammation, and stress, and higher AKI severity correlates with greater infection probability and stress response (37), leading to uneven WBC distribution. As a key acid-base buffer regulated by the kidneys, bicarbonate metabolism is disrupted in AKI: varying degrees of renal injury across AKI severities cause differential impairment in bicarbonate reabsorption/secretion, and bicarbonate levels inversely correlate with AKI incidence and prognosis (44).
Numerous studies have examined the impact of AKI on prognosis in patients with AMI. Skalsky et al. (45) found that AMI patients with stage 1 AKI who did not recover within 48 h, as well as those with stage 2–3 AKI without recovery within 96 h, had a significantly higher risk of mortality. However, their diagnosis and staging of AKI were based solely on serum creatinine levels, without considering the diagnostic significance of urine output. Similarly, Kanic et al. (35) reported that the incidence of AKI among AMI patients undergoing PCI was 8.5%. During an average follow-up of 4.2 ± 3.0 years, the mortality rates were 50.3% for stage 1 AKI, 56.9% for stage 2, and 87.2% for stage 3. The hazard ratios for all-cause mortality were 1.77, 1.85, and 6.30 for stages 1, 2, and 3, respectively, compared to patients without AKI. In another study, Sun et al. observed 1,371 AMI patients and found that the severity of AKI, as classified by the KDIGO criteria, was an independent risk factor for 30-day mortality. Stage 3 AKI was also identified as an independent predictor of mortality between 30 days and 5 years. However, like previous studies, their definition of AKI relied solely on serum creatinine levels, without incorporating assessments of GFR or urine output. A review by Kaltsas et al. (11) summarized key studies on AKI complicating AMI, emphasizing that all studies consistently showed AKI worsened patient outcomes, increasing mortality by two- to threefold both within the first 30 days and throughout the first year after the acute event. Furthermore, a study (46) identified serum calcium levels as a strong predictor of AKI in AMI patients. Consequently, we set our outcome measures at 28 days, 180 days, and 1 year to better evaluate the long-term prognostic characteristics of this patient population.
We utilized a more sensitive dynamic assessment method based on the KDIGO criteria to evaluate the occurrence of AKI in patients with AMI. As shown in Table 3, no significant differences in prognosis were observed between the Non-AKI and Mild-AKI groups in the multivariable-adjusted models. However, the Log-rank test indicated differences in 180-day and 1-year mortality rates, suggesting that these differences may have been driven by other covariates rather than AKI itself. Further subgroup analysis identified age as a potential contributing factor. Specifically, patients under 65 years with Mild-AKI had worse outcomes compared to those without AKI. As reported in a study (47), younger patients in the ICU are more sensitive to nephrotoxic drugs (e.g., vancomycin and calcineurin inhibitors), which significantly deteriorate the prognosis of young AKI patients. Our analysis also suggests that physicians might adopt more conservative treatment strategies for these patients, assuming that younger individuals have stronger renal compensatory capacity. This approach may lead to progression of mild AKI or delayed control of systemic effects, thereby impacting patient outcomes. Across all models for the three outcome events, patients with Moderate-to-severe AKI consistently had worse outcomes compared to those with Normal-or-mild AKI. Sensitivity analyses confirmed the robustness of these findings.
Patients with Mild-AKI, referred to as “subclinical AKI,” only reached stage 1 AKI during hospitalization and did not experience adverse prognostic effects. In contrast, many patients initially classified as stage 1 AKI progressed to stages 2–3, leading to significantly worse outcomes compared to those who remained at AKI stages 0–1. The dynamic evaluation of AKI stages is clinically significant, as it allows for early detection of “subclinical AKI,” enabling timely intervention to prevent “conversion” to Moderate-to-severe AKI. The risk factors for progression may include the use of vasopressors, loop diuretics, higher SOFA scores, lower eGFR, and elevated WBC counts.
The treatment of CRS remains challenging. Diuretics, a mainstay in managing fluid overload, have uncertain long-term benefits. High-dose intermittent furosemide seems safe and effective in acute heart failure, but its impact on severe kidney disease is unclear. Ultrafiltration shows promise in some aspects like weight loss, yet its overall efficacy is still debated. CARRESS-HF indicated that ultrafiltration might not be the best primary treatment for type 1 CRS. While neurohormonal modulation therapies such as vasopressin antagonists and nesiritide have not significantly improved clinical outcomes in large-scale trials (48), sacubitril/valsartan has demonstrated renal protective effects in patients with cardiorenal syndrome (49). Furthermore, both traditional vasopressin antagonists and sacubitril/valsartan have been proven to be safe in clinical use. RAAS inhibitors are beneficial for some patients with CRS, yet they carry risks like hyperkalemia. β-adrenergic blockers have shown efficacy in reducing mortality in heart failure, but their use in CRS patients needs more evidence (48). Meanwhile, research has shown (50) that psychological interventions for patients with AKI can help improve their clinical outcomes.
Future research should focus on identifying additional influential factors and developing machine learning and deep learning models to predict the risk of moderate-to-severe AKI in AMI patients, and develop targeted effective treatment strategies to improve patient outcomes.
Limitation
While MIMIC-IV's data provide detailed records of clinical information for critically ill patients, the single-center and retrospective design warrant careful consideration of generalizability. Clinical data from this center may differ from those in community, rural, or international settings, potentially influencing outcome estimates—particularly for subgroups underrepresented in the dataset.
The retrospective design introduces risks of selection bias and unmeasured confounding factors. Although rigorous statistical methods were used to mitigate these limitations, residual confounding from unrecorded variables (e.g., family medical history, socioeconomic status) cannot be fully eliminated, which may affect the robustness of our findings.
External validation in independent cohorts remains essential to confirm the stability of these results. Future research should prioritize prospective multicenter studies across diverse healthcare systems to evaluate consistency across populations with differing baseline risks and care environments. Such efforts will enhance the generalizability of this study's findings and help provide a reliable theoretical foundation for clinical practice, ensuring these insights can inform real-world medical decision-making.
Conclusion
Patients with Mild-AKI can be more accurately described as having “subclinical AKI,” as their prognosis is often comparable to that of Non-AKI patients. However, the prognosis for those with Moderate-to-severe AKI is significantly worse than for patients with Normal-or-mild AKI. This indicates that if “subclinical AKI” undergoes a “conversion” to Moderate-to-severe AKI during hospitalization, the patient's prognosis will deteriorate considerably. Therefore, the dynamic and sensitive early identification of “subclinical AKI” and its potential “conversion” to Moderate-to-severe AKI is of great importance for timely intervention and improved outcomes.
Data availability statement
The datasets presented in this article are not readily available due to MIMIC dataset requiring principal investigator approval for access, the data cannot be made publicly available. Requests to access the datasets should be directed to https://www.physionet.org/content/mimiciv/3.0/.
Ethics statement
The studies involving humans were approved by the project received approval from the Institutional Review Boards at both the Massachusetts Institute of Technology (MIT) and Beth Israel Deaconess Medical Center (BIDMC). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
PG: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Visualization, Writing – original draft, Writing – review & editing. FT: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. LD: Conceptualization, Data curation, Formal analysis, Investigation, Software, Visualization, Writing – original draft, Writing – review & editing. HYa: Conceptualization, Data curation, Methodology, Validation, Writing – original draft. WW: Data curation, Formal analysis, Methodology, Visualization, Writing – review & editing. CM: Formal analysis, Software, Visualization, Writing – review & editing. XB: Conceptualization, Methodology, Writing – review & editing. LR: Investigation, Visualization, Writing – review & editing. HYi: Investigation, Writing – review & editing. LM: Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Key Research and Development Program of Qinhuangdao [grant numbers (202301A099, 202301A267, and 20242067)].
Acknowledgments
We would like to express our gratitude to the Ascetic Practitioners in Critical Care (APCC) team and the easy Data Science for Medicine (easyDSM) team for generously sharing their expertise and code related to big data in critical care, as well as for providing access to the cross-platform Big Data Master of Critical Care (BDMCC) software (https://github.com/ningyile/BDMCC_APP). We are also deeply appreciative of the MIMIC-IV database official team for their efforts in making the database and accompanying codes open-source. Additionally, we extend our thanks to the Home for Researchers editorial team (http://www.home-for-researchers.com) for their assistance with language editing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1555478/full#supplementary-material
References
1. Tao F, Yang H, Wang W, Bi X, Dai Y, Zhu A, et al. Acute kidney injury prediction model utility in premature myocardial infarction. iScience. (2024) 27:109153. doi: 10.1016/j.isci.2024.109153
2. Lisa A, Carbone F, Liberale L, Montecucco F. The need to identify novel markers for early renal injury in cardiorenal syndrome. Cells. (2024) 13:1283. doi: 10.3390/cells13151283
3. McCallum W, Testani JM. Updates in cardiorenal syndrome. Med Clin North Am. (2023) 107:763–80. doi: 10.1016/j.mcna.2023.03.011
4. Mullens W, Verbrugge FH, Nijst P, Tang W. Renal sodium avidity in heart failure: from pathophysiology to treatment strategies. Eur Heart J. (2017) 38:1872–82. doi: 10.1093/eurheartj/ehx035
5. Mullens W, Abrahams Z, Skouri HN, Francis GS, Taylor DO, Starling RC, et al. Elevated intra-abdominal pressure in acute decompensated heart failure: a potential contributor to worsening renal function. J Am Coll Cardiol. (2008) 51:300–6. doi: 10.1016/j.jacc.2007.09.043
6. Peoples JN, Saraf A, Ghazal N, Pham TT, Kwong JQ. Mitochondrial dysfunction and oxidative stress in heart disease. Exp Mol Med. (2019) 51:1–13. doi: 10.1038/s12276-019-0355-7
7. Milkovic L, Cipak Gasparovic A, Cindric M, Mouthuy PA, Zarkovic N. Short overview of ROS as cell function regulators and their implications in therapy concepts. Cells. (2019) 8:793. doi: 10.3390/cells8080793
8. Gallo G, Lanza O, Savoia C. New insight in cardiorenal syndrome: from biomarkers to therapy. Int J Mol Sci. (2023) 24:5089. doi: 10.3390/ijms24065089
9. McWilliam SJ, Wright RD, Welsh GI, Tuffin J, Budge KL, Swan L, et al. The complex interplay between kidney injury and inflammation. Clin Kidney J. (2021) 14:780–8. doi: 10.1093/ckj/sfaa164
10. Akhter MS, Goodwin JE. Endothelial dysfunction in cardiorenal conditions: implications of endothelial glucocorticoid receptor-wnt signaling. Int J Mol Sci. (2023) 24:14261. doi: 10.3390/ijms241814261
11. Kaltsas E, Chalikias G, Tziakas D. The incidence and the prognostic impact of acute kidney injury in acute myocardial infarction patients: current preventive strategies. Cardiovasc Drugs Ther. (2018) 32:81–98. doi: 10.1007/s10557-017-6766-6
12. Bellomo R, Kellum JA, Ronco C. Defining and classifying acute renal failure: from advocacy to consensus and validation of the RIFLE criteria. Intensive Care Med. (2007) 33:409–13. doi: 10.1007/s00134-006-0478-x
13. Mehta RL, Kellum JA, Shah SV, Molitoris BA, Ronco C, Warnock DG, et al. Acute kidney injury network: report of an initiative to improve outcomes in acute kidney injury. Crit Care. (2007) 11:R31. doi: 10.1186/cc5713
14. Kdigo KD, Outcomes IG. Acute kidney injury work group: KDIGO clinical practice guideline for acute kidney injury. Kidney Int Suppl. (2012) 2:1–138. Available online at: https://www.researchgate.net/publication/284548602_Acute_Kidney_Injury_Work_Group_KDIGO_clinical_practice_guideline_for_acute_kidney_injury#citations
15. Rodrigues FB, Bruetto RG, Torres US, Otaviano AP, Zanetta DM, Burdmann EA. Incidence and mortality of acute kidney injury after myocardial infarction: a comparison between KDIGO and RIFLE criteria. PLoS ONE. (2013) 8:e69998. doi: 10.1371/journal.pone.0069998
16. Yaqub S, Hashmi S, Kazmi MK, Aziz Ali A, Dawood T, Sharif H, et al. Comparison of AKIN, KDIGO, and RIFLE definitions to diagnose acute kidney injury and predict the outcomes after cardiac surgery in a South Asian cohort. Cardiorenal Med. (2022) 12:29–38. doi: 10.1159/000523828
17. Johnson A, Bulgarelli L, Shen L, Gayles A, Shammout A, Horng S, et al. MIMIC-IV, a freely accessible electronic health record dataset. Sci Data. (2023) 10:1. doi: 10.1038/s41597-023-01945-2
18. Levey AS, Coresh J, Greene T, Stevens LA, Zhang YL, Hendriksen S, et al. Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med. (2006) 145:247–54. doi: 10.7326/0003-4819-145-4-200608150-00004
19. Funk MJ, Westreich D, Wiesen C, Stürmer T, Brookhart MA, Davidian M. Doubly robust estimation of causal effects. Am J Epidemiol. (2011) 173:761–7. doi: 10.1093/aje/kwq439
20. Zhao BR, Hu XR, Wang WD, Zhou Y. Cardiorenal syndrome: clinical diagnosis, molecular mechanisms and therapeutic strategies. Acta Pharmacol Sin. (2025) 46:1539–55. doi: 10.1038/s41401-025-01476-z
21. Feng M, McSparron JI, Kien DT, Stone DJ, Roberts DH, Schwartzstein RM, et al. Transthoracic echocardiography and mortality in sepsis: analysis of the MIMIC-III database. Intensive Care Med. (2018) 44:884–92. doi: 10.1007/s00134-018-5208-7
22. Landi A, Branca M, Leonardi S, Frigoli E, Vranckx P, Tebaldi M, et al. Transient vs in-hospital persistent acute kidney injury in patients with acute coronary syndrome. JACC Cardiovasc Interv. (2023) 16:193–205. doi: 10.1016/j.jcin.2022.10.009
23. Cho E, Kim M, Ko YS, Lee HY, Song M, Kim MG, et al. Role of inflammation in the pathogenesis of cardiorenal syndrome in a rat myocardial infarction model. Nephrol Dial Transplant. (2013) 28:2766–78. doi: 10.1093/ndt/gft376
24. Jin L, Li Q, Li J, Pan Y, Zou J, Wu X, et al. Apela inhibits systemic and renal inflammatory reactions in mice with type I cardiorenal syndrome. FASEB J. (2021) 35:e21907. doi: 10.1096/fj.202101030R
25. Deng T, Wei Z, Gael A, Deng X, Liu Y, Lai J, et al. Higenamine improves cardiac and renal fibrosis in rats with cardiorenal syndrome via ASK1 signaling pathway. J Cardiovasc Pharmacol. (2020) 75:535–44. doi: 10.1097/FJC.0000000000000822
26. Giam B, Kuruppu S, Chu PY, Smith AI, Marques FZ, Fiedler A, et al. N-acetylcysteine attenuates the development of renal fibrosis in transgenic mice with dilated cardiomyopathy. Sci Rep. (2017) 7:17718. doi: 10.1038/s41598-017-17927-5
27. Urbanek K, Cappetta D, Bellocchio G, Coppola MA, Imbrici P, Telesca M, et al. Dapagliflozin protects the kidney in a non-diabetic model of cardiorenal syndrome. Pharmacol Res. (2023) 188:106659. doi: 10.1016/j.phrs.2023.106659
28. Le Jemtel TH, Rajapreyar I, Selby MG, Payne B, Barnidge DR, Milic N, et al. Direct evidence of podocyte damage in cardiorenal syndrome type 2: preliminary evidence. Cardiorenal Med. (2015) 5:125–34. doi: 10.1159/000375130
29. Yang CC, Chen YT, Chen CH Li YC, Shao PL, Huang TH, et al. The therapeutic impact of entresto on protecting against cardiorenal syndrome-associated renal damage in rats on high protein diet. Biomed Pharmacother. (2019) 116:108954. doi: 10.1016/j.biopha.2019.108954
30. Zhao Y, Wang C, Hong X, Miao J, Liao Y, Hou FF, et al. Wnt/β-catenin signaling mediates both heart and kidney injury in type 2 cardiorenal syndrome. Kidney Int. (2019) 95:815–29
31. Nanto-Hara F, Kanemitsu Y, Fukuda S, Kikuchi K, Asaji K, Saigusa D, et al. The guanylate cyclase C agonist linaclotide ameliorates the gut-cardio-renal axis in an adenine-induced mouse model of chronic kidney disease. Nephrol Dial Transplant. (2020) 35:250–64
32. Zou D, Li Y, Sun G. Attenuation of circulating trimethylamine N-oxide prevents the progression of cardiac and renal dysfunction in a rat model of chronic cardiorenal syndrome. Front Pharmacol. (2021) 12:751380
33. Chuppa S, Liang M, Liu P, Liu Y, Casati MC, Cowley AW, et al. MicroRNA-21 regulates peroxisome proliferator-activated receptor alpha, a molecular mechanism of cardiac pathology in Cardiorenal Syndrome Type 4. Kidney Int. (2018) 93:375–89. doi: 10.1016/j.kint.2017.05.014
34. Xu Y, Cao L, Zou W, Yu R, Shen W. Panax notoginseng saponins inhibits NLRP3 inflammasome-mediated pyroptosis by downregulating lncRNA-ANRIL in cardiorenal syndrome type 4. Chin Med. (2023) 18:50. doi: 10.1186/s13020-023-00756-2
35. Kanic V, Kompara G, Suran D, Ekart R, Bevc S, Hojs R. Impact of KDIGO-defined acute kidney injury on mortality after percutaneous coronary intervention for acute myocardial infarction. Cardiorenal Med. (2018) 8:332–9. doi: 10.1159/000492287
36. Marenzi G, Assanelli E, Campodonico J, De Metrio M, Lauri G, Marana I, et al. Acute kidney injury in ST-segment elevation acute myocardial infarction complicated by cardiogenic shock at admission. Crit Care Med. (2010) 38:438–44. doi: 10.1097/CCM.0b013e3181b9eb3b
37. Gameiro J, Lopes JA. Complete blood count in acute kidney injury prediction: a narrative review. Ann Intensive Care. (2019) 9:87. doi: 10.1186/s13613-019-0561-4
38. Chen Q, Gu Q, Yin A, Cai D, Xiao T, Wang Y, et al. Neutrophil percentage as a potential biomarker of acute kidney injury risk and short-term prognosis in patients with acute myocardial infarction in the elderly. Clin Interv Aging. (2024) 19:503–15. doi: 10.2147/CIA.S455588
39. Kanabar K, Sharma YP, Krishnappa D, Santosh K, Dhudasia M. A study of the predictive role of multiple variables for the incidence of acute kidney injury and its outcomes in Indian patients with ST-elevation myocardial infarction and cardiogenic shock. Egypt Heart J. (2024) 76:123. doi: 10.1186/s43044-024-00557-4
40. Wei J, Cai D, Xiao T, Chen Q, Zhu W, Gu Q, et al. Artificial intelligence algorithms permits rapid acute kidney injury risk classification of patients with acute myocardial infarction. Heliyon. (2024) 10:e36051. doi: 10.1016/j.heliyon.2024.e36051
41. Chen DN, Du J, Xie Y, Li M, Wang RL, Tian R. Relationship between early serum sodium and potassium levels and AKI severity and prognosis in oliguric AKI patients. Int Urol Nephrol. (2021) 53:1171–87. doi: 10.1007/s11255-020-02724-3
42. Han SS, Baek SH, Ahn SY, Chin HJ, Na KY, Chae DW, et al. Anemia is a risk factor for acute kidney injury and long-term mortality in critically Ill patients. Tohoku J Exp Med. (2015) 237:287–95. doi: 10.1620/tjem.237.287
43. Feng LH, Su T, Huang L, Liao T, Lu Y, Wu L. Development and validation of a dynamic nomogram for acute kidney injury prediction in ICU patients with acute heart failure. Front Med. (2025) 12:1544024. doi: 10.3389/fmed.2025.1544024
44. Huang L, Wang F. Identification of L-shaped curve association between serum bicarbonate concentrations and short-term outcomes in patients with acute kidney injury: a retrospective cohort study. Ren Fail. (2025) 47:2462264. doi: 10.1080/0886022X.2025.2462264
45. Skalsky K, Shiyovich A, Shechter A, Gilutz H, Plakht Y. Recovery from acute kidney injury and long-term prognosis following acute myocardial infarction. Biomedicines. (2024) 12:1490. doi: 10.3390/biomedicines12071490
46. Wang Y, Lu Y, Liu C, Xiao J. Association between serum calcium level and the risk of acute kidney injury in patients with acute myocardial infarction: evidences from the MIMIC-IV database. Ren Fail. (2024) 46:2401137. doi: 10.1080/0886022X.2024.2401137
47. Fuhrman DY, Kane-Gill S, Goldstein SL, Priyanka P, Kellum JA. Acute kidney injury epidemiology, risk factors, and outcomes in critically ill patients 16-25 years of age treated in an adult intensive care unit. Ann Intensive Care. (2018) 8:26. doi: 10.1186/s13613-018-0373-y
48. Rangaswami J, Bhalla V, Blair J, Chang TI, Costa S, Lentine KL, et al. Cardiorenal syndrome: classification, pathophysiology, diagnosis, and treatment strategies: a scientific statement from the American Heart Association. Circulation. (2019) 139:e840–878. doi: 10.1161/CIR.0000000000000664
49. Rahman MM, Rahman A, Nishiyama A. Potential renoprotective effects and possible underlying mechanisms of angiotensin receptor-neprilysin inhibitors in cardiorenal syndrome. Front Med. (2024) 11:1451450. doi: 10.3389/fmed.2024.1451450
Keywords: myocardial infarction, acute kidney injury, prognosis, propensity score matching, doubly robust analysis
Citation: Guo P, Tao F, Du L, Yang H, Wang W, Ma C, Bi X, Ren L, Yin H and Ma L (2025) Acute kidney injury after myocardial infarction: prognostic implications via dual robust methods. Front. Med. 12:1555478. doi: 10.3389/fmed.2025.1555478
Received: 04 January 2025; Accepted: 03 July 2025;
Published: 22 July 2025.
Edited by:
Gregory Braden, University of Massachusetts Medical School, United StatesReviewed by:
Hua Miao, Northwest University, ChinaSuleyman Yıldız, Prof. Dr. Cemil Taşcioglu Sehir Hastanesi, Türkiye
Copyright © 2025 Guo, Tao, Du, Yang, Wang, Ma, Bi, Ren, Yin and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lixiang Ma, cWhkbWFsaXhpYW5nQDE2My5jb20=
†These authors have contributed equally to this work and share first authorship