Skip to main content
Download PDF

Development of a prognostic model to predict BLCA based on anoikis-related gene signature: preliminary findings

BMC Urology volume 23, Article number: 199 (2023) Cite this article

Abstract

Background

The prevalence of bladder urothelial carcinoma (BLCA) is significant on a global scale. Anoikis is a type of procedural cell death that has an important role in tumor invasion and metastasis. The advent of single-cell RNA sequencing (scRNA-seq) approaches has revolutionized the genomics field by providing unprecedented opportunities for elucidating cellular heterogeneity. Understanding the mechanisms associated with anoikis in BLCA is essential to improve its survival rate.

Methods

Data on BLCA and clinical information were acquired from the databases of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). ARGs were obtained from Genecards and Harmonizome databases. According to univariate Cox regression analysis, the least absolute shrinkage and selection operator (LASSO) algorithm was utilized to select the ARGs associated with the overall rate (OS). A multivariate Cox regression analysis was carried out to identify eight prognostic ARGs, leading to the establishment of a risk model. The OS rate of BLCA patients was evaluated using Kaplan–Meier survival analysis. To explore the molecular mechanism in low- and high-risk groups, we employed Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSVA). Immune infiltration landscape estimation was performed using ESTIMATE, CIBERSOT, and single sample gene set enrichment analysis (ssGSEA) algorithms. Patients were categorized into different subgroups through consensus clustering analysis.

We employed biological functional enrichment analysis and conducted immune infiltration analysis to examine the disparities in potential biological functions, infiltration of immune cells, immune activities, and responses to immunotherapy.

Results

We identified 647 ARGs and 37 survival-related genes. We further developed a risk scoring model to quantitatively assess the predictive capacity of ARGs. The high-risk score group exhibited an unfavorable prognosis, whereas the low-risk score group demonstrated a converse effect. We also found that the two groups of patients might respond differently to immune targets and anti-tumor drugs.

Conclusion

The nomogram with 8 ARGs may help guide treatment of BLCA. The systematic assessment of risk scores can help to design more individualized and precise treatment strategies for BLCA patients.

Peer Review reports

Introduction

Bladder urothelial carcinoma (BLCA) ranks as the ninth most prevalent cancer globally, exhibiting a significant fatality rate. It is estimated that around 550,000 individuals are diagnosed with this condition annually [1, 2]. Patients with non-muscle invasive BLCA have a high rate of recurrence after surgery [3]. Muscle-invasive BLCA has a poor prognosis, with a 5-year overall survival rate of 40–50% [4]. Because there is a deficiency in identifiable biomarkers, a portion of individuals suffering from severe BLCA are still unable to receive effective treatment. Therefore, it is particularly important to construct an effective BLCA risk prediction model.

Recently, it has been discovered that during the process of metastasis, cancer cells transport extracellular matrix (ECM) while migrating and proliferating. Anoikis is triggered when tumor cells detach from the ECM [5]. Anoikis is a specific form of apoptosis, and its occurrence is an important mechanism of tumor invasion and metastasis [6]. However, cancer cells have the ability to develop various strategies to evade anoikis, including the utilization of growth factors, manipulation of pH levels, and adaptation to oxidative stress [7,8,9]. This discovery suggests ARGs may be exploited for cancer therapy.

Revealing novel avenues for addressing biological and medical challenges, single-cell RNA sequencing (scRNA-seq) presents unprecedented opportunities [10]. ScRNA-seq approaches have uncovered novel biology in terms of tissue composition, transcriptional dynamics, and regulatory relationships between genes [11]. Currently, scRNA-seq technology has become the most advanced method to reveal the heterogeneity and complexity of RNA transcripts within a single-cell [12]. Researchers can use scRNA-seq technology to build high-resolution profiles, which may provide better solutions for disease treatment.

In this study, we explored to analyze the relationship between the differential expression of ARGs and BLCA. Subsequently, a risk assessment model was developed to forecast the prognosis and immunotherapy outcomes for patients with BLCA. Finally, we performed scRNA-seq to build the foundation for the design of treatment regimens and the selection of targeted drugs. Exploring ARGs' prognostic model can help develop more personalized and accurate treatment strategies for BLCA patients. This study is entirely based on publicly available data.

Materials and methods

Data collection

Gene expression and clinical data for patients with BLCA were obtained from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/) and Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/ geo/query/acc.cgi?acc = GSE19423) databases. Limma is a software package in R/Bioconductor that offers a comprehensive solution for the analysis of gene expression experiment data [13]. The SVA package has the effect of identifying and removing batches [14]. The integration of the two sets of gene expression data was performed using the "limma" and "sva" packages.

Identification of differential genes and characteristic differences of ARGs

ARGs were obtained from Genecards (https://www.genecards.org/) and Harmonizome (https://maayanlab.cloud/Harmonizome/) databases. Genecards and Harmonizome are public databases that can directly search for the required gene names. In the Genecards database, ARGs were selected with a relevance score > 0.4. All ARGs in the Harmonizome database were included. We utilized the limma package in R software to identify differentially expressed ARGs, using a threshold of |log2 fold change (FC)|> 1 and false discovery rate (FDR) < 0.05. The limma, survival and survminer packages were utilized for conducting survival analysis of AGRs. The RCircos package offers a convenient and versatile approach for generating R-based 2D track plots in the style of Circos [15]. CNV copy number data were obtained from UCSC Xena (https://xena.ucsc.edu/). We used RCircos package to further explore the alteration of these lost apoptosis-related genes on chromosomes and the location of each gene on chromosomes.

ARGs subtype identification

In order to identify distinct patterns of ARGs, a clustering analysis was conducted using the "consensusClusterPlus" algorithm. According to the optimal classification of cluster number = 2, the patients with BLCA were clustered into two distinct molecular subtypes for further analysis. The principal component analysis (PCA) is a technique used to reduce the dimensionality of datasets, enhancing interpretability while minimizing information loss [16]. t-Distributed Stochastic Neighbor Embedding (tSNE) is a well-suited technique for the visualization of high-dimensional data [17]. The Uniform Manifold Approximation and Projection (UMAP) is a viable approach that can effectively decrease the dimensionality of beta diversity distance matrices [18]. PCA, tSNE and UMAP were plotted by the Rtsne, limma, ggplot2 and UMAP packages in R. The infiltration level of 23 types of immune cells was assessed using the "GSVA" R packages through a single sample gene set enrichment analysis (ssGSEA) algorithm. The "GSVA" R package was used to analyze the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), enriching the ARGs in biological processes and signaling pathways. The distinct signal pathway between the two risk groups was described using gene set enrichment analysis (GSEA) to demonstrate the KEGG pathway enrichment.

Construction and evaluation of risk score prognostic model

The training set was utilized to perform least absolute shrinkage and selection operator (LASSO) Cox regression analysis using the 'glmnet' package in R, aiming to develop a prognostic signature for ARG. Finally, we conducted multivariate Cox regression analysis to identify highly correlated genes and construct a prognostic gene signature. The risk score was calculated using the following formula: risk score = gene Cox coefficient (Coef) × gene expression (Exp). Where Coef was the regression coefficient and Exp was expression level of ARGs. We randomly divided samples into train group and test group in 1:1 ratio, classified the training group into high- and low-risk groups according to the median risk score [19, 20]. To assess the model's ability to predict outcomes, we utilized the "Survival ROC" R package to generate Kaplan–Meier (KM) survival curves and time-dependent receiver operating characteristic (ROC) curves.

The "rms" and "survival" R packages were utilized to construct a prognostic nomogram incorporating all independent prognostic factors. By summing the points assigned to each factor, the total points accurately predicted the 1-, 3-, and 5-year survival probabilities for patients with BLCA. Calibration curves and C-Index values were generated to assess the reliability of the survival prediction.

Immune cell correlation analysis and drug sensitivity prediction

The ESTIMATE algorithm was utilized to calculate the stromal score, immune score, tumor purity, and ESTIMATE score in malignant tumor tissues for comparison between high-risk and low-risk groups [21]. The infiltration of immune cells and the presence of stroma in the tumor microenvironment (TME) were assessed using immune score and stromal score. The ESTIMATE score was calculated by summing up the stromal and immune scores to evaluate their combined effect. To determine the relative abundance of 22 different types of immune cells, cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) analysis was employed [22]. The "oncoPredict" R package was utilized to compute the half-maximal inhibitory concentration (IC50) values of chemotherapeutic and targeted drugs for each BLCA sample.

scRNA-seq and immunohistochemical (IHC) data extraction from public databases

The Tumor Immune Single-cell Hub 2 (TISCH2) is a scRNA-seq database (http://tisch.comp-genomics.org/) that aims to comprehensively characterize the TME at single-cell resolution [23]. To investigate the differences in ARG expression between normal and tumor tissues, we extracted the immunohistochemical results of 9 ARGs based on the public database of Human Protein Atlas (https://www.proteinatlas.org/). The flowchart of the present research is shown in Supplementary Figure S1.

Statistical analysis

In our study, all statistical analysis was performed by R 4.1.3. Related R packages, including “ggplot2”, “limma”, “heatmap”, and “RColorBrewer”, and other related R packages were used. A p < 0.05 indicated statistically significant for all analyses (*p < 0.05; ** p < 0.01; ***p < 0.001).

Results

Identification of ARGs in BLCA

First, we obtained 647 ARGs from the public databases (Supplementary Table S1). The univariate analysis and survival analysis were conducted for ARGs. The volcano plot and forest plot were drawn (Fig. 1A, B). This study incorporated 37 ARGs and utilized network diagrams to depict the intricate correlation between varying ARGs and their prognostic significance in BLCA (Fig. 1C). The chromosomal visualization of CNV alterations in ARGs is depicted in Fig. 1D. In addition, CNV-associated mutations were prevalent among the 37 ARGs. S100A7, RAD9A, MYC, RAC3, FASN, F10, SPHK1, SATB1, PDGFRA and ID2 showed significant CNV amplification, but ADAMTSL1, TAGLN, RPS6KA1, CRYAB, PBK and THBS1 exhibited notable CNV deletions (Fig. 1E).

Fig. 1
figure 1

Genetic variations and expression of ARGs in BLCA. A Volcano plot of 137 DEGs, red for up-regulated, green for down-regulated genes. B The forest plot shows the 37 ARGs via the univariate Cox regression analysis. C Network diagram showing the interaction of 37 ARGs in BLCA. D The localization of the 37 ARGs on 23 chromosomes. E CNV variation frequency of 37 ARGs in BLCA

Full size image

Characteristics of different ARGs

We categorized the BLCA patients into two distinct groups, labeled A and B, utilizing the NMF algorithm (Fig. 2A). The use of “PCA”, “tSNE”, and “UMAP” R packages could clearly distinguish between groups A and B (Fig. 2B-D). The Kaplan–Meier (KM) analysis revealed a significant disparity in patient survival rates between groups A and B (Fig. 2E). We also plotted a heatmap to visualize the detailed expression of ARGs in A and B groups (Fig. 2F). Regarding the differences in immune infiltration between the two groups, we observed a significant enrichment of the A group in Monocyte and CD56 cells (Fig. 2G). GO biopathway analysis showed that ARGs in the B group were obviously enriched in the following functional sets: regulation of microglia cell activation, cell adhesion mediated by integrin, astrocyte development and others (Fig. 2H). According to the KEGG analysis, B group exhibited significant activation of pathways related to cell adhesion molecules, chemokine signaling, cytokine-receptor interaction and systemic lupus erythematosus (Fig. 2I, J).

Fig. 2
figure 2

ARGs subtype identification. A Two ARG subtypes in CRC were identified by consensus clustering analysis. B-D PCA plot of risk score in ARGs cohort. E Survival analyses for the A and B groups. F-G Difference expression analysis and immune infiltration of ARGs in two groups. H-J GO and KEGG analyses for ARGs of the two groups

Full size image

The nomogram based on risk score in BLCA

To avoid overfitting, eight ARGs were further identified by the LASSO regression method (Fig. 3A, B). We constructed a risk score prognostic model by using 8 ARGs. Risk score = EXP CSPG4 × (0.164) + EXP CALR × (-0.321) + EXP NTF3 × (-0.424) + EXP RPS6KA1 × (-0.446) + EXP FASN × (0.316) + EXP F10 × (0.288) + EXP NTRK3 × (0.592) + EXP S100A7 × (0.073). Patients with BLCA were divided into low-risk and high-risk groups based on the median risk score. Finally, we built a nomogram to predicting the survival probabilities of BLCA (Fig. 3C-E).

Fig. 3
figure 3

Construction of risk score prognostic model. A, B The least absolute shrinkage and selection operator (LASSO) method of ARGs associated with prognosis. C Nomogram construction based on the ARGs prognostic signature and clinicopathological characteristics. D, E Cumulative frequency and calibration curves for risk score model in the high and low group

Full size image

Evaluation of the risk score model

Relevant findings from survival analysis and ROC curves indicated a less favorable prognosis in the high-risk group when compared to the low-risk group (Fig. 4A-D). To assess the validity of the risk prognostic model, we performed an independent prognostic analysis (Fig. 4E). The risk factor DCAs are presented in Fig. 4F-H. The heatmap revealed a higher level of ARG expression in the high-risk group as compared to the low-risk group (Fig. 4I). Moreover, we also found that group B had a higher risk score (Fig. 4J).

Fig. 4
figure 4

Validation of the risk score model. A, B Survival curve of the train and the test group. C, D Time-dependent receiver operating characteristics (ROC) of the train and the test group. E Forest plot of independent prognostic analysis. F–H The decision curve analysis. I Heatmap of ARGs expression between two groups. J Risk score of the two groups

Full size image

Correlation of risk scores with immunotherapy response

Immune cell differential analysis suggested that immune cell infiltration was different between the two groups, such as B cells naive, T cells CD4 memory resting and Dendritic cells resting (Fig. 5A, B). The correlation heatmap of immune cells was shown in Fig. 5C. We also assessed the value of risk scores in predicting response to immunotherapy. The correlation between risk scores and B cells naive and Macrophages M0 was found to be positive, while a negative correlation was observed with Dendritic cells resting, T cells CD4 memory activated, and T cells CD8 (Fig. 5D). TME difference analysis suggested significant differences in TME scores between the two groups in StromalScore and ESTIMATEScore (Fig. 5E).

Fig. 5
figure 5

Correlation of risk scores with immunotherapy response. A, B Heatmap and vioplot based on two groups. C The correlation of risk score and immune infiltration cells. D Immune-related heatmap. E Vioplot of TME difference analysis

Full size image

Drug sensitivity analysis

To search for the possibility of ARGs as prognostic markers in the individualized treatment of BLCA, we evaluated the relationship between drug risk scores and drug sensitivity in the treatment of BLCA. As shown in Fig. 6, the sensitivity of 34 anti-tumor drugs was significantly different in two groups (p < 0.05). This implies a potential role for these drugs in the future treatment of BLCA.

Fig. 6
figure 6

The relationship between drug risk scores and drug sensitivity

Full size image

ScRNA-seq and immunohistochemical (IHC) data extraction from public databases

We obtained the gene expression data for two single-cell RNA sequencing samples from the GEO database using accession number GSE130001 (Fig. 7A-C). FASN was widely available in endothelial, epithelial, fibroblasts and myofibroblasts (Fig. 7 D, E). F10 was widely available in fibroblasts and myofibroblasts (Fig. 7 F, G). CALR was significantly expressed in endothelial, epithelial, fibroblasts and myofibroblasts (Fig. 7 H, I). RPS6KA1 was widely available in endothelial, epithelial and fibroblasts (Fig. 7 J, K). NTF3 was widely available in fibroblasts and myofibroblasts, minimally expressed in epithelial (Fig. 7 L, M). CSPG4 was widely available in myofibroblasts, minimally expressed in endothelial, epithelial and fibroblasts (Fig. 7 N, O). IHC data revealed a notable upregulation of ARG expression levels in BLCA tissues compared to normal tissues (Fig. 8).

Fig. 7
figure 7

Single cell RNA sequencing (scRNA-seq) data extraction based on Tumor Immune Single-cell Hub 2 (TISCH2) database (http://tisch.comp-genomics.org/). A-C Basic data of BLCA_GSE130001. D-O ScRNA-seq data of FASN, F10, CALR, RPS6KA1, NTF3, CSPG4

Full size image
Fig. 8
figure 8

Immunohistochemical (IHC) data extraction based on public database (https:// www. proteinatlas.org/). Representative immunohistochemical staining for CALR, CSPG4, FASN, NTF3, NTRK3, RPS6KA1 and S100A7 in normal tissue and BLCA tissue

Full size image

Discussion

Anoikis, a programmed cell death caused by cell detachment from the extracellular matrix, is an important mechanism that prevents cell growth and attachment to inappropriate substrates independent of adhesions, thereby avoiding colonization of distant organs [24,25,26]. The scientific community has shown significant interest in the capacity of tumor cells to evade anoikis, given its crucial role in tumor progression and metastatic establishment [24, 27]. Intracellular signals such as cellular DNA damage and stress in the endoplasmic reticulum have the potential to induce anoikis, with mitochondria playing a pivotal role in regulating apoptosis [28]. This malfunction in the implementation of anoikis could potentially serve as an indicator for the invasion and movement of cancer cells, while also playing a role in the creation of metastases in distant organs and the emergence of resistance to drugs [29]. Hence, comprehending the mechanism underlying resistance to anoikis can aid in impeding tumor advancement and averting metastasis [6, 30].

In our research, we discovered 37 prognostic ARGs. Through LASSO analysis, out of which 8 were chosen as predictive models for estimating the overall survival (OS) in patients with bladder cancer (BLCA). The BLCA patients were classified using the consensus clustering algorithm, resulting in the identification of two potential subgroups. We conducted a comparative analysis of the data from both groups, revealing contrasting survival durations among BLCA patients in each group. Furthermore, there were notable disparities observed in terms of immune cell infiltration and immune targets among the two cohorts. Next, we constructed risk models based on 8 ARGs as a way to predict the survival and prognosis of BLCA patients. Based on the outcomes of ROC curves and an independent prognostic analysis, our constructed prognostic model demonstrated excellent predictive accuracy.

Growing evidence showed that the existence of anti-ARGs was strongly associated with tumor aggressiveness and tumor drug resistance [9, 31,32,33]. A study indicated that lncRNA APOC1P1-3 could promote breast cancer metastasis by specifically binding miRNA-188-3p to block Bcl-2 inhibition through anoikis-resistance [34]. After stromal deprivation, the lymphocyte spectrum specific Ets transcription factor SPIB was activated and directly enhanced SNAP47 transcription in certain lung cancer cells. Loss of adhesion-induced autophagy significantly increased anoikis-resistance through SPIB activation [26]. Here, we found ARGs alterations at the transcriptional level and showed a reciprocal correlation of BLCA. Down-regulation of NTF3 and RPS6KA1 genes had a positive effect on the survival of BLCA patients, and up-regulation of the other 6 genes had a negative effect on the survival of BLCA patients. Egan et al. found significantly elevated CSPG4 expression in aggressive thyroid cancer, which was strongly associated with poor prognosis [35]. CALR mutations had an important role in myeloproliferative neoplasms [36, 37]. NTF3 also had a significant correlation with prognosis and immune infiltration in hepatocellular carcinoma [38]. A research study demonstrated that the activation of the PI3K-Akt-mTOR pathway was facilitated by the interaction between the FLT3-ITD mutation and the CXCL12/CXCR4 axis. This activation, in turn, resulted in an increase in RPS6KA1 expression which regulated the expression of genes associated with resistance to multiple drugs, thereby contributing to drug indications [39]. Magali et al. suggested that FASN could be used for the de novo synthesis of fatty acids in human adipogenic enzymes, which were highly expressed in cancer cells [40]. In glioma patients, hypomethylation of the F10 promoter was responsible for the overexpression and aggressive biological behavior of the F10-encoded protein FX [41]. The NTRK3 gene was involved in regulating the malignant behavior of malignant melanoma and might be a new therapeutic target for malignant melanoma [42]. S100A7 could promote esophageal squamous cell carcinoma progression by activating oncogenic pathways and remodeling the TME [43]. Therefore, exploring ARGs is more useful for developing more personalized immunotherapy programs.

Finally, we performed a differential analysis of TME and drug sensitivity based on risk score prognostic model. We observed that BLCA patients with lower risk scores had significantly longer OS. When studying treatment outcomes based on risk scores in BLCA patients, we found correlation between drug sensitivity and risk score. This discovery validates the dependability of employing this risk score model for prognosticating the effectiveness of immunotherapy in patients with BLCA. In general, as a new prognostic risk model, the studies on these eight ARGs are superficial and further research is necessary in the future. To further validate the risk model, general patient information, including age, sex, tumor stage, and score of sensitive genes, should be collected. Our single-cell sequencing data indicated that ARGs were expressed to varying degrees in different cells of BLCA patients. Furthermore, we offer a comprehensive analysis of the contrasting expression patterns of BLCA marker genes in both tumor and normal tissue. IHC showed that the expression level of ARGs in the BLCA tissue was significantly higher than that in the normal tissue. However, there are few reports on F10. Therefore, it is necessary to further investigate the ARGs F10 in future research.

Needless to say, the study still has some limitations. First, the investigation exclusively utilized openly accessible datasets, potentially introducing selection bias. Secondly, all the findings in this manuscript were of speculation based on the transcriptional analyses using bioinformatics analysis, the exact mechanism of ARGs in BLCA need to be further investigated in vivo and in vitro. Additionally, there is a deficiency of crucial experimental evidence to authenticate the manifestation of predictive genes in BLCA. Thus, it is necessary to further explore the specific mechanism of action of ARGs through molecular and animal experiments.

Conclusion

In this study, we systematically investigated the prognostic model of BLCA risk score based on 8 ARGs and linked these models to TME, drug sensitivity, scRNA-seq, and IHC. The prognostic model based on risk scoring demonstrated its efficacy in accurately forecasting the prognosis of patients with BLCA and their response to immunotherapy. At the same time, the systematic assessment of risk scores could help to design more individualized and precise treatment strategies for BLCA patients.

Availability of data and materials

The study data are available from the TCGA database (https://portal.gdc.cancer.gov/) and the GEO (GSE19423) database (https://www.ncbi.nlm.nih.gov/geo/). ARGs were obtained from Genecards database (https://www.genecards.org/) and Harmonizome database (https://maayanlab.cloud/Harmonizome/). Single-cell sequencing based on Tumor Immune Single-cell Hub 2 database (http://tisch.comp-genomics.org/). IHC based on Human Protein Atlas (https://www.proteinatlas.org/). The anoikis-related genes can be found in the Supplementary Material. Further inquiries can be directed to the corresponding author.

Abbreviations

BLCA:

Bladder urothelial carcinoma

scRNA-seq:

Single-cell RNA sequencing

ARGs:

Anoikis-related genes

TCGA:

The Cancer Genome Atlas

GEO:

Gene Expression Omnibus

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GO:

Gene ontology

PCA:

Principal component analysis

IC50:

Half-maximal inhibitory concentration prediction

ECM:

Extracellular matrix

LASSO:

Least absolute shrinkage and selection operator

TME:

Tumor microenvironment

IHC:

Immunohistochemical

References

  1. Richters A, Aben KKH, Kiemeney L. The global burden of urinary bladder cancer: an update. World J Urol. 2020;38(8):1895–904. https://doi.org/10.1007/s00345-019-02984-4.

    Article  PubMed  Google Scholar 

  2. Sanli O, Dobruch J, Knowles MA, Burger M, Alemozaffar M, Nielsen ME, et al. Bladder cancer. Nat Rev Dis Primers. 2017;3:17022. https://doi.org/10.1038/nrdp.2017.22.

    Article  PubMed  Google Scholar 

  3. Cambier S, Sylvester RJ, Collette L, Gontero P, Brausi MA, van Andel G, et al. EORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive stage Ta-T1 urothelial bladder cancer patients treated with 1–3 years of maintenance bacillus calmette-guerin. Eur Urol. 2016;69(1):60–9. https://doi.org/10.1016/j.eururo.2015.06.045.

    Article  PubMed  Google Scholar 

  4. Wilson F, Joseph N, Choudhury A. Biomarkers in muscle invasive bladder cancer. Adv Clin Chem. 2022;107:265–97. https://doi.org/10.1016/bs.acc.2021.07.005.

    Article  PubMed  CAS  Google Scholar 

  5. Chen S, Gu J, Zhang Q, Hu Y, Ge Y. Development of biomarker signatures associated with anoikis to predict prognosis in endometrial carcinoma patients. J Oncol. 2021;2021:3375297. https://doi.org/10.1155/2021/3375297.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Adeshakin FO, Adeshakin AO, Afolabi LO, Yan D, Zhang G, Wan X. Mechanisms for modulating anoikis resistance in cancer and the relevance of metabolic reprogramming. Front Oncol. 2021;11:626577. https://doi.org/10.3389/fonc.2021.626577.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Du S, Yang Z, Lu X, Yousuf S, Zhao M, Li W, et al. Anoikis resistant gastric cancer cells promote angiogenesis and peritoneal metastasis through C/EBPbeta-mediated PDGFB autocrine and paracrine signaling. Oncogene. 2021;40(38):5764–79. https://doi.org/10.1038/s41388-021-01988-y.

    Article  PubMed  CAS  Google Scholar 

  8. Zhu Z, Fang C, Xu H, Yuan L, Du Y, Ni Y, et al. Anoikis resistance in diffuse glioma: the potential therapeutic targets in the future. Front Oncol. 2022;12:976557. https://doi.org/10.3389/fonc.2022.976557.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Lee HY, Son SW, Moeng S, Choi SY, Park JK. The role of noncoding RNAs in the regulation of anoikis and anchorage-independent growth in cancer. Int J Mol Sci. 2021;22(2). https://doi.org/10.3390/ijms22020627.

  10. Ziegenhain C, Vieth B, Parekh S, Reinius B, Guillaumet-Adkins A, Smets M, et al. Comparative analysis of single-cell RNA sequencing methods. Mol Cell. 2017;65(4):631-643 e634. https://doi.org/10.1016/j.molcel.2017.01.023.

    Article  PubMed  CAS  Google Scholar 

  11. Kolodziejczyk AA, Kim JK, Svensson V, Marioni JC, Teichmann SA. The technology and biology of single-cell RNA sequencing. Mol Cell. 2015;58(4):610–20. https://doi.org/10.1016/j.molcel.2015.04.005.

    Article  PubMed  CAS  Google Scholar 

  12. Jovic D, Liang X, Zeng H, Lin L, Xu F, Luo Y. Single-cell RNA sequencing technologies and applications: a brief overview. Clin Transl Med. 2022;12(3):e694. https://doi.org/10.1002/ctm2.694.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3. https://doi.org/10.1093/bioinformatics/bts034.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Zhang H, Meltzer P, Davis S. RCircos: an R package for Circos 2D track plots. BMC Bioinformatics. 2013;14:244. https://doi.org/10.1186/1471-2105-14-244.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Jolliffe IT, Cadima J. Principal component analysis: a review and recent developments. Philos Trans A Math Phys Eng Sci. 2016;374(2065):20150202. https://doi.org/10.1098/rsta.2015.0202.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Pezzotti N, Lelieveldt BPF, Van Der Maaten L, Hollt T, Eisemann E, Vilanova A. Approximated and user steerable tSNE for progressive visual analytics. IEEE Trans Vis Comput Graph. 2017;23(7):1739–52. https://doi.org/10.1109/TVCG.2016.2570755.

    Article  PubMed  Google Scholar 

  18. Armstrong G, Martino C, Rahman G, Gonzalez A, Vazquez-Baeza Y, Mishne G, et al. Uniform manifold approximation and projection (UMAP) reveals composite patterns and resolves visualization artifacts in microbiome data. mSystems. 2021;6(5):e0069121. https://doi.org/10.1128/mSystems.00691-21.

    Article  PubMed  Google Scholar 

  19. Zhao Z, Liu H, Zhou X, Fang D, Ou X, Ye J, et al. Necroptosis-related lncRNAs: predicting prognosis and the distinction between the cold and hot tumors in gastric cancer. J Oncol. 2021;2021:6718443. https://doi.org/10.1155/2021/6718443.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Meng T, Huang R, Zeng Z, Huang Z, Yin H, Jiao C, et al. Identification of prognostic and metastatic alternative splicing signatures in kidney renal clear cell carcinoma. Front Bioeng Biotechnol. 2019;7:270. https://doi.org/10.3389/fbioe.2019.00270.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. https://doi.org/10.1038/ncomms3612.

    Article  PubMed  CAS  Google Scholar 

  22. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7. https://doi.org/10.1038/nmeth.3337.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Han Y, Wang Y, Dong X, Sun D, Liu Z, Yue J, et al. TISCH2: expanded datasets and new tools for single-cell transcriptome analyses of the tumor microenvironment. Nucleic Acids Res. 2023;51(D1):D1425–31. https://doi.org/10.1093/nar/gkac959.

    Article  PubMed  Google Scholar 

  24. Paoli P, Giannoni E, Chiarugi P. Anoikis molecular pathways and its role in cancer progression. Biochim Biophys Acta. 2013;1833(12):3481–98. https://doi.org/10.1016/j.bbamcr.2013.06.026.

    Article  PubMed  CAS  Google Scholar 

  25. Taddei ML, Giannoni E, Fiaschi T, Chiarugi P. Anoikis: an emerging hallmark in health and diseases. J Pathol. 2012;226(2):380–93. https://doi.org/10.1002/path.3000.

    Article  PubMed  CAS  Google Scholar 

  26. Zhang H, Wang G, Zhou R, Li X, Sun Y, Li Y, et al. SPIB promotes anoikis resistance via elevated autolysosomal process in lung cancer cells. FEBS J. 2020;287(21):4696–709. https://doi.org/10.1111/febs.15272.

    Article  PubMed  CAS  Google Scholar 

  27. Valentijn AJ, Zouq N, Gilmore AP. Anoikis. Biochem Soc Trans. 2004;32(Pt3):421–5. https://doi.org/10.1042/BST0320421.

    Article  PubMed  CAS  Google Scholar 

  28. Zhi Z, Ouyang Z, Ren Y, Cheng Y, Liu P, Wen Y, et al. Non-canonical phosphorylation of Bmf by p38 MAPK promotes its apoptotic activity in anoikis. Cell Death Differ. 2022;29(2):323–36. https://doi.org/10.1038/s41418-021-00855-3.

    Article  PubMed  CAS  Google Scholar 

  29. Jin L, Chun J, Pan C, Alesi GN, Li D, Magliocca KR, et al. Phosphorylation-mediated activation of LDHA promotes cancer cell invasion and tumour metastasis. Oncogene. 2017;36(27):3797–806. https://doi.org/10.1038/onc.2017.6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Xu R, Yan Y, Zheng X, Zhang H, Chen W, Li H, et al. Aspirin suppresses breast cancer metastasis to lung by targeting anoikis resistance. Carcinogenesis. 2022;43(2):104–14. https://doi.org/10.1093/carcin/bgab117.

    Article  PubMed  CAS  Google Scholar 

  31. Frankel A, Rosen K, Filmus J, Kerbel RS. Induction of anoikis and suppression of human ovarian tumor growth in vivo by down-regulation of Bcl-X(L). Cancer Res. 2001;61(12):4837–41.

    PubMed  CAS  Google Scholar 

  32. Toricelli M, Melo FH, Peres GB, Silva DC, Jasiulionis MG. Timp1 interacts with beta-1 integrin and CD63 along melanoma genesis and confers anoikis resistance by activating PI3-K signaling pathway independently of Akt phosphorylation. Mol Cancer. 2013;12:22. https://doi.org/10.1186/1476-4598-12-22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Markouli M, Strepkos D, Basdra EK, Papavassiliou AG, Piperi C. Prominent Role of Histone Modifications in the Regulation of Tumor Metastasis. Int J Mol Sci. 2021;22(5). https://doi.org/10.3390/ijms22052778.

  34. Lu Q, Wang L, Gao Y, Zhu P, Li L, Wang X, et al. lncRNA APOC1P1–3 promoting anoikis-resistance of breast cancer cells. Cancer Cell Int. 2021;21(1):232. https://doi.org/10.1186/s12935-021-01916-w.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Egan CE, Stefanova D, Ahmed A, Raja VJ, Thiesmeyer JW, Chen KJ, et al. CSPG4 is a potential therapeutic target in anaplastic thyroid cancer. Thyroid. 2021;31(10):1481–93. https://doi.org/10.1089/thy.2021.0067.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Kim HY, Han Y, Jang JH, Jung CW, Kim SH, Kim HJ. Effects of CALR-Mutant type and burden on the phenotype of myeloproliferative neoplasms. Diagnostics (Basel).2022;12(11). https://doi.org/10.3390/diagnostics12112570.

  37. Olschok K, Han L, de Toledo MAS, Bohnke J, Grasshoff M, Costa IG, et al. CALR frameshift mutations in MPN patient-derived iPSCs accelerate maturation of megakaryocytes. Stem Cell Rep. 2021;16(11):2768–83. https://doi.org/10.1016/j.stemcr.2021.09.019.

    Article  CAS  Google Scholar 

  38. Liu R, Li R, Yu H, Liu J, Zheng S, Li Y, et al. NTF3 correlates with prognosis and immune infiltration in hepatocellular carcinoma. Front Med (Lausanne). 2021;8:795849. https://doi.org/10.3389/fmed.2021.795849.

    Article  PubMed  Google Scholar 

  39. Yu DH, Chen C, Liu XP, Yao J, Li S, Ruan XL. Dysregulation of miR-138–5p/RPS6KA1-AP2M1 is associated with poor prognosis in AML. Front Cell Dev Biol. 2021;9:641629. https://doi.org/10.3389/fcell.2021.641629.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Humbert M, Seiler K, Mosimann S, Rentsch V, Sharma K, Pandey AV, et al. Reducing FASN expression sensitizes acute myeloid leukemia cells to differentiation therapy. Cell Death Differ. 2021;28(8):2465–81. https://doi.org/10.1038/s41418-021-00768-1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Liu X, Tang H, Wang Z, Huang C, Zhang Z, She X, et al. F10 gene hypomethylation, a putative biomarker for glioma prognosis. J Neurooncol. 2012;107(3):479–85. https://doi.org/10.1007/s11060-011-0775-2.

    Article  PubMed  CAS  Google Scholar 

  42. Zhou X, He J, Wang Q, Ma T. MiRNA-128–3p restrains malignant melanoma cell malignancy by targeting NTRK3. Front Oncol. 2020;10:538894. https://doi.org/10.3389/fonc.2020.538894.

    Article  PubMed  Google Scholar 

  43. Lu Z, Zheng S, Liu C, Wang X, Zhang G, Wang F, et al. S100A7 as a potential diagnostic and prognostic biomarker of esophageal squamous cell carcinoma promotes M2 macrophage infiltration and angiogenesis. Clin Transl Med. 2021;11(7):e459. https://doi.org/10.1002/ctm2.459.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

The authors thank all investigators and participants from the open datasets.

Funding

This research was funded by Shandong Natural Science Foundation Innovation and Development Joint Fund (Grant No. ZR2022LSW008).

Author information

Authors and Affiliations

  1. Department of Urology, Jining No. 1 People’s Hospital, Jining, Shandong, China

    Shusheng Zhu & Qingsong Zhao

  2. Department of Urology, The First Hospital of Jilin University, Changchun, Jilin, China

    Yanpeng Fan

  3. Department of Urology, Affiliated Yantai Yuhuangding Hospital, Qingdao University, Yantai, 264000, Shandong, China

    Chao Tang

Authors
  1. Shusheng Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Qingsong Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yanpeng Fan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Chao Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CT contributed to the conception and design; SZ contributed to the drafting of the manuscript; QZ contributed to the article revision. YF contributed to the acquisition, analysis, and interpretation of data. SZ and QZ contributed equally to this work.

Corresponding author

Correspondence to Chao Tang.

Ethics declarations

Ethics approval and consent to participate

The summary statistics are publicly downloadable in the website. All of these data are de-identified, freely downloadable, and can be used without restriction. Due to such a re‐analysis of previously collected and published data, no additional ethics approval was needed.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure S1.

Flowchart of the present research.

Additional file 2: Supplementary Table S1.

647 anoikis-related genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, S., Zhao, Q., Fan, Y. et al. Development of a prognostic model to predict BLCA based on anoikis-related gene signature: preliminary findings. BMC Urol 23, 199 (2023). https://doi.org/10.1186/s12894-023-01382-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12894-023-01382-8

Keywords

BMC Urology

ISSN: 1471-2490

Contact us