Identification of tumor immunophenotypes associated with immunotherapy response in bladder cancer

OBJECTIVES

This study aims to reveal immunophenotypes associated with immunotherapy response in bladder cancer, identify the signature genes of immune subtypes, and provide new molecular targets for improving immunotherapy response.

METHODS

Bladder cancer immunophenotypes were characterized in the bulk RNA sequencing dataset GSE32894 and Imvigor210, and gene expression signatures were established to identify the immunophenotypes. Expression of gene signatures were validated in single-cell RNA sequencing dataset GSE145140 and human proteins expression data source. Investigation of Immunotherapy Response was performed in IMvigor210 dataset. Prognosis of tumor immunophenotypes was further analyzed.

RESULTS

Inflamed and immune-excluded immunophenotypes were characterized based on the tumor immune cell scores. Risk score models that were established rely on RNA sequencing profiles and overall survival of bladder cancer cohorts. The inflamed tumors had lower risk scores, and the low-risk tumors were more likely to respond to atezolizumab, receiving complete response/partial response (CR/PR). Patients who responded to atezolizumab had higher SRRM4 and lower NPHS1 and TMEM72 expression than the non-responders. SRRM4 expression was a protective factor for bladder cancer prognosis, while the NPHS1 and TMEM72 showed the opposite pattern.

CONCLUSION

This study provided a novel classification method for tumor immunophenotypes. Bladder cancer immunophenotypes can predict the response to immune checkpoint blockade. The immunophenotypes can be identified by the expression of signature genes.

Metadherin Promotes the Development of Bladder Cancer by Enhancing Cell Division

Background: Bladder cancer (BLCA) is a malignant tumor occurring in bladder mucosa. Metadherin (MTDH) has been implicated in tumor progression; however, its molecular biological mechanisms in BLCA remain unclear. Materials and Methods: Cell functions were tested after BLCA cells were transfected by both short hairpin RNAs and small interfering RNAs to silence MTDH. Furthermore, in-house RNA sequencing (RNA-seq) was performed with T24 cells after the knockdown of MTDH. In addition, MTDH-related pathways were explored. Finally, MTDH mRNA and protein expression levels were examined using multiple detection methods in BLCA tissues. Results: MTDH knockdown could largely inhibit cell proliferation, viability, and migration and induce apoptosis of BLCA cells. In-house RNA-seq showed that MTDH knockdown led to extracellular matrix organization and cell division. The integrated analysis showed that the comprehensive expression of MTDH at the mRNA level was 0.47 and that at the protein level was 0.54, based on 11 platforms, including 1485 BLCA and 180 non-BLCA samples. Conclusions: MTDH promotes the growth of BLCA cells through the pathway of cell division. This study provides new directions and biomarkers for future treatment.

Nitidine chloride regulates cell function of bladder cancer in vitro through downregulating Lymphocyte antigen 75

Nitidine chloride (NC) is effective on cancer in many tumors, but its effect on bladder cancer (BC) is unknown. We conducted cell function experiments to verify the antineoplastic effect of NC on BC cell lines (5637, T24, and UM-UC-3) in vitro. Then, mRNAs of NC-treated and NC-untreated BC cells were extracted for mRNA sequencing. Differentially expressed genes (DEGs), expression analysis, and drug molecular docking were conducted to discover the target gene of NC. Finally, functional enrichment was analyzed to explore the underlying mechanisms. NC dramatically inhibited proliferation, migration, and invasion, and it induced apoptosis and arrested the S and G2/M phases of BC cell lines. Lymphocyte antigen 75 (LY75) appeared to be the target of NC. LY75 was highly expressed and had the ability to distinguish BC tissue from non-cancerous tissue. Then, drug molecular docking confirmed the targeting relationship between NC and LY75. Gene enrichment analysis showed that the downregulated genes, after being treated with NC, were mainly enriched in pathways relevant to cell pathophysiological processes. NC inhibits BC cell proliferation, migration, and invasion, induces apoptosis, and arrests cell cycles by downregulating the expression of LY75. This study provides molecular and theoretical bases for NC treatment of BC.

A study comparing dusting to basketing for renal stones ≤ 2 cm during flexible ureteroscopy

OBJECTIVES

To compare the dusting efficiency and safety with basketing for treating renal stones ≤ 2 cm during flexible ureteroscopy (fURS).

MATERIALS AND METHODS

This study included 218 patients with renal stones ≤ 2 cm treated with fURS. Among them, 106 patients underwent dusting, and 112 patients underwent fragmentation with basket extraction. All patients were followed up for 3 months postoperatively. The operating time, lasing time, stone-free rate (SFR) and complication rate were compared.

RESULTS

The mean stone size in the dusting group was 1.3 cm, whereas 1.4 cm in the basketing group. The mean operative time was significantly lower in the dusting group than in the basketing group (43.1±11.7 minutes VS 60.5±13.4 minutes, P < 0.05), but the lasing time was significantly longer for the dusting group than for the basketing group (17.7±3.9 minutes VS 14.1±3.6 minutes, P < 0.05). SFR was significantly higher in the basketing group immediately after the operation and follow-up after 1 month (76.8% vs 55.7%, P= 0.001 and 88.4% vs 78.3%, P = 0.045). However, the SFR was similar for both groups (88.8% in the dusting group vs. 90.2% in the basketing group) after 3 months postoperatively. There was no statistical difference in the complication rates between the two groups.

CONCLUSIONS

Dusting has advantages in shortening the operation time and reducing the operation cost, but the lasing time was longer compared with the basketing. Although there is no difference in long-term effect, basketing is superior to dusting in terms of short-term SFR. Moreover, dusting should be avoided in some special cases and basketing a better choice. Both techniques are effective for the treatment of renal stones ≤ 2 cm and choice depends on patient demographic and stone characteristics.

Clinical significance and potential regulatory mechanism of overexpression of pituitary tumor-transforming gene transcription factor in bladder cancer

Background

Pituitary tumor transforming gene-1 (PTTG1) transcription factor is identified as carcinogenic and associated with tumor invasiveness, but its role in bladder cancer (BLCA) remains obscure. This research is intended to analyze the aberrant expression and clinical significance of PTTG1 in BLCA, explore the relationship between PTTG1 and tumor microenvironment characteristics and predict its potential transcriptional activity in BLCA tissue.

Methods

We compared the expression discrepancy of PTTG1 mRNA in BLCA and normal bladder tissue, using the BLCA transcriptomic datasets from GEO, ArrayExpress, TCGA, and GTEx. In-house immunohistochemical staining was implemented to determine the PTTG1 protein intensity. The prognostic value of PTTG1 was evaluated using the Kaplan-Meier Plotter. CRISPR screen data was utilized to estimate the effect PTTG1 interference has on BLCA cell lines. We predicted the abundance of the immune cells in the BLCA tumor microenvironment using the microenvironment cell populations-counter and ESTIMATE algorithms. Single-cell RNA sequencing data was applied to identify the major cell types in BLCA, and the dynamics of BLCA progression were revealed using pseudotime analysis. PTTG1 target genes were predicted by CistromeDB.

Results

The elevated expression level of PTTG1 was confirmed in 1037 BLCA samples compared with 127 non-BLCA samples, with a standardized mean difference value of 1.04. Higher PTTG1 expression status exhibited a poorer BLCA prognosis. Moreover, the PTTG1 Chronos genetic effect scores were negative, indicating that PTTG1 silence may inhibit the proliferation and survival of BLCA cells. With PTTG1 mRNA expression level increasing, higher natural killer, cytotoxic lymphocyte, and monocyte lineage cell infiltration levels were observed. A total of four candidate targets containing CHEK2, OCIAD2, UBE2L3, and ZNF367 were determined ultimately.

Conclusions

PTTG1 mRNA over-expression may become a potential biomarker for BLCA prognosis. Additionally, PTTG1 may correlate with the BLCA tumor microenvironment and exert transcriptional activity by targeting CHEK2, OCIAD2, UBE2L3, and ZNF367 in BLCA tissue.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09810-y.

Small Nucleolar RNAs (snoRNAs)-Based Risk Score Classifier Predicts Overall Survival in Bladder Carcinoma

Background

Bladder carcinoma (BLCA) is a leading cause of cancer-related deaths worldwide. The aim of this work was to develop an accurate stratification in predicting the prognosis and directing the treatment of BLCA patients based on small nucleolar RNAs (snoRNAs).

Material/Nethods

Expression profiles of snoRNAs were downloaded from the SNORic database. The expression profiles and clinical outcomes of BLCA patients were analyzed. Survival-associated snoRNAs were identified and used to develop a novel risk score classifier. Genes in the whole genome that were significantly correlated with the included prognostic snoRNAs were used for functional enrichment analysis.

Results

The results showed that age, American Joint Committee on Cancer (AJCC) stage, and tumor status were significantly correlated with overall survival (OS) of BLCA patients. We selected 12 survival-associated snoRNAs to build a prognostic signature. Patients were separated into high- and low-risk groups based on the median value of the risk score. Patients in the high-risk group and low-risk group have distinct clinical outcomes. The AJCC TNM stage showed moderate utility as a prognostic indicator for clinical outcome prediction. Then, clinical parameters and risk scores were entered in multivariate Cox analysis. Notably, the prognostic signature remained an independent significant prognostic risk factor. The pathway analysis suggested that these genes were enriched in several types of cancer and “Focal adhesion” pathways.

Conclusions

The prognostic signature defined by expression profiles of 12 survival-associated snoRNAs appears to be an excellent predictor of the clinical outcome of BLCA patients.

Development and validation of an immune prognostic classifier for clear cell renal cell carcinoma

BACKGROUND

Tumor-infiltrating immune cells are indispensable to the progression and prognosis of clear cell renal cell carcinoma (ccRCC).

OBJECTIVE

The aim of this study was to explore the clinical implications of immune cell infiltrates in ccRCC.

METHODS

The Cancer Genome Atlas (TCGA) database (N= 515) and E-MTAB-1980 cohort of patients (N= 101) were adopted to estimate the prognostic value of immune cell infiltration. Twenty-four types of immune cells were evaluated using single-sample gene set enrichment analysis. Cox regression analyses were conducted to develop an immune risk score.

RESULTS

Survival analyses revealed that 13 genes significantly associated with the overall survival (OS). Furthermore, multivariate Cox analysis identified an immune risk score on the basis of mast cells, natural killer CD56bright cells, T helper 17 (Th17) cells, and Th2 cells. The immune risk score was associated with OS, with hazard ratios of 2.72 (95% CI 2.17-3.40) and 3.24 (95% CI 1.64-6.44) in TCGA and E-MTAB-1980 datasets, respectively. This immune risk score was significantly correlated with some immunotherapy-related biomarkers.

CONCLUSIONS

We profiled a prognostic signature and established an immune risk score model for ccRCC, which could provide novel predictive markers for patients with ccRCC and an indicator for immunotherapy response measurement.

Expression levels and co‑targets of miRNA‑126‑3p and miRNA‑126‑5p in lung adenocarcinoma tissues: Αn exploration with RT‑qPCR, microarray and bioinformatic analyses

Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer. Previous studies have found that many microRNAs (miRNAs), including miRNA-126-3p, may play a critical role in the development of LUAD. However, no study of LUAD has researched the synergistic effects and co-targets of both miRNA-126-3p and miRNA-126-5p. The present study used real-time quantitative polymerase chain reaction (RT-qPCR) to explore the expression values of miRNA-126-3p and miRNA-126-5p in 101 LUAD and 101 normal lung tissues. Ten relevant microarray datasets were screened to further validate the expression levels of miRNA-126-3p and −5p in LUAD. Twelve prediction tools were employed to obtain potential targets of miRNA-126-3p and miRNA-126-5p. The results showed that both miRNA-126-3p and −5p were expressed significantly lower in LUAD. A significant positive correlation was also present between miRNA-126-3p and −5p expression in LUAD. In addition, lower expression of miRNA-126-3p and −5p was indicative of vascular invasion, lymph node metastasis (LNM), and a later tumor/node/metastasis (TNM) stage of LUAD. The authors obtained 167 targets of miRNA-126-3p and 212 targets of miRNA-126-5p; 44 targets were co-targets of both. Eight co-target genes (IGF2BP1, TRPM8, DUSP4, SOX11, PLOD2, LIN28A, LIN28B and SLC7A11) were initially identified as key genes in LUAD. The results of the present study indicated that the co-regulation of miRNA-126-3p and miRNA-126-5p plays a key role in the development of LUAD, which also suggests a fail-proof mode between miRNA-3p and miRNA-126-5p.

Evaluation of the HOXA11 level in patients with lung squamous cancer and insights into potential molecular pathways via bioinformatics analysis

Background

This study was carried out to discover the underlying role that HOXA11 plays in lung squamous cancer (LUSC) and uncover the potential corresponding molecular mechanisms and functions of HOXA11-related genes.

Methods

Twenty-three clinical paired LUSC and non-LUSC samples were utilized to examine the level of HOXA11 using quantitative real-time polymerase chain reaction (qRT-PCR). The clinical significance of HOXA11 was systematically analyzed based on 475 LUSC and 18 non-cancerous adjacent tissues from The Cancer Genome Atlas (TCGA) database. A total of 102 LUSC tissues and 121 non-cancerous tissues were available from Oncomine to explore the expressing profiles of HOXA11 in LUSC. A meta-analysis was carried out to further assess the differential expression of HOXA11 in LUSC, including in-house qRT-PCR data, expressing data extracted from TCGA and Oncomine databases. Moreover, the enrichment analysis and potential pathway annotations of HOXA11 in LUSC were accomplished via Gene Oncology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of hub genes and according correlations with HOXA11 were assessed to further explore the biological role of HOXA11 in LUSC.

Results

HOXA11 expression in LUSC had a tendency to be upregulated in comparison to adjacent non-cancerous tissues by qRT-PCR. TCGA data displayed that HOXA11 was remarkably over-expressed in LUSC compared with that in non-LUSC samples, and the area under curves (AUC) was 0.955 (P < 0.001). A total of 1523 co-expressed genes were sifted for further analysis. The most significant term enriched in the KEGG pathway was focal adhesion. Among the six hub genes of HOXA11, including PARVA, ILK, COL4A1, COL4A2, ITGB1, and ITGA5, five (with the exception of COL4A1) were significantly decreased compared with the normal lung tissues. Moreover, the expression of ILK was negatively related to HOXA11 (r = − 0.141, P = 0.002).

Conclusion

High HOXA11 expression may lead to carcinogenesis and the development of LUSC. Furthermore, co-expressed genes might affect the prognosis of LUSC.

Upregulation of HOXA13 as a potential tumorigenesis and progression promoter of LUSC based on qRT-PCR and bioinformatics

In this study, we investigated the levels of homeobox A13 (HOXA13) and the mechanisms underlying the co-expressed genes of HOXA13 in lung squamous cancer (LUSC), the signaling pathways in which the co-expressed genes of HOXA13 are involved and their functional roles in LUSC. The clinical significance of 23 paired LUSC tissues and adjacent non-tumor tissues were gathered. HOXA13 levels in LUSC were detected by quantitative real-time polymerase chain reaction (qRT-PCR). HOXA13 levels in LUSC from The Cancer Genome Atlas (TCGA) and Oncomine were analyzed. We performed receiver operator characteristic (ROC) curves of various clinicopathological features of LUSC. Co-expressed of HOXA13 were collected from MEM, cBioPortal and GEPIA. The functions and pathways of the most reliable overlapped genes were achieved from the Gene Otology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. The protein-protein interaction (PPI) networks were mapped using STRING. HOXA13 in LUSC were markedly upregulated compared with those in the non-cancerous controls as demonstrated by qRT-PCR (LUSC: 0.330±0.360; CONTROLS: 0.155±0.142; P=0.021). TCGA (LUSC: 6.388±2.097, CONTROLS: 1.157±0.719; P<0.001) and Hou's study from Oncomine (LUSC: 1.154±0.260; CONTROLS: 0.957±0.065; P=0.001) showed the same tendency. Meanwhile, the area under the curve (AUC) of TNM was calculated as 0.877 with P=0.002. Based on the HOXA13 expression data from TCGA, the ROC of the tissue types was calculated as AUC=0.971 (P<0.001). In addition, 506 genes were filtered as co-expression genes of HOXA13. The 3 most significant KEGG pathways were metabolic pathways (P=5.41E-15), the calcium signaling pathway (P=3.01E-11), and the cAMP signaling pathway (P=5.63E-11). MAPK1, GNG7, GNG12, PRKCA were selected as the hub genes. In conclusion, HOXA13 was upregulated and related to the TNM stage in LUSC. The expression of hub genes in LUSC might be deregulated by HOXA13. Moreover, the 4 co-expressed hub genes of HOXA13 might be crucial biomarkers for the diagnosis and prognosis of LUSC, as well as the development of novel therapeutic targets against LUSC.