Calculated identification of mutator-derived lncRNA signatures of genomic instability to predict the clinical outcome of muscle-invasive bladder cancer

Muscle-invasive bladder cancer (MIBC) is one of the most important type of bladder cancer, with a high morbidity and mortality rate. Studies have found that long non-coding RNA (lncRNA) plays a key role in maintaining genomic instability. However, Identification of lncRNAs related to genomic instability (GIlncRNAs) and their clinical significance in cancers have not been extensively studied yet. Here, we downloaded the lncRNA expression profiles, somatic mutation profiles and clinical related data in MIBC patients from The Cancer Genome Atlas (TCGA) database. A lncRNA computational framework was used to find differentially expressed GIlncRNAs. Multivariate Cox regression analysis was used to construct a genomic instability-related lncRNA signature (GIlncSig). Univariate and multivariate Cox analyses were used to assess the independent prognostic for the GIlncSig and other key clinical factors. We found 43 differentially expressed GIlncRNAs and constructed the GIlncSig with 6 GIlncRNAs in the training cohort. The patients were divided into two risk groups. The overall survival of patients in the high-risk group was lower than that in the low-risk group (P < 0.001), which were further verified in the testing cohort and the entire TCGA cohort. Univariate and multivariate Cox regression showed that the GIlncSig was an independent prognostic factor. In addition, the GIlncSig correlated with the genomic mutation rate of MIBC, indicating its potential as a measure of the degree of genomic instability. The GIlncSig was able to divide FGFR3 wild- and mutant-type patients into two risk groups, and effectively enhanced the prediction effect. Our study introduced an important reference for further research on the role of GIlncRNAs, and provided prognostic indicators and potential biological therapy targets for MIBC.

invasive bladder cancers progress to MIBC. The 5-year overall survival (OS) rate of MIBC patients is about 60% to 70% [2,3]. Males have a higher incidence, but female patients may already be in a more advanced stage at the time of diagnosis, resulting in a lower survival rate [4]. Advanced MIBC may spread to nearby lymph nodes or other organs, and even to distant organs, such as lungs and liver [5]. Although the current treatment methods include surgery, radiotherapy, chemotherapy, and immunotherapy, the recurrence rate of patients is relatively high and the survival rate of advanced patients is low [6][7][8]. Invasion and metastasis of MIBC are important factors for recurrence and poor prognosis [9,10]. Due to the complex molecular mechanism and clinical heterogeneity within the tumor, the therapeutic effect is different. Therefore, a new biomarker is urgently needed to more accurately assess the clinical prognosis of MIBC.
Studies have reported that genomic instability (GIN) and gene mutations are one of the hallmarks of cancer [11,12]. GIN can be observed in a variety of malignant tumors and pre-cancerous lesions [13]. The accumulation of GIN is closely related to tumor progression and survival, and is an important prognostic factor [14,15]. At the same time, the drug resistance of tumor is also related to GIN [16,17]. For example, chromosomal mutant of the epidermal growth factor receptor gene can make glioblastoma cells resistant to epidermal growth factor receptor inhibitors [17]. At the same time, studies have found that excessive GIN or certain gene mutation can also hinder cell growth and make cancer cells sensitive to drugs. This may be due to excessive genotoxicity and protein toxicity [18,19]. For example, cells with BRCA gene defects are sensitive to ionizing radiation [20,21]. These opposite effects, as well as the selective killing cancer cells of GIN, indicate that GIN is both a challenge and a potential opportunity for all kinds of cancer treatment [22,23].
Although we do not yet fully understand the molecular basis of GIN, abnormal transcription and post-transcriptional regulation are closely related to GIN [24]. Molecular markers have the potential to be able to quantitatively detect GIN. For example, Habermann et al. tested the gene expression profiles of 48 breast cancer patients and designed a GIN signature consisting of 12 genes [25]. Wang et al. constructed miRNA regulatory networks that is associated with DNA damage responses, and identified 10-miRNA signature related to outcome of ovarian cancer and GIN [26][27][28].
Long non-coding RNA (lncRNA) is a type of RNA that has no protein-coding ability and is usually more than 200 nt [29]. Studies have shown that it can regulate gene expression during or after transcription [30]. Many lncR-NAs have specific expression profiles in specific cancers, and the abnormal expression of lncRNAs may play a key role in the occurrence, proliferation, progression or metastasis of tumor [31]. More and more evidences show that lncRNA plays a vital role in maintaining GIN [32]. For example, Mendell et al. showed that DNA damage can activate specific lncRNAs, interact with proteins involved in DNA repair and replication, and then maintain genomic stability [33]. Betts et al. confirmed that two lncRNAs, CUPID1 and CUPID2, can regulate the expression of genes related to DNA repair [34]. LncRNA DDSR1 can maintain genome stability by regulating the expression of DNA damage genes or binding to DNA damage proteins [35,36]. Although existing studies have confirmed that several lncRNAs are related to genome stability, lncRNAs related to genomic instability (GIlncR-NAs) and their clinical guidance value in cancer are still largely unexplored.
In this study, we tried to construct a genomic instability-derived lncRNA signature (GIlncSig) combining lncRNA expression profiles and somatic mutation profiles, and to evaluate and verify its prognostic ability and independent prognostic value, so as to guide clinical application of the GIlncSig in MIBC.

Data collection
The clinical characteristics, RNA-seq expression data and somatic mutation information of MIBC patients were collected from the Cancer Genome Atlas (TCGA) database (https:// portal. gdc. cancer. gov/). A total of 399 cases with paired lncRNA and mRNA expression profiles, somatic mutation information and clinicopathological parameters including individual survival status, survival time, T stage, N stage, M stage, age and gender were obtained. All MIBC patients used in this study were randomly divided into two groups, named training cohort and testing cohort respectively. A total of 200 patients in the training cohort were used to establish a prognosticrelated risk model. The testing cohort (n = 199) and the entire TCGA cohort (n = 399) were used to verify the prognostic-related risk model. Table 1 shows brief information on clinical and pathological characteristics. We checked the relevant data of FGFR3 mutation patients on the cBioPortal.

Identification of genomic instability-associated lncRNAs
In order to identify lncRNAs related to genomic instability (GIlncRNAs), we used a calculation framework derived from the mutation hypothesis designed by Bao et al. [37]. It is roughly as follows, we calculated the number of somatic mutations in each patient and ranked the patients in descending order of the number of somatic mutations. The bottom 25% of patients were defined as the genomic stable (GS) group, and the top 25% of patients belonged to the genomic unstable (GU) group. The lncRNA expression profile of the two groups was compared using the 'limma' package of R version 4.0.2. Differentially expressed lncRNAs (fold change > 1, false discovery rate (FDR) adjustment P < 0.05) were defined as GIlncRNAs.

Identification of a genomic instability-derived lncRNA signature
Univariate Cox regression proportional hazard analysis was performed on the expression level of GIlncR-NAs to find the lncRNA related to prognosis in MIBC patients. Then these GIlncRNAs were included in the multivariate Cox regression to obtain the correlation coefficient. The GIlncSig was constructed to predict the outcome. The risk score was calculated by the The median value of patients in the training cohort was used as a risk cutoff point, patients were divided into high-risk groups and low-risk groups.

Pathway enrichment analysis
We used the 'limma' package in R to find mRNAs paired with GIlncRNAs, considered the top 10 mRNAs as lncRNA-related partners, and constructed a co-expression network. In order to predict the potential function of lncRNA, we used the 'clusterProfiler' and 'ggplot2' packages in R to perform Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on the mRNAs co-expressed with lncRNAs.

Cell culture and siRNA knockdown
We obtained bladder cancer (T24 and 5637) cell lines from the Type Culture Collection in Chinese Academy of Sciences (Shanghai, China). T24 cells were cultured in Dulbecco's modified Eagle medium (Gibco, Grand Island, NY, USA) and 5637 cells were cultured in 1640 medium (Gibco), supplemented with 10% foetal bovine serum (FBS, Gibco) in a humidified incubator at 37 °C with 5% CO2. We transfected si-CFAP58-DT, and respective negative controls into cultured T24 and 5637 cells using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) following standard procedure. Following transfection for 48 h, the follow-up experiments were carried out. The sequence for si-CFAP58-DT was as follows: CGC CTC TAA TTC TCC AAT ACATC.

Cell proliferation and apoptosis assays
We used the Cell Counting Kit 8 (Gibco) assay to assess cell proliferation. We placed 2000 transfected cells per well into 96-well plates in triplicate wells for 5 days and detected cell proliferation at 1, 2, 3, 4, and 5 d after seeding according to the standard manufacturer's instructions.
We used the Annexin V-APC/7-AAD Apoptosis Kit (MultiSciences; Hangzhou China) to detect apoptosis according to the manufacturer's instructions. In short, after washing in PBS, cells were incubated in 7-AAD staining solutions and APC Annexin-V in the dark at room temperature for 5 min. After incubation, the cells were collected by a C6 flow cytometer (BD, NY, USA) and FlowJo 10.4 software was used for analysis. Apoptotic cells population include 7-AAD-and APC Annexin-V + (undergoing apoptosis) cells and 7-AAD + and APC Annexin-V + (end of apoptosis or death) cells.

Statistical analysis
The 'survfit' and 'survdiff ' functions in the 'survival' package of R were used to compare the survival differences between the high-and low-risk groups. The Kaplan-Meier method was used to present the survival rate of each prognostic risk group. The performance of the GIlncSig was evaluated by a time-dependent receiver operating characteristic (ROC) curve. Univariate and multivariate Cox analyses were used to assess the independent prognostic for the GIlncSig and other key clinical factors. All statistical analyses were handled with R version 4.0.2.

Identification of genomic instability-related lncRNAs in MIBC patients
The Workflow in Fig. 1 is a brief summary of this research. In order to extract lncRNAs related to GIN, we divided the patients into GS-like group and GU-like group according to the cumulative number of somatic mutations in each patient. Then the lncRNA expression profiles between the two groups were compared. There were a total of 43 significantly differentially expressed lncRNAs (Additional file 1: Fig. S1a). All 399 samples were divided into two groups based on the expression levels of 43 differentially expressed lncRNAs (Fig. 2a). The group with low accumulation of somatic mutations was called GS-like group, and the other group was GUlike group. The median value of somatic accumulated mutations in the GS-like group was significantly lower than that in the GU-like group (P < 0.001, Fig. 2b). We then compared the expression levels of the UBQLN4 gene (a driver of GIN) in the GS-like and GU-like group. The expression of UBQLN4 in the GS-like group was significantly lower than that in the GU-like group (P < 0.001, Fig. 2c).

Enrichment analysis
In order to determine whether these lncRNAs were related to GIN, we performed a functional enrichment analysis of its related genes to predict potential functions. First, the top 10 mRNAs of 43 differentially expressed lncRNAs were obtained, and the lncRNAs-mRNA coexpression network was constructed (Fig. 2d). The KEGG pathway enrichment analysis was performed on the mRNA co-expressed by lncRNA. Seven significantly enriched pathways were found, most of which are related to GIN, including Alcoholism and Systemic lupus erythematosus (Fig. 2e). These results indicated that 43 differentially expressed lncRNAs were related to GIN, and these differentially expressed lncRNAs were defined as candidate GIlncRNAs.

Development of a GIlncSig for outcome prediction in the training cohort
To further research the prognostic effects of these GIl-ncRNAs, we randomly divided 399 MIBC samples into a training cohort (n = 200) and a testing cohort (n = 199). Univariate Cox proportional hazard regression was used to analyze the relationship between the expression levels of 43 GIlncRNAs and OS in the training cohort. It was found that the expression levels of 6 GIlncRNAs were significantly correlated with the prognosis of MIBC patients (P < 0.05, Fig. 3a). Then, multivariate Cox analysis was used to construct the GIlncSig. Table 2 summarized the detailed information of 6 lncRNAs including their Ensembel ID, coefficients and the results of multivariate Cox regression analysis. Among the GIlncSig, the coefficients of lncRNA LINC01419, AC087392.1, LINC02762 and CFAP58-DT were positive, indicating that they may be risk factors, and their high expression was related to a poor prognosis. While lncRNA LINC02446 and AC078880.3 may be protective factors, and their high expression showed a better prognosis. We used the above formula to obtain the risk score of each patient in the training cohort, and then divided these patients into a high-risk group and a low-risk group based on the threshold of the median risk score. Kaplan-Meier analysis indicated that the survival results of patients in the high-risk group were significantly lower than those in the low-risk group (P < 0.001, Fig. 3b). The time-dependent ROC curve analysis of the GIlncSig showed that the area under the curve (AUC) was 0.764 (Fig. 3c). We sorted the patients according to the risk score, and observed the expression level of the GIlncSig, the count of somatic mutations, and the risk score in the patients. For patients with high risk scores, the expression level of risk lncR-NAs were up-regulated, while the expression level of protective lncRNAs were down-regulated. Conversely, the patients with low scores showed the opposite result (Fig. 3d). Comparison analysis indicated that there were significant differences in the somatic mutation patterns of patients in the low-risk group and the high-risk group.  The number of somatic mutations in the high-risk group was significantly lower than that in the low-risk group (P = 0.008, Fig. 3e). The risk score, survival time and survival status of each patient in the training cohort were shown in Fig. 3f.

Validation of the GIlncSig in the testing cohort and the entire TCGA MIBC cohort
In order to examine the predictive ability and stability of the GIlncSig, we verified it in the testing cohort (n = 199) and the entire TCGA MIBC cohort (n = 399). The same GIlncSig and risk thresholds derived from the training cohort was applied to the testing cohort and the entire TCGA MIBC cohort. The 199 patients in the testing cohort were divided into a high-risk group (n = 100) and a low-risk group (n = 99). The overall survival of patients in the high-risk group were significantly lower than those in the low-risk group (P = 0.020, Fig. 4a). The time-dependent ROC curve analysis of the GIlncSig in the testing cohort showed that the AUC was 0.630 (Additional file 1: Fig. S1b). The distribution of somatic mutation counts and the expression level of the GIlncSig in the testing cohort samples were shown in Fig. 4c. The number of somatic mutations in the high-risk group and the low-risk group were similar (P = 0.97, Fig. 4e). The risk score, survival b Boxplots of somatic mutations in the GS-like group and GU-like group. The accumulation of somatic mutations in the GS-like group is significantly lower than that in the GU-like group. c Boxplots of the expression level of UBQLN4 in the GU-like group and GS-like group. The expression level of UBQLN4 in GS-like group was significantly lower than that in GU-like group. d Co-expression network of differential lncRNA related to genomic instability and their corresponding mRNAs. The blue circles represent lncRNAs, and the red circles represent mRNAs. e KEGG pathway analysis the main signalling pathways for mRNAs related to differential lncRNAs time and survival status of each patient in the testing cohort were shown in Fig. 4g. The prognostic and prediction performance of the GIl-ncSig in the entire TCGA MIBC cohort was similar to the previous results. Patients in the entire TCGA MIBC cohort were divided into a high-risk group (n = 200) and a low-risk group (n = 199). The survival time of patients in the high-risk group were significantly lower than those in the low-risk group (P < 0.001, Fig. 4b). The time-dependent ROC curve analysis of the GIlncSig in the entire TCGA MIBC cohort showed that the AUC was 0.686 (Additional file 1: Fig. S1c). The distribution of somatic mutation counts and the expression level of the GIlncSig in the entire TCGA MIBC cohort samples were shown in Fig. 4d. The number of somatic mutations in the high-risk group was lower than that in the low-risk group (P = 0.069, Fig. 4f ). The risk score, survival time and survival status of each patient in the entire TCGA MIBC cohort were shown in Fig. 4h.

Independent prognostic analysis of the GIlncSig and clinical relevance analysis
To assess the independent prognostic value of the GIlnc-Sig and its relationship with common clinical variables, multivariate Cox analyses were applied on age, gender, pathological grade, pathological stage, and the GIlnc-Sig risk score model. The results showed that the GIl-ncSig was significantly related to the OS in each cohort (Table 3). Therefore, it is determined that the GIlnc-Sig was an independent prognostic indicator for MIBC patients. Except for the GIlncSig, other indicators also have important prognostic values. In order to further   . 4 Evaluating the predictive power of the GIlncSig in the testing cohort and entire TCGA cohort. Kaplan-Meier survival curve estimates of overall survival of high-risk and low-risk groups in the testing cohort (a) and entire TCGA cohort (b). Expression heatmap of six lncRNAs and the distribution of somatic mutation in the testing cohort (c) and entire TCGA cohort (d). Boxplots of the distribution of somatic mutations in the testing cohort (e) and entire TCGA cohort (f). Distribution of risk score, survival outcome of patients based on data in the testing cohort (g) and entire TCGA cohort (h) verify, we carried out a stratified analysis of each index. Briefly, we divided the patients into young patients and elderly patients according to age, with 65 years old as the demarcation point. Similarily, the patients were divided into male and female groups by gender, divided into highgrade and low-grade groups according to pathological classification. Subsequently, the patients were divided into early group (stage T2) and advanced group (stage T3-4), without lymph node metastasis (stage N0) group and lymph node metastasis (stage N1-3) group, without distant metastasis (stage M0) group, and distant metastasis (stage M1) group. Then the GIlncSig was used to divide the patients in each group into high-risk or lowrisk groups. Finally, it was found that in the elderly patient group, male group, high-grade group, T3-4 group, N1-3 group and M0 group, the OS rate of the high-risk groups were significantly lower than those of the low-risk groups (P < 0.05; Fig. 5a-f ).

Combination of the GIlncSig and FGFR3 mutation status
FGFR3 mutations have an important value in bladder cancer [38]. From cBioPortal, we found that the OS of patients with FGFR3 mutations was better than that of wild-type ones (p < 0.05, Fig. 6a). Further analysis showed that in the training, testing and entire TCGA MIBC cohort, the proportion of FGFR3 mutation patients in the high-risk group was significantly lower than that in the low-risk group (Fig. 6b). The results demonstrated that the GIlncSig was also related to the FGFR3 mutation status and may be a mutation marker of the FGFR3. Therefore, we further tested whether the GIlncSig combined with FGFR3 mutation status could better predict the outcome. We used the GIlncSig and FGFR3 mutation status to divide patients into four risk groups, which were FGFR3 Mutation/high-risk group, FGFR3 Mutation/ low-risk group, FGFR3 Wide/high-risk group and FGFR3 Wide/low-risk group. The survival curves of the four groups were significantly different (P < 0.001, Fig. 6c). The FGFR3 Mutation/low-risk group had the best survival results. The results showed that the GIlncSig can effectively distinguish the survival results of FGFR3 Mutation and FGFR3 Wide patients respectively. At the same time, we tested the relationship between FGFR3 mRNA expression and the GIlncSig. It was found that the higher the expression level of FGFR3 mRNA, the lower the risk of patients (P = 0.019, Fig. 6d).

Discussion
GIN plays an important role in the occurrence, progression and recurrence of cancer. The pattern and degree of GIN are of great significance to the diagnosis and prognosis of tumors [39,40]. However, the quantitative measurement of the degree of GIN has always been a challenge. Evidence shows that the abnormal gene transcription and replication and inactivation of DNA repair can lead to GIN [41]. Although a large number of studies have reported that GIN plays an important role in the occurrence and development of tumors, studies have found that excessive GIN may hinder cancer cell growth and reduce the risk of disease [11]. This study found that the GIlncSig was correlated with the genomic mutation rate of MIBC, indicating its potential as a measure of GIN. Moreover, the number of mutations in the low-risk group was higher than that in the high-risk group, which provided an important reference for further judgment of the disease prognosis.
In recent years, lncRNA is an important part of tumor biology. It is closely related to the occurrence and progression of tumors, and can be used as a new type of prognostic marker [42,43]. For example, LncRNA GClnc1 is associated with the progression of bladder cancer [44]. In addition, lncRNAs are closely related to genome stability. Studies have found that lncRNA NORAD [45] and GUARDIN [46] play an important role in maintaining the stability of the genome. However, there are few studies on the relationship between lncRNA and cancer related to GIN. Therefore, we combined lncRNA expression profiles and somatic mutation profiles, found out GIlncRNAs, and then constructed a GIlncSig to predict the survival prognosis in MIBC. In the training cohort, the GIlncSig divides patients into high-and low-risk groups, which was verified on the testing and entire TCGA MIBC cohorts. This may be used as a promising prognostic indicator of MIBC, and may find specific GIlncRNAs in MIBC, which can provide an important reference for the occurrence, development and prognosis of the disease.
First, we screened 43 GIlncRNAs and performed a pathway enrichment analysis on their co-expressed mRNAs, which were mainly enriched in Alcoholism, Systemic lupus erythematosus. These pathways were related to gene instability [47][48][49][50]. The GIlncSig in this study was composed of 6 lncRNAs. LINC01419, AC087392.1, Fig. 6 Relationship between the GIlncSig and FGFR3 somatic mutation. a Kaplan-Meier survival curve for BCa patients stratified by the FGFR3 mutation. b The proportion of FGFR3 mutation in high-risk and low-risk groups in the training cohort, entire TCGA cohort and testing cohort. c Kaplan-Meier curve of overall survival is shown for MIBC patients classified according to FGFR3 mutation status and the GIlncSig. d FGFR3 expression in patients of low-risk and high-risk groups in the entire TCGA cohort LINC02762 and CFAP58-DT were risk factors and were related to poor prognosis, while LINC02446 and AC078880.3 were protective factors and were related to better prognosis. According to the results of Kaplan-Meier survival analysis, the GIlncSig had a strong prognostic performance and was significantly associated with the prognostic survival of MIBC patients. Using univariate and multivariate Cox regression analysis, we conducted independent prognostic analysis of risk factors such as age, gender, stage and the GIlncSig. It was found that risk score was an independent prognostic factor of MIBC. It indicates that it may be superior to other clinicopathological features and can be used as a better prognostic biomarker. Further stratified analysis showed that in the elderly patient group, male group, high-grade group, T3-4 group, N1-3 group and M0 group, the OS rate of the high-risk groups were significantly lower than those of the low-risk groups. The GIlncSig can predict the survival prognosis of patients in each subgroup more accurately.
Among the 6 lncRNAs in the GIlncSig, AC087392.1, AC078880.3 and LINC02762 have not been reported. LINC01419 can promote the growth and metastasis of hepatocellular carcinoma [51], gastric cancer cells [52] and osteosarcoma cells [53]. Xie et al. found that CFAP58-DT can enhance the innate immune response to viral infection [54], but it has not been reported in tumor treatment. Research by Zhang et al. showed that LINC02446 inhibits the proliferation and metastasis of bladder cancer cells [55]. We further analyzed the effect of CFAP58-DT on bladder cancer cells. The CCK-8 proliferation test showed that compared with si-control T24 and 5637 cells, the proliferation ability of CFAP58-DT-knockdown T24 and 5637 cells was significantly reduced. Flow cytometry found that the apoptosis rate of CFAP58-DT-knockdown T24 and 5637 cells was significantly increased. This finding indicates that CFAP58-DT may play a carcinogenic effect in bladder cancer, but the specific mechanism needs to be further studied. This model not only plays an important role in predicting the prognosis of MIBC, but also provides a new biomarker to guide the clinical treatment of MIBC patients.
We added a correlation analysis between the FGFR3 mutation and risk score to enhance the prognostic performance of our GIlncSig. This is not available in previous studies. A large number of studies have confirmed that FGFR3 mutation is a common genetic activity in bladder cancer and has become a promising molecular target for recurrence, prognosis, and treatment targets in bladder cancer [56,57]. FGFR3 mutations in bladder cancer may have better OS [38]. Consistent with previous studies, the FGFR3 mutation rate of patients in the high-risk group was significantly lower than that in the low-risk group, suggesting that the GIlncSig may be able to predict FGFR3 mutation status. In addition, we also observed that the GIlncSig can significantly distinguish the different clinical outcomes of FGFR3 wild-type and mutant-type patients and the FGFR3 Mutation/low-risk group has the best survival results. These results indicated that the combination of the GIlncSig may have better prognostic significance than FGFR3 mutation status alone. At the same time, we observed that the higher the expression level of FGFR3 mRNA, the lower the risk of patients, which may be due to the increased transcription of FGFR3 caused by FGFR3 mutation. The specific mechanism needs further study.
Although our research provides an effective method to better evaluate GIN and an important reference for the prognosis of MIBC patients. There are still some limitations. Firstly, our GIlncSig was constructed and verified in the TCGA data cohort, more independent data cohorts and biological experiments are needed to verify its robustness and reproducibility. Secondly, it is necessary to use prospective clinical data to verify its clinical applicability, to ensure its stability and reproducibility, and to explore its potential mechanism of action.

Conclusions
In conclusion, the GIlncSig was constructed by combining the lncRNA expression profiles, the somatic mutation profiles and clinical information of MIBC, which can be used as an independent prognostic marker for MIBC patients and successfully stratified the risk subgroups of MIBC patients. The combination of the GIlncSig and FGFR3 mutation status can improve the predictive effect of patients. Understanding the underlying mechanism and significance of the GIlncSig in MIBC can provide advice for determining the therapeutic target of MIBC.