Analysis of clinical baseline data
A total of 106 healthy individuals and 114 children with severe pneumonia were enrolled in this study. There were no statistically significant differences in gender distribution, age, or body mass index (BMI) between the two groups (all P > 0.05), indicating balanced baseline characteristics and good comparability (Table 1). This reduces the potential interference of these confounding factors on subsequent indicator differences.
However, compared to healthy individuals, children with severe pneumonia showed significantly higher levels of WBC, neutrophil count (NEU), PCT, and CRP (all P < 0.001), while LYM was notably decreased (P < 0.001). These findings demonstrate that children with severe pneumonia undergo a strong systemic inflammatory response, which may be accompanied by immune dysfunction. Additionally, the significant increases in brain natriuretic peptide (BNP) and D-dimer (D-D) levels in the patient group (both P < 0.001) suggest that severe pneumonia can affect cardiac function, indicating that the disease extends beyond pulmonary inflammation alone.
Clinical value of miR-136-5p in children with severe pneumonia
Results showed that miR‑136‑5p expression was significantly upregulated in children with severe pneumonia compared with healthy individuals, with a Hodges‑Lehmann median difference of 0.290 (95% CI: 0.233–0.332, P < 0.001) (Figure 1A). Subsequently, ROC curve analysis was performed to evaluate its diagnostic value, which revealed that miR-136-5p exhibited favorable discriminatory efficacy for childhood severe pneumonia compared with healthy controls [AUC = 0.845, 95% CI = 0.795–0.895, optimal cut‑off = 1.245, Youden index = 0.536, sensitivity = 64.9%, specificity = 88.7%, positive likelihood ratio (LR+) = 5.73, negative likelihood ratio (LR−) = 0.4] (Figure 1B).
To clarify the association between miR-136-5p and disease-related indicators, further correlation analysis was conducted between miR-136-5p expression and clinical parameters in children with severe pneumonia. As shown in Table 2, miR-136-5p expression in the peripheral blood of patients was significantly correlated with core clinical inflammatory and functional indicators, specifically exhibiting a positive correlation with WBC (r = 0.803), NEU (r = 0.832), PCT (r = 0.718), CRP (r = 0.829), BNP (r = 0.661), and D‑dimer (r = 0.710) levels, as well as a negative correlation with LYM (r = -0.826) (all P < 0.05).
Patients were stratified into low-expression and high-expression groups based on the median miR-136-5p level. Kaplan‑Meier time‑to‑event analysis demonstrated that individuals in the high‑expression group experienced significantly earlier composite adverse events compared with those in the low‑expression group (Log Rank P = 0.0369, Figure 1C). Univariable Cox regression was performed for all candidate variables, and variables with P < 0.05 were included in the multivariable Cox model. As shown in Table 3, univariable analysis revealed that miR‑136‑5p, NEU, PCT, BNP, D‑dimer, and CRP were significantly associated with composite adverse events. After multivariable adjustment, Z‑score‑standardized miR‑136‑5p, interpreted as a 1‑SD increase in expression, remained an independent prognostic predictor of composite adverse events among children with severe pneumonia (HR = 4.198, 95% CI = 1.657–10.640, P = 0.002) (Figure 1D).
Collectively, these findings suggest that miR‑136‑5p may have discriminatory potential for severe pneumonia compared with healthy controls. Elevated miR‑136‑5p expression is independently associated with composite adverse events, and its expression provides information relevant to clinical prognosis.
Effects of miR-136-5p on inflammation and cell viability of THP-1 macrophages
An in vitro inflammatory model was established using THP-1 macrophages, with inflammation induced by LPS. Compared with the control group, LPS stimulation significantly increased the concentrations of IL-6, TNF-α, and IL-1β in the cell supernatant (Figure 2A), confirming the successful establishment of the LPS-induced inflammatory model. To investigate the biological function of miR-136-5p in this inflammatory model, the efficacy of the cell transfection system was first verified. Results showed that miR-136-5p expression was already significantly upregulated in cells after LPS stimulation. Accordingly, transfection with miR-136-5p mimics increased its expression, whereas inhibitors decreased it (Figure 2B), confirming the efficacy and reliability of the transfection system.
Cell viability assay results demonstrated that LPS stimulation significantly reduced the viability of THP-1 cells. Transfection with miR-136-5p mimics further reduced cell viability, while transfection with inhibitors partially reversed the LPS-induced decrease in cell viability (Figure 2C), suggesting that miR-136-5p can enhance the LPS-associated reduction in cell viability on THP-1 cells. Further detection of inflammatory factor secretion levels confirmed that compared with the LPS group, transfection with miR-136-5p mimics significantly promoted the release of TNF-α, IL-6, and IL-1β from LPS-treated cells. In contrast, transfection with miR-136-5p inhibitors effectively suppressed the secretion of these pro-inflammatory factors (Figure 2D–F). These results indicate that miR-136-5p exerts a pro-inflammatory regulatory role in the LPS-induced inflammatory response of THP-1 macrophages.
miR-136-5p regulates cell viability via targeted inhibition of AGO1 in THP-1 cells
To systematically identify the downstream targets of miR-136-5p, potential target genes were predicted using four independent databases. The intersection of these predictions yielded 39 common candidate targets (Figure 3A). GO enrichment analysis was performed on these in silico predicted candidate genes to characterize their putative functional distribution (Figure 3B). Core functional terms closely associated with inflammatory regulation and miRNA action mechanisms served as the functional criteria for subsequent candidate prioritization: kinase binding, ubiquitin‑like protein ligase activity, mRNA metabolic processes, and intracellular receptor signaling pathway.
Based on the functional annotation combined with TargetScan-assessed targeting reliability, AGO1 was selected as the priority target for experimental validation. In silico prediction from TargetScan revealed that the 3’ untranslated region of AGO1 harbors two evolutionarily conserved binding sites for miR‑136‑5p, including one canonical 8mer site and one 7mer‑m8 site, with a total context++ score of −0.10, indicating a specific and evolutionarily conserved targeting interaction between the two. A dual-luciferase reporter assay confirmed that miR-136-5p directly binds to the 3'-UTR of AGO1, as miR-136-5p mimics significantly reduced luciferase activity of the wild-type AGO1 reporter but not the mutant construct (Figure 4A). The regulatory relationship was further validated at the cellular level. Detection of target gene expression revealed that AGO1 expression in THP-1 cells was significantly decreased following LPS stimulation. Transfection with miR-136-5p mimics further downregulated AGO1 expression, whereas transfection with miR-136-5p inhibitors reversed this effect (Figure 4B). These results indicate that miR-136-5p can target and suppress AGO1 expression.
To clarify the clinical relevance of this regulatory axis, analysis of clinical samples demonstrated that AGO1 expression was significantly lower in children with severe pneumonia compared to healthy controls, with a Hodges‑Lehmann median difference of −0.390 (95% CI: −0.418 to −0.361, P < 0.001) (Figure 4C). Moreover, AGO1 expression showed a significant negative correlation with miR-136-5p levels (r = -0.864, P < 0.001; Figure 4D), providing further evidence for the targeted regulation of AGO1 by miR-136-5p in the clinical context. To confirm that miR-136-5p exerts its biological function by targeting AGO1, a rescue experiment was conducted. Results showed that cell viability was significantly reduced in the LPS + miR-136-5p mimic group, while overexpression of AGO1 effectively reversed this cell viability effect (Figure 4E). In line with these findings, the secretion levels of pro-inflammatory factors (TNF-α, IL-6, and IL-1β) were markedly elevated in the LPS + miR-136-5p mimic group, and this pro-inflammatory phenotype was significantly attenuated by AGO1 overexpression (Figures 4F–H). We next examined whether silencing endogenous AGO1 could fully recapitulate the cellular phenotypes induced by miR-136-5p overexpression. Under LPS challenge, AGO1 knockdown significantly decreased THP-1 cell viability and increased the production of TNF-α, IL-6, and IL-1β compared with the si-NC control (Figures 4I–L). Importantly, when comparing the two independent transfection systems, the magnitude of cell viability and inflammatory response triggered by AGO1 silencing was consistently less pronounced than that observed in miR-136-5p-mimic-transfected cells, suggesting that AGO1 downregulation partially, but not entirely, accounts for the pro-inflammatory effects of miR-136-5p.
Collectively, the results of this investigation indicate that miR‑136‑5p exhibits discriminatory capacity to distinguish children with severe pneumonia from healthy controls, and its expression is informative for clinical prognosis. Its upregulated expression serves as an independent prognostic predictor of composite adverse events in children with severe pneumonia. Mechanistically, the findings suggest that miR‑136‑5p promotes inflammatory responses and reduced cell viability in LPS‑challenged THP‑1‑derived macrophages at least partially by targeting AGO1. These findings highlight the miR‑136‑5p/AGO1 regulatory pathway as one candidate for further therapeutic‑oriented investigation.
DATA AVAILABILITY:
All raw data derived from clinical samples and cell experiments supporting the findings of this study are included within Supplementary File 1. The public databases utilized in this work are available at the following URLs: TargetScanHuman v8.0: https://www.targetscan.org/vert_80/; miRDB: https://mirdb.org/index.html; miRWalk: http://mirwalk.umm.uni‑heidelberg.de/; ENCORI: https://rnasysu.com/encori/index.php; Metascape: https://metascape.org/gp/#/main/step1.

Figure 1: Expression, diagnostic efficacy, and prognostic value of miR-136-5p in childhood severe pneumonia. (A) Relative expression levels of miR-136-5p in peripheral samples from healthy individuals (n = 106) and children with severe pneumonia (n = 114). Two‑tailed Mann‑Whitney U test, ***P < 0.001. (B) ROC curve evaluating the discriminatory performance of miR-136-5p in children with severe pneumonia compared with healthy controls. The embedded table below the curve summarizes key diagnostic indices, including optimal cut‑off value, AUC, Youden index, sensitivity, specificity, their corresponding 95% Wilson‑score confidence intervals, and likelihood ratios. (C) Kaplan‑Meier time‑to‑event analysis for composite adverse events stratified by miR‑136‑5p expression levels (low-expression group, n = 58; high-expression group, n = 56). Event‑free probability denotes the probability of remaining free from composite adverse events (clinical deterioration or death). (D) Forest plot of multivariate Cox proportional‑hazards regression for composite adverse events. Hazard ratios (HRs) and corresponding 95% confidence intervals are displayed. An HR > 1 indicates an increased risk of composite adverse events. A P value < 0.05 was considered statistically significant. Please click here to view a larger version of this figure.

Figure 2: Regulatory effects of miR-136-5p on inflammation and cell viability in LPS-induced THP-1 macrophages. (A) Concentrations of IL-6, TNF-α, and IL-1β in the supernatant of THP-1 macrophages from the control group and LPS-stimulated group. (B) Relative expression levels of miR-136-5p in different treatment groups measured by qRT-PCR. (C) Viability of THP-1 macrophages in different treatment groups assessed by CCK-8 assay. (D–F) Concentrations of TNF-α, IL-6, and IL-1β in the supernatant of THP-1 macrophages across different treatment groups were measured by ELISA. Data are presented as mean ± SD. An independent-samples t-test was used for two-group comparisons, whereas a one‑way ANOVA followed by Tukey’s post‑hoc test was adopted for multi‑group comparisons. n ≥ 3; ***P < 0.001, **P < 0.01, *P < 0.05. Please click here to view a larger version of this figure.

Figure 3: Validation analysis of miR-136-5p targeting AGO1. (A) Venn diagram illustrating the overlap of miR-136-5p target genes predicted by the StarBase, miRDB, miRWalk, and TargetScan databases. (B) Bubble plot depicting GO enrichment analysis of the 39 target genes. Please click here to view a larger version of this figure.

Figure 4: Functional analysis of miR-136-5p targeting AGO1. (A) Validation of the binding interaction between miR-136-5p and AGO1 using a dual-luciferase reporter assay. (B) Relative expression levels of AGO1 in different treatment groups measured by qRT-PCR. (C) Relative expression levels of AGO1 in healthy individuals and children with severe pneumonia, determined by qRT-PCR. (D) Correlation scatter plot of miR-136-5p and AGO1 expression in children with severe pneumonia. (E) Functional rescue experiment assessing the effect of AGO1 overexpression on miR-136-5p-induced changes in THP-1 cell viability. (F–H) Functional rescue experiments assessing the effects of AGO1 overexpression on miR‑136‑5p‑induced alterations in the secretion levels of inflammatory factors (TNF‑α, IL‑6, and IL‑1β) in THP‑1 cell culture supernatants, as determined by ELISA. (I) Experiment assessing the effect of AGO1 knockdown on THP‑1 cell viability under LPS stimulation. (J–L) Experiments assessing the effects of AGO1 knockdown on the secretion levels of inflammatory factors (TNF‑α, IL‑6, and IL‑1β) in LPS‑treated THP‑1 cell culture supernatants, as determined by ELISA. Data are presented as mean ± SD; n ≥ 3; ***P < 0.001, *P < 0.05, ns, no significant difference. Two‑way ANOVA was used for panel A; two‑tailed non‑parametric Mann‑Whitney U test for panel C; one‑way ANOVA followed by Tukey’s post‑hoc test for panels B, E–L. Spearman’s correlation analysis was performed for panel D. Please click here to view a larger version of this figure.
Table 1: Comparison of clinical data between children with severe pneumonia and healthy individuals. BMI, body mass index; WBC, white blood cell; NEU, neutrophil count; LYM, lymphocyte count; PCT, procalcitonin; CRP, C-reactive protein; BNP, brain natriuretic peptide. Categorical variables (gender, age) were compared using the χ2 test. Continuous variables were presented as median (25th, 75th percentile) and compared by the Mann‑Whitney U test. The Hodges‑Lehmann median difference with 95% confidence interval (CI) was reported only for variables with statistically significant between‑group differences (P < 0.05). Please click here to download this file.
Table 2: Correlation between miR-136-5p and clinical characteristics. Normality was tested using the Kolmogorov–Smirnov (K‑S) test. BNP and D‑dimer followed a normal distribution and were analyzed by Pearson correlation. All other variables were analyzed by Spearman's rank correlation. Separate columns present the lower and upper bounds of 95% confidence intervals for correlation coefficients. Abbreviations: BMI = body mass index; WBC = white blood cell; NEU = neutrophil count; LYM = lymphocyte; PCT = procalcitonin; BNP = brain natriuretic peptide; D-D = D-dimer; CRP = C-reactive protein. Please click here to download this file.
Table 3: Univariable and multivariable Cox proportional‑hazards regression analyses for composite adverse events. HR, hazard ratio; CI, confidence interval. Event: composite adverse events defined as clinical deterioration or death within the 28‑day follow‑up window. miR-136-5p expression was Z-score standardized and modeled as a continuous variable; HR corresponds to the risk change per 1-standard-deviation increase in miR-136-5p expression. Variables with P < 0.05 in univariable Cox regression were entered into the multivariable Cox model. The proportional-hazards assumption was verified using the time-dependent covariate test. Total participants n = 114, total events = 39, censored observations = 75. A total of 9 candidate variables were subjected to univariable Cox regression screening, and 6 variables were included in the final multivariable model, yielding an events-per-variable (EPV) ratio of 6.5. Please click here to download this file.
Supplementary File 1: Raw data supporting the clinical and experimental analyses. Please click here to download this file.