- 1Laboratory of Cellular Immunity, Institute of Clinical Immunology, Shuguang Hospital Affiliated to Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 2Department of Liver Intensive Care Unit, Shanghai Public Health Clinical Center, Fudan University, Shanghai, China
- 3Institute of Infectious Diseases of Integrated Traditional Chinese and Western Medicine, Shanghai, China
Background: Acute kidney injury (AKI) is a frequent complication in patients with hepatitis B virus-related acute-on-chronic liver failure (HBV-ACLF) and is associated with high rates of mortality. We aimed to estimate serum high mobility group protein 1 (HMGB1) levels in hepatitis B virus-related acute-on-chronic liver failure patients and analyze their clinical value in the development and outcomes of Acute kidney injury.
Methods: A total of 251 consecutive patients with hepatitis B virus-related acute-on-chronic liver failure were enrolled in this retrospective study. Using the International Club of Ascites staging criteria of Acute kidney injury, 153 patients developed Acute kidney injury. The clinical data of patients were collected and serum levels of high mobility group protein 1 were measured by ELISA. All patients were followed up until death or for a minimum of 3 months. Early prediction and prognostic implications of high mobility group protein 1 in Hepatitis B Virus-Related Acute-on-Chronic Liver Failure Patients with Acute Kidney Injury were investigated in different cohorts, including a propensity score-matched ACLF cohort.
Results: Among all individuals with hepatitis B virus-related acute-on-chronic liver failure, the incidence of Acute kidney injury was 61.0% (153/251). The patients who developed stage 2/3 Acute kidney injury showed the highest high mobility group protein 1 levels, followed by those who developed stage 1 Acute kidney injury, and those without Acute kidney injury showed the lowest high mobility group protein 1 levels. Moreover, high mobility group protein 1 levels were significantly higher in non-survivors than in survivors among hepatitis B virus-related acute-on-chronic liver failure patients with Acute kidney injury. Furthermore, analysis of the area under the receiver operating characteristic curve (AUROC) indicated that serum high mobility group protein 1 levels (pre-matching: AUC = 0.740; post-matching: AUC = 0.661) may be a potential predictive factor for Acute kidney injury development and that high mobility group protein 1 (AUC = 0.727) might be a reliable biomarker for prognosis in patients with Acute kidney injury.
Conclusion: In patients with hepatitis B virus-related acute-on-chronic liver failure, Acute kidney injury is universal. Acute kidney injury and its stages negatively influence the 90-day transplant-free mortality rate. Serum high mobility group protein 1 levels can serve as a positive predictor of Acute kidney injury development, and high mobility group protein 1 might also be a prognostic biomarker for Acute kidney injury among hepatitis B virus-related acute-on-chronic liver failure patients.
Highlights
1) Elevated HMGB1 levels correlated with AKI stage in patients with HBV-ACLF.
2) HMGB1 can be used as an early AKI predictor in individuals with HBV-ACLF.
3) HMGB1 can serve as a prognostic biomarker for AKI in individuals with HBV-ACLF.
Introduction
Acute kidney injury (AKI) is a major clinical concern of acute-on-chronic liver failure (ACLF). In China, patients with hepatitis B virus (HBV)-related acute-on-chronic liver failure (HBV-ACLF) are characterized by high-level systemic inflammation, organ failure and high 90-day mortality (50%–70%) (Li et al., 2016; Wu et al., 2018). The occurrence of AKI in ACLF usually aggravates the disease and frequently leads to a negative prognostic impact. A recent meta-analysis concluded that approximately 40% of ACLF patients had AKI complications and that AKI is substantially associated with increased short-term mortality (Jiang et al., 2020). The identification of AKI is therefore of great importance for prevention and early treatment. Early prediction of AKI development provides a critical therapeutic time window and will provide the most appropriate and timely interventions for patients, which may contribute to better outcomes for patients with AKI. However, the current clinical detection based on creatinine and urine output (UO) is not effective in the early detection of AKI. AKI diagnosis is currently based on an increase in serum creatinine (sCr) and weight-adjusted hourly UO. Unfortunately, Cr is unreliable for predicting acute changes in kidney function. Over the past decades, researchers have evaluated more than 26 new biomarkers for AKI, including cystatin C (CysC) (Jha et al., 2021), neutrophil gelatinase-associated lipocalin (Nickolas et al., 2012), angiotensinogen (Chen C et al., 2016), and liver-type fatty acid-binding protein (Ferguson et al., 2010). However, the high cost, instability and limitations of these markers hinder their wide application in clinical and biochemical settings, and most of these markers need continuous systemic evaluation, which is not cost-effective. To date, only CysC has small-scale clinical applications.
High mobility group box 1 (HMGB1), a proinflammatory cytokine, was discovered several decades ago (Goodwin et al., 1973) and is one of the most abundant and highly conserved proteins in eukaryotic cells. HMGB1 regulates cellular proliferation, differentiation, tissue regeneration and inflammation (Gu et al., 2016), which contributes to kidney disease pathogenesis. The pathogenic mechanism for AKI in patients with ACLF may be related to bacterial endotoxin, inflammatory cytokines, and reduction of effective renal blood volume caused by substantial ascites. A cross-sectional study showed that serum HMGB1 levels increased in AKI patients and that elevated HMGB1 levels were positively associated with leukocyte count and negatively associated with proteinuria (Zakiyanov et al., 2013). Infection accompanied by the release of HMGB1 contributes to circulatory dysfunction and is one of the major contributing factors of AKI. Experimental studies have shown that HMGB1 release from podocytes is induced by lipopolysaccharide (LPS), which subsequently exacerbates AKI (Gao et al., 2021). The release of HMGB1 worsened circulatory dysfunction. Recently, a study showed that serum HMGB1 levels had a high correlation with infection and a moderate correlation with AKI and death (Vilela et al., 2018). Another study including 25 cirrhotic patients with AKI demonstrated that the presence of high HMGB1 levels was associated with decreased survival (de Oliveira Gomes et al., 2020). Therefore, serum HMGB1 levels may be associated with the development or progression of AKI. Moreover, the aberrant expression of HMGB1 may also be related to worse prognosis in AKI patients.
AKI is a common and serious complication of HBV-ACLF. However, comparatively little is known about the impact of HMGB1 on ACLF-associated AKI. In this study, we measured the serum levels of HMGB1 in all enrolled patients, analyzed the clinical relevance of HMGB1 levels to ACLF-associated AKI and further explored the early predictive and prognostic values of HMGB1 in ACLF-associated AKI.
Methods
Study patients and design
From January 2019 through January 2022 at Shanghai Public Health Clinical Center, 410 hospitalized patients with HBV-ACLF were reviewed in this retrospective study. The flow chart is shown in Figure 1. In light of the unanimous suggestions of the Asian Pacific Association for the Study of the Liver (APASL) (Sarin et al., 2019), patients matching the following criteria were recruited for this research: 1) a history of HBV-related cirrhosis or chronic liver disease that was either previously recognized or gone undetected; 2) 18–80 years old; 3) coagulation abnormalities [prothrombin activity ≤40% or international normalized ratio (INR) ≥ 1.5] with serum total bilirubin (TBIL) ≥ 5 mg/dl (85 µmol/L); and 4) ascites and hepatic encephalopathy (HE) clinical manifestations within 4 weeks. Patients were excluded if they had one of exclusion criteria: 1) inconsistent with APASL criteria; 2) co-infection with the human immunodeficiency virus or the hepatitis A/C/E virus; 3) patients who had other chronic liver diseases; 4) patients who had severe extra-hepatic diseases; 5) patients who had previously used hepatotoxic, immunosuppressive, or anticancer medications; 6) patients who had severe disease complications which may affect survival time; 7) patients who had AKI at admission; 8) pregnancy. Finally, this research included 251 individuals with HBV-ACLF. The baseline data and demographic information for participants at admission were collected, including demographic data, complication rate, and laboratory measures. Peripheral blood was collected at baseline to measure HMGB1 levels. In addition, detailed methods of score calculation are specified in the Supplemental Methods, including the model for end-stage liver disease (MELD), Child-Turcotte-Pugh (Child-Pugh), CLIF Consortium Organ Failure (CLIF-OF) and CLIF-Consortium ACLF (CLIF-C ACLF) scores. Based on the established guidelines (Sarin et al., 2009; European Association for the Study of the Liver, 2010; Sarin et al., 2014; Vilstrup et al., 2014), patients with HBV-ACLF were managed.
FIGURE 1. Flow chart of all enrolled HBV-ACLF patients with AKI after applying the inclusion and exclusion criteria.
Definition of AKI or renal response to treatment
AKI diagnoses and classification were performed in accordance with the International Club of Ascites (ICA) (Angeli et al., 2015). We evaluated renal function at baseline based on the most recent sCr level within 3 months prior to admission or the sCr value measured first after hospitalization. AKI was defined as follows: AKI stage 1 was defined as sCr increased by ≥ 0.3 mg/dl (26.5 µmol/L) in 48 h or to 1.5–2.0 times baseline during the first 7 days; AKI stage 2 was defined as sCr increased 2.0–3.0 times from baseline; and AKI stage 3 was defined as sCr increased 3.0 times from baseline or an increase by ≥ 4.0 mg/dl (353.6 μmol/L) or initiation of renal replacement treatment (RRT).
To compare the renal outcome among different stages of AKI, we made the following definition. We defined Stage 1 (0.3) as an increase of ≥0.3 mg/dl in the sCr level in 48 h only; Stage 1 (50%) as a 50% increase above baseline sCr; Stage 3 as not requiring RRT; and Stage 3 (RRT) as a requirement of RRT.
The peak sCr was defined by the highest sCr value during hospital admission. We defined stage 2–3 AKI as severe AKI. The maximum AKI stage (AKImax) was abstracted from medical records during hospitalization.
According to the ICA, renal responses to treatment were classified as complete response (CR), partial response (PR), or absence of response (AR). CR represents an sCr return to within 0.3 mg/dl at hospital discharge. PR represents a decrease in AKI stage compared with their maximal AKI stage or a sCr reduction to ≥0.3 mg/dl. AR represents without recovery compared with their maximal AKI stage.
Study end points
All patients were followed up until death or at least 3 months. Our primary outcome was 90-day survival. Liver transplantation (LT) or death was considered the adverse endpoint. Our secondary outcome measure was AKI outcome (CR, PR, AR) defined according to ICA criteria.
Measurement of serum HMGB1
Quantification of HMGB1 in plasma was performed using ELISA kits (Cat. ST51011, Shino-test Corporation, Kanagawa, Japan). The assay was conducted in accordance with the manufacturers’ instructions. The specimens were stored in a -80°C freezer until further handling.
Statistical analysis
In this study, normally distributed continuous variables are expressed as the mean ± SD, and categorical variables are expressed as frequencies and percentages. Quantitative variables were compared using Student’s t test, and categorical variables were compared using the chi-square test or Fisher’s exact test as appropriate.
Continuous variables that were not normally distributed are presented as the median (interquartile range [IQR]) and were compared by using Mann-Whitney U tests or Kruskal–Wallis tests. Continuous variables of three or more groups were analyzed using the Kruskal–Wallis test. Propensity score–matching analysis was used to reduce the effect of selection bias and potential confounding between the no-AKI and AKI groups. For propensity score matching, a nearest-neighbor 1:1 matching scheme with a caliper size of 0.02 was used. Survival analyses were performed using the Kaplan-Meier method and the log-rank test. Multivariate analysis was carried out using Cox regression models and logistic regression. The HMGB1 levels were divided into four groups according to the HMGB1 quartile levels (negative: <25%; low: 25–49.9%; moderate: 50–74.9%; high: ≥75%), and the p-value for trend was obtained by regression analysis. The performance of HMGB1 was assessed by the area under the receiver operating characteristic curve (AUC). To identify a differential effect of HMGB1 on outcomes per subgroup, interaction analysis using the Cox model was conducted, where an interaction term was composed of HMGB1 and a subgroup variable and presented with forest plots.
All statistical analyses were performed with SPSS 25.0 (IBM, Seattle, WA, USA) software or in R (version 4.0.4, http://www.r-project.org). Two-sided p < 0.05 was considered statistically significant.
Results
Demographic and clinical characteristics of the participants
A total of 251 patients with HBV-ACLF were enrolled. A total of 153 (61.0%) of 251 patients with HBV-ACLF developed AKI after hospital admission. Among those with AKI, 53 patients (34.6%) had stage 1 AKI, 61 patients (39.9%) had stage 2 AKI, and 39 patients (25.5%) had stage 3 AKI. The clinical baseline characteristics of the no-AKI and AKI stage 1–3 groups in the present study are shown in Table 1. The patients with AKI were older than patients without AKI (p = 0.001). The patients with AKI had higher TBIL (p < 0.001), INR (p = 0.035), HBVDNA (p = 0.046), sCr (p < 0.001), peak sCr (p < 0.001), CysC (p < 0.001), WBC (p = 0.012), and serum lactate (p = 0.023) levels, and lower Na (p = 0.008) levels than the patients without AKI. Patients with AKI more frequently presented with HE (p = 0.001), infection (p = 0.001) and ascites (p = 0.024) at admission than patients without AKI. Compared with individuals without AKI, individuals with AKI had higher MELD (p < 0.001), Child–Pugh (p < 0.001) and CLIF-C ACLF (p < 0.001) scores.
TABLE 1. Baseline demographic and clinical parameters of HBV-ACLF patients with and without AKI according to ICA criteria.
The serum levels of HMGB1 in patients with stage 3 AKI [6.49 [3.32–9.84] mg/dl] and stage 2 AKI [4.62 [2.74–6.6] mg/dl] were meaningfully higher than those with stage 1 AKI [3.46 [2.54–4.9] mg/dl], and serum HMGB1 levels in patients with stage 1 AKI were relatively higher than those without AKI [2.67 [1.74–3.63] mg/dl] (Table 1 and Supplementary Figure S1A). However, no significant difference was observed between the AKI stage 3 and AKI stage 2 groups (p = 0.193, Supplementary Figure S1A). To elucidate the clinical significance of HMGB1 for ACLF, we compared serum HMGB1 levels among HBV-ACLF patients between different groups of organ failure according to CLIF-C OF scores. A positive correlation was found between serum HMGB1 levels and the degrees of brain failure or coagulation failure but not liver failure (Supplementary Figure S1B, S1C, S1D).
To evaluate whether HMGB1 levels constitute a risk factor independent of early AKI prediction in HBV-ACLF patients, as confounding covariates, we included age, ascites, CysC, sCr, and CLIF-C ACLF scores and performed 1:1 propensity score matching (n = 138, 69 patients with AKI were matched to 69 patients without AKI). Baseline characteristics of pre- and post-matching cohorts are displayed in Table 2. Hospitalized HBV-ACLF patients who experienced AKI had significantly higher baseline HMGB1 levels than well-matched patients hospitalized without AKI (p = 0.001).
TABLE 2. Characteristics of 251 pre-matching or 138 post-matching ACLF patients (1:1 propensity score matching (No-AKI: AKI), covariables age, ascites, CysC, sCr, and CLIF-C ACLF score).
Clinical characteristics of non-survivors and survivors in AKI individuals with HBV-ACLF
According to the 90-day prognosis, the AKI group was divided into a survival group (n = 64) and a non-survival group (n = 89). We found that non-survivors were older than survivors (p = 0.044). Meanwhile, non-survivors had a higher frequency of HE (p<0.001), infection (p = 0.001) and ascites (p = 0.001). More importantly, the levels of TBIL (p < 0.001), INR (p = 0.010), peak sCr (p < 0.001), serum lactate (p = 0.015) and HMGB1 (p < 0.001) in non-survivors were also significantly higher than those in survivors. Likewise, the MELD (p = 0.001), Child-Pugh (p < 0.001) and CLIF-C ACLF (p < 0.001) scores were higher at baseline for the non-survival group than for the survival group. The above characteristics in non-survivors and survivors are listed in Table 3.
TABLE 3. Baseline characteristics of HBV-ACLF patients with AKI based on different prognostic groups at admission.
Longitudinal follow-up evaluation and kidney outcomes
The percentages of outcome events (death, survival, LT) between the no-AKI and AKI stage 1/2/3 groups during the 90-day follow-up period in all HBV-ACLF patients are shown in Supplementary Figure S2. AKI patients had a higher LT-free mortality or incidence of adverse outcomes than those without AKI at 90 days (45.8% vs. 10.2%, p < 0.001; 58.2% vs. 15.3%, p < 0.001).
Subsequently, the relationship between different AKI stages and the survival rate was analyzed by using Kaplan-Meier survival curves plotted using 90-day prognosis. Transplant-free survival rates were significantly decreased in patients who developed AKI compared with those who did not develop AKI at 90 days (pre-matching: 41.8% vs. 84.7%, p < 0.0001, Figure 2A; post-matching: 50.7% vs. 89.9%, p < 0.0001, Figure 2C). The transplant-free survival rate of patients without AKI was significantly higher than that of patients with Stage 1AKI (pre-matching: 84.7% vs. 64.2%, p < 0.01, Figure 2B; post-matching: 89.9% vs. 66.7%, p < 0.01, Figure 2D). The transplant-free survival rate was substantially higher among patients in the Stage 1 AKI group than among patients in the Stage 2 AKI group (64.2% vs. 37.7%, p < 0.01, Figure 2B) and there was no significant difference in the transplant-free survival rates between patients with Stage 2 AKI and those with Stage 3 AKI in the pre-matching cohort (p > 0.05, Figure 2B). The transplant-free survival rate was substantially higher among patients in the Stage 2 AKI group than among patients in the Stage 3 AKI group (52.4% vs. 22.2%, p < 0.05, Figure 2D) and there was no significant difference in the transplant-free survival rates between patients with Stage 1 AKI and those with Stage 2 AKI in the post-matching cohort (p > 0.05, Figure 2D).
FIGURE 2. Survival curves for patients with adverse outcomes are shown. (A) The Kaplan–Meier curves are classified according to the presence or absence of AKI during the 90-day follow-up period in the pre-matching cohort. (B) The Kaplan–Meier curves are classified according to the presence or absence of AKI and its stage during the 90-day follow-up period in the pre-matching cohort. (C) The Kaplan–Meier curves are classified according to the presence or absence of AKI during the 90-day follow-up period in the post-matching cohort. (D) The Kaplan–Meier curves are classified according to the presence or absence of AKI and its stage during the 90-day follow-up period in the post-matching cohort.
As expected, there was a correlation between kidney outcomes and AKI stages in the pre-matching cohort. Among patients who developed AKI, 17 (11.1%) patients developed Stage 1 (0.3), 36 (23.5%) patients developed Stage 1 (50%), 61 (39.9%) patients developed Stage 2, 30 patients (19.6%) developed Stage 3, and nine patients (5.9%) developed Stage 3 (RRT) (Supplementary Table S1). The CR was only 37.3% in the entire AKI population, 20.3% of patients still had stage 1 AKI, 23.5% of patients still had stage 2 AKI, and 19% of patients had stage 3 AKI. These results are shown in Supplementary Table S1. We found that the CR rate decreased with increasing AKI stage in the entire AKI cohort [i.e., from 76.5% in Stage 1 (0.3) to 22.2% in Stage 3 (RRT) (p < 0.001)]. The renal outcomes of surviving patients are listed in Supplementary Table S2. The CR was 64.1%, PR was 12.5% and AR was 23.4% among survivors. The incidence of CR was markedly higher in the surviving patients than in the total AKI population (64.1% vs. 37.3%, p < 0.001).
Performance of serum HMGB1 levels for predicting AKI development
First, an AUROC analysis was used to evaluate the performance of HMGB1 and CysC in predicting the development of AKI in the pre-matching cohort as shown in Figure 3A. Serum HMGB1 showed strong predictive ability, with a high AUROC value [0.740 (95% CI: 0.680–0.800)], which was comparable to the AUROC value obtained with CysC [0.728 (95% CI: 0.664–0.791)]. The joint indicator (HMGB1+ CysC) showed stronger predictive ability [0.772 (95% CI: 0.715–0.830)] than serum HMGB1 (p = 0.012). HMGB1 as a potential biomarker predicted AKI development with a sensitivity of 93.9% and specificity of 44.4%, in which the cut-off value was set at 4.850 ng/ml (Supplementary Table S3). Likewise, serum HMGB1 showed strong predictive ability, with a high AUROC value [0.661 (95% CI: 0.570–0.752)] in the post-matching cohort, as shown in Figure 4.
FIGURE 3. AUROC analysis was performed to compare the efficacy of serum HMGB1 and CysC levels in predicting the development of AKI (A) or severe AKI (B) in pre-matching HBV-ACLF patients during the 90-day follow-up period (AUROC, area under the receiver operating characteristic curve; HBV-ACLF, hepatitis B virus-related acute-on-chronic liver failure; AKI, acute kidney injury; CysC, cystatin C; HMGB1, high-mobility group box 1).
FIGURE 4. AUROC analysis of serum HMGB1 levels in predicting the development of AKI in the post-matching HBV-ACLF cohort during the 90-day follow-up period (AUROC, area under the receiver operating characteristic curve; HBV-ACLF, hepatitis B virus-related acute-on-chronic liver failure; AKI, acute kidney injury; HMGB1, high-mobility group box 1).
Next, we further explored the potential value of HMGB1 in predicting the development of severe AKI in all enrolled HBV-ACLF patients as shown in Figure 3B. Serum HMGB1 showed strong predictive ability, with a high AUROC value [0.746 (95% CI: 0.680–0.812)], which was comparable to the AUROC value obtained with CysC [0.744 (95% CI: 0.682–0.806)]. The joint indicator (HMGB1+ CysC) showed stronger predictive ability [0.785 (95% CI: 0.724–0.845)] than serum HMGB1 (p = 0.007). Furthermore, we further explored the potential predictive value of HMGB1 in the non-cirrhotic (n = 121) and cirrhotic subgroups (n = 130) of all HBV-ACLF patients as shown in Supplementary Figure S3. In the non-cirrhotic or cirrhotic subgroups, HMGB1 and CysC had certain early predictive value for AKI development in HBV-ACLF patients, and there was no statistically significant difference between the two single indicators. In the non-cirrhotic subgroup, HMGB1 and CysC had greater predictive value as joint indicators than serum HMGB1 (p = 0.004). Concurrently, in the cirrhotic subgroup, HMGB1 and CysC had greater predictive value as joint indicators than serum CysC (p = 0.039). Partial ROC-related information and data for each predictor variable are shown in Supplementary Table S3 and Supplementary Table S4.
Finally, given the essential role of HMGB1 in the development of AKI, we sought to analyze the possible relationship through three-step adjustment by building unadjusted or adjusted models. Model I unadjusted; Model II adjusted for age; Model III adjusted for the covariates in Model II plus HE, infection, and ascites; Model IV adjusted for the covariates in Model III plus Na, TBIL, and INR. These selected covariates were identified through the univariate linear regression analysis in Table 1. After adjusting for the effects of several potential confounders, the role of HMGB1 was not meaningfully affected, with each increase in serum HMGB1 levels, there was an increased risk of AKI development (p-value for trend = 0.005) as shown in Table 4.
TABLE 4. Adjusted effects of HMGB1 level as a predictor of AKI among individuals with HBV-ACLF during the 90-day follow-up period (n = 251).
As shown in Supplementary Figure S4, to further assess the role of HMGB1 levels in AKI development in HBV-ACLF individuals, the subgroup results were analyzed. Interactions were observed in all subgroups, including age, cirrhosis, HE, infection, ascites, TBIL, INR, sCr, CysC, lactate, MELD, Child-Pugh and CLIF-C ACLF scores except Na.
Potential prognostic value of HMGB1 for AKI in individuals with ACLF
Next, we explored the prognostic value of HMGB1 in AKI patients, and an AUROC analysis was used to compare the performance of HMGB1, MELD, Child-Pugh and CLIF-C ACLF scores separately according to 90-day prognosis. From Figure 5A, discrimination was relatively higher for HMGB1 [0.727 (95% CI: 0.648–0.807) than the AUROC of the MELD [0.663 (95% CI: 0.574–0.751)], Child-Pugh [0.722 (95% CI: 0.645–0.800)], and CLIF-C ACLF [0.708 (95% CI: 0.627–0.789)] scores, although these differences did not reach statistical significance. Partial ROC-related information and data for each predictor variable are listed in Figure 5B. Furthermore, we further explored the potential prognostic value of HMGB1 in the non-cirrhotic (n = 68) and cirrhotic subgroups (n = 85) of all AKI patients. It is interesting that the AUROC of HMGB1 was above 0.7 in the cirrhotic subgroup [0.800 (95% CI: 0.706–0.894)] and below 0.7 in the non-cirrhotic subgroup [0.644 (95% CI: 0.512–0.776)] (Supplementary Figure S5 and Supplementary Table S5).
FIGURE 5. Performance of serum HMGB1 level, MELD, Child-Pugh and CLIF-C ACLF scores were compared in predicting 90-day adverse outcomes in AKI patients using AUROC (A). (B) Partial ROC-related information and data for each predictor variable (AUROC, area under the receiver operating characteristic curve; AKI, acute kidney injury; HMGB1, high-mobility group box 1; MELD, model for end-stage liver disease; Child-Pugh, Child-Turcotte-Pugh; CLIF-C ACLF; CLIF-Consortium ACLF).
Furthermore, multivariate regression analyses were performed to identify independent risk factors related to 90-day poor outcomes of all AKI individuals. In unadjusted or adjusted models, we found an association between HMGB1 level as a continuous or categorical variable and 90-day prognosis. Importantly, the risk of 90-day adverse outcomes increased with higher HMGB1 values (p-value for trend = 0.004) in Model Ⅳ adjusted for age, infection, and CLIF-C ACLF (Supplementary Table S6).
To further assess the role of HMGB1 levels on 90-day outcomes of AKI individuals, the subgroup results were analyzed. According to the 90-day prognosis, interactions were observed in all subgroups, including age, cirrhosis, HE, infection, ascites, AKI stage, TBIL, INR, sCr, lactate, MELD, Child-Pugh and CLIF-C ACLF scores (p < 0.05 Supplementary Figure S6).
Discussion
AKI is a common and life-threatening complication in end-stage liver disease. The 30-day mortality in AKI patients with ACLF remains very high (approximately 50%) (Maiwall et al., 2017). The early diagnosis and timely intervention of AKI in ACLF is extremely important to improve prognosis. Circulating HMGB1 has pleiotropic physiologic effects in several diseases and has potential as a biomarker for relevant diseases. Our study investigated the role of HMGB1 in AKI in individuals with HBV-ACLF. There were three major findings in our study. First, AKI is quite common among patients with HBV-ACLF. Patients with AKI have a significantly worse prognosis than those without AKI. Elevated HMGB1 levels are associated with AKI stage in patients with HBV-ACLF. Second, HMGB1 can be used as an early AKI predictor in patients with ACLF. Finally, serum HMGB1 can serve as a potential prognostic biomarker for AKI in individuals with HBV-ACLF. HMGB1 is a highly conserved nucleoprotein implicated in several bioactivities. During the past decade, many researchers have proposed that HMGB1 acts as a key signaling molecule in kidney and liver diseases (Chen et al., 2014; Chen Q et al., 2016; Zhao et al., 2020). Among patients with HBV-ACLF, high levels of HMGB1 correlate closely with coagulation and brain injury. A recent study also found that HMGB1 levels were significantly increased in patients with HE after transjugular intrahepatic portosystemic shunt (Chen et al., 2020). However, we found no association between HMGB1 levels and liver injury, which may be due to inherent differences in the populations studied. However, Cai J et al. found higher HMGB1 expression in HBV-ACLF patients (Cai et al., 2016).
Previous studies have suggested that CysC can be used for early prediction of AKI (Hou et al., 2021; Ye et al., 2021; Galic et al., 2022), including in patients with HBV-ACLF complicated by AKI (Wan et al., 2013). A recent study demonstrated that HMGB1 is a potential preprocedural predictor of contrast-induced AKI after percutaneous coronary artery intervention (Mo et al., 2022). When compared to CysC, serum HMGB1 had a stronger predictive potential for AKI or severe AKI (relatively higher AUC values). The most meaningful AUC value for the joint indicator (HMGB1+CysC) had excellent sensitivity (0.724) and specificity (0.706) for predicting AKI and sensitivity (0.947) and specificity (0.500) for predicting severe AKI (Supplementary Table S3). On multivariable analysis, baseline HMGB1 was an independent predictor of AKI in HBV-ACLF patients after adjusting for all confounding factors. Multiple studies have shown that serum HMGB1 is predictive of the outcomes of end-stage liver and renal disease (Leelahavanichkul et al., 2011; Zhang et al., 2014; Dai et al., 2021; Dos Santos Pinheiro et al., 2022). Here, we showed that serum HMGB1 might be a prognostic biomarker for AKI among HBV-ACLF patients, while it was superior in the cirrhotic subgroup. HMGB1, infection, and Child-Pugh score are considered to be independent risk factors for 90-day adverse clinical outcomes of AKI among HBV-ACLF patients. A large retrospective multicenter study of 1032 patients showed that the MELD score, ascites, sepsis/infection, and acute variceal bleeding were independent risk factors predicting the development of AKI in ACLF patients with underlying cirrhosis in China (Zang et al., 2016). Concurrently, elevated serum HMGB1 levels appear to be clearly associated with sterile inflammation and infection (Kang et al., 2017). In summary, we believe that the serum HMGB1 level shows great value in facilitating clinical decision-making for the clinical management of patients with HBV-ACLF complicated with AKI. In our study, patients with AKI had a greater risk of complications such as HE, infection, and ascites than those without AKI. These findings were similar to our previous study (Yuan et al., 2017), which suggested that spontaneous bacterial peritonitis had some prognostic value for patients with AKI (Yuan et al., 2017). Therefore, early proactive intervention for complications can significantly improve the prognosis of ACLF-associated AKI patients.
In addition, ACLF-associated AKI may involve not only tissue hypoperfusion but, more importantly, capillary leukocyte infiltration, vascular microthrombosis, and cell apoptosis (Zaccherini et al., 2021). HMGB1 is a non-histone chromosomal protein known to regulate gene transcription and DNA repair by binding to DNA or chromatin via receptor for advanced glycation end-products (RAGE) and Toll-like receptor 4 (TLR4). One clinical study showed that HMGB1, RAGE and interleukin (IL)-17 expression is increased in the liver tissue of HBV-ACLF patients, and the HMGB1 and RAGE interaction may contribute to inflammation of the liver, enhancing the expression of IL-17 (Jhun et al., 2015). Intriguingly, a recent experimental study showed that compared with wild-type control mice, mice lacking both RAGE and soluble isoforms (sRAGE) had more severe AKI, with enhanced renal tubular damage, macrophage infiltration, and fibrosis (Miyagawa et al., 2022). However, additional mechanisms or pathways are certainly involved, for example, autophagy-mediated HMGB1 release or the novel epigenetic role of HMGB1. The precise molecular mechanisms remain to be elucidated in future work. Increasing evidence indicates that HMGB1 targeting is a potential AKI therapy. A more recent animal study found that ulinastatin combined with thrombomodulin significantly improved LPS-induced liver and kidney pathological structure and functional injury; additionally, these improvements in survival of endotoxic rats by inhibiting liver and kidney cell apoptosis, promoting proliferation, and inhibiting inflammation and oxidative injury (Zhang et al., 2022). To date, numerous scholars have concentrated on HMGB1 neutralizing antibodies, and cytokines released by immune cells can be reduced by inhibiting HMGB1-TLR4 binding in vitro (Venereau et al., 2012). In drug-induced acute liver failure and ACLF models, there were significant decreases in serum HMGB1 levels and a reduction in liver injury after intervention with the anti-HMGB1 antibody, which ultimately contributed to the significantly improved animal survival (Takano et al., 2010; Li et al., 2013). Nevertheless, translation of these promising findings into the clinical realm remains to be done.
Nevertheless, despite many strengths, the current study has some possible limitations. First, the data of our cohort were collected in a consistent fashion. However, this is a single-center study; therefore, the results should be validated in multicenter studies. In addition, we analyzed HMGB1 levels only at baseline, and we were unable to evaluate dynamic changes in HMGB1 levels over time. Finally, the design of our study does not allow us to draw mechanistic conclusions about the role of HMGB1 in the underlying pathophysiology of ACLF. Further studies are required to elucidate the pathophysiological and molecular regulatory mechanism of HMGB1 in ACLF-associated AKI.
Conclusion
In conclusion, our study suggests that elevated HMGB1 levels are implicated in AKI stage in HBV-ACLF patients and that the determination of serum HMGB1 levels may be of important early prediction and prognostic implications in HBV-ACLF patients with AKI. To confirm our findings, further investigation in larger real-world studies is warranted to confirm the clinical significance and prognostic value of HMGB1 in ACLF-associated AKI.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Shanghai Public Health Clinical Center. The patients/participants provided their written informed consent to participate in this study.
Author contributions
YL, WY, and MF contributed equally and share first authorship. study concept, YG, ML, ZQ, YL, WY, and MF; study design, YG, ML, ZQ, YL, WY, and MF; enrollment of patients and data collection, HG, XZ, XM, LJ, YZ, YG, and JW; analysis and interpretation of data, YL, WY, and MF; critical revision of the manuscript for important intellectual content, YG, ML, and ZQ; Drafting of the manuscript, YL, WY, and MF. All authors have read, corrected, and approved the final version of the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81874436, 82074155, 81973773, China); the “Shuguang Program” supported by the Shanghai Education Development Foundation and Shanghai Municipal Education Commission (No.18SG39, China); the Program of Shanghai Academic/Technology Research Leader (No.20XD1423500, China); the Clinical Research Plan of SHDC (No. SHDC2020CR3089B, SHDC2020CR1037B, China); and the Shanghai Key Clinical Specialty Construction Project (No.shslczdzk01201, China).
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.
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/fphar.2022.1031790/full#supplementary-material
References
Angeli, P., Ginès, P., Wong, F., Bernardi, M., Boyer, T. D., Gerbes, A., et al. (2015). Diagnosis and management of acute kidney injury in patients with cirrhosis: revised consensus recommendations of the international Club of ascites. J. hepatology 62 (4), 968–974. doi:10.1016/j.jhep.2014.12.029
Cai, J., Han, T., Nie, C., Jia, X., Liu, Y., Zhu, Z., et al. (2016). Biomarkers of oxidation stress, inflammation, necrosis and apoptosis are associated with Hepatitis B-related acute-on-chronic liver failure. Clin. Res. hepatology gastroenterology 40 (1), 41–50. doi:10.1016/j.clinre.2015.06.009
Chen, R., Yi, P., Zhou, R., Xiao, M. F., Huang, Z. B., Tang, D. L., et al. (2014). The role of HMGB1-RAGE axis in migration and invasion of hepatocellular carcinoma cell lines. Mol. Cell. Biochem. 390, 271–280. doi:10.1007/s11010-014-1978-6
Chen, C., Yang, X., Lei, Y., Zha, Y., Liu, H., Ma, C., et al. (2016). Urinary biomarkers at the time of AKI diagnosis as predictors of progression of AKI among patients with acute cardiorenal syndrome. Clin. J. Am. Soc. Nephrol. CJASN 11 (9), 1536–1544. doi:10.2215/CJN.00910116
Chen, Q., Guan, X., Zuo, X., Wang, J., and Yin, W. (2016). The role of high mobility group box 1 (HMGB1) in the pathogenesis of kidney diseases. Acta Pharm. Sin. B 6 (3), 183–188. doi:10.1016/j.apsb.2016.02.004
Chen, Q., Zhang, Y., Yue, Z. D., Zhao, H. W., Wang, L., Fan, Z. H., et al. (2020). High-mobility group protein B1: A predictive biomarker for hepatic encephalopathy after transjugular intrahepatic portosystemic shunt. J. Hepatobiliary Pancreat. Sci. 27 (8), 522–530. doi:10.1002/jhbp.770
Dai, Q., Xie, D., Zhang, C., Zhu, L., Xu, Y., Li, K., et al. (2021). Osthole blocks HMGB1 release from the nucleus and confers protective effects against renal ischemia-reperfusion injury. Front. Physiol. 12, 735425. doi:10.3389/fphys.2021.735425
de Oliveira Gomes, C., de Andrade, M., Guedes, L., Rocha, H. C., Guimaraes, R. G., Carvalho, F. A. C., et al. (2020). Evaluation of the biomarkers HMGB1 and IL-6 as predictors of mortality in cirrhotic patients with acute kidney injury. Mediat. Inflamm. 2020, 2867241. doi:10.1155/2020/2867241
Dos Santos Pinheiro, C., de Oliveira Gomes, C. G., Ribeiro Lima Machado, C., Guedes, L. R., Rocha, H. C., Guimaraes, R. G., et al. (2022). Performance of high mobility protein group 1 and interleukin-6 as predictors of outcomes resulting from variceal bleeding in patients with advanced chronic liver disease. Inflammation 45 (2), 544–553. doi:10.1007/s10753-021-01565-1
European Association for the Study of the Liver (2010). EASL clinical practice guidelines on the management of ascites, spontaneous bacterial peritonitis, and hepatorenal syndrome in cirrhosis. J. Hepatol. 53 (3), 397–417. doi:10.1016/j.jhep.2010.05.004
Ferguson, M., Vaidya, V., Waikar, S., Collings, F. B., Sunderland, K. E., Gioules, C. J., et al. (2010). Urinary liver-type fatty acid-binding protein predicts adverse outcomes in acute kidney injury. Kidney Int. 77 (8), 708–714. doi:10.1038/ki.2009.422
Galic, S., Milosevic, D., Filipovic-Grcic, B., Rogić, D., Vogrinc, Ž., Ivančan, V., et al. (2022). Early biochemical markers in the assessment of acute kidney injury in children after cardiac surgery. Ther. Apher. Dial. 26 (3), 583–593. doi:10.1111/1744-9987.13736
Gao, Z., Lu, L., and Chen, X. (2021). Release of HMGB1 in podocytes exacerbates lipopolysaccharide-induced acute kidney injury. Mediat. Inflamm. 2021, 5220226. doi:10.1155/2021/5220226
Goodwin, G., Sanders, C., and Johns, E. (1973). A new group of chromatin-associated proteins with a high content of acidic and basic amino acids. Eur. J. Biochem. 38 (1), 14–19. doi:10.1111/j.1432-1033.1973.tb03026.x
Gu, J., Xu, R., Li, Y., Zhang, J., and Wang, S. (2016). MicroRNA-218 modulates activities of glioma cells by targeting HMGB1. Am. J. Transl. Res. 8 (9), 3780–3790. doi:10.3892/ol.2020.11714
Hou, Y., Deng, Y., Hu, L., He, L., Yao, F., Wang, Y., et al. (2021). Assessment of 17 clinically available renal biomarkers to predict acute kidney injury in critically ill patients. J. Transl. Int. Med. 9 (4), 273–284. doi:10.2478/jtim-2021-0047
Jha, P., Jha, A., Dayal, V., Jha, S., and Kumar, A. (2021). Baseline serum cystatin C as a marker of acute kidney injury in patients with acute-on-chronic liver failure. Indian J. gastroenterology 40 (6), 563–571. doi:10.1007/s12664-021-01201-8
Jhun, J., Lee, S., Kim, H., Her, Y. M., Byun, J. K., Choi, J. Y., et al. (2015). HMGB1/RAGE induces IL-17 expression to exaggerate inflammation in peripheral blood cells of Hepatitis B patients. J. Transl. Med. 13, 310. doi:10.1186/s12967-015-0663-1
Jiang, W., Hu, Y., Sun, Y., Shen, Y., and Xun, Y. (2020). Prevalence and short-term outcome of acute kidney injury in patients with acute-on-chronic liver failure: A meta-analysis. J. viral Hepat. 27 (8), 810–817. doi:10.1111/jvh.13287
Kang, R. X. Y., Zhang, Q., Hou, W., Jiang, Q., Zhu, S., Liu, J., et al. (2017). Intracellular HMGB1 as a novel tumor suppressor of pancreatic cancer. Cell Res. 27 (7), 916–932. doi:10.1038/cr.2017.51
Leelahavanichkul, A., Huang, Y., Hu, X., Zhou, H., Tsuji, T., Chen, R., et al. (2011). Chronic kidney disease worsens sepsis and sepsis-induced acute kidney injury by releasing High Mobility Group Box Protein-1. Kidney Int. 80 (11), 1198–1211. doi:10.1038/ki.2011.261
Li, X., Wang, L., Wang, L., Han, X., Yang, F., Gong, Z., et al. (2013). Blockade of high-mobility group box-1 ameliorates acute on chronic liver failure in rats. Inflamm. Res. 62 (7), 703–709. doi:10.1007/s00011-013-0624-1
Li, H., Chen, L., Zhang, N., Li, S. T., Zeng, B., Pavesi, M., et al. (2016). Characteristics, diagnosis and prognosis of acute-on-chronic liver failure in cirrhosis associated to hepatitis B. Sci. Rep. 6, 25487. doi:10.1038/srep25487
Maiwall, R., Sarin, S., Kumar, S., Jain, P., Kumar, G., Bhadoria, A. S., et al. (2017). Development of predisposition, injury, response, organ failure model for predicting acute kidney injury in acute on chronic liver failure. Liver Int. 37 (10), 1497–1507. doi:10.1111/liv.13443
Miyagawa, T., Iwata, Y., Oshima, M., Ogura, H., Sato, K., Nakagawa, S., et al. (2022). Soluble receptor for advanced glycation end products protects from ischemia- and reperfusion-induced acute kidney injury. Biol. open 11 (1), bio058852. doi:10.1242/bio.058852
Mo, C., Ma, X., Jian, W., Huang, Q., Zheng, W., Yang, Z., et al. (2022). High mobility group box 1 and homocysteine as preprocedural predictors for contrast-induced acute kidney injury after percutaneous coronary artery intervention. Int. Urol. Nephrol. 54 (7), 1663–1671. doi:10.1007/s11255-021-03050-y
Nickolas, T., Schmidt-Ott, K., Canetta, P., Forster, C., Singer, E., Sise, M., et al. (2012). Diagnostic and prognostic stratification in the emergency department using urinary biomarkers of nephron damage: a multicenter prospective cohort study. J. Am. Coll. Cardiol. 59 (3), 246–255. doi:10.1016/j.jacc.2011.10.854
Sarin, S. K., Kumar, A., Almeida, J. A., Chawla, Y. K., Fan, S. T., Garg, H., et al. (2009). Acute-on-chronic liver failure: consensus recommendations of the asian pacific association for the study of the liver (APASL). Hepatol. Int. 3 (1), 269–282. doi:10.1007/s12072-008-9106-x
Sarin, S. K., Kedarisetty, C. K., Abbas, Z., Amarapurkar, D., Bihari, C., Chan, A. C., et al. (2014). Acute-on-chronic liver failure: consensus recommendations of the asian pacific association for the study of the liver (APASL) 2014. Hepatol. Int. 8 (4), 453–471. doi:10.1007/s12072-014-9580-2
Sarin, S., Choudhury, A., Sharma, M., Maiwall, R., Al Mahtab, M., Rahman, S., et al. (2019). Correction to: Acute-on-chronic liver failure: consensus recommendations of the asian pacific association for the study of the liver (APASL): an update. Hepatol. Int. 13 (6), 826–828. doi:10.1007/s12072-019-09980-1
Takano, K., Shinoda, M., Tanabe, M., Miyasho, T., Yamada, S., Ono, S., et al. (2010). Protective effect of high-mobility group box 1 blockade on acute liver failure in rats. Shock (Augusta, Ga) 34 (6), 573–579. doi:10.1097/SHK.0b013e3181df0433
Venereau, E., Casalgrandi, M., Schiraldi, M., Antoine, D. J., Cattaneo, A., De Marchis, F., et al. (2012). Mutually exclusive redox forms of HMGB1 promote cell recruitment or proinflammatory cytokine release. J. Exp. Med. 209 (9), 1519–1528. doi:10.1084/jem.20120189
Vilela, E., Pinheiro, C., Saturnino, S., Gomes, C., Nascimento, V., and Andrade, M. (2018). Evaluation of the behavior of levels of Hmgb1 and Il6 as predictors of infection, acute kidney injury and mortality in cirrhotic patients with variceal bleeding. Arq. Gastroenterol. 55 (4), 338–342. doi:10.1590/S0004-2803.201800000-75
Vilstrup, H., Amodio, P., Bajaj, J., Cordoba, J., Ferenci, P., Mullen, K. D., et al. (2014). Hepatic encephalopathy in chronic liver disease: 2014 practice guideline by the American association for the study of liver diseases and the European association for the study of the liver. Hepatology 60 (2), 715–735. doi:10.1002/hep.27210
Wan, Z. H., Wang, J. J., You, S. L., Liu, H. L., Zhu, B., Zang, H., et al. (2013). Cystatin C is a biomarker for predicting acute kidney injury in patients with acute-on-chronic liver failure. World J. Gastroenterol. 19 (48), 9432–9438. doi:10.3748/wjg.v19.i48.9432
Wu, T., Li, J., Shao, L., Xin, J., Jiang, L., Zhou, Q., et al. (2018). Development of diagnostic criteria and a prognostic score for Hepatitis B virus-related acute-on-chronic liver failure. Gut 67 (12), 2181–2191. doi:10.1136/gutjnl-2017-314641
Ye, Z., Liu, H., Zhao, B., Fu, H., Li, Y., and Chen, L. (2021). Correlation and diagnostic value of serum cys-C, RBP4, and NGAL with the condition of patients with traumatic acute kidney injury. Evid. Based Complement. Altern. Med. 2021, 4990941. doi:10.1155/2021/4990941
Yuan, W., Zhang, Y. Y., Zhang, Z. G., Zou, Y., Lu, H. Z., and Qian, Z. P. (2017). Risk factors and outcomes of acute kidney injury in patients with hepatitis B virus-related acute-on-chronic liver failure. Am. J. Med. Sci. 353 (5), 452–458. doi:10.1016/j.amjms.2017.03.005
Zaccherini, G., Weiss, E., and Moreau, R. (2021). Acute-on-chronic liver failure: Definitions, pathophysiology and principles of treatment. JHEP Rep. innovation hepatology 3 (1), 100176. doi:10.1016/j.jhepr.2020.100176
Zakiyanov, O., Kriha, V., Vachek, J., Zima, T., Tesar, V., and Kalousova, M. (2013). Placental growth factor, pregnancy-associated plasma protein-A, soluble receptor for advanced glycation end products, extracellular newly identified receptor for receptor for advanced glycation end products binding protein and high mobility group box 1 levels in patients with acute kidney injury: a cross sectional study. BMC Nephrol. 14, 245. doi:10.1186/1471-2369-14-245
Zang, H., Liu, F., Liu, H., You, S., Zhu, B., Wan, Z., et al. (2016). Incidence, risk factors and outcomes of acute kidney injury (AKI) in patients with acute-on-chronic liver failure (ACLF) of underlying cirrhosis. Hepatol. Int. 10 (5), 807–818. doi:10.1007/s12072-016-9756-z
Zhang, L., Han, J., Wu, H., Liang, X., Zhang, J., Li, J., et al. (2014). The association of HMGB1 expression with clinicopathological significance and prognosis in hepatocellular carcinoma: a meta-analysis and literature review. PLoS One 9 (10), e110626. doi:10.1371/journal.pone.0110626
Zhang, X., Su, C., Zhao, S., Li, J., and Yu, F. (2022). Combination therapy of Ulinastatin with Thrombomodulin alleviates endotoxin (LPS) - induced liver and kidney injury via inhibiting apoptosis, oxidative stress and HMGB1/TLR4/NF-κB pathway. Bioengineered 13 (2), 2951–2970. doi:10.1080/21655979.2021.2024686
Keywords: hepatitis B virus-related acute-on-chronic liver failure, acute kidney injury, high mobility group protein 1, prediction, prognosis
Citation: Liu Y, Yuan W, Fang M, Guo H, Zhang X, Mei X, Zhang Y, Ji L, Gao Y, Wang J, Qian Z, Li M and Gao Y (2023) Determination of HMGB1 in hepatitis B virus-related acute-on-chronic liver failure patients with acute kidney injury: Early prediction and prognostic implications. Front. Pharmacol. 13:1031790. doi: 10.3389/fphar.2022.1031790
Received: 13 September 2022; Accepted: 13 December 2022;
Published: 13 January 2023.
Edited by:
Runping Liu, Beijing University of Chinese Medicine, ChinaReviewed by:
Xinhua Li, Third Affiliated Hospital of Sun Yat-sen University, ChinaTianyan Chen, The First Affiliated Hospital of Xi’an Jiaotong University, China
Copyright © 2023 Liu, Yuan, Fang, Guo, Zhang, Mei, Zhang, Ji, Gao, Wang, Qian, Li and Gao. 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: Zhiping Qian, cWlhbnpoaXBpbmdAc2hhcGhjLm9yZw==; Man Li, bGltYW4xMjEwMDBAc2h1dGNtLmVkdS5jbg==; Yueqiu Gao, Z2FveXVlcWl1QHNodXRjbS5lZHUuY24=
†These authors have contributed equally to this work