CD4+ conventional T cells-related genes signature is a prognostic indicator for ovarian cancer.

It is believed that ovarian cancer (OC) is the most deadly form of gynecological cancer despite its infrequent occurrence, which makes it one of the most salient public health concerns. Clinical and preclinical studies have revealed that intratumoral CD4+ T cells possess cytotoxic capabilities and were capable of directly killing cancer cells. This study aimed to identify the CD4+ conventional T cells-related genes (CD4TGs) with respect to the prognosis in OC. We obtained the transcriptome and clinical data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CD4TGs were first identified from single-cell datasets, then univariate Cox regression was used to screen prognosis-related genes, LASSO was conducted to remove genes with coefficient zero, and multivariate Cox regression was used to calculate riskscore and to construct the CD4TGs risk signature. Kaplan-Meier analysis, univariate Cox regression, multivariate Cox regression, time-dependent receiver operating characteristics (ROC), decision curve analysis (DCA), nomogram, and calibration were made to verify and evaluate the risk signature. Gene set enrichment analyses (GSEA) in risk groups were conducted to explore the tightly correlated pathways with the risk group. The role of riskscore has been further explored in the tumor microenvironment (TME), immunotherapy, and chemotherapy. A risk signature with 11 CD4TGs in OC was finally established in the TCGA database and furtherly validated in several GEO cohorts. High riskscore was significantly associated with a poorer prognosis and proven to be an independent prognostic biomarker by multivariate Cox regression. The 1-, 3-, and 5-year ROC values, DCA curve, nomogram, and calibration results confirmed the excellent prediction power of this model. Compared with the reported risk models, our model showed better performance. The patients were grouped into high-risk and low-risk subgroups according to the riskscore by the median value. The low-risk group patients tended to exhibit a higher immune infiltration, immune-related gene expression and were more sensitive to immunotherapy and chemotherapy. Collectively, our findings of the prognostic value of CD4TGs in prognosis and immune response, provided valuable insights into the molecular mechanisms and clinical management of OC.

Hua T ，Liu DX ，Zhang XC ，Li ST ，Yan P ，Zhao Q ，Chen SB ... -  《Frontiers in Immunology》

Establishment of an ovarian cancer exhausted CD8+T cells-related genes model by integrated analysis of scRNA-seq and bulk RNA-seq.

Ovarian cancer (OC) was the fifth leading cause of cancer death and the deadliest gynecological cancer in women. This was largely attributed to its late diagnosis, high therapeutic resistance, and a dearth of effective treatments. Clinical and preclinical studies have revealed that tumor-infiltrating CD8+T cells often lost their effector function, the dysfunctional state of CD8+T cells was known as exhaustion. Our objective was to identify genes associated with exhausted CD8+T cells (CD8TEXGs) and their prognostic significance in OC. We downloaded the RNA-seq and clinical data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CD8TEXGs were initially identified from single-cell RNA-seq (scRNA-seq) datasets, then univariate Cox regression, the least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression were utilized to calculate risk score and to develop the CD8TEXGs risk signature. Kaplan-Meier analysis, univariate Cox regression, multivariate Cox regression, time-dependent receiver operating characteristics (ROC), nomogram, and calibration were conducted to verify and evaluate the risk signature. Gene set enrichment analyses (GSEA) in the risk groups were used to figure out the closely correlated pathways with the risk group. The role of risk score has been further explored in the homologous recombination repair deficiency (HRD), BRAC1/2 gene mutations and tumor mutation burden (TMB). A risk signature with 4 CD8TEXGs in OC was finally built in the TCGA database and further validated in large GEO cohorts. The signature also demonstrated broad applicability across various types of cancer in the pan-cancer analysis. The high-risk score was significantly associated with a worse prognosis and the risk score was proven to be an independent prognostic biomarker. The 1-, 3-, and 5-years ROC values, nomogram, calibration, and comparison with the previously published models confirmed the excellent prediction power of this model. The low-risk group patients tended to exhibit a higher HRD score, BRCA1/2 gene mutation ratio and TMB. The low-risk group patients were more sensitive to Poly-ADP-ribose polymerase inhibitors (PARPi). Our findings of the prognostic value of CD8TEXGs in prognosis and drug response provided valuable insights into the molecular mechanisms and clinical management of OC.

Hua T ，Liu DX ，Zhang XC ，Li ST ，Wu JL ，Zhao Q ，Chen SB ... -  《-》

Identification and validation of a seven m6A-related lncRNAs signature predicting prognosis of ovarian cancer.

Long non-coding RNAs (lncRNAs) play an important role in angiogenesis, immune response, inflammatory response and tumor development and metastasis. m6 A (N6-methyladenosine) is one of the most common RNA modifications in eukaryotes. The aim of our research was to investigate the potential prognostic value of m6A-related lncRNAs in ovarian cancer (OC). The data we need for our research was downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. Pearson correlation analysis between 21 m6A regulators and lncRNAs was performed to identify m6A-related lncRNAs. Univariate Cox regression analysis was implemented to screen for lncRNAs with prognostic value. A least absolute shrinkage and selection operator (LASSO) Cox regression and multivariate Cox regression analyses was used to further reduct the lncRNAs with prognostic value and construct a m6A-related lncRNAs signature for predicting the prognosis of OC patients. Two hundred seventy-five m6A-related lncRNAs were obtained using pearson correlation analysis. 29 m6A-related lncRNAs with prognostic value was selected through univariate Cox regression analysis. Then, a seven m6A-related lncRNAs signature was identified by LASSO Cox regression. Each patient obtained a riskscore through multivariate Cox regression analyses and the patients were classified into high-and low-risk group using the median riskscore as a cutoff. Kaplan-Meier curve revealed that the patients in high-risk group have poor outcome. The receiver operating characteristic curve revealed that the predictive potential of the m6A-related lncRNAs signature for OC was powerful. The predictive potential of the m6A-related lncRNAs signature was successfully validated in the GSE9891, GSE26193 datasets and our clinical specimens. Multivariate analyses suggested that the m6A-related lncRNAs signature was an independent prognostic factor for OC patients. Moreover, a nomogram based on the expression level of the seven m6A-related lncRNAs was established to predict survival rate of patients with OC. Finally, a competing endogenous RNA (ceRNA) network associated with the seven m6A-related lncRNAs was constructed to understand the possible mechanisms of the m6A-related lncRNAs involed in the progression of OC. In conclusion, our research revealed that the m6A-related lncRNAs may affect the prognosis of OC patients and identified a seven m6A-related lncRNAs signature to predict the prognosis of OC patients.

Song Y ，Qu H 《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》

Construction of a novel mRNA-signature prediction model for prognosis of bladder cancer based on a statistical analysis.

Bladder cancer (BC) is a common malignancy neoplasm diagnosed in advanced stages in most cases. It is crucial to screen ideal biomarkers and construct a more accurate prognostic model than conventional clinical parameters. The aim of this research was to develop and validate an mRNA-based signature for predicting the prognosis of patients with bladder cancer. The RNA-seq data was downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were screened in three datasets, and prognostic genes were identified from the training set of TCGA dataset. The common genes between DEGs and prognostic genes were narrowed down to six genes via Least Absolute Shrinkage and Selection Operator (LASSO) regression, and stepwise multivariate Cox regression. Then the gene-based risk score was calculated via Cox coefficient. Time-dependent receiver operating characteristic (ROC) and Kaplan-Meier (KM) survival analysis were used to assess the prognostic power of risk score. Multivariate Cox regression analysis was applied to construct a nomogram. Decision curve analysis (DCA), calibration curves, and time-dependent ROC were performed to assess the nomogram. Finally, functional enrichment of candidate genes was conducted to explore the potential biological pathways of candidate genes. SORBS2, GPC2, SETBP1, FGF11, APOL1, and H1-2 were screened to be correlated with the prognosis of BC patients. A nomogram was constructed based on the risk score, pathological stage, and age. Then, the calibration plots for the 1-, 3-, 5-year OS were predicted well in entire TCGA-BLCA patients. Decision curve analysis (DCA) indicated that the clinical value of the nomogram was higher than the stage model and TNM model in predicting overall survival analysis. The time-dependent ROC curves indicated that the nomogram had higher predictive accuracy than the stage model and risk score model. The AUC of nomogram time-dependent ROC was 0.763, 0.805, and 0.806 for 1-year, 3-year, and 5-year, respectively. Functional enrichment analysis of candidate genes suggested several pathways and mechanisms related to cancer. In this research, we developed an mRNA-based signature that incorporated clinical prognostic parameters to predict BC patient prognosis well, which may provide a novel prognosis assessment tool for clinical practice and explore several potential novel biomarkers related to the prognosis of patients with BC.

Li J ，Cao J ，Li P ，Yao Z ，Deng R ，Ying L ，Tian J ... -  《BMC CANCER》