Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 12 January 2023
Sec. Pharmacoepidemiology
This article is part of the Research Topic Application of Data Mining in Pharmaceutical Research View all 12 articles

Association between glucocorticoid use and all-cause mortality in critically ill patients with heart failure: A cohort study based on the MIMIC-III database

Jia-Liang Zhu,&#x;Jia-Liang Zhu1,2Liang Hong&#x;Liang Hong3Shi-Qi YuanShi-Qi Yuan4Xiao-Mei Xu,Xiao-Mei Xu1,2Jian-Rui Wei
Jian-Rui Wei2*Hai-Yan Yin
Hai-Yan Yin1*
  • 1Department of Intensive Care Unit, The First Affiliated Hospital of Jinan University, Guangzhou, Guangdong, China
  • 2Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, Guangdong, China
  • 3Department of Intensive Care Unit, Nanjing First Hospital, Nanjing Medical University, Nanjing, China
  • 4Department of Neurology, The First Affiliated Hospital of Jinan University, Guangzhou, Guangdong, China

Background: Heart failure (HF) is the terminal stage of various heart diseases. Conventional treatments have poor efficacy, and diuretic resistance can present. Previous studies have found that the use of glucocorticoids can enhance the diuretic effect of patients with heart failure and reduce heart failure symptoms. However, the relationship between glucocorticoid use and mortality in patients with heart failure in intensive care units is unclear.

Objectives: The aim of this study was to determine the association between glucocorticoid use and all-cause mortality in critically ill patients with heart failure. Methods: The information on patients with heart failure in this study was extracted from the MIMIC-III (Medical Information Mart for Intensive Care-III) database. Patients in the glucocorticoid and non-glucocorticoid groups were matched using propensity scores. The Kaplan-Meier method was used to explore the difference in survival probability between the two groups. A Cox proportional-hazards regression model was used to analyze the hazard ratios (HRs) for the two patient groups. Subgroup analyses were performed with prespecified stratification variables to demonstrate the robustness of the results.

Results: The study included 9,482 patients: 2,099 in the glucocorticoid group and 7,383 in the non-glucocorticoid group. There were 2,055 patients in each group after propensity-score matching. The results indicated that the non-glucocorticoid group was not significantly associated with reduced mortality in patients with heart failure during the 14-day follow-up period [HRs = .901, 95% confidence interval (CI) = .767–1.059]. During the follow-up periods of 15–30 and 15–90 days, the mortality risk was significantly lower in the non-glucocorticoid group than in the glucocorticoid group (HRs = .497 and 95% CI = .370–.668, and HRs = .400 and 95% CI = .310–.517, respectively). Subgroup analyses indicated no interaction among each stratification variable and glucocorticoid use.

Conclusion: Glucocorticoid use was associated with an increased mortality risk in critically ill patients with heart failure.

Introduction

Heart failure (HF) is caused by abnormalities in heart structure and function, is characterized by water and sodium retention, often manifests as congestive pulmonary edema or vena cava congestion, and has typical symptoms of dyspnea, edema, and fatigue (Ponikowski et al., 20162016). It has become one of the most common causes of death worldwide, and has gradually become a serious public health problem due to the aging population and its increasing prevalence in the elderly (Vergaro et al., 2019). However, in addition to the aggravation of physiological symptoms, severe HF can also be accompanied by oliguria, electrolyte imbalance, and renal function deterioration. Conventional treatments often fail to significantly relieve symptoms, which increases the difficulty of treatment and the mortality risk.

The relationship between glucocorticoids and the cardiovascular system is well known to be complex. Studies have found that in the treatment of patients with acute decompensated heart failure (ADHF) and diuretic resistance, glucocorticoid use can enhance the effect of diuresis, increase urine output, and improve renal function (Massari et al., 2012). One study assessed a case of a diuretic-resistant patient with ADHF treated using methylprednisolone as a conventional treatment. After 3 days of glucocorticoid treatment the HF symptoms were significantly relieved, and the level of brain natriuretic peptide had decreased by 46% (Massari et al., 2012). A study by Bayliss found that glucocorticoids can promote urinary sodium excretion in patients with congestive HF and induce serum sodium to return to normal levels (Bayliss, 1966). The mechanism of the aforementioned beneficial effects may be related to a 30%–50% increase in renal blood flow and improved renal function (Raisz et al., 1957; Pechet et al., 1959). Studies of HF by Liu et al. found that patients in the glucocorticoid group had an increased urine output and glomerular filtration rate (GFR) compared with those in the placebo group, and these effects were observed after 3 days of glucocorticoid use (Liu et al., 2006a). There are a few studies that have suggested that glucocorticoids can be used safely in the short term and can lead to improvements in the clinical treatment of patients with ADHF (Liu and Liu, 2014). The results of a large, randomized, controlled clinical trial of methylprednisolone for treating acute myocardial infarction indicated that glucocorticoid use reduced the mortality risk in patients with new myocardial infarction by 50% (Stubbs, 1986). However, using glucocorticoids for HF also has certain adverse effects, such as increasing the risk of nosocomial infection, which makes this intervention controversial (Dec, 2007; Mosterd and Hoes, 2007). Glucocorticoids are the drugs most frequently associated with adverse drug events in hospitalized patients in the United States (AHRQ, 2006). Metabolic and cardiovascular adverse events are among the most common and serious adverse events caused by glucocorticoids (Fardet and Fève, 2014).

Research on the association between glucocorticoid use and all-cause mortality in patients with HF in intensive care units has produced unclear results. The aim of this study was to determine the associations of glucocorticoid use with 30- and 90-day all-cause mortality risk in patients with HF.

Materials and methods

Data collection

The data analyzed in this study were extracted from the Medical Information Mart for Intensive Care-III (MIMIC-III) database. This large database was established in 2003 by the Beth Israel Deaconess Medical Center and Massachusetts Institute of Technology. The latest version (version 1.4) of the MIMIC-III database collects admission data on more than 40,000 patients during 2001–2012. The database contains detailed patient information, including basic patient demographics, laboratory test results, imaging data, comorbidities, and other medical information (Yang et al., 2020). The database also includes times of admission, discharge, and death. Information on patients who died after discharge can be obtained from the social security database. We completed the relevant courses and passed the examination to obtain a certificate to access the database (certificate number: 45848365).

Patients are fully deidentified in the MIMIC-III database, and so this study did not require approval from the ethics committees of the hospitals. We used Structured Query Language to extract information on variables from the MIMIC-III database. Variables for which the proportion of missing data exceeded 10% were excluded, and missing values for the remaining variables were calculated using the multiple imputation method. The variables in this study included basic patient demographics, laboratory test results, comorbidities, and vital-sign variables (Wu et al., 2021). Demographic information included age and sex. Laboratory test results included potassium, sodium, hemoglobin, lymphocyte, platelet count, activated partial thromboplastin time (APTT), white blood cell (WBC) count, oxygen saturation (SpO2), and glucose. Vital signs included heart rate (HR), systolic blood pressure (SBP), diastolic blood pressure (DBP), respiratory rate (Rr), and body temperature (T). Comorbidities included congestive HF, hypertension, chronic pulmonary disease, diabetes, renal failure, liver disease, peptic ulcer, obesity, and anemia.

Inclusion and exclusion criteria

Patients diagnosed with HF according to the disease codes of the ninth revision of the International Classification of Diseases were included in the study. Exclusion criteria were as follows: 1) age < 18 years, 2) age > 90 years, and 3) not the first hospitalization.

Patient outcomes

The outcomes of this study were 30- and 90-day all-cause mortality in patients with HF after admission.

Statistical analysis

The study population was divided into glucocorticoid and non-glucocorticoid groups. Normally distributed continuous variables are expressed as mean ± standard-deviation values, and non-normally distributed continuous variables are expressed as quartiles. Student’s t-test and Kruskal–Wallis test were used to assess the significance of differences between the glucocorticoid and non-glucocorticoid groups. Categorical variables are expressed as frequencies and percentages, and the chi-square test was used to compare differences between groups. Propensity-score matching (PSM) was used to address the imbalance of baseline characteristics between the glucocorticoid and non-glucocorticoid groups. Patients were matched at a 1:1 ratio using estimated propensity scores with a .05 caliper width. The standardized mean difference (SMD) was used to assess whether baseline characteristics were balanced, with SMD < .1 indicating that they were. Multicollinearity among covariates was tested using the variance inflation factor (VIF), with VIF < 5 indicating that there was no multicollinearity.

We counted the number of outcome events at different follow-up periods and compared the differences between the two groups. The Kaplan-Meier (K-M) method was used to estimate the 30- and 90-day survival probabilities of patients in the glucocorticoid and non-glucocorticoid groups, and the log-rank test was used to compare differences in survival probabilities between the two groups. According to the intersection of the survival curves between the glucocorticoid and non-glucocorticoid groups, combined with our clinical experience, we divided the 30-day follow-up period into 1–14 and 15–30 days, and the 90-day period into 1–14 and 15–90 days. We performed landmark analyses for endpoints at different follow-up periods, and differences between the glucocorticoid and non-glucocorticoid groups were compared using the log-rank test.

Cox proportional-hazards regression was used to determine the differences in the risk of death between patients in the glucocorticoid and non-glucocorticoid groups, and results were expressed using hazard ratios (HRs) and 95% confidence intervals (95% CIs). We selected the population of the glucocorticoid group as a reference and established three Cox proportional-hazards regression models: 1) Model I, unadjusted model; 2) Model II, included age, sex, HR, SBP, DBP, Rr, T, and SpO2; and 3) Model III, Model II plus potassium, sodium, hemoglobin, lymphocyte, platelet count, APTT, WBC count, glucose, congestive HF, hypertension, chronic pulmonary disease, diabetes, renal failure, liver disease, peptic ulcer, obesity, and anemia. The performance of each model was evaluated by calculating the C statistic.

We also identified prespecified age, sex, congestive HF, diabetes, hypertension, chronic lung disease, renal failure, liver disease, peptic ulcer, obesity, and anemia. Patients were divided into two groups based on age, with the cutoff of 65 years. In each subgroup, Cox proportional-hazards regression analysis was performed after establishing Model III, and the results were presented using forest plots. We also assessed whether each variable interacted with glucocorticoid exposure to impact all-cause mortality in each follow-up period.

A probability value of p-value < .05 was considered significant. We used R software (version 4.2.0) for statistical analyses. The R packages used included gtsummary, tableone, dplyr, remotes, jskm, survival, MatchIt, foreign, forestplot, tidyverse, ISwR, car, lattice, MASS, nnet, mice, and cobalt.

Results

As shown in Figure 1, 10,372 patients with HF were extracted from the MIMIC-III database. After excluding those aged < 18 and > 90 years and who had previously been hospitalized, 9,482 patients were included in the final analysis. Before PSM, 2,099 patients who had used glucocorticoids were assigned to the glucocorticoid group and 7,383 patients who had not used glucocorticoids were assigned to the non-glucocorticoid group. Table 1 lists the baseline characteristics of the glucocorticoid and non-glucocorticoid groups before matching. Supplementary Table S1 lists the VIFs for each covariate in the three follow-up periods, all of which were < 5, indicating no multicollinearity among the variables. The mean ages of the glucocorticoid and non-glucocorticoid groups were 71.28 and 73.01 years, respectively. The numbers of males and females were 979 and 1,120 in the glucocorticoid group, respectively, and 4,035 and 3,348 in the non-glucocorticoid group. The numbers of patients with congestive HF in the glucocorticoid and non-glucocorticoid groups were 2,048 (97.6%) and 7,066 (95.7%), respectively. The WBC and lymphocyte counts were significantly higher in the glucocorticoid group than in the non-glucocorticoid group. Overall, the baseline characteristics of the two groups of patients were unbalanced for most variables.

FIGURE 1
www.frontiersin.org

FIGURE 1. Selection of study population from MIMIC-Ⅲ database.

TABLE 1
www.frontiersin.org

TABLE 1. Baseline characteristics of population before PSM.

In PSM, 2,055 patients with glucocorticoids were matched with the same number of patients in the non-glucocorticoid group. The baseline characteristics of patients in the glucocorticoid and non-glucocorticoid groups after PSM are listed in Supplementary Table S2. Figure 2 presents the SMD of each covariate before and after PSM. All covariates had SMD < .1 after matching, with p-value > .05, indicating that the covariates in the two groups had been balanced.

FIGURE 2
www.frontiersin.org

FIGURE 2. Standardized mean difference (SMD) of variables before and after propensity score matching.

We assessed differences in all-cause mortality between the glucocorticoid and non-glucocorticoid groups during different follow-up periods. Before PSM, in all follow-up periods (14, 15–30, and 15–90 days), the mortality rate was significantly higher in the glucocorticoid group than in the non-glucocorticoid group (p-value < .05). Similarly, the mortality rate was significantly higher in the glucocorticoid group than in the non-glucocorticoid group after PSM; however, there was no significant difference between the two groups during the 14-day follow-up period after PSM (p-value > .05). The results are listed in Table 2.

TABLE 2
www.frontiersin.org

TABLE 2. Mortality rates in patients with heart failure in the glucocorticoids group and non-glucocorticoids groups.

K-M curves were used to plot the 30- and 90-day survival probabilities of the glucocorticoid and non-glucocorticoid groups after PSM, which are shown in Figure 3, respectively. Figure 4 present the results of the landmark analysis for 30- and 90-day survival probabilities, respectively. The 30- and 90-day survival probabilities differed significantly between the two patient groups (p-value < .05). The survival curves of the glucocorticoid and non-glucocorticoid groups intersected at 14 days. After 14 days, the mortality risk was slightly higher in the non-glucocorticoid group than in the glucocorticoid group, but the difference was not significant (p-value > .05). After 15–30 and 15–90 days of follow-up, the mortality risk was significantly lower in the non-glucocorticoid group than in the glucocorticoid group (p-value < .05).

FIGURE 3
www.frontiersin.org

FIGURE 3. Kaplan-Meier survival curves of non-glucocorticoids group and non-glucocorticoids group after PSM. (A) 30-day mortality; (B) 90-day mortality.

FIGURE 4
www.frontiersin.org

FIGURE 4. Landmark analysis discriminating between survival probability before and after 14 day of follow-up. (A) Landmark analysis during 30-day follow-up; (B) Landmark analysis during 90-day follow-up.

In the Cox proportional-hazards regression analysis, we calculated the C statistic for Model III. Before and after PSM, the C statistics for Model III were .739 and .746 for 14 days, respectively, .769 and .760 for 15–30 days, and .775 and .758 for 15–90 day, indicating that the adopted model had good accuracy; the results are listed in Supplementary Table S3. After PSM, within the 14-day follow-up period and with the glucocorticoid group as a reference, the HRs (95% CI) values of the non-glucocorticoid group were .972 (.830–1.139) for Model I, .921 (.785–1.081) for Model II, and .901 (.767–1.059) for Model III. The difference between the two groups for Model III was not significant (p-value > .05). During the 15–30 day follow-up period, the risk of death was significantly lower in the non-glucocorticoid group than in the glucocorticoid group (p-value < .05); the HRs (95% CI) values of the non-glucocorticoid group were .496 (.370–.665) for Model I, .502 (.374–.672) for Model II, and .497 (.370–.668) for Model III. Similarly, during the 15–90 day follow-up period, the HRs (95% CI) values of the non-glucocorticoid groups were .398 (.309–.513) for Model I, .401 (.311–.517) for Model II, and .400 (.310–.517) for Model III, respectively, and the differences between the risks of death were significant (p-value < .05). The results of the Cox proportional-hazards regression analysis are presented in Table 3.

TABLE 3
www.frontiersin.org

TABLE 3. Association between the non-glucocorticoid group and all-cause mortality.

The subgroup analysis results indicated that stratification variables and glucocorticoid exposure did not interact (P for interaction > .05). Supplementary Figures S1–S3 show the subgroup analysis results. Within 14 days, there were no significant differences in the risk of death between the glucocorticoid and non-glucocorticoid groups for each stratification variable (p-value > .05), which was consistent with the overall risk of death, indicating that the results were robust. In the 15–30 and 15–90 day follow-up periods, the risk of death differed significantly between the glucocorticoid and non-glucocorticoid groups (p-value < .05) for all stratification variables except for obesity and anemia, which was consistent with the overall risk of death, indicating that the results were robust.

Discussion

HF is a complex clinical syndrome caused by abnormal changes in cardiac structure and/or function via various mechanisms, which results in ventricular systolic and/or diastolic dysfunction. It is a severe manifestation or terminal stage of various heart diseases. In the early stage of HF, the stress response significantly increases the levels of epinephrine, glucagon, and glucocorticoids to above the physiological ranges (Luger et al., 1987). In the later stage, due to low cardiac output, the adrenal cortex is ischemic and hypoxic due to hypoperfusion, so its function of synthesizing corticosteroids is impaired, resulting in decreased blood cortisol.

Glucocorticoids are steroid compounds secreted by the zona fasciae in the adrenal cortex and is mostly regulated by the hypothalamic-pituitary-adrenal axis (Liu et al., 2019a). Physiological doses of glucocorticoids can regulate substance metabolism in the body and maintain its life activities. Supraphysiological doses of glucocorticoids have anti-inflammatory, anti-infective, antishock, and immunosuppressive effects, so they are widely used in clinical practice. The effects of standard treatment drugs are often less strong in the late stage of HF than in the early stage. The adrenal glands are often in an excited state due to repeated acute HF attacks, which eventually leads to a decrease in glucocorticoid secretion, especially in elderly patients with HF accompanied by adrenal insufficiency and infection.

Various animal experiments have found that glucocorticoids can enhance myocardial contractility, which occurs via numerous mechanisms. First, the increased release or enhanced effect of endogenous catecholamines indirectly leads to increased myocardial contractility (Tecklenberg et al., 1973; Jain et al., 2004). Second, it may also be caused by β-adrenergic receptor upregulation in the myocardium (Nishimura et al., 1997). Glucocorticoids can also protect kidney function. Previous studies have found that glucocorticoids dilate renal blood vessels, enhance the sensitivity of nephrons to diuretics, and increase the GFR (Liu et al., 2007; Smets et al., 2012; Liu et al., 2016), via the following mechanisms: First, glucocorticoids upregulate natriuretic peptide receptor A expression in the inner medullary collecting duct cells to enhance the sensitivity of the kidneys to natriuretic peptides and promote the diuretic effect in patients with HF (Gardner et al., 1986; Damjancic and Vierhapper, 1990a; Liu et al., 2010; Liu et al., 2011). Second, glucocorticoids reduce the production and secretion of vasopressin, and downregulate arginine vasopressin receptor expression (Erkut et al., 1998; Liu et al., 2006b; Greenwood et al., 2015; Zhu et al., 2020); the production of prostaglandins, nitric oxide, and dopamine is also increased, resulting in improved renal blood flow and GFR (Gong et al., 2008; Tokudome et al., 2009; Butts and Phillips, 2013; Liu et al., 2019b). Third, an animal study found that glucocorticoids somewhat affected the synthesis and release of atrial natriuretic peptide (ANP) (Garcia et al., 1985), up-regulating ANP receptors in vascular endothelial cells (Lanier-Smith and Currie, 1990). Likewise, studies have found that glucocorticoids may modulate ANP-mediated natriuresis and diuresis in humans (Damjancic and Vierhapper, 1990b). Glucocorticoids also enhance pulmonary edema clearance. Glucocorticoids can adjust the activity of alveolar epithelial Na+ channel and Na+-K+ ATPase, increase the tension in pulmonary blood vessels, and reduce the permeability of capillaries, so it can increase alveolar fluid clearance and accelerate the curing of pulmonary edema (Barquin et al., 1997; Ingbar et al., 1997).

However, excessive cortisol secretion and the use of various synthetic glucocorticoids increase the incidence rates of diabetes and cardiovascular disease (Plotz et al., 1952). Glucocorticoid use in non-diabetic patients increases their risk of developing diabetes by two-to fourfold (Conn and Poynard, 1994; Gurwitz et al., 1994; Blackburn et al., 2002; Gulliford et al., 2006). The mechanism of glucocorticoid-induced diabetes involves increasing insulin resistance, thereby affecting the glucose metabolism process, which is similar to the mechanism of type 2 diabetes. In the liver, increased insulin resistance can increase basal glucose production. Insulin also promotes intracellular glucose utilization by stimulating glucose transporter type 4, but glucocorticoids interfere with these signaling pathways and with glycogen synthesis. The above two mechanisms may lead to abnormal glucose metabolism and increase blood glucose levels (Ruzzin et al., 2005; Burén et al., 2008). An observational study found that low-density lipoprotein (LDL) was elevated after glucocorticoid therapy in patients with asthma, heart transplant, kidney transplant, and rheumatoid arthritis, which led to abnormal lipid metabolism (Sholter and Armstrong, 2000). An in vitro study found that hydrocortisone and dexamethasone alter LDL degradation in fibroblasts and macrophages, which may contribute to atherosclerosis (Henze et al., 1983; Hirsch and Mazzone, 1986).

Dexamethasone is also known to cause cardiotoxicity (de Salvi Guimarães et al., 2017a). Studies have found that long-term dexamethasone use may lead to cardiac fibrosis, increased norepinephrine-induced vasoconstriction, increased apoptosis, and decreased angiogenesis (Qi and Rodrigues, 2007; de Salvi Guimarães et al., 2017a). Dexamethasone can also cause cphy, which is an adaptive ardiac hypertrochange in response to dexamethasone-induced hypertension (Walker, 2007; Stahn and Buttgereit, 2008; Bal et al., 2009). Bézard et al. (2021) assessed the effect of dexamethasone on cardiac function using the Langendorff perfused heart model in healthy rats. The model is free from other confounding factors, allowing for precise studies of how the heart responds to dexamethasone (Sutherland and Hearse, 2000). Their study found that dexamethasone significantly reduced overall left ventricle (LV) function, as well as being caused LV systolic and diastolic dysfunction to various degrees. It is particularly interesting that the above-mentioned changes were also observed in healthy rats (Bézard et al., 2021). Similarly, de Salvi Guimarães et al. (2017b) found that dexamethasone-treated rats also developed multiple circulatory complications, including increased blood pressure, myocardial fibrosis, and cardiomyocyte apoptosis, which led to myocardial remodeling and diastolic dysfunction. That study also confirmed that the mechanisms underlying these changes included impaired calcium handling and activation of the calcineurin signaling pathway. It has also been found that systemic glucocorticoid use may increase the risk of deep-vein thrombosis, especially in patients with pulmonary embolism (Johannesdottir et al., 2013). Despite the effects of glucocorticoids on sodium excretion, diuresis, and increased GFR, glucocorticoid use increased the likelihood of non-acute complications such as diabetes, dyslipidemia, cardiotoxicity, and deep-vein thrombosis, which will be detrimental to the prognosis of patients with HF.

In this study, the WBC count was higher in the glucocorticoid group than in the non-glucocorticoid group (13.06 ± 20.53 vs. 11.76 ± 9.66 k/μl), as was the proportion of patients with chronic respiratory disease (45.7% vs. 28.8%). This suggests that the higher risk of death in the glucocorticoid group is caused by glucocorticoid-use-related diseases, such as chronic respiratory disease and infection. Some studies also found that serum inflammatory markers were associated with the incidence of cardiovascular events, and disease activity was directly related to the occurrence of cardiovascular events (Ross, 1999). We found that the association between glucocorticoid use and increased mortality in patients with HF persisted after adjusting for chronic respiratory disease and WBC counts.

We performed subgroup analyses of age, sex, congestive HF, diabetes, hypertension, chronic lung disease, kidney failure, liver disease, peptic ulcer, obesity, and anemia. The results indicated that there was no interaction between stratification variables and glucocorticoid exposure. The results of each subgroup analysis and the overall results also remained consistent at different follow-up times, indicating that our results were robust.

The highlight of our study was the division of follow-up into two time periods to more precisely account for the occurrence of outcome events. This study had some limitations. First, database limitations meant that we did not perform a subgroup analysis of total glucocorticoid dosage and administration duration. Second, the left ventricle ejection fraction was not analyzed in the study because it had more than 40% missing values. In addition, our study is a single-center study, and we hope that in future studies, multi-center studies can be conducted to obtain more information.

In conclusion, this retrospective study of a large database found that glucocorticoid use was associated with increased all-cause mortality in critically ill patients with HF. Our findings suggest that the use of glucocorticoids should not be recommended for critically ill patients with HF, but this requires further validation in prospective studies or randomized controlled trials.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://mimic.mit.edu/.

Author contributions

J-LZ and LH contributed equally to the work. J-RW and H-YY contributed to the study concept and study design and contributed equally. J-LZ performed statistical analysis and data interpretation. J-LZ were responsible for the quality control of data and algorithms. S-QY, LH, and X-MX performed literature research and data extraction. All authors contributed to writing of the manuscript and approved the final version.

Funding

This work was supported by National Natural Science Foundation of China (Nos 82072232; 81871585); Science and Technology Program of Guangzhou, China (No. 202201020028); Special Projects in Key Areas of General Colleges and Universities in Guangdong Province (No. 2022ZDZX 2003); The 2021 Annual Medical Teaching and Education Management Reform Research Project of Jinan University (No. 2021YXJG029).

Acknowledgments

The authors wish to thank the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center for the MIMIC project.

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.2023.1118551/full#supplementary-material

References

Ahrq, , Healthcare cost and utilization project (HCUP) statistical briefs [Internet]. Rockville (MD): Agency for Healthcare Research and Quality US; 2006. PMID: 21413206.

Google Scholar

Bal, M. P., de Vries, W. B., Steendijk, P., Homoet-van der Kraak, P., van der Leij, F. R., Baan, J., et al. (2009). Histopathological changes of the heart after neonatal dexamethasone treatment: Studies in 4-8-and 50-week-old rats. Pediatr. Res. 66 (1), 74–79. doi:10.1203/PDR.0b013e3181a283a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Barquin, N., Ciccolella, D. E., Ridge, K. M., and Sznajder, J. I. (1997). Dexamethasone upregulates the Na-K-ATPase in rat alveolar epithelial cells. Am. J. Physiol. 273 (4), L825–L830. doi:10.1152/ajplung.1997.273.4.L825

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayliss, R. I. (1966). Corticosteroids in heart disease. Br. Med. J. 2 (5516), 721–724. PMID: 5917388; PMCID: PMC1943945. doi:10.1136/bmj.2.5516.721

PubMed Abstract | CrossRef Full Text | Google Scholar

Bézard, M., Oghina, S., Vitiello, D., Kharoubi, M., Kordeli, E., Galat, A., et al. (2021). Dexamethasone is associated with early deaths in light chain amyloidosis patients with severe cardiac involvement. PLoS One 16 (9), e0257189. PMID: 34525116; PMCID: PMC8443042. doi:10.1371/journal.pone.0257189

PubMed Abstract | CrossRef Full Text | Google Scholar

Blackburn, D., Hux, J., and Mamdani, M. (2002). Quantification of the risk of corticosteroid-induced diabetes mellitus among the elderly. J. Gen. Intern Med. 17 (9), 717–720. PMID: 12220369; PMCID: PMC1495107. doi:10.1046/j.1525-1497.2002.10649.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Burén, J., Lai, Y. C., Lundgren, M., Eriksson, J. W., and Jensen, J. (2008). Insulin action and signalling in fat and muscle from dexamethasone-treated rats. Arch. Biochem. Biophys. 474 (1), 91–101. Epub 2008 Feb 29. PMID: 18328801. doi:10.1016/j.abb.2008.02.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Butts, K. A., and Phillips, A. G. (2013). Glucocorticoid receptors in the prefrontal cortex regulate dopamine efflux to stress via descending glutamatergic feedback to the ventral tegmental area. Int. J. Neuropsychopharmacol. 16 (8), 1799–1807. PMID: 23590841. doi:10.1017/S1461145713000187

PubMed Abstract | CrossRef Full Text | Google Scholar

Conn, H. O., and Poynard, T. (1994). Corticosteroids and peptic ulcer: meta-analysis of adverse events during steroid therapy. J. Intern Med. 236 (6), 619–632. doi:10.1111/j.1365-2796.1994.tb00855.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Damjancic, P., and Vierhapper, H. (1990). Permissive action of glucocorticoid substitution therapy on the effects of atrial natriuretic peptide (hANP) in patients with adrenocortical insufficiency. Exp. Clin. Endocrinol. 95 (3), 315–321. doi:10.1055/s-0029-1210971

PubMed Abstract | CrossRef Full Text | Google Scholar

Damjancic, P., and Vierhapper, H. (1990). Permissive action of glucocorticoid substitution therapy on the effects of atrial natriuretic peptide (hANP) in patients with adrenocortical insufficiency. Exp. Clin. Endocrinol. 95 (3), 315–321. doi:10.1055/s-0029-1210971

PubMed Abstract | CrossRef Full Text | Google Scholar

de Salvi Guimarães, F., de Moraes, W. M., Bozi, L. H., Souza, P. R., Antonio, E. L., Bocalini, D. S., et al. (2017). Dexamethasone-induced cardiac deterioration is associated with both calcium handling abnormalities and calcineurin signaling pathway activation. Mol. Cell Biochem. 424 (1-2), 87–98. Epub 2016 Oct 19. PMID: 27761848. doi:10.1007/s11010-016-2846-3

PubMed Abstract | CrossRef Full Text | Google Scholar

de Salvi Guimarães, F., de Moraes, W. M., Bozi, L. H., Souza, P. R., Antonio, E. L., Bocalini, D. S., et al. (2017). Dexamethasone-induced cardiac deterioration is associated with both calcium handling abnormalities and calcineurin signaling pathway activation. Mol. Cell Biochem. 424 (1-2), 87–98. Epub 2016 Oct 19. PMID: 27761848. doi:10.1007/s11010-016-2846-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Dec, G. W. (2007). Management of acute decompensated heart failure. Curr. Probl. Cardiol. 32 (6), 321–366. PMID: 17531903. doi:10.1016/j.cpcardiol.2007.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Erkut, Z. A., Pool, C., and Swaab, D. F. (1998). Glucocorticoids suppress corticotropin-releasing hormone and vasopressin expression in human hypothalamic neurons. J. Clin. Endocrinol. Metab. 83 (6), 2066–2073. doi:10.1210/jcem.83.6.4881

PubMed Abstract | CrossRef Full Text | Google Scholar

Fardet, L., and Fève, B. (2014). Systemic glucocorticoid therapy: A review of its metabolic and cardiovascular adverse events. Drugs 74 (15), 1731–1745. PMID: 25204470. doi:10.1007/s40265-014-0282-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, R., Debinski, W., Gutkowska, J., Kuchel, O., Thibault, G., Genest, J., et al. (1985). Gluco- and mineralocorticoids may regulate the natriuretic effect and the synthesis and release of atrial natriuretic factor by the rat atria in vivo. Biochem. Biophys. Res. Commun. 131 (2), 806–814. PMID: 2932106. doi:10.1016/0006-291x(85)91311-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardner, D. G., Hane, S., Trachewsky, D., Schenk, D., and Baxter, J. D. (1986). Atrial natriuretic peptide mRNA is regulated by glucocorticoids in vivo. Biochem. Biophys. Res. Commun. 139 (3), 1047–1054. PMID: 2945558. doi:10.1016/s0006-291x(86)80283-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, R., Latif, S., Morris, D. J., and Brem, A. S. (2008). Co-localization of glucocorticoid metabolizing and prostaglandin synthesizing enzymes in rat kidney and liver. Life Sci. 83 (21-22), 725–731. Epub 2008 Sep 30. PMID: 18930067. doi:10.1016/j.lfs.2008.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Greenwood, M., Greenwood, M. P., Mecawi, A. S., Loh, S. Y., Rodrigues, J. A., Paton, J. F., et al. (2015). Transcription factor CREB3L1 mediates cAMP and glucocorticoid regulation of arginine vasopressin gene transcription in the rat hypothalamus. Mol. Brain 8 (1), 68. PMID: 26503226; PMCID: PMC4624382. doi:10.1186/s13041-015-0159-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulliford, M. C., Charlton, J., and Latinovic, R. (2006). Risk of diabetes associated with prescribed glucocorticoids in a large population. Diabetes Care 29 (12), 2728–2729. doi:10.2337/dc06-1499

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurwitz, J. H., Bohn, R. L., Glynn, R. J., Monane, M., Mogun, H., and Avorn, J. (1994). Glucocorticoids and the risk for initiation of hypoglycemic therapy. Arch. Intern Med. 154 (1), 97–101. PMID: 8267494. doi:10.1001/archinte.154.1.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Henze, K., Chait, A., Albers, J. J., and Bierman, E. L. (1983). Hydrocortisone decreases the internalization of low density lipoprotein in cultured human fibroblasts and arterial smooth muscle cells. Eur. J. Clin. Invest. 13 (2), 171–177. doi:10.1111/j.1365-2362.1983.tb00083.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, L. J., and Mazzone, T. (1986). Dexamethasone modulates lipoprotein metabolism in cultured human monocyte-derived macrophages. Stimulation of scavenger receptor activity. J. Clin. Invest. 77 (2), 485–490. PMID: 3944266; PMCID: PMC423369. doi:10.1172/JCI112327

PubMed Abstract | CrossRef Full Text | Google Scholar

Ingbar, D. H., Duvick, S., Savick, S. K., Schellhase, D. E., Detterding, R., Jamieson, J. D., et al. (1997). Developmental changes of fetal rat lung Na-K-ATPase after maternal treatment with dexamethasone. Am. J. Physiol. 272 (4), L665–L672. doi:10.1152/ajplung.1997.272.4.L665

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, P., Armando, I., Juorio, A. V., Barden, N., Benicky, J., and Saavedra, J. M. (2004). Decreased hypothalamic and adrenal angiotensin II receptor expression and adrenomedullary catecholamines in transgenic mice with impaired glucocorticoid receptor function. Neuroendocrinology 80 (3), 171–180. Epub 2004 Dec 2. PMID: 15583474. doi:10.1159/000082358

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesdottir, S. A., Horváth-Puhó, E., Dekkers, O. M., Cannegieter, S. C., Jørgensen, J. O., Ehrenstein, V., et al. (2013). Use of glucocorticoids and risk of venous thromboembolism: A nationwide population-based case-control study. JAMA Intern Med. 173 (9), 743–752. doi:10.1001/jamainternmed.2013.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Lanier-Smith, K. L., and Currie, M. G. (1990). Effect of glucocorticoids on the binding of atrial natriuretic peptide to endothelial cells. Eur. J. Pharmacol. 178 (1), 105–109. doi:10.1016/0014-2999(90)94800-d

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., Zhang, T. N., Knight, J. K., and Goodwin, J. E. (2019). The glucocorticoid receptor in cardiovascular health and disease. Cells 8 (10), 1227. PMID: 31601045; PMCID: PMC6829609. doi:10.3390/cells8101227

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., Zhang, T. N., Knight, J. K., and Goodwin, J. E. (2019). The glucocorticoid receptor in cardiovascular health and disease. Cells 8 (10), 1227. PMID: 31601045; PMCID: PMC6829609. doi:10.3390/cells8101227

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Chen, H., Zhou, C., Ji, Z., Liu, G., Gao, Y., et al. (2006). Potent potentiating diuretic effects of prednisone in congestive heart failure. J. Cardiovasc Pharmacol. 48 (4), 173–176. doi:10.1097/01.fjc.0000245242.57088.5b

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Chen, Y., Kang, Y., Ni, Z., Xiu, H., Guan, J., et al. (2011). Glucocorticoids improve renal responsiveness to atrial natriuretic peptide by up-regulating natriuretic peptide receptor-A expression in the renal inner medullary collecting duct in decompensated heart failure. J. Pharmacol. Exp. Ther. 339 (1), 203–209. Epub 2011 Jul 7. PMID: 21737535. doi:10.1124/jpet.111.184796

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Ge, N., Zhai, J. L., and Zhang, J. X. (2016). Dexamethasone-induced diuresis is associated with inhibition of the renin-angiotensin-aldosterone system in rats. Kaohsiung J. Med. Sci. 32 (12), 614–619. Epub 2016 Oct 27. PMID: 27914612. doi:10.1016/j.kjms.2016.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Guan, J., Kang, Y., Xiu, H., Chen, Y., Deng, B., et al. (2010). Inhibition of dehydration-induced water intake by glucocorticoids is associated with activation of hypothalamic natriuretic peptide receptor-A in rat. PLoS One 5 (12), e15607. PMID: 21187974; PMCID: PMC3004933. doi:10.1371/journal.pone.0015607

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Liu, G., Zhou, C., Ji, Z., Zhen, Y., and Liu, K. (2007). Potent diuretic effects of prednisone in heart failure patients with refractory diuretic resistance. Can. J. Cardiol. 23 (11), 865–868. PMID: 17876376; PMCID: PMC2651362. doi:10.1016/s0828-282x(07)70840-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., and Liu, K. (2014). Cardiac outcome prevention effectiveness of glucocorticoids in acute decompensated heart failure: COPE-ADHF study. J Cardiovasc Pharmacol. 63 (4), 333–338. PMID: 24710470. doi:10.1097/FJC.0000000000000048

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, R. Y., Unmehopa, U. A., Zhou, J. N., and Swaab, D. F. (2006). Glucocorticoids suppress vasopressin gene expression in human suprachiasmatic nucleus. J Steroid Biochem Mol Biol 98 (4), 248–253. Epub 2006 Feb 17. PMID: 16483767. doi:10.1016/j.jsbmb.2005.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Luger, A., Deuster, P. A., Kyle, S. B., Gallucci, W. T., Montgomery, L. C., Gold, P. W., et al. (1987). Acute hypothalamic-pituitary-adrenal responses to the stress of treadmill exercise. Physiologic adaptations to physical training. N. Engl. J. Med. 316 (21), 1309–1315. doi:10.1056/NEJM198705213162105

PubMed Abstract | CrossRef Full Text | Google Scholar

Massari, F., Mastropasqua, F., Iacoviello, M., Nuzzolese, V., Torres, D., and Parrinello, G. (2012). The glucocorticoid in acute decompensated heart failure: Dr Jekyll or Mr Hyde? Am. J. Emerg. Med. 30 (3), 517.e5–10. doi:10.1016/j.ajem.2011.01.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Mosterd, A., and Hoes, A. W. (2007). Clinical epidemiology of heart failure. Heart 93 (9), 1137–1146. PMID: 17699180; PMCID: PMC1955040. doi:10.1136/hrt.2003.025270

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimura, H., Yoshikawa, T., Kobayashi, N., Anzai, T., Nagami, K., Handa, S., et al. (1997). Effects of methylprednisolone on hemodynamics and beta-adrenergic receptor signaling in rabbits with acute left ventricular failure. Heart Vessels 12 (2), 84–91. doi:10.1007/BF02820871

PubMed Abstract | CrossRef Full Text | Google Scholar

Pechet, M. M., Bowers, B., and Bartter, F. C. (1959). Metabolic studies with a new series of 1, 4-diene steroids. II. Effects in normal subjects of prednisone, prednisolone, and 9 alpha-fluoroprednisolone. J. Clin. Invest. 38 (4), 691–701. PMID: 13641421; PMCID: PMC293209. doi:10.1172/JCI103848

PubMed Abstract | CrossRef Full Text | Google Scholar

Plotz, C. M., Knowlton, A. I., and Ragan, C. (1952). The natural history of Cushing's syndrome. Am. J. Med. 13 (5), 597–614. PMID: 12996538. doi:10.1016/0002-9343(52)90027-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ponikowski, P., Voors, A. A., Anker, S. D., Bueno, H., Cleland, J. G. F., Coats, A. J. S., et al. (2016). 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC)Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur. Heart J. 37 (27), 2129–2200. 20. Erratum in: Eur Heart J. 2016 Dec 30;: PMID: 27206819. doi:10.1093/eurheartj/ehw128

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, D., and Rodrigues, B. (2007). Glucocorticoids produce whole body insulin resistance with changes in cardiac metabolism. Am. J. Physiol. Endocrinol. Metab. 292 (3), E654–E667. Epub 2006 Oct 31. PMID: 17077342. doi:10.1152/ajpendo.00453.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Raisz, L. G., Mcneely, W. F., Saxon, L., and Rosenbaum, J. D. (1957). The effects of cortisone and hydrocortisone on water diuresis and renal function in man. J. Clin. Invest. 36 (6), 767–779. PMID: 13439015; PMCID: PMC441751. doi:10.1172/JCI103481

PubMed Abstract | CrossRef Full Text | Google Scholar

Ross, R. (1999). Atherosclerosis-an inflammatory disease. N. Engl. J. Med. 340 (2), 115–126. doi:10.1056/NEJM199901143400207

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruzzin, J., Wagman, A. S., and Jensen, J. (2005). Glucocorticoid-induced insulin resistance in skeletal muscles: Defects in insulin signalling and the effects of a selective glycogen synthase kinase-3 inhibitor. Diabetologia 48 (10), 2119–2130. Epub 2005 Aug 3. PMID: 16078016. doi:10.1007/s00125-005-1886-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sholter, D. E., and Armstrong, P. W. (2000). Adverse effects of corticosteroids on the cardiovascular system. Can. J. Cardiol. 16 (4), 505–511. PMID: 10787466.

PubMed Abstract | Google Scholar

Smets, P. M., Lefebvre, H. P., Aresu, L., Croubels, S., Haers, H., Piron, K., et al. (2012). Renal function and morphology in aged Beagle dogs before and after hydrocortisone administration. PLoS One 7 (2), e31702. Epub 2012 Feb 29. PMID: 22393368; PMCID: PMC3290534. doi:10.1371/journal.pone.0031702

PubMed Abstract | CrossRef Full Text | Google Scholar

Stahn, C., and Buttgereit, F. (2008). Genomic and nongenomic effects of glucocorticoids. Nat. Clin. Pract. Rheumatol. 4 (10), 525–533. PMID: 18762788. doi:10.1038/ncprheum0898

PubMed Abstract | CrossRef Full Text | Google Scholar

Stubbs, D. F. (1986). Positive synergism between diuretics and methylprednisolone following acute myocardial infarction. J. Int. Med. Res. 14 (1), 21–24. doi:10.1177/03000605860140S104

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutherland, F. J., and Hearse, D. J. (2000). The isolated blood and perfusion fluid perfused heart. Pharmacol. Res. 41 (6), 613–627. doi:10.1006/phrs.1999.0653

PubMed Abstract | CrossRef Full Text | Google Scholar

Tecklenberg, P. L., Mullin, E. M., Stinson, E. B., and Morrow, A. G. (1973). The effects of massive doses of methylprednisolone on myocardial contractility and peripheral vascular resistance. Am. Heart J. 85 (2), 216–226. PMID: 4687122. doi:10.1016/0002-8703(73)90463-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokudome, S., Sano, M., Shinmura, K., Matsuhashi, T., Morizane, S., Moriyama, H., et al. (2009). Glucocorticoid protects rodent hearts from ischemia/reperfusion injury by activating lipocalin-type prostaglandin D synthase-derived PGD2 biosynthesis. J. Clin. Invest. 119 (6), 1477–1488. PMID: 19451694; PMCID: PMC2689117. doi:10.1172/JCI37413

PubMed Abstract | CrossRef Full Text | Google Scholar

Vergaro, G., Aimo, A., Prontera, C., Ghionzoli, N., Arzilli, C., Zyw, L., et al. (2019). Sympathetic and renin-angiotensin-aldosterone system activation in heart failure with preserved, mid-range and reduced ejection fraction. Int. J. Cardiol. 296, 91–97. Epub 2019 Aug 17. PMID: 31443984. doi:10.1016/j.ijcard.2019.08.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, B. R. (2007). Glucocorticoids and cardiovascular disease. Eur. J. Endocrinol. 157 (5), 545–559. PMID: 17984234. doi:10.1530/EJE-07-0455

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W. T., Li, Y. J., Feng, A. Z., Li, L., Huang, T., Xu, A. D., et al. (2021). Data mining in clinical big data: The frequently used databases, steps, and methodological models. Mil. Med. Res. 8 (1), 44. PMID: 34380547; PMCID: PMC8356424. doi:10.1186/s40779-021-00338-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Li, Y., Liu, Q., Li, L., Feng, A., Wang, T., et al. (2020). Brief introduction of medical database and data mining technology in big data era. J. Evid. Based Med. 13 (1), 57–69. Epub 2020 Feb 22. PMID: 32086994; PMCID: PMC7065247. doi:10.1111/jebm.12373

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Huang, Y., Li, S., Ge, N., Li, T., Wang, Y., et al. (2020). Glucocorticoids reverse diluted hyponatremia through inhibiting arginine vasopressin pathway in heart failure rats. J. Am. Heart Assoc. 9 (10), e014950. Epub 2020 May 9. PMID: 32390535; PMCID: PMC7660850. doi:10.1161/JAHA.119.014950

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: glucocorticoids, heart failure, all-cause mortality, retrospective cohort study, MIMIC-III database

Citation: Zhu J-L, Hong L, Yuan S-Q, Xu X-M, Wei J-R and Yin H-Y (2023) Association between glucocorticoid use and all-cause mortality in critically ill patients with heart failure: A cohort study based on the MIMIC-III database. Front. Pharmacol. 14:1118551. doi: 10.3389/fphar.2023.1118551

Received: 07 December 2022; Accepted: 02 January 2023;
Published: 12 January 2023.

Edited by:

Zhenyu Pan, Xi’an Children’s Hospital, China

Reviewed by:

Hua Jiang, Sichuan Academy of Medical Sciences and Sichuan Provincial People’s Hospital, China
Yuli Huang Huang, Southern Medical University, China

Copyright © 2023 Zhu, Hong, Yuan, Xu, Wei and Yin. 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: Jian-Rui Wei, Jianruiw2021@163.com; Hai-Yan Yin, yinhaiyan1867@126.com

These authors have contributed equally to this work

Disclaimer: 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.