Neutrophil dysregulation differentiates pediatric septic shock biomarker-based mortality-risk strata: insights from weighted gene co-expression network and transcriptomic analyses.

Background: Pediatric sepsis is a leading cause of global mortality, particularly among children, with limited therapeutic options beyond antibiotics and organ support. The Pediatric Sepsis Biomarker Risk Model (PERSEVERE-II) stratifies mortality risk in pediatric septic shock, yet the molecular mechanisms underlying high mortality risk remain incompletely understood. Methods: We analyzed whole blood transcriptomes collected from 81 children with septic shock on day 1 of meeting study criteria. Patients were stratified into high- and low-mortality risk groups according to the PERSEVERE-II biomarker risk model. Using weighted gene co-expression network analysis (WGCNA) and differential gene expression analyses, we identified molecular pathways and transcription factors (TFs) associated with mortality risk. Cell type differences were inferred using CIBERSORTx and using a reference single-cell dataset inclusive of neutrophils and their subsets. Findings: We identified distinct molecular profiles with high-risk patients displaying significant overexpression of genes related to neutrophil degranulation and innate immunity, alongside suppressed adaptive immune responses. The predominance of developing neutrophils underscored a major role of emergency granulopoiesis. Key TFs identified, including Interpretation: Our results suggest that neutrophil dysregulation underpins the high mortality risk conferred by the PERSEVERE-II model. The identified transcriptional regulators may provide potential targets to mitigate neutrophil dysregulation and improve outcomes among high-risk patients.