Abstract

Background. T-helper 17 (Th17) and CD4+CD25+ T-regulatory (Treg) cells play important roles in the pathogenesis of hepatitis B virus-related acute-on-chronic liver failure (HBV-ACLF). This study is aimed at investigating shifts in Treg/Th17 balance in the peripheral blood of HBV-ACLF patients at different disease stages. Methods. Sixty HBV-ACLF patients, admitted to the First Hospital of Hunan University of Chinese Medicine, China, including early-stage (), middle-stage (), and late-stage patients (), were enrolled in the study. In addition, 20 patients with chronic hepatitis B and 20 healthy volunteers were also included in the study as controls. Flow cytometry, cytometric bead array, and quantitative real-time PCR protocols were used to evaluate the expression of Treg and Th17 cells as well as of related cytokines. Results. The levels of Th17 cells and their effectors interleukin- (IL-) 17A, IL-23, and tumor necrosis factor-α increased with disease progression. Similarly, Treg cells and their effector cytokines transforming growth factor-β and IL-10 also increased. Although Treg and Th17 levels were positively correlated, the latter were always at higher numbers. Noteworthy, the Treg/Th17 ratio gradually decreased and was negatively correlated with ACLF severity. FoxP3 levels in the peripheral blood gradually decreased with ACLF progression, whereas ROR-γt gradually increased. Serum c-reactive protein, procalcitonin, and lipopolysaccharide were also upregulated with disease progression and positively correlated with Th17 abundance. Further, Th17, IL-17A, and IL-23 were independent risk factors for ACLF. A prognostic model for HBV-ACLF was established, with a correct prediction rate of 90.00% (54/60). Conclusion. Treg/Th17 imbalance occurs throughout the pathogenic course of HBV-ACLF, with an imbalance shift toward Th17. Hence, the Th17-mediated inflammatory response drives HBV-ACLF-associated inflammation and supports the pathological mechanisms of liver failure.

1. Introduction

Acute-on-chronic liver failure (ACLF) is a clinical condition characterized by an exacerbation of acute jaundice and coagulation dysfunction caused by various factors on the basis of chronic liver disease [1, 2]. ACLF progresses rapidly and has a high short-term mortality rate. Based on the severity of clinical symptoms, the disease course can be divided into early, middle, and late stages. In China, hepatitis B virus- (HBV-) related ACLF is the most common type of liver failure, accounting for over 80% of all liver failure cases [3]. The pathogenesis of HBV-ACLF is a complex immunopathological process, which remains poorly known. A growing body of evidence suggests that persistent inflammatory response and immune dysregulation are the core mechanisms underlying HBV-ACLF [4, 5]. An antiviral inflammatory response is observed in the early and advanced stages of the disease, followed by endotoxin-driven inflammation during the recovery stage, which drives an imbalance of anti- and proinflammatory factors that are closely related to prognosis [6].

CD4+ T cells play a central role in the adaptive immune response during liver failure. In particular, T-helper 17 (Th17) and T-regulatory (Treg) cells modulate the inflammatory response in liver failure; thus, these cell populations have become core research targets [7]. Tregs express the transcription factor forkhead box protein 3 (FoxP3) and mainly secrete interleukin- (IL-) 10, transforming growth factor-β (TGF-β), and IL-35, which are anti-inflammatory cytokines. Th17 cells express the transcription factor retinoid orphan receptor γt (ROR-γt) and mainly secrete IL-17A, IL-23, tumor necrosis factor-α (TNF-α), and IL-6, which are proinflammatory [8]. Thus, Treg and Th17 cells have opposing functions. However, they are also interdependent and regulate each other. Several clinical studies have shown that the Treg/Th17 balance is associated with the severity of HBV-ACLF and that Treg/Th17 dysregulation can reflect aberrant immune function in these patients. Moreover, increased or decreased levels of associated cytokines are related to the degree of liver injury, thus having a guiding role in clinical prognosis [915]. However, there is no consensus on the specific Treg/Th17 ratio changes observed in the peripheral blood of HBV-ACLF patients.

Therefore, in the present study, we determined the frequency of Treg/Th17 cells and the expression of associated transcription factors and cytokines in the peripheral blood of patients with HBV-ACLF in early, middle, and late disease stages. The overall aim of the study was to investigate the regulatory role of the Treg/Th17 axis at different stages of viral- and endotoxin-mediated inflammation during HBV-ACLF, thereby providing new directions for the clinical treatment of liver failure.

2. Materials and Methods

2.1. Patients

A total of 60 patients with HBV-ACLF, admitted to the First Hospital of Hunan University of Chinese Medicine, China, between February 2018 and December 2018, were enrolled, including early-stage (), middle-stage (), and advanced-stage patients (). Twenty patients with chronic hepatitis B (CHB) on the immunoactive phase (including HBeAg-negative CHB or HBeAg-positive CHB) and 20 healthy individuals were recruited as controls. All patients were aged between 18 and 65 years and were negative for hepatitis A, C, D, E, or HIV infections. Enrolled patients did not have hepatocellular carcinoma or any extrahepatic malignancy. None of the patients had history of alcohol abuse, hepatotoxic drug use, chronic autoimmune, metabolic, renal, cardiac, or pulmonary diseases, nor had received antiviral treatment or immunosuppressive therapy during the previous 3 months. The study protocol was approved by the Ethics Committee of the First Hospital of Hunan University of Chinese Medicine (No. PJ20170720), and written informed consent was obtained from all participants.

2.2. Diagnostic Criteria of ACLF Staging

Early-stage ACLF was described as follows: , or international normalized ratio (INR) ranging from 1.5 to 1.9, no hepatic encephalopathy (HE), or other complications. Middle-stage ACLF was described as follows: , PTA within 20–30% or international normalized ratio (INR) ranging from 1.9 to 2.6, I-II HE and/or ascites, and infection. Late-stage ACLF was described as follows: , or , and III–IV HE (or hepatorenal syndrome, gastrointestinal bleeding, and severe infection).

2.3. Prognostic Criteria for HBV-ACLF Patients

Patients with HBV-ACLF () were divided into improvement () and deterioration () groups based on their condition at discharge or occurrence of death at the hospital. In the improvement group, symptoms and signs of liver failure disappeared or improved significantly, with a considerable improvement in liver function indexes (serum bilirubin and alanine aminotransferase decreased by over 50%). In the deterioration group, symptoms and signs did not improve, and deterioration was observed, along with organ failure, worsening of other vital signs, or death after ineffective treatment.

2.4. Clinical Data Collection

Demographic data, laboratory measurements (e.g., WBC count, lymphocyte (LYM) count, serum albumin (ALB), total bilirubin (TBil), alanine aminotransferase (ALT), aspartate aminotransferase (AST), PTA, prothrombin time (PT), INR, and creatinine levels), and information on HBV infection biomarker levels, virus infection time, and HBV DNA levels were collected. Laboratory data, HBV infection biomarkers, and HBV DNA levels were obtained at week 0.

2.5. Flow Cytometry

Peripheral blood mononuclear cells (PBMCs) were isolated from 2 mL of fresh heparinized blood via Ficoll-Hypaque (4A Biotech, Beijing, China) density gradient centrifugation. RPMI1640 (150 μL) and cell stimulation cocktail (0.5 μL) were added to the isolated PBMCs, evenly mixed, and incubated in a 37°C water bath for 6 h. The activated PBMCs were then incubated with anti-CD4, anti-CD8, anti-CD3, and anti-CD25 antibodies (Becton, Dickinson and Company, San Jose, CA, USA) at room temperature for 20 min. Anti-IL-17A (5 μL) and anti-FoxP3 (5 μL) antibodies were then added, followed by incubation at room temperature for 30 min. After staining, the cells were washed twice with PBS and fixed in 1% paraformaldehyde. The stained cells were subsequently processed on an Accuri C6 FCM (Becton, Dickinson and Company) for analysis. All data were analyzed using the Accuri C6 FCM software.

2.6. Cytometric Bead Array (CBA)

Blood samples for the measurement of cytokines were collected early in the morning. The CBA assay (CBA immunoassay kit; Becton, Dickinson and Company) was used for quantitative determination of serum cytokines (Human Inflammatory Kit), according to the manufacturer’s instructions. The Inflammatory CBA Kit comprises microbeads coupled to mAbs against IL-6, IL-10, TGF-β, IL-17A, and IL-23. A secondary phycoerythrin-labeled anti-cytokine antibody was used, and the concentration of individual cytokines was determined based on fluorescence intensity. Data were acquired on a FACSVerse flow cytometer (Becton, Dickinson and Company). Sample analysis was performed using BD FCAP Array 3.0 software (Becton, Dickinson and Company).

2.7. Quantitative Real-Time PCR (qRT-PCR)

Total RNA was extracted from PBMCs using the Redzol reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s protocol. The extracted RNA was reverse-transcribed to first-strand cDNA using the Transcriptor First Strand cDNA Synthesis Kit (Promega, Madison, WI, USA). Gene expression was quantified on an ABI 7500 thermal cycler (Applied Biosystems, Waltham, MA, USA) using the SYBR Green PCR Master Mix Kit (Promega). The amplification program was as follows: initial denaturation at 94°C for 5 min, followed by 40 cycles at 94°C for 1 min, 58°C for 1 min, 72°C for 1 min, and 83°C for 1 min. β-Actin was used as internal control. The primer sequences used are shown in Table 1. For each sample, PCR was performed twice with triplicates for each sample, and the data were analyzed using the thermal cycler software to calculate the ΔCt value.

2.8. Model for End-Stage Liver Disease Score (MELD Score)

The MELD score was determined as per the following equation:

.

2.9. Statistical Analysis

Statistical analyses were performed using the SPSS software (ver. 25.0; IBM Corp., Armonk, NY, USA), and the figures were produced using Prism 8.0 (GraphPad Software, San Diego, CA, USA). All data are expressed as the or median (minimum and maximum) for continuous variables and as a number (percentage) for categorical variables. Differences between three or more groups were analyzed using one-way ANOVA, and LSD multiple comparison was used for comparison between groups. The correlations of CRP, PCT, and LPS with Th17 and Treg frequency were determined via Pearson’s or Spearman’s correlation analyses. A two-sided value < 0.05 was considered to indicate statistical significance.

3. Results

3.1. Participant Characteristics

In total, 60 HBV-ACLF patients, 20 CHB patients, and 20 healthy controls (NC) were included in the study as per the inclusion and exclusion criteria. The number of male participants in each group was greater. There were no significant differences in the percentage of hepatitis B, antigen (HBeAg) positivity, and HBV DNA quantity among the early-stage ACLF (ACLF-E), mid-stage ACLF (ACLF-M), late-stage ACLF (ACLF-L), and CHB patients. There were significant differences in WBC and LYM counts, and serum ALB levels between all groups. WBC count in the ACLF-L group was significantly higher than in the NC and CHB groups (). LYM count and ALB levels in HBV-ACLF patients were significantly lower than those in the NC and CHB groups (). Furthermore, LYM count and ALB levels in HBV-ACLF patients decreased with disease progression. There was no significant difference in the levels of ALT, AST, and TBil between the ACLF-E, ACLF-M, and ACLF-L groups. In contrast, significant differences in PT, INR, and MELD scores were observed between these groups. The demographic and clinical characteristics of the participants are shown in Table 2.

3.2. Treg and Th17 Frequencies Increase with the Progression of ACLF

We assessed the frequency of Treg and Th17 cells by flow cytometry and calculated the Treg/Th17 ratio. The frequency of Th17 cells in the ACLF-E, ACLF-M, and ACLF-L groups was significantly higher than in patients with CHB or normal controls. The frequencies of both cell populations gradually increased with the progression of ACLF, from early through late-stage disease. A significant difference between the late- and early-stage groups (, ) was observed. Moreover, a significant difference was also observed between the late- and middle-stage groups (). Frequency was positively correlated with the severity of ACLF. Overall, Th17 cells were more abundant than Treg cells, resulting in a decrease of the Treg/Th17 ratio with disease progression (Figures 13).

3.3. Serum IL-10, TGF-β, TNF-α, IL-17A, and IL-23 Levels Increase with the Progression of ACLF

CBA was performed to quantify the serum levels of cytokines IL-10, TGF-β, TNF-α, IL-17A, and IL-23 in each participant. Compared with the NC group, the levels of IL-10 and TGF-β were lower in the ACLF-E group, whereas TNF-α and IL-17A were higher. The serum levels of all cytokines were significantly higher in the ACLF-M and ACLF-L groups when compared to those in the NC group. Moreover, the levels of IL-10, TGF-β, TNF-α, IL-17A, and IL-23 were significantly lower in the ACLF-E and ACLF-M groups than those in the CHB group (, ). The serum levels of Treg-secreted cytokines IL-10 and TGF-β gradually increased with ACLF progression from early- to mid- and late-stage disease (, ). Furthermore, the levels of TNF-α, IL-17A, and IL-23 secreted by Th17 cells gradually increased, with a significant difference between the ACLF-L and ACLF-E groups (, ) (Figure 4).

3.4. FoxP3 Expression Decreases with the Progression of ACLF, Whereas That of ROR-γt Increases

qRT-PCR was performed to determine the expression Treg-specific transcription factor FoxP3 and Th17-specific transcription factor ROR-γt in peripheral blood. Compared with the NC group, FoxP3 and ROR-γt were significantly upregulated in the CHB, ACLF-E, ACLF-M, and ACLF-L groups. FoxP3 expression was highest in the CHB group, whereas that of ROR-γt was highest in the ACLF-L group. The expression of FoxP3 gradually decreased with ACLF progression, with the difference between the ACLF-L and ACLF-E groups being statistically significant (). In contrast, the levels of ROR-γt gradually increased with disease progression, with significant differences among the three groups () (Figure 5).

3.5. CRP, PCT, and LPS Levels Increase with ACLF Progression

The dysregulated inflammatory response is a major pathological mechanism in ACLF. Inflammatory mediators are strongly associated with the severity of ACLF [16]. The serum levels of CRP, PCT, and LPS were determined in all subject groups, revealing that all increased with the progression of ACLF. Significant differences were observed between the ACLF-E and ACLF-L groups as well as between the ACLF-M and ACLF-L groups () (Figure 6).

3.6. CRP, PCT, and LPS Serum Levels Are Positively Correlated with Th17, but Not with Treg Frequency

We analyzed the correlation of CRP, PCT, and LPS levels with Th17 and Treg frequency. The serum levels of all three inflammatory biomarkers were positively correlated with Th17 frequency in all HBV-ACLF patients (Figure 7). In contrast, serum CRP, PCT, and LPS exhibited no correlation with Treg levels (Figure 8).

3.7. Binary Logistic Regression Analysis of ACLF Risk Factors

Among the 60 patients with HBV-ACLF included in the study, 43 showed symptom improvement and 17 experienced disease aggravation. The correlation between the frequency of Treg and Th17 cells, as well as IL-10, TGF-β, TNF-α, IL-17A, IL-23, FoxP3, and ROR-γt levels with ACLF prognosis, was analyzed using binary logistic regression (Table 3). The results indicated that Th17 (, , 95% confidence ), IL-17A (, , ), and IL-23 (, , ) levels were independent risk factors for ACLF. The logistic regression model was established as follows: . The Hosmer-Lemeshow test was used to assess goodness of the fit in the logistic regression model (, ). A good fit was observed, and the logistic regression model was statistically significant. Receiver operating characteristic curve analysis revealed an area under the curve of 0.956, which was higher than that of any independent risk factor. The maximum Youden index value was 0.825, with a corresponding sensitivity and specificity of 88.4% and 94.1%, respectively (Figure 9). The 60 patients with HBV-ACLF were selected as a validation cohort. With a prediction probability of 0.5 as threshold, the total correct predication rate was 90.00% (54/60), the positive predictive value was 93.00% (43/40), and the negative predictive value was 82.4% (14/17).

4. Discussion

Recent studies have shown that immune activation plays a major role in the occurrence and development of liver failure. Persistent inflammatory response and immune dysregulation are considered core mechanisms in HBV-ACLF pathophysiology [4, 5] that determine the severity of tissue damage and patient prognosis. Therefore, the identification of novel targets for immunosuppression and to restore the immune balance based on the immune profile of patients at different disease stages is a major goal within HBV-ACLF research. In the current study, we confirmed that the frequencies of Treg and Th17 cells in peripheral blood increase with HBV-ACLF progression. These cell populations were previously described as closely associated with the occurrence and development of HBV-ACLF [1719]. When compared with Tregs, Th17 cells increased to a greater extent in ACLF patients, and the two CD4+ T cell subsets were positively correlated with each other.

The close relationship between Th17 and Treg cells has been acknowledged in recent years as they share a common developmental pathway, with TGF-β being necessary for the differentiation of both. Th17 cells are activated by TGF-β, whereas Treg cells are activated by TGF-β as well as IL-10 [2022]. Th17 cells are potent proinflammatory effectors, participating in liver injury and viral clearance after HBV infection. Studies have reported that IL-17A levels and Th17 cell frequency in the peripheral blood of patients with CHB are positively correlated with liver injury. Circulating Th17 cells mainly accumulate within the liver of CHB patients, and their frequency gradually increases during progression from CHB to ACLF [23]. Treg cells play an important role in inducing and maintaining immune tolerance and terminating the immune response. Herein, a Treg/Th17 imbalance was observed in patients with CHB and HBV-ACLF. Consistent with previous studies, we confirmed that Th17 and Tregs in the peripheral blood of CHB and HBV-ACLF patients were significantly increased compared with healthy controls. Furthermore, the highest Treg frequency was observed in patients with active CHB, whereas Th17 cells were most abundant in HBV-ACLF patients. Previous studies have shown that Treg cell-mediated immunosuppression mainly occurs during the CHB stage, whereas Th17-mediated inflammatory damage is sustained throughout the HBV-ACLF stage [23]. Although Th17 and Treg cells are known to be significantly increased in ACLF patients, our understanding of their frequency and interaction throughout the different disease stages is limited.

In the current study, we observed that the increase of Th17 cell frequency in CHB and HBV-ACLF patients occurred concurrently to an increase in Treg cell frequency. Noteworthy, Th17 frequency in ACLF patients was higher than that of Tregs. More importantly, the Treg/Th17 ratio gradually decreased with the progression of disease. These data highlight the involvement of both immune cell populations in HBV-ALCF. Furthermore, dynamic shifts in both are important for the maintenance of immune balance and for avoiding tissue damage. During HBV-ACLF pathogenesis, inflammation gradually transitions from an antiviral response during the early and middle stages to an endotoxin-mediated inflammatory response in the late disease stages. During HBV-mediated immune injury, immune hyperactivation is observed in ACLF patients. Moreover, the Treg/Th17 balance is dysregulated [24], and aberrant release of Th17 cytokines (IL-17, IL-23, and TNF-α) is observed, resulting in a proinflammatory state. In addition, Treg cells and their effector cytokines TGF-β and IL-10 are also upregulated via compensatory induction. However, their anti-inflammatory counterbalance effect is insufficient to effectively inhibit the aberrant inflammatory response. LPS-driven inflammation plays a major role during the late stage of ACLF. Intestinal mucosal permeability is increased, allowing the entry of gut microbes into circulation, which in turn leads to LPS-induced intestinal endotoxemia [25, 26]. Furthermore, the Treg/Th17 imbalance worsens, with a significant increase in Th17 frequency, and consequent release of effector cytokines IL-17, IL-23, and TNF-α will in turn exacerbate the immune response. Herein, CRP, PCT, and LPS levels were also found to be increased and positively correlated with Th17 frequency. The risk of infection is also considerably greater in these patients; thus, infection may be a major factor aggravating the course of ACLF [27]. Taken together, the observed increase in Treg levels during ACLF may be a negative feedback effect for counterbalancing the Th17 response. However, the greater increase in Th17 cells overwhelms this homeostatic response, driving uncontrolled inflammation and leading to end-stage liver disease.

FoxP3 and ROR-γt are transcription factors specific to Treg and Th17 cells, respectively, and drive their differentiation. Moreover, high concentrations of TGF-β induce FoxP3 expression in immature CD4+ T cells and promote their differentiation into Treg cells [28]. In contrast, TGF-β along with IL-6 or IL-21 can induce ROR-γt expression and promote subsequent differentiation of Th17 cells [29]. In turn, IL-21 inhibits FoxP3 expression and promotes Th17 cell differentiation by modulating TGF-β signal transduction [30]. Some studies have shown that increased FOXP3 expression is involved in the pathogenesis of CHB [31]. In addition, FoxP3 and ROR-γt levels were herein found to be significantly higher in HBV-ACLF patients compared with healthy and CHB subjects, suggesting that the expression of both transcription factors is related to the severity of disease [14, 17]. Hence, Treg and Th17 cells in patients with HBV-ACLF increase with the progression of the disease through the early, middle, and late disease stages. Furthermore, the imbalanced Treg/Th17 was found to shift toward Th17, which agrees with previous findings. In the peripheral blood of patients with HBV-ACLF, ROR-γt expression changed concurrently with the Th17 cell frequency, whereas FoxP3 expression did not change in concert with the Treg levels. Interestingly, FoxP3 expression began to gradually decrease with an increase in Tregs, and no definite correlation was observed between the two. Recent studies have shown that Treg cells can downregulate FoxP3 expression under certain conditions, thus losing their immunosuppressive function [32]. The key factor leading to FoxP3 downregulation is the high level of inflammatory cytokines within the surrounding environment, significantly suppressing Treg function and preventing the Treg-mediated inhibition of inflammation. In contrast, here, we found that ROR-γt expression was positively associated with Th17 dominance. Thus, it is possible that FoxP3 expression is influenced by other genes, as well as differences between the liver and the blood microenvironment.

Finally, we designed a prognostic prediction model for HBV-ACLF based on the herein collected Treg- and Th17-related data. Currently, only few studies have provided HBV-ACLF prognostic models. Binary logistic regression analysis indicated that Th17, IL-17A, and IL-23 are independent risk factors for ACLF. However, Treg cells and related cytokines were not significantly correlated to disease progression, and the relationship between Treg cells and prognosis of ACLF remains unclear. To the best of our knowledge, there is no consensus on whether Treg cells can effectively control the Th17 response in human disease. Importantly, the Treg/Th17 ratio was associated with ACLF severity. The prognostic model established herein furthers our understanding on the relationship between Th17 and Tregs in HBV-ACLF, as well as their association with disease prognosis. These results suggest that Th17 cells and the Treg/Th17 axis may represent candidates for further study in ACLF, especially in orthotopic liver transplantation.

This study is aimed at investigating shifts in Treg/Th17 balance and both pro/anti-inflammatory cytokine level expression in the peripheral blood of HBV-ACLF patients at different disease stages. However, there are some limitations that compromise this study. First, this study only detected the peripheral blood lymphocyte frequency of the subjects, excluding the same lymphocyte population frequency in concomitant biopsy samples. This may not be sufficient to reflect what is seen in the liver. In addition, this is a single-center study. Due to the small sample size, the conclusion may have limited generalizability. Therefore, these conclusions need more cases to be verified.

5. Conclusion

In conclusion, the Treg/Th17 cell imbalance occurs throughout all stages of HBV-ACLF pathogenesis. The frequency of Th17 cells and the levels of their effectors (IL-17A, IL-23, and TNF-α) were upregulated with disease progression, resulting in a proinflammatory reaction. In parallel, a compensatory induction of Treg cells and their effector cytokines (TGF-β and IL-10) was also observed, yet an imbalance toward Th17 signals remained. During late-stage endotoxin-driven inflammation, the Treg/Th17 imbalance was further aggravated. Th17 cells, IL-17A, and IL-23 may represent valuable independent risk factors for the prognosis of HBV-ACLF. The current findings provide new insight into the pathogenesis of HBV-ACLF and its clinical prognosis.

Data Availability

The data used to support the findings of this study are available from the first author or corresponding author upon request.

Ethical Approval

The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All procedures performed in this study involving human participants were conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by the ethics committee of the First Affiliated Hospital of Hunan University of Chinese Medicine (No. PJ20170720).

Informed consent was provided by all participants.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Chen B, Tan NH, and Du S designed and performed research; Peng J provided serum samples; Tan NH and Peng J collected the data, performed the statistical analysis, and drafted the manuscript; Tan NH and Chen B revised the manuscript.

Acknowledgments

We are grateful to the experimental platform of the Key Laboratory of Liver Pathology of the First Affiliated Hospital of Hunan University of Chinese Medicine. We are grateful to the Hepatology Department of the First Affiliated Hospital of Hunan University of Chinese Medicine for providing the facility for specimen collection. We would like to thank Editage (http://www.editage.cn/) for English language editing. The research is funded by the National Science Foundation of China (No. 81673959), the Hunan Provincial Natural Science Foundation of China (No. 2020JJ4475), and the Scientific Research Project of Hunan Provincial Health Commission of China (No. 202103081452).