Preclinical analysis of novel prognostic transcription factors and immune-related gene signatures for bladder cancer via TCGA-based bioinformatic analysis

Gene Expression Monotonicity across Bladder Cancer Stages Informs on the Molecular Pathogenesis and Identifies a Prognostic Eight-Gene Signature

Despite advancements in molecular classification, tumor stage and grade still remain the most relevant prognosticators used by clinicians to decide on patient management. Here, we leverage publicly available data to characterize bladder cancer (BLCA)’s stage biology based on increased sample sizes, identify potential therapeutic targets, and extract putative biomarkers. A total of 1135 primary BLCA transcriptomes from 12 microarray studies were compiled in a meta-cohort and analyzed for monotonal alterations in pathway activities, gene expression, and co-expression patterns with increasing stage (Ta–T1–T2–T3–T4), starting from the non-malignant tumor-adjacent urothelium. The TCGA-2017 and IMvigor-210 RNA-Seq data were used to validate our findings. Wnt, MTORC1 signaling, and MYC activity were monotonically increased with increasing stage, while an opposite trend was detected for the catabolism of fatty acids, circadian clock genes, and the metabolism of heme. Co-expression network analysis highlighted stage- and cell-type-specific genes of potentially synergistic therapeutic value. An eight-gene signature, consisting of the genes AKAP7, ANLN, CBX7, CDC14B, ENO1, GTPBP4, MED19, and ZFP2, had independent prognostic value in both the discovery and validation sets. This novel eight-gene signature may increase the granularity of current risk-to-progression estimators.

A Transcription Factor-Based Risk Model for Predicting the Prognosis of Prostate Cancer and Potential Therapeutic Drugs

Background

Prostate cancer (PC) is one of the most critical cancers affecting men's health worldwide. The development of many cancers involves dysregulation or mutations in key transcription factors. This study established a transcription factor-based risk model to predict the prognosis of PC and potential therapeutic drugs.

Materials and Methods

In this study, RNA-sequencing data were downloaded and analyzed using The Cancer Genome Atlas dataset. A total of 145 genes related to the overall survival rate of PC patients were screened using the univariate Cox analysis. The Kdmist clustering method was used to classify prostate adenocarcinoma (PRAD), thereby determining the cluster related to the transcription factors. The support vector machine-recursive feature elimination method was used to identify genes related to the types of transcription factors and the key genes specifically upregulated or downregulated were screened. These genes were further analyzed using Lasso to establish a model. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used for the functional analysis. The TIMER algorithm was used to quantify the abundance of immune cells in PRAD samples. The chemotherapy response of each GBM patient was predicted based on the public pharmacogenomic database, Genomics of Drug Sensitivity in Cancer (GDSC, http://www.cancerrxgene.org). The R package "pRRophetic" was applied to drug sensitivity (IC50) value prediction.

Results

We screened 10 genes related to prognosis, including eight low-risk genes and two high-risk genes. The receiver operating characteristic (ROC) curve was 0.946. Patients in the high-risk score group had a poorer prognosis than those in the low-risk score group. The average area under the curve value of the model at different times was higher than 0.8. The risk score was an independent prognostic factor. Compared with the low-risk score group, early growth response-1 (EGR1), CACNA2D1, AC005831.1, SLC52A3, TMEM79, IL20RA, CRACR2A, and FAM189A2 expressions in the high-risk score group were decreased, while AC012181.1 and TRAPPC8 expressions were increased. GO and KEGG analyses showed that prognosis was related to various cancer signaling pathways. The proportion of B_cell, T_cell_CD4, and macrophages in the high-risk score group was significantly higher than that in the low-risk score group. A total of 25 classic immune checkpoint genes were screened out to express abnormally high-risk scores, and there were significant differences. Thirty mutant genes were identified; in the high- and low-risk score groups, SPOP, TP53, and TTN had the highest mutation frequency, and their mutations were mainly missense mutations. A total of 36 potential drug candidates for the treatment of PC were screened and identified.

Conclusions

Ten genes of both high-and low-risk scores were associated with the prognosis of PC. PC prognosis may be related to immune disorders. SPOP, TP53, and TTN may be potential targets for the prognosis of PC.

A novel immune-related gene pair prognostic signature for predicting overall survival in bladder cancer

Background

Bladder cancer (BC) is the ninth most common malignant tumor. We constructed a risk signature using immune-related gene pairs (IRGPs) to predict the prognosis of BC patients.

Methods

The mRNA transcriptome, simple nucleotide variation and clinical data of BC patients were downloaded from The Cancer Genome Atlas (TCGA) database (TCGA-BLCA). The mRNA transcriptome and clinical data were also extracted from Gene Expression Omnibus (GEO) datasets (GSE31684). A risk signature was built based on the IRGPs. The ability of the signature to predict prognosis was analyzed with survival curves and Cox regression. The relationships between immunological parameters [immune cell infiltration, immune checkpoints, tumor microenvironment (TME) and tumor mutation burden (TMB)] and the risk score were investigated. Finally, gene set enrichment analysis (GSEA) was used to explore molecular mechanisms underlying the risk score.

Results

The risk signature utilized 30 selected IRGPs. The prognosis of the high-risk group was significantly worse than that of the low-risk group. We used the GSE31684 dataset to validate the signature. Close relationships were found between the risk score and immunological parameters. Finally, GSEA showed that gene sets related to the extracellular matrix (ECM), stromal cells and epithelial-mesenchymal transition (EMT) were enriched in the high-risk group. In the low-risk group, we found a number of immune-related pathways in the enriched pathways and biofunctions.

Conclusions

We used a new tool, IRGPs, to build a risk signature to predict the prognosis of BC. By evaluating immune parameters and molecular mechanisms, we gained a better understanding of the mechanisms underlying the risk signature. This signature can also be used as a tool to predict the effect of immunotherapy in patients with BC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08486-0.