Prognostic stratification of sepsis through DNA damage response based RiskScore system: insights from single-cell RNA-sequencing and transcriptomic profiling

Background A novel risk scoring system, predicated on DNA damage response (DDR), was developed to enhance prognostic predictions and potentially inform the creation of more effective therapeutic protocols for sepsis. Methods To thoroughly delineate the expression profiles of DDR markers within the context of sepsis, an analytical approach utilizing single-cell RNA-sequencing (scRNA-seq) was implemented. Our study utilized single-cell analysis techniques alongside weighted gene co-expression network analysis (WGCNA) to pinpoint the genes that exhibit the most substantial associations with DNA damage response (DDR). Through Cox proportional hazards LASSO regression, we distinguished DDR-associated genes and established a risk model, enabling the stratification of patients into high- and low-risk groups. Subsequently, we carried out an analysis to determine our model’s predictive accuracy regarding patient survival. Moreover, we examined the distinct biological characteristics, various signal transduction routes, and immune system responses in sepsis patients, considering different risk categories and outcomes related to survival. Lastly, we conducted experimental validation of the identified genes through in vivo and in vitro assays, employing RT-PCR, ELISA, and flow cytometry. Results Both single-cell RNA sequencing (scRNA-seq) and bulk transcriptomic analyses have demonstrated a strong correlation between DNA damage response (DDR) levels and sepsis prognosis. Specific cell subtypes, including monocytes, megakaryocytes, CD4+ T cells, and neutrophils, have shown elevated DDR activity. Cells with increased DDR scores exhibited more robust and numerous interactions with other cell populations. The weighted gene co-expression network analysis (WGCNA) and single-cell analyses revealed 71 DDR-associated genes. We developed a four-gene risk scoring system using ARL4C, CD247, RPL7, and RPL31, identified through univariate COX, LASSO COX regression, and log-rank (Mantel-Cox) tests. Nomograms, calibration plots, and decision curve analyses (DCA) regarding these specific genes have provided significant clinical benefits for individuals diagnosed with sepsis. The study suggested that individuals categorized as lower-risk demonstrated enhanced infiltration of immune cells, upregulated expression of immune regulators, and a more prolific presence of immune-associated functionalities and pathways. RT-qPCR analyses on a sepsis rat model revealed differential gene expression predominantly in the four targeted genes. Furthermore, ARL4C knockdown in sepsis model in vivo and vitro caused increased inflammatory response and a worse prognosis. Conclusion The delineated DDR expression landscape offers insights into sepsis pathogenesis, whilst our riskScore model, based on a robust four-gene signature, could underpin personalized sepsis treatment strategies.