Viral infection, typically disregarded, has a significant role in burns. However, there is still a lack of biomarkers and immunotherapy targets related to viral infections in burns.
Virus-related genes (VRGs) that were extracted from Gene Oncology (GO) database were included as hallmarks. Through unsupervised consensus clustering, we divided patients into two VRGs molecular patterns (VRGMPs). Weighted gene co-expression network analysis (WGCNA) was performed to study the relationship between burns and VRGs. Random forest (RF), least absolute shrinkage and selection operator (LASSO) regression, and logistic regression were used to select key genes, which were utilized to construct prognostic signatures by multivariate logistic regression. The risk score of the nomogram defined high- and low-risk groups. We compared immune cells, immune checkpoint-related genes, and prognosis between the two groups. Finally, we used network analysis and molecular docking to predict drugs targeting
We established two VRGMPs, which differed in monocytes, neutrophils, dendritic cells, and T cells. In WGCNA, genes were divided into 14 modules, and the black module was correlated with VRGMPs. A total of 65 genes were selected by WGCNA, STRING, and differential expression analysis. The results of GO enrichment analysis were enriched in Th1 and Th2 cell differentiation, B cell receptor signaling pathway, alpha-beta T cell activation, and alpha-beta T cell differentiation. Then the 2-gene signature was constructed by RF, LASSO, and LOGISTIC regression. The signature was an independent prognostic factor and performed well in ROC, calibration, and decision curves. Further, the expression of immune cells and checkpoint genes differed between high- and low-risk groups.
This is the first VRG-based signature (including 2 key genes validated by qPCR) for predicting survival, and it could provide vital guidance to achieve optimized immunotherapy for immunosuppression in burns.