A seven-gene prognosis model to predict biochemical recurrence for prostate cancer based on the TCGA database

Analysis of risk factors affecting the prognosis of angiosarcoma patients: a retrospective study

This study aimed to identify prognostic factors influencing the survival of angiosarcoma patients and to explore the relationship between peripheral blood indicators and patient prognosis. A retrospective analysis was conducted on the clinical data collected from 105 angiosarcoma patients treated at China-Japan Union Hospital of Jilin University from January 2004 to April 2019, with an additional 50 patients included as external validation cohort. The median survival time for the study cohort was 1395 days, with 66.7% of patients (n=70) dying during the follow-up period. Significant differences were observed between the survival and death groups in age (P=0.022), primary tumor site (P=0.013), tumor size (P=0.008), and metastasis (P=0.018). Analysis of peripheral blood indicators showed that white blood cell (WBC) (P=0.006), platelet (PLT) (P=0.019), platelet-to-lymphocyte ratio (PLR) (P<0.001), and systemic immune-inflammation index (SII) (P=0.036) were significantly lower in the survival group, while lymphocyte (LYM) (P<0.001), albumin (ALB) (P<0.001), and prognostic nutritional index (PIN) (P<0.001) were significantly higher in the survival group. Multivariate Cox regression analysis identified SII (P=0.049, HR=0.551, 95% CI: 0.304-0.998), primary tumor site (P=0.001, HR=0.405, 95% CI: 0.235-0.699), metastasis (P=0.029, HR=1.864, 95% CI: 1.066-3.26), and chemotherapy (P=0.004, HR=0.434, 95% CI: 0.245-0.768) as independent prognostic factors affecting patients' 5-year survival. A nomogram model constructed based on these factors demonstrated high accuracy and stability in predicting 1-year, 3-year, and 5-year survival rates, with area under the curve (AUC) values of 0.836, 0.837, and 0.803, respectively, as validated by calibration curves and receiver operating characteristic (ROC) analysis. External validation further confirmed the model's reliability. Additionally, significant interactions were found between SII and primary tumor site (P=0.005) as well as chemotherapy (P=0.045). In conclusion, SII, primary tumor site, metastasis, and chemotherapy are crucial prognostic factors for angiosarcoma, and the developed nomogram provides a reliable tool for predicting survival outcomes.

A glycolysis-related signature to improve the current treatment and prognostic evaluation for breast cancer

Background

As a heterogeneous malignancy, breast cancer (BRCA) shows high incidence and mortality. Discovering novel molecular markers and developing reliable prognostic models may improve the survival of BCRA.

Methods

The RNA-seq data of BRCA patients were collected from the training set The Cancer Genome Atlas (TCGA)-BRCA and validation set GSE20685 in the Gene Expression Omnibus (GEO) databases. The "GSVA" R package was used to calculate the glycolysis score for each patient, based on which all the patients were divided into different glycolysis groups. The "limma" package was employed to perform differentially expression genes (DEGs) analysis. Key signature genes were selected by performing un/multivariate and least absolute shrinkage and selection operator (LASSO) C regression and used to develop a RiskScore model. The ESTIMATE and MCP-Counter algorithms were used for quantifying immune infiltration level. The functions of the genes were validated using Western blot, colony formation, transwell and wound-healing assay.

Results

The glycolysis score and prognostic analysis showed that high glycolysis score was related to tumorigenesis pathway and a poor prognosis in BRCA as overactive glycolysis inhibited the normal functions of immune cells. Subsequently, we screened five key prognostic genes using the LASSO Cox regression analysis and used them to establish a RiskScore with a high classification efficiency. Based on the results of the RiskScore, it was found that patients in the high-risk group had significantly unfavorable immune infiltration and prognostic outcomes. A nomogram integrating the RiskScore could well predict the prognosis for BRCA patients. Knockdown of PSCA suppressed cell proliferation, invasion and migration of BRCA cells.

Conclusion

This study developed a glycolysis-related signature with five genes to distinguish between high-risk and low-risk BRCA patients. A nomogram developed on the basis of the RiskScore was reliable to predict BRCA survival. Our model provided clinical guidance for the treatment of BRCA patients.

Glycolysis‑related lncRNA may be associated with prognosis and immune activity in grade II‑III glioma

Glucose metabolism, as a novel theory to explain tumor cell behavior, has been intensively studied in various tumors. The present study explored the long non-coding RNAs (lncRNAs) related to glycolysis in grade II-III glioma, aiming to provide a promising target for further research. Pearson correlation analysis was used to identify glycolysis-related lncRNAs. Univariate/multivariate Cox regression analysis and the Least Absolute Shrinkage and Selection Operator algorithm were applied to identify glycolysis-related lncRNAs to construct a prognosis prediction model. Subsequently, multi-dimensional evaluations were used to verify whether the risk model could predict the prognosis and survival rate of patients with grade II-III glioma. Finally, it was verified by functional experiments. The present study finally identified seven glycolysis-related lncRNAs (CRNDE, AC022034.1, RHOQ-AS1, AL159169.2, AL133215.2, AC007098.1 and LINC02587) to construct a prognosis prediction model. The present study further investigated the underlying immune microenvironment, somatic landscape and functional enrichment pathways. Additionally, individualized immunotherapeutic strategies and candidate compounds were identified to guide clinical treatment. The experimental results demonstrated that CRNDE could increase the proliferation of SHG-44 cells. In conclusion, a large sample of human grade II-III glioma in The Cancer Genome Atlas database was used to construct a risk model using glycolysis-related lncRNAs to predict the prognosis of patients with grade II-III glioma.

Comprehensive analyses and experimental verification of NETs and an EMT gene signature for prognostic prediction, immunotherapy, and chemotherapy in pancreatic adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is an aggressive malignancy with high mortality and poor prognosis. Neutrophil extracellular traps (NETs) and the epithelial-mesenchymal transition (EMT) significantly influence on the progression of various cancers. However, the underlying relevance of NETs- and EMT-associated genes on the outcomes of patients with PAAD remains to be elucidated. Transcriptome RNA sequencing data, together with clinical information and single-cell sequencing data of PAAD were collected from public databases. In the TCGA-PAAD cohort, ssGSEA was used to calculate NET and EMT scores. WGCNA was used to determine the key gene modules. A risk model with eight NET- and EMT-related genes (NERGs) was established using LASSO and multivariate Cox regression analysis. Patients in the reduced risk (RR) group showed better prognostic values compared with those in the elevated risk (ER) group. The prognostic model exhibited reliable and robust prediction when validated using an external database. The distributions of risk genes were explored in a single-cell sequencing data set. Immune infiltration, immune cycle, and immune checkpoints were compared between the RR and ER groups. Moreover, potential chemotherapeutic drugs were examined. DCBLD2 was identified as a key gene in PAAD cell lines by qRT-PCR, and was highly expressed in PAAD tissues. GSEA demonstrated that DCBLD2 induced the EMT. Transwell assays and western blotting showed that cell invasion and EMT induction were significantly reduced after DCBLD2 knockdown. Collectively, we constructed a prognosis model based on a NET and EMT gene signature, providing a valuable perspective for the prognostic evaluation and management of PAAD patient.

Prognostic prediction using a gene signature developed based on exhausted T cells for liver cancer patients

Background

Liver hepatocellular carcinoma (LIHC) is a solid primary malignancy with poor prognosis. This study discovered key prognostic genes based on T cell exhaustion and used them to develop a prognostic prediction model for LIHC.

Methods

SingleR's annotations combined with Seurat was used to automatically annotate the single-cell clustering results of the LIHC dataset GSE166635 downloaded from the Gene Expression Omnibus (GEO) database and to identify clusters related to exhausted T cells. Patients were classified using ConsensusClusterPlus package. Next, weighted gene co-expression network analysis (WGCNA) package was employed to distinguish key gene module, based on which least absolute shrinkage and selection operator (Lasso) and multi/univariate cox analysis were performed to construct a RiskScore system. Kaplan-Meier (KM) analysis and receiver operating characteristic curve (ROC) were employed to evaluate the efficacy of the model. To further optimize the risk model, a nomogram capable of predicting immune infiltration and immunotherapy sensitivity in different risk groups was developed. Expressions of genes were measured by quantitative real-time polymerase chain reaction (qRT-PCR), and immunofluorescence and Cell Counting Kit-8 (CCK-8) were performed for analyzing cell functions.

Results

We obtained 18,413 cells and clustered them into 7 immune and non-immune cell subpopulations. Based on highly variable genes among T cell exhaustion clusters, 3 molecular subtypes (C1, C2 and C3) of LIHC were defined, with C3 subtype showing the highest score of exhausted T cells and a poor prognosis. The Lasso and multivariate cox analysis selected 7 risk genes from the green module, which were closely associated with the C3 subtype. All the patients were divided into low- and high-risk groups based on the medium value of RiskScore, and we found that high-risk patients had higher immune infiltration and immune escape and poorer prognosis. The nomogram exhibited a strong performance for predicting long-term LIHC prognosis. In vitro experiments revealed that the 7 risk genes all had a higher expression in HCC cells, and that both liver HCC cell numbers and cell viability were reduced by knocking down MMP-9.

Conclusion

We developed a RiskScore model for predicting LIHC prognosis based on the scRNA-seq and RNA-seq data. The RiskScore as an independent prognostic factor could improve the clinical treatment for LIHC patients.