Development and validation a prognostic model based on natural killer T cells marker genes for predicting prognosis and characterizing immune status in glioblastoma through integrated analysis of single-cell and bulk RNA sequencing.

Glioblastoma (GBM) is an aggressive and unstoppable malignancy. Natural killer T (NKT) cells, characterized by specific markers, play pivotal roles in many tumor-associated pathophysiological processes. Therefore, investigating the functions and complex interactions of NKT cells is great interest for exploring GBM. We acquired a single-cell RNA-sequencing (scRNA-seq) dataset of GBM from Gene Expression Omnibus (GEO) database. The weighted correlation network analysis (WGCNA) was employed to further screen genes subpopulations. Subsequently, we integrated the GBM cohorts from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases to describe different subtypes by consensus clustering and developed a prognostic model by least absolute selection and shrinkage operator (LASSO) and multivariate Cox regression analysis. We further investigated differences in survival rates and clinical characteristics among different risk groups. Furthermore, a nomogram was developed by combining riskscore with the clinical characteristics. We investigated the abundance of immune cells in the tumor microenvironment (TME) by CIBERSORT and single sample gene set enrichment analysis (ssGSEA) algorithms. Immunotherapy efficacy assessment was done with the assistance of Tumor Immune Dysfunction and Exclusion (TIDE) and The Cancer Immunome Atlas (TCIA) databases. Real-time quantitative polymerase chain reaction (RT-qPCR) experiments and immunohistochemical profiles of tissues were utilized to validate model genes. We identified 945 NKT cells marker genes from scRNA-seq data. Through further screening, 107 genes were accurately identified, of which 15 were significantly correlated with prognosis. We distinguished GBM samples into two distinct subtypes and successfully developed a robust prognostic prediction model. Survival analysis indicated that high expression of NKT cell marker genes was significantly associated with poor prognosis in GBM patients. Riskscore can be used as an independent prognostic factor. The nomogram was demonstrated remarkable utility in aiding clinical decision making. Tumor immune microenvironment analysis revealed significant differences of immune infiltration characteristics between different risk groups. In addition, the expression levels of immune checkpoint-associated genes were consistently elevated in the high-risk group, suggesting more prominent immune escape but also a stronger response to immune checkpoint inhibitors. By integrating scRNA-seq and bulk RNA-seq data analysis, we successfully developed a prognostic prediction model that incorporates two pivotal NKT cells marker genes, namely, CD44 and TNFSF14. This model has exhibited outstanding performance in assessing the prognosis of GBM patients. Furthermore, we conducted a preliminary investigation into the immune microenvironment across various risk groups that contributes to uncover promising immunotherapeutic targets specific to GBM.

Hu J ，Xu L ，Fu W ，Sun Y ，Wang N ，Zhang J ，Yang C ，Zhang X ，Zhou Y ，Wang R ，Zhang H ，Mou R ，Du X ，Li X ，Hu S ，Xie R ... -  《-》

Comprehensive analysis of scRNA-Seq and bulk RNA-Seq reveals dynamic changes in the tumor immune microenvironment of bladder cancer and establishes a prognostic model.

The prognostic management of bladder cancer (BLCA) remains a great challenge for clinicians. Recently, bulk RNA-seq sequencing data have been used as a prognostic marker for many cancers but do not accurately detect core cellular and molecular functions in tumor cells. In the current study, bulk RNA-seq and single-cell RNA sequencing (scRNA-seq) data were combined to construct a prognostic model of BLCA. BLCA scRNA-seq data were downloaded from Gene Expression Omnibus (GEO) database. Bulk RNA-seq data were obtained from the UCSC Xena. The R package "Seurat" was used for scRNA-seq data processing, and the uniform manifold approximation and projection (UMAP) were utilized for downscaling and cluster identification. The FindAllMarkers function was used to identify marker genes for each cluster. The limma package was used to obtain differentially expressed genes (DEGs) affecting overall survival (OS) in BLCA patients. Weighted gene correlation network analysis (WGCNA) was used to identify BLCA key modules. The intersection of marker genes of core cells and genes of BLCA key modules and DEGs was used to construct a prognostic model by univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO) analyses. Differences in clinicopathological characteristics, immune microenvironment, immune checkpoints, and chemotherapeutic drug sensitivity between the high and low-risk groups were also investigated. scRNA-seq data were analyzed to identify 19 cell subpopulations and 7 core cell types. The ssGSEA showed that all 7 core cell types were significantly downregulated in tumor samples of BLCA. We identified 474 marker genes from the scRNA-seq dataset, 1556 DEGs from the Bulk RNA-seq dataset, and 2334 genes associated with a key module identified by WGCNA. After performing intersection, univariate Cox, and LASSO analysis, we obtained a prognostic model based on the expression levels of 3 signature genes, namely MAP1B, PCOLCE2, and ELN. The feasibility of the model was validated by an internal training set and two external validation sets. Moreover, patients with high-risk scores are predisposed to experience poor OS, a larger prevalence of stage III-IV, a greater TMB, a higher infiltration of immune cells, and a lesser likelihood of responding favorably to immunotherapy. By integrating scRNA-seq and bulk RNA-seq data, we constructed a novel prognostic model to predict the survival of BLCA patients. The risk score is a promising independent prognostic factor that is closely correlated with the immune microenvironment and clinicopathological characteristics.

Tan Z ，Chen X ，Zuo J ，Fu S ，Wang H ，Wang J ... -  《Journal of Translational Medicine》

Expression of hub genes of endothelial cells in glioblastoma-A prognostic model for GBM patients integrating single-cell RNA sequencing and bulk RNA sequencing.

This study aimed to use single-cell RNA-seq (scRNA-seq) to discover marker genes in endothelial cells (ECs) and construct a prognostic model for glioblastoma multiforme (GBM) patients in combination with traditional high-throughput RNA sequencing (bulk RNA-seq). Bulk RNA-seq data was downloaded from The Cancer Genome Atlas (TCGA) and The China Glioma Genome Atlas (CGGA) databases. 10x scRNA-seq data for GBM were obtained from the Gene Expression Omnibus (GEO) database. The uniform manifold approximation and projection (UMAP) were used for downscaling and cluster identification. Key modules and differentially expressed genes (DEGs) were identified by weighted gene correlation network analysis (WGCNA). A non-negative matrix decomposition (NMF) algorithm was used to identify the different subtypes based on DEGs, and multivariate cox regression analysis to model the prognosis. Finally, differences in mutational landscape, immune cell abundance, immune checkpoint inhibitors (ICIs)-associated genes, immunotherapy effects, and enriched pathways were investigated between different risk groups. The analysis of scRNA-seq data from eight samples revealed 13 clusters and four cell types. After applying Fisher's exact test, ECs were identified as the most important cell type. The NMF algorithm identified two clusters with different prognostic and immunological features based on DEGs. We finally built a prognostic model based on the expression levels of four key genes. Higher risk scores were significantly associated with poorer survival outcomes, low mutation rates in IDH genes, and upregulation of immune checkpoints such as PD-L1 and CD276. We built and validated a 4-gene signature for GBM using 10 scRNA-seq and bulk RNA-seq data in this work.

Zhao S ，Ji W ，Shen Y ，Fan Y ，Huang H ，Huang J ，Lai G ，Yuan K ，Cheng C ... -  《BMC CANCER》

Establishment of an ovarian cancer omentum metastasis-related prognostic model by integrated analysis of scRNA-seq and bulk RNA-seq.

Ovarian cancer has the highest mortality rate among gynecological malignant tumors, and it preferentially metastasizes to omental tissue, leading to intestinal obstruction and death. scRNA-seq is a powerful technique to reveal tumor heterogeneity. Analyzing omentum metastasis of ovarian cancer at the single-cell level may be more conducive to exploring and understanding omentum metastasis and prognosis of ovarian cancer at the cellular function and genetic levels. The omentum metastasis site scRNA-seq data of GSE147082 were acquired from the GEO (Gene Expression Omnibus) database, and single cells were clustered by the Seruat package and annotated by the SingleR package. Cell differentiation trajectories were reconstructed through the monocle package. The ovarian cancer microarray data of GSE132342 were downloaded from GEO and were clustered by using the ConsensusClusterPlus package into omentum metastasis-associated clusters according to the marker genes gained from single-cell differentiation trajectory analysis. The tumor microenvironment (TME) and immune infiltration differences between clusters were analyzed by the estimate and CIBERSORT packages. The expression matrix of genes used to cluster GSE132342 patients was extracted from bulk RNA-seq data of TCGA-OV (The Cancer Genome Atlas ovarian cancer), and least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were performed to establish an omentum metastasis-associated gene (OMAG) signature. The signature was then tested by GSE132342 data. Finally, the clinicopathological characteristics of TCGA-OV were screened by univariate and multivariate Cox regression analysis to draw the nomogram. A total of 9885 cells from 6 patients were clustered into 18 cell clusters and annotated into 14 cell types. Reconstruction of differentiation trajectories divided the cells into 5 branches, and a total of 781 cell trajectory-related characteristic genes were obtained. A total of 3769 patients in GSE132342 were subtyped into 3 clusters by 74 cell trajectory-related characteristic genes. Kaplan-Meier (K-M) survival analysis showed that the prognosis of cluster 2 was the worst, P < 0.001. The TME analysis showed that the ESTIMATE score and stromal score in cluster 2 were significantly higher than those in the other two clusters, P < 0.001. The immune infiltration analysis showed differences in the fraction of 8 immune cells among the 3 clusters, P < 0.05. The expression data of 74 genes used for GEO clustering were extracted from 379 patients in TCGA-OV, and combined with survival information, 10 candidates for OMAGs were filtered by LASSO. By using multivariate Cox regression, the 6-OMAGs signature was established as RiskScore = 0.307*TIMP3 + 3.516*FBN1-0.109*IGKC + 0.209*RPL21 + 0.870*UCHL1 + 0.365*RARRES1. Taking TCGA-OV as the training set and GSE132342 as the test set, receiver operating characteristic (ROC) curves were drawn to verify the prognostic value of 6-OMAGs. Screened by univariate and multivariate Cox regression analysis, 3 (age, cancer status, primary therapy outcome) of 5 clinicopathological characteristics were used to construct the nomogram combined with risk score. We constructed an ovarian cancer prognostic model related to omentum metastasis composed of 6-OMAGs and 3 clinicopathological features and analyzed the potential mechanism of these 6-OMAGs in ovarian cancer omental metastasis.

Zhang D ，Lu W ，Cui S ，Mei H ，Wu X ，Zhuo Z ... -  《Journal of Ovarian Research》

Systematic identification, development, and validation of prognostic biomarkers involving the tumor-immune microenvironment for glioblastoma.

Gliomas are infiltrative neoplasms with a highly invasive nature. Due to its distinct genomic, genetic and epigenetic features, the immune prognostic signature (IPS) and immune microenvironment of glioblastoma (GBM) merit further research. We aimed to explore prognosis-related immune genes and develop an IPS model for predicting prognosis in GBM. RNA-sequencing data, as well as clinical information, from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) public cohorts were analyzed. To develop the IPS, least absolute shrinkage and selection operator (LASSO) Cox analysis was performed for immune-related genes that were differentially expressed between GBM and normal tissues. Then, interaction effects of the IPS on the immune microenvironment were systematically analyzed; the precise prognostic model was developed based on the IPS and clinical data and was then further validated. A total of 21 immune prognostic genes were identified based on GBM microenvironment status. An 8-gene IPS was established, and the GBM patients were effectively stratified into low- and high-risk groups in the TCGA cohort as a training set. Univariate and multivariate Cox analyses revealed that IPS was an independent prognostic factor, and the prognostic performance of individual IPS genes was systematically illustrated. In addition, a comprehensive and novel nomogram model was initially established to estimate overall survival in TCGA-GBM patients, and high-risk patients had higher levels of dendritic cell and neutrophil infiltration. Furthermore, the nomogram model was developed and validated in the CGGA validation set. The low-risk IPS was linked to a stronger response to anti-PD-L1 immunotherapy and clinical advantages in the IMvigor210 cohort. This novel IPS with promising biomarkers classifies GBM patients into subgroups with distinct clinical outcomes and immunophenotypes. Our findings and this resource may help to characterize the immune microenvironment, inform cancer immunotherapy and facilitate the development of precision immuno-oncology.

Zhao B ，Wang Y ，Wang Y ，Chen W ，Liu PH ，Kong Z ，Dai C ，Wang Y ，Ma W ... -  《-》