Tumor Expression Profile Analysis Developed and Validated a Prognostic Model Based on Immune-Related Genes in Bladder Cancer.

Background: Bladder cancer (BLCA) ranks 10th in incidence among malignant tumors and 6th in incidence among malignant tumors in males. With the application of immune therapy, the overall survival (OS) rate of BLCA patients has greatly improved, but the 5-year survival rate of BLCA patients is still low. Furthermore, not every BLCA patient benefits from immunotherapy, and there are a limited number of biomarkers for predicting the immunotherapy response. Therefore, novel biomarkers for predicting the immunotherapy response and prognosis of BLCA are urgently needed. Methods: The RNA sequencing (RNA-seq) data, clinical data and gene annotation files for The Cancer Genome Atlas (TCGA) BLCA cohort were extracted from the University of California, Santa Cruz (UCSC) Xena Browser. The BLCA datasets GSE31684 and GSE32894 from the Gene Expression Omnibus (GEO) database were extracted for external validation. Immune-related genes were extracted from InnateDB. Significant differentially expressed genes (DEGs) were identified using the R package "limma," and Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for the DEGs were performed using R package "clusterProfiler." Least absolute shrinkage and selection operator (LASSO) regression analysis were used to construct the signature model. The infiltration level of each immune cell type was estimated using the single-sample gene set enrichment analysis (ssGSEA) algorithm. The performance of the model was evaluated with receiver operating characteristic (ROC) curves and calibration curves. Results: In total, 1,040 immune-related DEGs were identified, and eight signature genes were selected to construct a model using LASSO regression analysis. The risk score of BLCA patients based on the signature model was negatively correlated with OS and the immunotherapy response. The ROC curve for OS revealed that the model had good accuracy. The calibration curve showed good agreement between the predictions and actual observations. Conclusions: Herein, we constructed an immune-related eight-gene signature that could be a potential biomarker to predict the immunotherapy response and prognosis of BLCA patients.

Dong B ，Liang J ，Li D ，Song W ，Zhao S ，Ma Y ，Song J ，Zhu M ，Yang T ... -  《Frontiers in Genetics》

Identification of a Prognostic Signature Associated With the Homeobox Gene Family for Bladder Cancer.

Background: Bladder cancer (BLCA) is a common malignant tumor of the genitourinary system, and there is a lack of specific, reliable, and non-invasive tumor biomarker tests for diagnosis and prognosis evaluation. Homeobox genes play a vital role in BLCA tumorigenesis and development, but few studies have focused on the prognostic value of homeobox genes in BLCA. In this study, we aim to develop a prognostic signature associated with the homeobox gene family for BLCA. Methods: The RNA sequencing data, clinical data, and probe annotation files of BLCA patients were downloaded from the Gene Expression Omnibus database and the University of California, Santa Cruz (UCSC), Xena Browser. First, differentially expressed homeobox gene screening between tumor and normal samples was performed using the "limma" and robust rank aggregation (RRA) methods. The mutation data were obtained with the "TCGAmutation" package and visualized with the "maftools" package. Kaplan-Meier curves were plotted with the "survminer" package. Then, a signature was constructed by logistic regression analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using "clusterProfiler." Furthermore, the infiltration level of each immune cell type was estimated using the single-sample gene set enrichment analysis (ssGSEA) algorithm. Finally, the performance of the signature was evaluated by receiver-operating characteristic (ROC) curve and calibration curve analyses. Results: Six genes were selected to construct this prognostic model: TSHZ3, ZFHX4, ZEB2, MEIS1, ISL1, and HOXC4. We divided the BLCA cohort into high- and low-risk groups based on the median risk score calculated with the novel signature. The overall survival (OS) rate of the high-risk group was significantly lower than that of the low-risk group. The infiltration levels of almost all immune cells were significantly higher in the high-risk group than in the low-risk group. The average risk score for the group that responded to immunotherapy was significantly lower than that of the group that did not. Conclusion: We constructed a risk prediction signature with six homeobox genes, which showed good accuracy and consistency in predicting the patient's prognosis and response to immunotherapy. Therefore, this signature can be a potential biomarker and treatment target for BLCA patients.

Dong B ，Liang J ，Li D ，Song W ，Song J ，Zhu M ，Zhao S ，Ma Y ，Yang T ... -  《-》

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》

Development and validation of a model based on immunogenic cell death related genes to predict the prognosis and immune response to bladder urothelial carcinoma.

Immunogenic cell death (ICD) has been categorized as a variant of regulated cell death that is capable of inducing an adaptive immune response. A growing body of evidence has indicated that ICD can modify the tumor immune microenvironment by releasing danger signals or damage-associated molecular patterns (DAMPs), potentially enhancing the efficacy of immunotherapy. Consequently, the identification of biomarkers associated with ICD that can classify patients based on their potential response to ICD immunotherapy would be highly advantageous. Therefore the goal of the study is to better understand and identify what patients with bladder urothelial carcinoma (BLCA) will respond to immunotherapy by analyzing ICD signatures and investigate ICD-related prognostic factors in the context of BLCA. The data obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases regarding BLCA and normal samples was categorized based on ICD-related genes (IRGs). Specifically, we conducted an immunohistochemical (IHC) experiment to validate the expression levels of Calreticulin (CALR) in both tumor and adjacent tissues, and evaluated its prognostic significance using the Kaplan-Meier (KM) curve. Subsequently, the samples from TCGA were divided into two subtypes using consensus clustering. To obtain a more comprehensive comprehension of the biological functions, we utilized Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA). The calculation of immune landscape between two subtypes was performed through ESTIMATE and CIBERSORT. Risk models were constructed using Cox and Lasso regression and their prognosis predictive ability was evaluated using nomogram, receiver operating characteristic (ROC), and calibration curves. Finally, Tumor Immune Dysfunction and Exclusion (TIDE) algorithms was utilized to predict the response to immunotherapy. A total of 34 IRGs were identified, with most of them exhibiting upregulation in BLCA samples. The expression of CALR was notably higher in BLCA compared to the adjacent tissue, and this increase was associated with an unfavorable prognosis. The differentially expressed genes (DEGs) associated with ICD were linked to various immune-related pathways. The ICD-high subtypes exhibited an immune-activated tumor microenvironment (TME) compared to the ICD-low subtypes. Utilizing three IRGs including CALR, IFNB1, and IFNG, a risk model was developed to categorize BLCA patients into high- and low-risk groups. The overall survival (OS) was considerably greater in the low-risk group compared to the high-risk group, as evidenced by both the TCGA and GEO cohorts. The risk score was identified as an independent prognostic parameter (all p < 0.001). Our model demonstrated good predictive ability (The area under the ROC curve (AUC), AUC1-year= 0.632, AUC3-year= 0.637, and AUC5-year =0.653). Ultimately, the lower risk score was associated with a more responsive immunotherapy group. The potential of the ICD-based risk signature to function as a marker for evaluating the prognosis and immune landscape in BLCA suggests its usefulness in identifying the suitable population for effective immunotherapy against BLCA.

Chen L ，Lin J ，Wen Y ，Chen Y ，Chen CB ... -  《Frontiers in Oncology》

Identification of a novel defined inflammation-related long noncoding RNA signature contributes to predicting prognosis and distinction between the cold and hot tumors in bladder cancer.

Bladder cancer (BLCA) is one of the most frequently diagnosed urological malignancies and is the 4th most common cancer in men worldwide. Molecular targets expressed in bladder cancer (BLCA) are usually used for developing targeted drug treatments. However, poor prognosis and poor immunotherapy efficacy remain major challenges for BLCA. Numerous studies have shown that long non-coding RNAs (LncRNAs) play an important role in the development of cancer. However, the role of lncRNAs related to inflammation in BLCA and their prognostic value remain unclear. Therefore, this study is aimed to explore new potential biomarkers that can predict cancer prognosis. We downloaded BLCA-related RNA sequencing data from The Cancer Genome Atlas (TCGA) and searched for inflammation-related prognostic long non-coding RNAs (lncRNAs) by univariate Cox (uniCox) regression and co-expression analysis. We used the least absolute shrinkage and selection operator (LASSO) analysis to construct an inflammation-related lncRNA prognosis risk model. Samples were divided into high-risk score (HRS) group and low-risk score (LRS) group based on the median value of risk scores. The independent variable factors were identified by univariate Cox (uni-Cox) and multivariate Cox (multi-Cox) regression analyses, and receiver operating characteristic (ROC) curves were used to compare the role of different factors in predicting outcomes. Nomogram and Calibration Plot were generated by the R package rms to analyze whether the prediction results are correct and show good consistency. Correlation coefficients were calculated by Pearson analysis. The Kaplan-Meier method was used to assess the prognostic value. The expression of 7 lncRNAs related with inflammation was also confirmed by qRT-PCR in BLCA cell lines. Kyoto Encyclopedia of Gene and Genome (KEGG) pathways that were significantly enriched (P < 0.05) in each risk group were identified by the GSEA software. The R package pRRophetic was used to predict the IC50 of common chemotherapeutic agents. TIMER, XCELL, QUANTISEQ, MCPCOUNTER, EPIC and CIBERSORT were applied to quantify the relative proportions of infiltrating immune cells. We also used package ggpubr to evaluate TME scores and immune checkpoint activation in LRS and HRS populations. R package GSEABase was used to analyze the activity of immune cells or immune function. Different clusters of principal component analysis (PCA), t-distribution random neighborhood embedding (t-SNE), and Kaplan-Meier survival were analyzed using R package Rtsne's. The R package ConsensesClusterPlus was used to class the inflammation-related lncRNAs. In this study, a model containing 7 inflammation-related lncRNAs was constructed. The calibration plot of the model was consistent with the prognosis prediction outcomes. The 1-, 3-, and 5-year ROC curve (AUC) were 0.699, 0.689, and 0.699, respectively. High-risk patients were enriched in lncRNAs related with tumor invasion and immunity, and had higher levels of immune cell infiltration and immune checkpoint activation. Hot tumors and cold tumors were effectively distinguished by clusters 2 and 3 and cluster 1, respectively, which indicated that hot tumors are more susceptible to immunotherapy. Our study showed that inflammation-related LncRNAs are closely related with BLCA, and inflammation-related lncRNA can accurately predict patient prognosis and effectively differentiate between hot and cold tumors, thus improving individualized immunotherapy for BLCA patients. Therefore, this study provides an effective predictive model and a new therapeutic target for the prognosis and clinical treatment of BLCA, thus facilitating the development of individualized tumor therapy.

Xiong X ，Chen C ，Li X ，Yang J ，Zhang W ，Wang X ，Zhang H ，Peng M ，Li L ，Luo P ... -  《Frontiers in Oncology》