- 1Key Laboratory of Molecular Biology for Infectious Diseases, Ministry of Education, Institute for Viral Hepatitis, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of Respiratory and Critical Care Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 3Department of Intensive Care Unit, The Second Affiliated Hospital, Chongqing Medical University, Chongqing, China
- 4Department of Hepatobiliary Surgery, The Second Affiliated Hospital, Chongqing Medical University, Chongqing, China
Background: Sepsis is a life-threatening disease with challenges in clinical management due to delayed diagnosis and immunosuppression. Lymphopenia is a key prognostic indicator in sepsis. microRNAs (miRNAs) are recognized as key immune modulators affecting all stages of inflammation. In this study, we investigated the expression of miR-146b-5p and lymphocyte subsets in sepsis patients and analyzed their correlation with the aim of establishing their combined predictive value for clinical outcomes and advancing personalized treatment strategies.
Methods: From January 2020 to July 2024, we enrolled 191 patients diagnosed with sepsis at the ICU of the Second Affiliated Hospital of Chongqing Medical University and collected basic clinical data. These patients were categorized into two groups: nonsurvivors (n = 117) and survivors (n = 74). Correlation analysis was employed to analyze the correlation between miR-146b-5p and lymphocyte subsets with disease severity. Binary logistic regression and Cox regression analyses were employed to identify independent risk factors influencing the prognosis of sepsis. The predictive value of miR-146b-5p and lymphocyte subsets for sepsis prognosis was assessed using receiver operating characteristic (ROC) curves.
Results: miR-146b-5p expression was significantly lower in sepsis patients compared to HD group, with levels in the nonsurvivor group being lower than those in the survivor group. Survival curves for miR-146b-5p indicated that lower levels of miR-146b-5p (<0.272) were associated with a higher mortality rate (HR 3.063). The absolute counts of lymphocytes (Lym), CD3+ T cells, CD4+ T cells, and CD8+ T cells were significantly lower in the sepsis group compared to the HD group. Testing lymphocyte counts at different time points revealed that absolute counts of CD3+ T cells, CD4+ T cells, and CD8+ T cells were consistently lower in the sepsis group across all time intervals. miR-146b-5p levels were predictive of patient prognosis, with the combination of miR-146b-5p and APACHE II scores yielding the highest AUC.
Conclusion: The early-stage sepsis-associated downregulation of miR-146b-5p serves as a promising biomarker for severity stratification and prognostic evaluation. The combination of miR-146b-5p with APACHE II scores enhances diagnostic accuracy. Additionally, dynamic monitoring of lymphocyte subsets may facilitate the evaluation of immune status and guide personalized treatment strategies.
1 Introduction
Sepsis is a life-threatening organ dysfunction caused by the dysregulated host immune response to infection and is an important global health problem (1). In 2017, an estimated 48.9 million incident cases of sepsis and 11 million sepsis-related deaths (representing 19.7% of all global deaths) were reported worldwide (2). Despite advances in clinical management of sepsis, sepsis remains the leading cause of death in the intensive care unit (ICU), with long-term mortality rates ranging from 40% to 80% (3–5). Moreover, sepsis treatment imposes a heavy burden on healthcare resources and significantly affects patients’ quality of life. Early diagnosis of sepsis is challenging and immune dysregulation during the course of the disease often leads to immunosuppression, which increases the risk of secondary infections and death (6–10). Therefore, early and accurate prediction of the patient’s condition and prognosis is of paramount clinical importance (11, 12).
The immune system in sepsis patients is severely compromised, marked by a loss of innate immune function and subsequent suppression of adaptive immunity (13, 14). Lymphopenia, characterized by a reduction in lymphocyte counts, is commonly observed in sepsis, with a significant depletion of CD4+ T cells (15). This depletion is strongly associated with increased mortality, underscoring the crucial role of CD4+ T cells as predictive biomarkers and potential therapeutic targets (16).
With the in-depth study of immune mechanisms in sepsis, microRNAs (miRNAs) have emerged as key regulators of immune responses. MiRNAs are a class of small, noncoding RNAs that usually negatively regulate gene expression by inducing messenger RNA (mRNA) cleavage or repressing translation (17). Aberrant miRNA expression has been observed in sepsis, indicating their potential utility in diagnosis and prognosis (18). In recent years, many studies have shown that miRNAs can be used as sepsis biomarkers, offering significant value in disease diagnosis, severity assessment and prognostic evaluation (19–24). MiRNAs are recognized as crucial immune regulators that modulate various phases of inflammation (23, 25), and their dysregulation has been reported in various autoimmune and autoinflammatory diseases, including rheumatoid arthritis, systemic lupus erythematosus (SLE) and inflammatory bowel disease (IBD) (26–28). Using miRNA gene chip technology, we identified miR-146b-5p as a key regulator of inflammation, capable of downregulating pro-inflammatory factors such as TNF-α, IL-6, and HMGB1 in primary peritoneal macrophages of mice induced by Candida albicans (29). However, the role of miR-146b-5p in human sepsis remains unexplored.
This study aimed to investigate the expression and changes of miR-146b-5p and lymphocyte subsets in patients with sepsis. We sought to analyze the correlation between miR-146b-5p and lymphocyte subsets and to assess the predictive value for sepsis prognosis. By elucidating the roles of miR-146b-5p and lymphocyte subsets in sepsis, we aim to identify novel biomarkers for early diagnosis and prognostic assessment, ultimately facilitating the development of individualized therapeutic strategies for sepsis.
2 Materials and methods
2.1 Patients and healthy donors
This prospective cohort study included patients and healthy controls from the ICU of the Second Affiliated Hospital of Chongqing Medical University, China. The study was conducted in accordance with the Declaration of Helsinki and approved by the Human Research Ethics Committee of Chongqing Medical University (Scientific Ethics Pre-Review No. (2018) 110). All procedures adhered to the relevant ethical guidelines and regulations. We rigorously enforced a three-level consent mechanism. First, informed consent was obtained from participants if they were able to provide written informed consent. Second, we rigorously sought informed consent from guardians (e.g., spouses, children) before conducting the study in patients who were unconscious or unconscious. Third, if the patients had no guardians, we would not take such patients as research subjects, and our study strictly complied with the requirements of clinical research medical ethics. A total of 191 consecutive patients aged ≥18 years old were admitted to the ICU between January 2020 and July 2024. Sepsis was diagnosed according to the Sepsis-3 criteria, which defines sepsis as life-threatening organ dysfunction caused by a dysregulated host response to infection. Organ dysfunction was determined by an acute increase of 2 or more points in the Sequential Organ Failure Assessment (SOFA) score. The exclusion criteria included: malignant tumor, autoimmune disease, hematopoietic system disease and immunodeficiency disease; use of radiotherapy, chemotherapy and immunosuppressant within half a year. Healthy donors (HD) were randomly selected from outpatients at the Second Affiliated Hospital of Chongqing Medical University to serve as controls. In order to ensure that the healthy control group and the case group were matched on key characteristics, we took the following approach: (1) Age matching: Individuals with a similar age distribution to the case group were screened to reduce the potential effect of age differences. (2) Gender matching: the gender composition of the healthy controls was determined according to the sex ratio of the case group to ensure that the sex ratio was consistent. (3) Basic health status matching: healthy donors were strictly screened and those with a history of chronic diseases were excluded to ensure the comparability of health status with the case group. Informed consent was obtained from all participants providing samples. All sepsis patients were treated in accordance with the guidelines of the Surviving Sepsis Campaign.
For the two indicators of main interest in this study, miR146b-5p and lymphocyte subsets, there was no missing data for miR146b-5p, and 12 (6%) were missing data for lymphocyte subsets. For the baseline indicators, 76 (39%) were missing data for IL-1β, 114 (59%) for IL-2, 70 (36%) for IL-6, 146 (76%) for IL-8, and 70 (36%) for IL-10. There were 71 (37%) missing data for TNF-α, 15 (7%) missing data for CRP, 21 (10%) missing data for PCT, 97 (50%) missing data for SOFA score, and 46 (18%) missing data for APACHE II score. There were 16 cases (8%) missing data for white blood cells, 15 cases (7%) missing data for neutrophils, 15 cases (7%) missing data for lymphocytes, and 15 cases (7%) missing data for monocytes. For continuous variables with a proportion of missing data (5% - 10%) that were missing at random, we used mean imputation to preserve the completeness of the data. Multiple imputation was used for variables with a high proportion of missing data (>10%) and a complex missing mechanism. The patient selection flow chart, grouping, and data analysis process are presented in Figure 1.
2.2 Sample collection
Blood samples were collected at four time points during patient admission: days 1-3 (D1-3), days 4-7 (D4-7), days 8-12 (D8-12), and days 13-21 (D13-21). These samples were used for the analysis of miR-146b-5p expression and lymphocyte subsets. At admission, the following baseline demographic and clinical information was recorded: age, gender, Acute Physiology and Chronic Health Evaluation (APACHE II) scores, SOFA scores, complete blood count results, and cytokine levels. The primary endpoint was 28-day mortality.
Additional clinical data were collected throughout the study, including vital signs, medical history, laboratory test results, diagnoses, and clinical outcomes. The SOFA and APACHE II scores were calculated within the first 24 hours following the diagnosis of sepsis. In addition, 2 mL of venous blood was collected in EDTA vacutainers within 24 hours of sepsis diagnosis. Peripheral blood mononuclear cells (PBMCs) were isolated from these blood samples using density gradient centrifugation with Ficoll-Hypaque (TBD Science; Tianjin, China), following the manufacturer’s instructions.
2.3 Quantitative real-time polymerase chain reaction analysis
Total RNA was isolated from cells using TRIzol Reagent (Thermo Fisher Scientific). One milliliter of TRIzol and 200 microliters of chloroform were added to the sample and mixed vigorously for 15 seconds, followed by incubation on ice for 10 minutes. The samples were centrifuged at 4°C for 10 minutes at 12,000 rpm to achieve complete phase separation. The upper aqueous phase was collected, and isopropanol was added in a 1:1 ratio to a new centrifuge tube, which was mixed gently and incubated on ice for an additional 10 minutes. After centrifugation at 4°C for 10 minutes at 12,000 rpm, the resulting pellets were washed once with 75% ethanol. The precipitated RNA was resuspended in 30 microliters of RNase-free water. To assess the quantity and quality of RNA, the samples were measured using a NanoDrop spectrophotometer (NanoDrop, Thermo) and analyzed by agarose gel electrophoresis. A total of 1 microgram of RNA was reverse transcribed into cDNA using the PrimeScript RT reagent kit (Takara, Dalian, China). Quantitative real-time PCR analysis was conducted using the SYBR Premix Ex Taq™ II kit (Takara, Dalian, China). Primers were purchased from Sangon, Shanghai, China. The primers for human miR-146b-5p were: 5’- CGCGTGAGAACTGAATTCCAT -3’(sense) and 5’-AGTGCAGGGTCCGAGGTATT -3’ (antisense), and U6 primers were 5’-AGAGAAGATTAGCATGGCCCCTG-3’ (sense) and 5’-AGTGCAGGGTCCGAGGTATT-3’ (antisense). The thermal cycling conditions for qPCR were as follows: 95°C for 30 s (step 1), 40 cycles of 95°C for 5 s, 60°C for 30 s (step 2); 95°C for 15 s, 65°C for 30 s, 95°C for 15 s (step 3), and then maintained at 4°C. The reactions were performed on a CFX Connect Sequence Detection System (Bio-Rad). The mRNA expression levels of target genes were normalized to U6 using the 2-ΔΔCtmethod where ΔCt = target gene Ct-U6 Ct and ΔΔCt = ΔCt treatment-ΔCt control. Three independent experiments were conducted in triplicate.
2.4 Determination of lymphocyte subsets
2 mL of whole blood was collected in ethylenediaminetetraacetic acid (EDTA)-containing tubes and analyzed within 6 hours of collection. PBMCs were stained with fluorochrome-conjugated antibodies specific for CD3, CD4, CD8, CD19, and CD56 (eBioscience, San Diego, USA). After a 30-minute incubation in the dark, the cells were washed, centrifuged, and analyzed using a FACSCanto II flow cytometer (BD, Heidelberg, Germany). FlowJo software (Treestar Inc., Ashland, USA) was used for subsequent analysis.
Cell populations were identified based on antibody positivity and doublet exclusion as follows: T cells were defined as CD3+CD56−, with further analysis for CD4+ and CD8+ T cell subsets. Natural killer (NK) cells were defined as CD3−CD19−CD56+, and B cells were defined as CD19+. The assay was conducted in accordance with the manufacturer’s instructions, with each sample being stained and analyzed following standardized procedures.
2.5 Routine blood and cytokine tests
Routine blood tests were performed in the hospital laboratory using an automated hematology analyzer (Sysmex XN-2000, Sysmex Corporation). Peripheral blood samples (2-3 mL) were collected in tubes containing anticoagulant (e.g., EDTA), and processed according to the manufacturer’s instructions. Blood samples were analyzed using the XN-2000 fully automated hematology analyzer kit (model: XN-02, Sysmex Corporation). The test parameters included white blood cell count (WBC), neutrophil count (NEU), lymphocyte count (LYM), and mononuclear cell count (MONO). Absolute cell counts were automatically calculated using the instrument’s differential white blood cell analysis.
Cytokine analysis was conducted using the Cytometric Bead Array (CBA) assay (BD, Heidelberg, Germany) to perform multiplex analysis of plasma cytokines. The cytokines analyzed included Interleukin (IL)-1β, IL-2, IL-6, IL-8, IL-10, and Tumor Necrosis Factor-α (TNF-α). The assay was conducted according to the manufacturer’s instructions.
2.6 Statistical analysis
Statistical analyses were conducted using SPSS version 25.0 (IBM, Armonk, USA) and GraphPad Prism software version 9.0 (GraphPad Software, La Jolla, USA). Two-group comparisons were performed using the independent samples t tests or chi-square tests. Continuous variables that did not follow a normal distribution were analyzed using the Mann-Whitney U test for two-group comparisons and the Kruskal-Wallis test for multiple-group comparisons. Spearman correlation analysis was employed to assess the correlation between the two variables. Binary logistic regression was utilized to analyze the risk factors influencing prognosis, and Cox regression was applied to examine the independent risk factors affecting the 28-day mortality of patients with sepsis.
Survival curves were generated using the Kaplan-Meier method, with grouping criteria established based on cut-off values determined from receiver operating characteristic (ROC) curve analysis, and differences between groups were assessed using the log-rank test. The area under curve (AUC) was calculated to evaluate the diagnostic efficacy of miR-146b-5p and lymphocyte subsets. The indicators used to establish survival curves included the expression levels of miR-146b-5p and the absolute counts of Lym, CD3+ T cells, CD4+ T cells, and CD8+ T cells, as well as CD4+ T cell percentage (CD3%), CD4+ T cell percentage (CD4%), CD8+ T cell percentage (CD8%), B cell percentage (B%), and NK cell percentage (NK%), and the CD4+/CD8+ ratio. The corresponding cut-off values for these indicators were as follows: miR-146b-5p (0.272), Lym (1.24), the absolute counts of CD3+ T cells (0.145), CD4+ T cells (0.065), CD8+ T cells (0.059), CD3% (53.98), CD4% (23.18), CD8% (27.76), B% (5.91), NK% (8.6), and CD4+/CD8+ ratio (2.54).
3 Results
3.1 Clinical characteristics of 28-day survivors and nonsurvivors
Patients were divided into two groups based on 28-day mortality. A total of 191 sepsis patients were included in the study, with 117 in the non-survivor group and 74 in the survivor group. The mean age in the non-survivor group was 70.30 ± 13.14 years, with 69.23% of patients being male. In contrast, the survivor group had a mean age of 68.46 ± 17.60 years, with 75.68% of patients being male. No significant differences in age or gender were observed between the survivor and non-survivor groups. Laboratory examinations revealed no significant differences in leukocyte, lymphocyte, neutrophil, and monocyte counts between the two groups. However, nonsurvivor patients exhibited significantly higher APACHE II scores of 24.12 ± 8.14 and SOFA scores of 9.89 ± 4.16 than the survivor group. Additionally, levels of IL-6 and IL-10 were markedly elevated in nonsurvivor patients (IL-6, 1307.91 ± 3254.43 pg/ml; IL-10, 123.94 ± 278.73 pg/ml) compared to survivor patients (IL-6, 379.28 ± 902.33 pg/ml, p = 0.018; IL-10, 30.00 ± 73.85 pg/ml, p = 0.005) (Table 1).
3.2 Differential expression of miR-146b-5p in patients with sepsis and its association with patient survival
We analyzed the differences in the relative expression of miR-146b-5p between sepsis patients and HD. The relative expression of miR-146b-5p was significantly lower in the sepsis group than in the HD group. Furthermore, within the sepsis cohort, nonsurvivors exhibited significantly lower miR-146b-5p expression compared to survivors (Figure 2A).

Figure 2. Differential expression of miR-146b-5p in sepsis patients and its relationship with survival outcomes. (A) Box plots displaying miR-146b-5p expression levels in PBMCs from healthy donors (HD) and sepsis patients (nonsurvivors and survivors). (B) Trend of miR-146b-5p expression levels in PBMCs of sepsis patients at different time intervals(D1-3, D4-7, D8-12, D13-21) Time intervals represent days after admission. (C) Comparison of miR-146b-5p expression between nonsurvivors and survivors at multiple time points. (D) Kaplan-Meier survival curves showing the 28-day survival of sepsis patients classified into high (miR-146b-5p > 0.272) and low (miR-146b-5p ≤ 0.272) expression groups. Each dot represents an individual patient’s expression level. The boxes indicate the lower and upper quartile values with the line indicating the median value, the whiskers indicate the lowest and highest values. The Mann-Whitney U test was used for comparisons between two groups and the Kruskal-Wallis test was applied for analysis of three or more groups. Survival differences were analyzed with log-rank tests. The p value indicates the level of statistical significance.
Peripheral blood samples from sepsis patients were collected at four distinct time intervals: days 1-3 (D1-3), days 4-7 (D4-7), days 8-12 (D8-12), and days 13-21 (D13-21) to assess the dynamic changes in miR-146b-5p expression. A comparison of miR-146b-5p expression across these time points with baseline levels from the HD group revealed significantly lower expression in the D1-3, D4-7, and D8-12 groups compared to the HD group (all P < 0.0001). We observed an increasing trend in the expression of miR-146b-5p over time, with levels in the D4-7 group significantly higher than those in the D1-3 group (p = 0.0073), the expression in the D13-21 group was greater than that in the D4-7 group (p = 0.0200). However, no significant difference was found between the D13-21 group and the HD group (p = 0.2763) (Figure 2B).
The relative expression of miR-146b-5p was compared between the nonsurvivor and survivor groups across the time intervals (D 1-3, D 4-7, D 8-12, and D 13-21). At D1-3, miR-146b-5p levels were significantly lower in the nonsurvivor group compared to the survivor group. There is no significant difference in other time intervals. This observation suggests that lower miR-146b-5p expression during the early stages of sepsis is associated with an increased risk of mortality Figure 2C).
The 28-day prognosis of patients with varying levels of miR-146b-5p was evaluated using an optimal cut-off value as the classification criterion. All patients were divided into a high miR-146b-5p group (miR-146b-5p > 0.272) and a low miR-146b-5p group (miR-146b-5p ≤ 0.272). Notably, the Kaplan-Meier survival curve indicated that patients with miR-146b-5p expression below 0.272 had a significantly higher risk of death compared to those with higher expression (log-rank χ² = 39.46, p < 0.001) (Figure 2D). Taken together, these findings suggest that lower relative expression of miR-146b-5p in PBMCs is associated with increased mortality in septic patients.
3.3 Differential expression of lymphocyte subsets in patients with sepsis and its relation to patient survival
We compared lymphocyte subsets between sepsis patients and HD group. The absolute counts of lymphocytes (Lym), CD3+ T cells, CD4+ T cells, and CD8+ T cells were significantly lower in sepsis patients compared to HD. However, B% and NK% showed an increasing trend in sepsis patients compared to HD. We further examined the changes in T cell subsets between nonsurvivors and survivors of sepsis patients. The CD3% and CD4% were significantly lower in nonsurvivors compared to survivors. Notably, no significant difference was observed in the CD4+/CD8+ T cell ratio between sepsis patients and HD, which may be attributed to the concurrent reduction of both CD4+ and CD8+ T cells in sepsis patients (Figure 3).

Figure 3. Box plots displaying lymphocyte subsets levels in PBMCs from healthy donors (HD) and sepsis patients (nonsurvivors and survivors). Each dot represents an individual patient’s expression level. The boxes indicate the lower and upper quartile values with the line indicating the median value, the whiskers indicate the lowest and highest values. Group comparisons were performed by Mann-Whitney U test. The p value indicates the level of statistical significance.
We further compared the changes in the four time intervals. The results indicated that the absolute counts of lymphocytes, CD3+ T cells, CD4+ T cells, and CD8+ T cells were significantly lower than those in the HD group during all periods (p value were less than 0.0001). B% was reduced in the D1-3 and D4-7 groups compared to HD (p = 0.0152 and p = 0.0005, respectively). Additionally, CD8% and NK% were lower in the D4-7 group compared to HD (p = 0.0077 and p = 0.0003, respectively). No significant differences were observed in the CD4%, CD3%, and CD4+/CD8+ T cell ratio compared to HD. Furthermore, NK% was significantly lower in the D4-7 group than in the D1-3 group (p = 0.0175). Furthermore, dynamic analysis revealed that, compared to HD group, used as the baseline, the total lymphocyte count, as well as the absolute counts of CD3+ T cells, CD4+ T cells, and CD8+ T cells, were significantly lower across all four time points (D1-3, D4-7, D8-12, and D13-21). The percentage of B cells increased significantly during the early stage of sepsis, within the first week, and subsequently declined. In contrast, the percentage of NK cells showed the opposite pattern, with a significant decrease in the first week after admission, followed by a gradual increase over time. No significant changes were observed in the proportions of other lymphocyte subsets, including CD3+ T cells, CD4+ T cells, and CD8+ T cells, across the different time intervals (Figure 4).

Figure 4. Trend of lymphocyte subsets in PBMCs of sepsis patients at different time intervals (D1-3, D4-7, D8-12, and D13-21). Time intervals represent days after admission. Each dot represents an individual patient’s expression level. The boxes indicate the lower and upper quartile values with the line indicating the median value, the whiskers indicate the lowest and highest values. Differences between groups were calculated using Kruskal-Wallis test. The p value indicates the level of statistical significance. indicates the level of statistical significance.
The expression of lymphocyte subsets in both the nonsurvivor and survivor groups was compared across the intervals D1-3, D4-7, D8-12, and D13-21. During the D1-3 period, the absolute counts of CD4+ T cells and CD3% were significantly higher in the survivor group compared to the nonsurvivor group. However, the differences in absolute counts of Lym, CD4+ T cells, CD8+ T cells, and CD3+ T cells and CD4%, CD8%, B cell%, NK%, and CD4+/CD8+ T cell ratio were not statistically significant. This finding suggests that in the early stage (within one week of admission), the levels of lymphocyte subsets in the two groups were similar; while in the latter two time periods, the survival group had a higher trend than the nonsurvivor group, although this difference was not statistically significant (Figure 5).

Figure 5. Comparison of Lymphocyte subsets between nonsurvivors and survivors at multiple time points. Time intervals represent days after admission. Each dot represents an individual patient’s expression level. The boxes indicate the lower and upper quartile values with the line indicating the median value, the whiskers indicate the lowest and highest values. Differences between groups were calculated using Mann-Whitney U test. The p value indicates the level of statistical significance.
Kaplan-Meier (K-M) curves showed that among lymphocyte subsets, patients in the group with high expression of CD3%, CD4%, and CD8%, as well as in the group with high expression of the absolute counts of CD3+ T cells, CD4+ T cells, and CD8+ T cells had a significantly lower mortality rate than patients in the group with low expression of these lymphocytes after 28 days (χ2 = 13.16, p = 0.0003; χ2 = 16.09, p < 0.0001; χ2 = 5.676, p = 0.0172; χ2 = 9.894, p = 0.0017; χ2 = 20.06, p < 0.0001; χ2 = 5.642, p = 0.0175). No significant differences were observed between the two groups regarding B%, NK%, CD4+/CD8+ T cell ratio, and Lym (Figure 6).

Figure 6. Kaplan-Meier survival curves for two groups of patients with different lymphocyte subsets levels. Kaplan-Meier survival curves demonstrating lymphocyte subsets (Lym, CD3+ T cell counts, CD4+ T cell counts, CD8+ T cell counts, CD3+ T cell (%), CD4+ T cell (%), CD8+ T cell (%), CD4+/CD8+, B cell (%), NK cell (%) in sepsis patients in the high and low levels groups for 28-day survival. HR indicates hazard ratio; 95% CI, 95% confidence interval. Survival differences were analyzed with log-rank tests. The p value indicates the level of statistical significance.
3.4 Correlation study of relative expression of miR-146b-5p and lymphocyte subsets and disease severity in patients with sepsis
Spearman correlation coefficients were employed to examine the relationship between miR-146b-5p expression levels, lymphocyte subsets, and APACHE II and SOFA scores to assess their correlation with disease severity. Notably, CD3% and CD4% exhibited a negative correlation with the APACHE II scores in sepsis patients (r = -0.214, p = 0.0103; r = -0.237, p = 0.0043), while CD3% showed a negative correlation with the SOFA score (r = -0.242, p = 0.0287). When analyzing the relative expression of miR-146b-5p, no significant correlation was observed with either the APACHE II or SOFA scores in the overall assessment. During days 1 to 3 (D1-3), the relative expression of miR-146b-5p, CD3%, and CD4% demonstrated negative correlations with the APACHE II scores in sepsis patients (r = -0.224, p = 0.0241; r = -0.280, p = 0.0046; r = -0.240, p = 0.0156). From days 4 to 7 (D4-7), a negative correlation was identified between miR-146b-5p expression and the CD4+/CD8+ T cell ratio with the APACHE II scores (r = -0.435, p = 0.0338; r = -0.592, p = 0.0426), while positive correlations were observed between CD8% and NK% with the SOFA scores (r = 0.677, p = 0.0155; r = 0.681, p = 0.0211). During days 8 to 12 (D8-12), no associations were found between miR-146b-5p expression levels and lymphocyte subsets with the APACHE II and SOFA scores. From days 13 to 21 (D13-21), absolute counts of CD4+ T cells were positively correlated with the SOFA score (r = 0.900, p = 0.0374; r = 1.000, p < 0.0001) and negatively correlated with NK% (r = -1.000, p < 0.0001). These findings suggest that lower expression levels of miR-146b-5p, CD3%, and CD4% are associated with higher APACHE II scores, which indicate a more severe condition in sepsis patients. However, it is important to note that the strength of these correlations is relatively weak (Table 2).

Table 2. Correlation between miR-146b-5p expression and lymphocyte subsets in relation to disease severity.
3.5 Logistic regression analysis of miR-146b-5p and lymphocyte subsets as independent predictors of sepsis
Univariate logistic regression results indicated that the expression level of miR-146b-5p, CD3%, CD4%, APACHE II scores, and SOFA scores were independent risk factors for sepsis prognosis (p = 0.003, 0.002, 0.037, 0.001, 0.001). Additionally, multivariate logistic regression results demonstrated that the expression level of miR-146b-5p and CD3% were significant risk factors influencing the prognosis of sepsis patients. The above findings suggest that miR-14bb-5p expression and CD3% may serve as reliable biomarkers for predicting disease outcomes (p = 0.027, 0.032) (Table 3).

Table 3. Logistic regression analysis (univariate and multivariate) of miR-146b-5p and lymphocyte subsets as independent predictors of sepsis.
3.6 Cox proportional hazard model of relative expression of miR-146b-5p and 28-day mortality
The Cox proportional hazard model with time-dependent covariates indicated that, in the univariate analysis, the relative expression of miR-146b-5p in PBMCs, APACHE II scores, SOFA scores, CD3%, CD4%, procalcitonin (PCT) and IL-6 were significantly associated with the 28-day mortality. The hazard ratios (HR) for these variables were as follows: relative expression of miR-146b-5p (HR = 0.452, 95% CI 0.297-0.690), APACHE II scores (HR =1.064, 95% CI 1.041-1.087), SOFA scores (HR = 1.155, 95% CI 1.077-1.239), CD3% ((HR = 0.976, 95% CI 0.964-0.988), CD4% with hazard ratios of (HR =0.981, 95% CI 0.966-0.996), PCT (HR = 1.003, 95% CI 1.001-1.006), and IL-6 (HR = 1.000, 95% CI 1.000-1.000) were significantly associated with 28-day mortality in the univariate analysis. In the multivariate analysis, the relative expression of miR-146b-5p and the APACHE II scores remained significant, with hazard ratios of 0.145 (95% CI, 0.044-0.476) and 1.103 (95% CI, 1.027-1.184). These findings suggest that miR-146b-5p expression and APACHE II scores may serve as more reliable predictors of disease outcomes (Table 4).

Table 4. Cox regression analysis (univariate and multivariate) of factors associated with 28-day mortality in patients with sepsis.
3.7 Correlation analysis of relative expression of miR-146b-5p with lymphocyte subsets
We used Spearman correlation analysis to examine the relationship between miR-146b-5p and lymphocyte subsets. Overall, miR-146b-5p exhibited a positive correlation with CD3% and CD4% (r = 0.242, p = 0.001; r = 0.155, p = 0.038). During the D1-3 period, miR-146b-5p showed positive correlations with CD3%, CD4%, and absolute counts of CD3+ T and CD4+ T cells (r = 0.312, p = 0.001; r = 0.204, p = 0.039; r = 0.233, p = 0.028; r = 0.233, p = 0.018). In the D4-7 period, miR-146b-5p was positively correlated with CD3% and CD4% (r = 0.453, p = 0.045; r = 0.522, p = 0.018). Conversely, during D8-12, miR-146b-5p exhibited negative correlations with CD4%, NK%, and CD4+/CD8+ T cell ratios (r = -0.566, p = 0.001; r = -0.658, p = 0.014; r = -0.658, p = 0.014). Moreover, at D13-21, miR-146b-5p was not correlated with any lymphocyte subsets. These findings indicate that miR-146b-5p is correlated with certain lymphocyte subsets, although the correlation is relatively weak, suggesting that changes in miR-146b-5p expression may affect the immune status of patients (Table 5).
3.8 ROC curves for variables predicting 28-day mortality risk in sepsis
The logistic regression analysis and Cox regression analysis were used to evaluate the predictive value of serum miR-146b-5p, APACHE II scores, CD3%, and CD4% levels for the prognosis of patients with sepsis (Figure 7A). The AUC values for miR-146b-5p, APACHE II scores, CD3%, and CD4% in predicting the prognosis of sepsis patients were 0.691, 0.761, 0.633, and 0.611, respectively (Table 6). To further improve predictive accuracy, logistic regression analysis was used to calculate the combined prediction probability of miR-146b-5p, APACHE II scores, CD3%, and CD4%. ROC curve analysis was then conducted to evaluate the prognostic value of these combined indicators in sepsis patients (Figure 7B). The AUC results for the three indicators were 0.785, 0.765, and 0.646, respectively, all of which were superior to those predicted by any single index (Table 6). Notably, the combination of miR-146b-5p and the APACHE II scores yielded the highest AUC.

Figure 7. Receiver operating characteristic (ROC) curve for predicting 28-day mortality in patients with sepsis. (A) ROC curve of miR-146b-5p expression, APACHE II scores, CD3%, and CD4% levels. (B) ROC curve of miR-146b-5p expression, in combined with the APACHE II scores, CD3%, and CD4%. ROC curves were plotted using the true positive rate (sensitivity) against the false positive rate (1- specificity) at various cut-off points for miR-146b-5p expression, APACHE II scores, CD3%, and CD4% in patients with sepsis. The areas under the ROC curves were calculated to show the values of AUC. The p value indicates the level of statistical significance.
4 Discussion
The pathogenesis of sepsis involves complex systemic inflammatory network effects, genetic polymorphisms, immune dysfunction, coagulation abnormalities, tissue damage, and abnormal host responses to different infectious pathogenic microorganisms and their toxins, which are closely related to multi-system and multi-organ pathophysiological alterations in the organism (30, 31). As a complex systemic inflammatory response syndrome triggered by infection, its pathological process undergoes a dynamic change from hyperinflammatory response to immunosuppression, the latter being one of the key factors leading to high morbidity and mortality (32–37).
Lymphocyte subsets reflect the immunity of the organism, T-lymphocytes, play a central role in mediating cellular immunity and are key indicators of immune status in sepsis (13). CD3+ T cells represent the total number of mature T-lymphocytes, are essential for adaptive immune responses. Double-negative T cells (DNs, CD-CD8-), which are predominantly localized in the skin and mucosal tissues, participating in innate immunity and forming the first line of defense against infections. CD4+ T cells (helper T cells) are responsible for relaying pathogenic signals, while CD8+ T cells (cytotoxic T cells) directly eliminate infected or abnormal cells through cytotoxic mechanisms (38). A reduction in CD3+ T cells, CD4+ T cells, and CD8+ T cells indicate suppressed immune function, and abnormalities in immune function (e.g., lymphopenia) alongside an imbalance of inflammatory factors have been suggested to play a significant role in the development of sepsis. However, it is important to note that the predictive value of individual lymphocyte subsets for clinical outcomes remains limited. A comprehensive evaluation of immune markers, including lymphocyte subsets, cytokine levels, and other inflammatory mediators, is necessary to accurately assess immune status and predict prognosis in sepsis patients.
In recent years, miRNAs have emerged as critical regulators of inflammation and immune responses, making them a promising area of research in sepsis (20, 39–41). Interestingly, the same miRNA can exhibit different functions in different diseases. Studies have found that miR-146b-5p is associated with a variety of diseases, including sepsis, cervical cancer, and hepatocellular carcinoma (HCC) (42–44). In HCC, miRNAs such as miR-15a, miR-125b, and miR-122 have been proposed as potential biomarkers and therapeutic targets (45, 46). Recent studies have highlighted the diagnostic and prognostic potential of miRNA signatures in sepsis. Yao et al. reported that circulating levels of miR-25 were significantly lower in sepsis patients compared to those with systemic inflammatory response syndrome (SIRS), and lower miR-25 levels were associated with worse clinical outcomes (47). Formosa A et al. found that miR-182, miR-143, miR-145, miR-146a, miR-150, and miR-155 were dysregulated in sepsis patients, and dysregulated miRNAs have immunological associations with clinical disease in sepsis, especially, the decrease in miR-146a correlated with elevated IL-6 levels and monocyte proliferation, indicating immunological associations with clinical disease in sepsis (20, 39, 48–50). Martin Jouza et al. identified several miRNAs independently validated by different teams in particular, which are miR-16a, miR-16, miR-96-5p, miR-141, miR-181a, and miR-1184 as having strong potential for early diagnosis of neonatal sepsis (51). This study focuses on miR-146b-5p, and through time course analysis, it reveals its dynamic changes during the early and subsequent stages of sepsis. It was found that miR-146b-5p remains at a consistently lower level in nonsurvivors and shows a certain correlation with the dynamic changes of lymphocyte subsets, providing a new perspective for understanding immune dysregulation in sepsis.
In this study, we analyzed ICU patients with sepsis and found no significant differences in age, gender, or baseline laboratory parameters (e.g., leukocyte, lymphocyte, neutrophil, and monocyte counts) between survivors and nonsurvivors. However, nonsurvivors exhibited higher APACHE II and SOFA scores, along with elevated levels of IL-6 and IL-10, highlighting the utility of these scores as indicators of disease severity and prognosis in critically ill patients (52–55).
We further investigated the differential expression and dynamic trends of miR-146b-5p and lymphocyte subset levels in patients with sepsis, comparing nonsurvivors to survivors. Our results revealed that miR-146b-5p levels were significantly lower in sepsis patients compared to healthy controls during the early stages of the disease, with nonsurvivors exhibiting even lower levels than survivors. Notably, miR-146b-5p levels began to recover after two weeks, suggesting a potential role in the immune response to sepsis. We also observed a weak negative correlation between miR-146b-5p levels and APACHE II scores during the early stages of sepsis, further supporting its potential as a prognostic marker. Survival analysis indicated that patients with miR-146b-5p levels below 0.272 had a significantly higher mortality rate (HR = 3.063), underscoring the clinical relevance of this miRNA in sepsis prognosis.
Our findings also revealed significant abnormalities in immune function among sepsis patients, particularly in the quantity and distribution of lymphocyte subsets. Compared to HD, sepsis patients exhibited markedly reduced absolute counts of Lym, CD3+ T cells, CD4+ T cells, and CD8+ T cells, as well as decreased B% and NK%. These findings are consistent with previous studies demonstrating lymphopenia and immune suppression in sepsis patients. Longitudinal analysis further revealed that CD3+, CD4+, and CD8+ T cell counts remained lower in sepsis patients than in healthy controls at all time points during the early phase of the disease. Additionally, differences in CD4% and CD3%between survivors and nonsurvivors may reflect the severity of immune dysfunction and its impact on clinical outcomes. Although the correlations between CD3%, CD4%, and APACHE II scores were weak, they suggest that immune cell percentages may serve as indirect markers of disease severity in sepsis patients.
Importantly, we found that peripheral blood levels of miR-146b-5p were significantly lower in nonsurviving sepsis patients compared to HD, and these levels were positively correlated with CD3% and CD4%. This suggests a potential link between miR-146b-5p expression and immune function in sepsis. Patients with reduced lymphocyte subsets also exhibited lower miR-146b-5p levels compared to those with normal lymphocyte counts, further supporting the role of miR-146b-5p in immune regulation. These findings highlight the potential of miR-146b-5p as a novel therapeutic target for sepsis. The precise mechanism by which miR-146b-5p regulates T cell function in sepsis has yet to be fully elucidated. Our research group has currently validated in vitro that miR-146b-5p can directly and positively regulate the expression of the target gene NF-κB in macrophages co-induced by LPS and ATP. This regulation occurs through the promotion of pyroptosis-related proteins such as GSDMD, caspase-11, and ASC, which further enhances the release of inflammatory factors IL-6 and IL-1β, thereby participating in cellular pyroptosis. These findings indicate that miRNA can regulate cellular pyroptosis in sepsis by targeting proteins involved in pyroptosis-related signaling pathways (56). In our study, we observed a lower relative expression of miR-146b-5p and higher levels of IL-6 in non-survivors, suggesting that miR-146b-5p may influence sepsis prognosis by modulating inflammatory responses. The NF-κB pathway plays a central role in T cell activation and the production of inflammatory factors, indicating that miR-146b-5p may regulate T cell function through the potential NF-κB/IL-6 pathway. Future studies are needed to validate these hypotheses and identify specific targets of miR-146b-5p in the regulation of T cell function in sepsis, which would provide new therapeutic targets for improving immune function and clinical outcomes in sepsis patients.
To identify independent risk factors influencing sepsis prognosis, we employed binary logistic regression (including univariate and multivariate analyses) and Cox regression (also incorporating univariate and multivariate analyses). These analyses provide a comprehensive evaluation of clinical prognostic factors and offer insights for personalized treatment strategies. Finally, we utilized ROC curves to assess the predictive value of miR-146b-5p, lymphocyte subsets, and APACHE II scores in sepsis prognosis. Our results demonstrated that miR-146b-5p levels effectively predicted patient prognosis, exhibiting the largest AUC in the ROC analysis. The diagnostic efficacy was further enhanced when miR-146b-5p was combined with the APACHE II scores, suggesting its potential utility in clinical practice. However, the predictive value of lymphocyte subsets requires further validation.
The limitations of our research should be noted. Firstly, our research was conducted at a single center in Chongqing, China, which may limit the applicability of our findings to other regions of China or different countries. To mitigate this limitation, we strictly adhered to the inclusion and exclusion criteria, ensuring a certain degree of universality in our data. Additionally, as the Second Affiliated Hospital of Chongqing Medical University is a teaching hospital, the cases we collected can, to some extent, reflect the situation of sepsis patients in other affiliated hospitals of Chongqing Medical University. Secondly, due to resource constraints, the sample size we collected was limited. Therefore, it is necessary to conduct large-scale, multi-center, randomized controlled trials to further validate our findings and comprehensively assess their universality and specificity in diverse populations.
5 Conclusion
In conclusion, our study highlights the significant downregulation of miR-146b-5p during the early stages of sepsis and its potential as a biomarker for assessing disease severity and prognostic risk. The combination of miR-146b-5p with APACHE II scores improved diagnostic accuracy, underscoring its clinical relevance. Additionally, dynamic monitoring of lymphocyte subsets may aid in evaluating immune status and guiding personalized treatment strategies. The positive correlation between miR-146b-5p levels and certain lymphocyte subsets suggests that miR-146b-5p may modulate immune responses by influencing T-cell function. However, the precise mechanisms underlying this relationship remain unclear and warrant further investigation through cellular experiments and animal models.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Human Research Ethics Committee of Chongqing Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
SC: Writing – original draft, Formal analysis, Software, Visualization, Data curation, Conceptualization, Methodology, Writing – review & editing, Validation, Investigation. JJW: Visualization, Writing – original draft, Data curation, Formal analysis. PFW: Data curation, Writing – review & editing, Project administration. XLL: Data curation, Writing – review & editing. HLL: Writing – review & editing, Methodology. XYL: Visualization, Writing – review & editing. XWZ: Writing – review & editing, Visualization. XMZ: Software, Writing – review & editing. CXW: Writing – review & editing, Supervision, Formal analysis. HS: Conceptualization, Writing – review & editing, Supervision, Resources, Investigation, Funding acquisition.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the National Natural Science Foundation of China (No. 81871608).
Acknowledgments
We sincerely thank the Human Research Ethics Committee of Chongqing Medical University for its important support in the design, initiation, conduct, and completion of this study. We would like to thank all the medical and nursing staff of the Department of Intensive Care Unit, the Second Affiliated Hospital, Chongqing Medical University for their professional support in case data collection and sample processing and Xiaoyan Wang, PhD, for his review and input on statistical analyses. At the same time, we keep in mind the cooperation and trust of all the patients and families who participated in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Glossary
APACHE: Acute physiology and chronic health evaluation
AUC: Area under the curve
CI: Confidence interval
CRP: C-reactive protein
EDTA: Ethylenediaminetetraacetic acid
HCC: Hepatocellular carcinoma
HD: Healthy donors
HMGB1: High mobility group box 1 protein
IBD: Inflammatory bowel disease
ICU: Intensive care unit
IL-1β: Interleukin-1β
IL-2: Interleukin-2
IL-6: Interleukin-6
IL-8: Interleukin-8
IL-10: Interleukin-10
K-M: Kaplan-meier
LYM: Lymphocyte
MONO: Mononuclear
mRNA: Messenger RNA
NEU: Neutrophil
NK: Natural killer
PBMCs: Peripheral blood mononuclear cells
PCT: Procalcitonin
ROC: Receiver operating characteristic
SIRS: Systemic inflammatory response syndrome
SLE: Systemic lupus erythematosus
SOFA: Sequential organ failure assessment
TNF-α: Tumor necrosis factor alpha
WBC: White blood cell count
References
1. Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. (2016) 315:801–10. doi: 10.1001/jama.2016.0287
2. Rudd KE, Johnson SC, Agesa KM, Shackelford KA, Tsoi D, Kievlan DR, et al. Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study. Lancet. (2020) 395:200–11. doi: 10.1016/S0140-6736(19)32989-7
3. Mayr FB, Yende S, and Angus DC. Epidemiology of severe sepsis. Virulence. (2014) 5:4–11. doi: 10.4161/viru.27372
4. Venet F, Pachot A, Debard AL, Bohe J, Bienvenu J, Lepape A, et al. Increased percentage of CD4+CD25+ regulatory T cells during septic shock is due to the decrease of CD4+CD25- lymphocytes. Crit Care Med. (2004) 32:2329–31. doi: 10.1097/01.ccm.0000145999.42971.4b
5. Fleischmann C, Scherag A, Adhikari NK, Hartog CS, Tsaganos T, Schlattmann P, et al. Assessment of global incidence and mortality of hospital-treated sepsis. Current estimates and limitations. Am J Respir Crit Care Med. (2016) 193:259–72. doi: 10.1164/rccm.201504-0781OC
6. Goh KH, Wang L, Yeow AYK, Poh H, Li K, Yeow JJL, et al. Artificial intelligence in sepsis early prediction and diagnosis using unstructured data in healthcare. Nat Commun. (2021) 12:711. doi: 10.1038/s41467-021-20910-4
7. Jaurila H, Koivukangas V, Koskela M, Gäddnäs F, Myllymaa S, Kullaa A, et al. NMR based metabolomics in human sepsis and healthy serum. Metabolites. (2020) 10. doi: 10.3390/metabo10020070
8. Choi YJ, Kim SB, Kim JH, Park SH, Park MS, Kim JM, et al. Impaired polyfunctionality of CD8(+) T cells in severe sepsis patients with human cytomegalovirus reactivation. Exp Mol Med. (2017) 49:e382. doi: 10.1038/emm.2017.146
9. Carvelli J, Piperoglou C, Bourenne J, Farnarier C, Banzet N, Demerlé C, et al. Imbalance of circulating innate lymphoid cell subpopulations in patients with septic shock. Front Immunol. (2019) 10:2179. doi: 10.3389/fimmu.2019.02179
10. Walton AH, Muenzer JT, Rasche D, Boomer JS, Sato B, Brownstein BH, et al. Reactivation of multiple viruses in patients with sepsis. PloS One. (2014) 9:e98819. doi: 10.1371/journal.pone.0098819
11. Ameri A, Ahmed HM, Pecho RDC, Arabnozari H, Sarabadani H, Esbati R, et al. Diverse activity of miR-150 in Tumor development: shedding light on the potential mechanisms. Cancer Cell Int. (2023) 23:261. doi: 10.1186/s12935-023-03105-3
12. Angus DC, Linde-Zwirble WT, Lidicker J, Clermont G, Carcillo J, and Pinsky MR. Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Crit Care Med. (2001) 29:1303–10. doi: 10.1097/00003246-200107000-00002
13. Ryan T, Coakley JD, and Martin-Loeches I. Defects in innate and adaptive immunity in patients with sepsis and health care associated infection. Ann Transl Med. (2017) 5:447. doi: 10.21037/atm.2017.09.21
14. Luperto M and Zafrani L. T cell dysregulation in inflammatory diseases in ICU. Intensive Care Med Exp. (2022) 10:43. doi: 10.1186/s40635-022-00471-6
15. Cabrera-Perez J, Condotta SA, James BR, Kashem SW, Brincks EL, Rai D, et al. Alterations in antigen-specific naive CD4 T cell precursors after sepsis impairs their responsiveness to pathogen challenge. J Immunol. (2015) 194:1609–20. doi: 10.4049/jimmunol.1401711
16. Hohlstein P, Gussen H, Bartneck M, Warzecha KT, Roderburg C, Buendgens L, et al. Prognostic relevance of altered lymphocyte subpopulations in critical illness and sepsis. J Clin Med. (2019) 8. doi: 10.3390/jcm8030353
17. Baek D, Villen J, Shin C, Camargo FD, Gygi SP, and Bartel DP. The impact of microRNAs on protein output. Nature. (2008) 455:64–71. doi: 10.1038/nature07242
18. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. (2004) 116:281–97. doi: 10.1016/s0092-8674(04)00045-5
19. Szilagyi B, Fejes Z, Poliska S, Pocsi M, Czimmerer Z, Patsalos A, et al. Reduced miR-26b expression in megakaryocytes and platelets contributes to elevated level of platelet activation status in sepsis. Int J Mol Sci. (2020) 21. doi: 10.3390/ijms21030866
20. Formosa A, Turgeon P, and Dos Santos CC. Role of miRNA dysregulation in sepsis. Mol Med. (2022) 28:99. doi: 10.1186/s10020-022-00527-z
21. Fabri-Faja N, Calvo-Lozano O, Dey P, Terborg RA, Estevez MC, Belushkin A, et al. Early sepsis diagnosis via protein and miRNA biomarkers using a novel point-of-care photonic biosensor. Anal Chim Acta. (2019) 1077:232–42. doi: 10.1016/j.aca.2019.05.038
22. Sandquist M and Wong HR. Biomarkers of sepsis and their potential value in diagnosis, prognosis and treatment. Expert Rev Clin Immunol. (2014) 10:1349–56. doi: 10.1586/1744666X.2014.949675
23. Testa U, Pelosi E, Castelli G, and Labbaye C. miR-146 and miR-155: two key modulators of immune response and tumor development. Noncoding RNA. (2017) 3. doi: 10.3390/ncrna3030022
24. Iro MA and Soundara Pandi SP. Clinical application of non-coding RNAs in sepsis. Curr Opin Infect Dis. (2020) 33:530–9. doi: 10.1097/QCO.0000000000000688
25. Boyd JH, Russell JA, and Fjell CD. The meta-genome of sepsis: host genetics, pathogens and the acute immune response. J Innate Immun. (2014) 6:272–83. doi: 10.1159/000358835
26. Ouboussad L, Hunt L, Hensor EMA, Nam JL, Barnes NA, Emery P, et al. Profiling microRNAs in individuals at risk of progression to rheumatoid arthritis. Arthritis Res Ther. (2017) 19:288. doi: 10.1186/s13075-017-1492-9
27. Hosokawa K, Muranski P, Feng X, Keyvanfar K, Townsley DM, Dumitriu B, et al. Identification of novel microRNA signatures linked to acquired aplastic anemia. Haematologica. (2015) 100:1534–45. doi: 10.3324/haematol.2015.126128
28. Kalla R, Ventham NT, Kennedy NA, Quintana JF, Nimmo ER, Buck AH, et al. MicroRNAs: new players in IBD. Gut. (2015) 64:504–17. doi: 10.1136/gutjnl-2014-307891
29. Cheng J, Tian X, Wu C, Wang J, Liu H, Cheng S, et al. MiR- 146b-5p inhibits Candida albicans-induced inflammatory response through targeting HMGB1 in mouse primary peritoneal macrophages. Heliyon. (2025) 11:e41464. doi: 10.1016/j.heliyon.2024.e41464
30. Tauber SC, Eiffert H, Bruck W, and Nau R. Septic encephalopathy and septic encephalitis. Expert Rev Anti Infect Ther. (2017) 15:121–32. doi: 10.1080/14787210.2017.1265448
31. Cheng Z, Abrams ST, Toh J, Wang SS, Wang Z, Yu Q, et al. The critical roles and mechanisms of immune cell death in sepsis. Front Immunol. (2020) 11:1918. doi: 10.3389/fimmu.2020.01918
32. Torres LK, Pickkers P, and van der Poll T. 1Sepsis-induced immunosuppression. Annu Rev Physiol. (2022) 84:157–81. doi: 10.1146/annurev-physiol-061121-040214
33. Padovani CM and Yin K. Immunosuppression in sepsis: biomarkers and specialized pro-resolving mediators. Biomedicines. (2024) 12. doi: 10.3390/biomedicines12010175
34. van der Poll T, van de Veerdonk FL, Scicluna BP, and Netea MG. The immunopathology of sepsis and potential therapeutic targets. Nat Rev Immunol. (2017) 17:407–20. doi: 10.1038/nri.2017.36
35. Venet F and Monneret G. Advances in the understanding and treatment of sepsis-induced immunosuppression. Nat Rev Nephrol. (2018) 14:121–37. doi: 10.1038/nrneph.2017.165
36. van der Poll T, Shankar-Hari M, and Wiersinga WJ. The immunology of sepsis. Immunity. (2021) 54:2450–64. doi: 10.1016/j.immuni.2021.10.012
37. Pei F, Yao RQ, Ren C, Bahrami S, Billiar TR, Chaudry IH, et al. Expert consensus on the monitoring and treatment of sepsis-induced immunosuppression. Mil Med Res. (2022) 9:74. doi: 10.1186/s40779-022-00430-y
38. Abd El-Kader SM and Al-Shreef FM. Inflammatory cytokines and immune system modulation by aerobic versus resisted exercise training for elderly. Afr Health Sci. (2018) 18:120–31. doi: 10.4314/ahs.v18i1.16
39. Antonakos N, Gilbert C, Théroude C, Schrijver IT, and Roger T. Modes of action and diagnostic value of miRNAs in sepsis. Front Immunol. (2022) 13:951798. doi: 10.3389/fimmu.2022.951798
40. Ma C, Liu K, Wang F, Fei X, Niu C, Li T, et al. Neutrophil membrane-engineered Panax ginseng root-derived exosomes loaded miRNA 182-5p targets NOX4/Drp-1/NLRP3 signal pathway to alleviate acute lung injury in sepsis: experimental studies. Int J Surg. (2024) 110:72–86. doi: 10.1097/JS9.0000000000000789
41. Zhao X, Xie J, Duan C, Wang L, Si Y, Liu S, et al. ADAR1 protects pulmonary macrophages from sepsis-induced pyroptosis and lung injury through miR-21/A20 signaling. Int J Biol Sci. (2024) 20:464–85. doi: 10.7150/ijbs.86424
42. Lehto P, Skarp S, Saukko T, Säkkinen H, Syrjälä H, Kerkelä R, et al. Postmortem analyses of myocardial microRNA expression in sepsis. Sci Rep. (2024) 14:29476. doi: 10.1038/s41598-024-81114-6
43. Kurata JS. MicroRNA-focused CRISPR-Cas9 library screen reveals fitness-associated miRNAs. RNA(New York, N.Y.) (2018) 24:966–81. doi: 10.1261/rna.066282.118
44. Ding W, Yang H, Gong S, Shi W, Xiao J, Gu J, et al. Candidate miRNAs and pathogenesis investigation for hepatocellular carcinoma based on bioinformatics analysis. Oncol Lett. (2017) 13:3409–14. doi: 10.3892/ol.2017.5913
45. Amr KS, Elmawgoud Atia HA, Elazeem Elbnhawy RA, and Ezzat WM. Early diagnostic evaluation of miR-122 and miR-224 as biomarkers for hepatocellular carcinoma. Genes Dis. (2017) 4:215–21. doi: 10.1016/j.gendis.2017.10.003
46. Kim HS, Kim JS, Park NR, Nam H, Sung PS, Bae SH, et al. Exosomal miR-125b exerts anti-metastatic properties and predicts early metastasis of hepatocellular carcinoma. Front Oncol. (2021) 11:637247. doi: 10.3389/fonc.2021.637247
47. Yao L LZ, Zhu J, Li B, Chai C, and Tian Y. Clinical evaluation of circulating microRNA-25 level change in sepsis and its potential relationship with oxidative stress. Int J Clin Exp Pathol. (2015) 8(7):7675–84.
48. Maiese A, Scatena A, Costantino A, Chiti E, Occhipinti C, La Russa R, et al. Expression of microRNAs in sepsis-related organ dysfunction: A systematic review. Int J Mol Sci. (2022) 23. doi: 10.3390/ijms23169354
49. Bindayna K. MicroRNA as sepsis biomarkers: A comprehensive review. Int J Mol Sci. (2024) 25. doi: 10.3390/ijms25126476
50. Zhou J, Chaudhry H, Zhong Y, Ali MM, Perkins LA, Owens WB, et al. Dysregulation in microRNA expression in peripheral blood mononuclear cells of sepsis patients is associated with immunopathology. Cytokine. (2015) 71:89–100. doi: 10.1016/j.cyto.2014.09.003
51. Jouza M, Bohosova J, Stanikova A, Pecl J, Slaby O, and Jabandziev P. MicroRNA as an early biomarker of neonatal sepsis. Front Pediatr. (2022) 10:854324. doi: 10.3389/fped.2022.854324
52. Tian Y, Yao Y, Zhou J, Diao X, Chen H, Cai K, et al. Dynamic APACHE II score to predict the outcome of intensive care unit patients. Front Med (Lausanne). (2021) 8:744907. doi: 10.3389/fmed.2021.744907
53. Knaus WA, Draper EA, Wagner DP, and Zimmerman JE. APACHE II: a severity of disease classification system. Crit Care Med. (1985) 13:818–29. doi: 10.1097/00003246-198510000-00009
54. Raith EP, Udy AA, Bailey M, McGloughlin S, MacIsaac C, Bellomo R, et al. Prognostic Accuracy of the SOFA Score, SIRS Criteria, and qSOFA Score for In-Hospital Mortality among Adults with Suspected Infection Admitted to the Intensive Care Unit. JAMA. (2017) 317:290–300. doi: 10.1001/jama.2016.20328
55. Thakur R, Naga Rohith V, and Arora JK. Mean SOFA score in comparison with APACHE II score in predicting mortality in surgical patients with sepsis. Cureus. (2023) 15:e36653. doi: 10.7759/cureus.36653
Keywords: diagnosis and prognosis of sepsis, immunomodulation, lymphocyte subsets, miR-146b-5p, sepsis
Citation: Cheng S, Wang J, Wang P, Long X, Liu H, Li X, Zhang X, Zhao X, Sun H and Wu C (2025) Relationship between peripheral blood miR-146b-5p levels and lymphocyte subsets in patients with sepsis and its predictive value for prognosis: a single-center study. Front. Immunol. 16:1608576. doi: 10.3389/fimmu.2025.1608576
Received: 09 April 2025; Accepted: 06 June 2025;
Published: 26 June 2025.
Edited by:
Juan Bautista De Sanctis, Palacký University Olomouc, CzechiaReviewed by:
Georgia Damoraki, National and Kapodistrian University of Athens, GreeceAlexis Hipólito García, Central University of Venezuela, Venezuela
Copyright © 2025 Cheng, Wang, Wang, Long, Liu, Li, Zhang, Zhao, Sun and Wu. 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: Hang Sun, MzAwNjEzQGNxbXUuZWR1LmNu; Chuanxin Wu, MzAwMzk1QGNxbXUuZWR1LmNu
†These authors have contributed equally to this work