Bioinformatic identification of prognostic indicators in bladder cancer

KRT13-expressing epithelial cell population predicts better response to chemotherapy and immunotherapy in bladder cancer: Comprehensive evidences based on BCa database

Neoadjuvant chemotherapy (NAC) prior to surgery and immune checkpoint therapy (ICT) has revolutionized bladder cancer (BCa) treatment. Patients likely to benefit from these therapies need to be accurately stratified; however, this remains a major clinical challenge. In the present study, single-cell RNA sequencing was used to evaluate the predictive ability of an epithelial cell population highly expressing keratin 13 (KRT13) to assess therapeutic response in BCa. The presence of KRT13-enriched tumors indicated favorable outcomes after NAC and superior response to ICT in patients with BCa. Furthermore, KRT13 population characteristics appeared to be closely related to changes in the immune microenvironment in the vicinity of this cell population. We constructed a prognostic model using an artificial neural network based on the gene signatures in the KRT13 population; the model demonstrated strong robustness and superiority. Additionally, a user-friendly and open-access web application named BCa database was developed for researchers to study BCa by mining the connective map database.

A single-sample mRNA molecular classification of bladder cancer predicting prognosis and response to immunotherapy

Background

As an immunogenic cancer, crosstalk between cancer cells and immune cells has been gradually recognized in bladder cancer (BC). Several studies have emphasized the clinical significance of the molecular stratification of BC without highlighting the role of the immune microenvironment. Although immunotherapy acted as a prospective treatment, more precise molecular stratification should be established to select those sensitive to immunotherapy.

Methods

To select specific immune genes forming subtypes indicating disparate prognoses, we performed bioinformatic analysis using BC transcriptomic profiles from six published datasets, with 408 BC samples in The Cancer Genome Atlas (TCGA) database and 295 individuals in International Cancer Genome Consortium (ICGC) database. Survival analyses were conducted using Kaplan-Meier curves, while Kruskal-Wallis tests were applied to test the differences among groups. Except for unsupervised clustering based on the differential expression of genes, we additionally performed binomial logistic regression, focusing on the mRNA level of a single sample.

Results

Unsupervised clustering showed that 4 clusters captured the best segmentation. After validation with survival data and simplification using binomial logistic regression, we found that cluster B and cluster D showed worse survival outcomes (P=0.012). Considering the similar survival outcomes of these two clusters, we recombined and performed another survival analysis, which also showed significant survival differences (P=0.0041). Bonding with clinical data, a greater proportion of risk factors were assigned to the worse prognosis subtype, especially showing higher grades in the subtype (P<0.001). In addition, immune cell infiltration, single nucleotide polymorphism (SNP) and copy number variation (CNV) all showed differences between clusters, indicating changes in the immune microenvironment and mutation burden. Through phenotypical analysis, we found metabolism and proliferation phenotypes associated with the immune clusters and mutually exclusive in BC, of which proliferation contributed to worse outcomes. Using the tumor immune dysfunction and exclusion (TIDE) score, a worse immunotherapy benefit was predicted in clusters B&D, defined as the worse prognosis subtype.

Conclusions

With this novel clustering criterion based on immune-related genes, we provide a better understanding of the immune microenvironment, further guiding the use of immunotherapy.

Identification of an immune gene-associated prognostic signature in patients with bladder cancer

A deeper understanding of the interaction between tumor cell and the immune microenvironment in bladder cancer may help select predictive and prognostic biomarkers. The current study aims to construct a prognostic signature for bladder cancer by analysis of molecular characteristics, as well as tumor-immune interactions. RNA-sequencing and clinical information from bladder cancer patients were downloaded from the TCGA database. The single sample Gene Sets Enrichment Analysis (ssGSEA) and Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) were employed to separate the samples into two clusters. Lasso Cox regression was performed to construct an immune gene signature for bladder cancer. The correlation between key target genes of immune checkpoint blockade and the prognostic signature was also analyzed. Dataset from Gene Expression Omnibus (GEO) was retrieved for validation. Two immunophenotypes and immunological characteristics were identified, and a 17-immune gene signature was constructed to provide an independent prognostic signature for bladder cancer. The signature was verified through external validation and correlated with genomic characteristics and clinicopathologic features. Finally, a nomogram was generated from the clinical characteristics and immune signature. Our study reveals a tumor-immune microenvironment signature useful for prognosis in bladder cancer. The results provide information on the potential development of treatment strategies for bladder cancer patients. Prospective studies are warranted to validate the prognostic capability of this model, but these data highlight the role of the microenvironment in the clinical outcome of patients.

Accounting for B-cell Behavior and Sampling Bias Predicts Anti-PD-L1 Response in Bladder Cancer

Cancer immunotherapy is predominantly based on T cell-centric approaches. At the same time, the adaptive immune response in the tumor environment also includes clonally produced immunoglobulins and clonal effector/memory B cells that participate in antigen-specific decisions through their interactions with T cells. Here, we investigated the role of infiltrating B cells in bladder cancer via patient dataset analysis of intratumoral immunoglobulin repertoires. We showed that the IgG1/IgA ratio is a prognostic indicator for several subtypes of bladder cancer and for the whole IMVigor210 anti-PD-L1 immunotherapy study cohort. A high IgG1/IgA ratio associated with the prominence of a cytotoxic gene signature, T-cell receptor signaling, and IL21-mediated signaling. Immunoglobulin repertoire analysis indicated that effector B-cell function, rather than clonally produced antibodies, was involved in antitumor responses. From the T-cell side, we normalized a cytotoxic signature against the extent of immune cell infiltration to neutralize the artificial sampling-based variability in immune gene expression. Resulting metrics reflected proportion of cytotoxic cells among tumor-infiltrating immune cells and improved prediction of anti-PD-L1 responses. At the same time, the IgG1/IgA ratio remained an independent prognostic factor. Integration of the B-cell, natural killer cell, and T-cell signatures allowed for the most accurate prediction of anti-PD-L1 therapy responses. On the basis of these findings, we developed a predictor called PRedIctive MolecUlar Signature (PRIMUS), which outperformed PD-L1 expression scores and known gene signatures. Overall, PRIMUS allows for reliable identification of responders among patients with muscle-invasive urothelial carcinoma, including the subcohort with the low-infiltrated "desert" tumor phenotype.

Bioinformatic analysis identifying FGF1 gene as a new prognostic indicator in clear cell Renal Cell Carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) has been the commonest renal cell carcinoma (RCC). Although the disease classification, diagnosis and targeted therapy of RCC has been increasingly evolving attributing to the rapid development of current molecular pathology, the current clinical treatment situation is still challenging considering the comprehensive and progressively developing nature of malignant cancer. The study is to identify more potential responsible genes during the development of ccRCC using bioinformatic analysis, thus aiding more precise interpretation of the disease METHODS: Firstly, different cDNA expression profiles from Gene Expression Omnibus (GEO) online database were used to screen the abnormal differently expressed genes (DEGs) between ccRCC and normal renal tissues. Then, based on the protein-protein interaction network (PPI) of all DEGs, the module analysis was performed to scale down the potential genes, and further survival analysis assisted our proceeding to the next step for selecting a credible key gene. Thirdly, immunohistochemistry (IHC) and quantitative real-time PCR (QPCR) were conducted to validate the expression change of the key gene in ccRCC comparing to normal tissues, meanwhile the prognostic value was verified using TCGA clinical data. Lastly, the potential biological function of the gene and signaling mechanism of gene regulating ccRCC development was preliminary explored.

RESULTS

Four cDNA expression profiles were picked from GEO database based on the number of containing sample cases, and a total of 192 DEGs, including 39 up-regulated and 153 down-regulated genes were shared in four profiles. Based on the DEGs PPI network, four function modules were identified highlighting a FGF1 gene involving PI3K-AKT signaling pathway which was shared in 3/4 modules. Further, both the IHC performed with ccRCC tissue microarray which contained 104 local samples and QPCR conducted using 30 different samples confirmed that FGF1 was aberrant lost in ccRCC. And Kaplan-Meier overall survival analysis revealed that FGF1 gene loss was related to worse ccRCC patients survival. Lastly, the pathological clinical features of FGF1 gene and the probable biological functions and signaling pathways it involved were analyzed using TCGA clinical data.

CONCLUSIONS

Using bioinformatic analysis, we revealed that FGF1 expression was aberrant lost in ccRCC which statistical significantly correlated with patients overall survival, and the gene's clinical features and potential biological functions were also explored. However, more detailed experiments and clinical trials are needed to support its potential drug-target role in clinical medical use.

STAT1 is a key gene in a gene regulatory network related to immune phenotypes in bladder cancer: An integrative analysis of multi-omics data

The immunophenotype of bladder cancer plays a pivotal role in the prognosis of cancer, but the effect of different epigenetic factors on different immunophenotypes in bladder tumours remains unclear. This study used multi-omics data analysis to provide molecular basis support for different immune phenotypes. Unsupervised cluster analysis revealed distinct subclusters with higher (subcluster B2) or lower cytotoxic immune phenotypes (subcluster A1) related to PD-L1 and IFNG expression. Mutational landscape analyses showed that the mutation level of TP53 in subcluster B1 was highest than other subclusters, and subcluster B1 had a lower frequency of concurrent mutation than subcluster A2. A total of 2364 differentially expressed genes were identified between subclusters A2 and B1, and the main functions of the up-regulated genes in subcluster B1 were enriched in the activation of T cells and other related pathways. We found that STAT1 was a key gene in a gene regulatory network related to immune phenotypes in bladder cancer. Finally, we constructed a prognostic prediction model by LASSO Cox regression which could distinguish high-risk and low-risk cases significantly. In conclusion, the present study addressed a field synopsis between genetic and epigenetic events in immune phenotypes of bladder cancer.