- Open Access
Prognostic value of immune-related lncRNA pairs in patients with bladder cancer
World Journal of Surgical Oncology volume 19, Article number: 304 (2021)
The characteristics of immune-related long non-coding ribonucleic acids (ir-lncRNAs), regardless of their specific levels, have important implications for the prognosis of patients with bladder cancer.
Based on The Cancer Genome Atlas database, original transcript data were analyzed. The ir-lncRNAs were obtained using a coexpression method, and their differentially expressed pairs (DE-ir-lncRNAs) were identified by univariate analysis. The lncRNA pairs were verified using a Lasso regression test. Thereafter, receiver operating characteristic curves were generated, and an optimal risk model was established. The clinical value of the model was verified through the analysis of patient survival rates, clinicopathological characteristics, presence of tumor-infiltrating immune cells, and chemotherapy efficacy evaluation.
In total, 49 pairs of DE-ir-lncRNAs were identified, of which 21 were included in the Cox regression model. A risk regression model was established on the premise of not involving the specific expression value of the transcripts.
The method and model used in this study have important clinical predictive value for bladder cancer and other malignant tumors.
The incidence and mortality of bladder cancer (BLCA) were approximately 500,000 and 200,000 worldwide, respectively, in 2020 . Muscle-invasive bladder cancer (MIBC) accounts for approximately 25% of patients with BLCA . Bacillus Calmette–Guérin, a type of Mycobacterium, reportedly prevents recurrence in patients with non-muscle–invasive bladder cancer, with the majority progressing to the MIBC subtype. With the development of immune checkpoint inhibitors (ICIs), patients who have been treated with pembrolizumab as second-line therapy during the KEYNOTE045 trial [3, 4] reportedly achieve approximately 10.3-month survival with an anti-tumor response (objective response rate) of 21.1%, which is greater than that in the chemotherapy group (11%). In addition, ICIs (atezolizumab and pembrolizumab) have been confirmed to be effective as first-line therapy, based on the results of NCT02108652  and KEYNOTE052  phase II clinical trials. Therefore, the European Medicines Agency (EMA) and the U.S. Food and Drug Agency have approved atezolizumab and pembrolizumab as first-line treatments for metastatic cisplatin-ineligible MIBC, restricted to cisplatin-unfit patients with PD-L1-high tumors. Although PD-L1 is a predictor of efficacy , other useful biomarkers related to ICIs for patients with BLCA need to be further explored to guide clinical practice.
Long non-coding RNAs (lncRNAs), with a transcript length of more than 200 nucleotides, are abundant, occupying more than 80 human transcripts . Recently, lncRNAs have been considered significant regulators of organic biological processes, including normal development and tumorigenesis. For example, the urothelial carcinoma-associated lncRNA (UCA1) , which is overexpressed in BLCA compared with healthy tissues, was reportedly associated with cisplatin sensitivity by modulating miR-196a-5p via the regulation of CREB. In addition, some lncRNAs have been reported to regulate the tumor microenvironment by targeting genes implicated in the function of immune cells [9,10,11]. Moreover, some immune-related lncRNA (ir-lncRNA) signatures have been recently identified in BLCA, and their expression is associated with the survival of patients with BLCA [12,13,14]. However, all these prognostic models were established based on the expression of lncRNA. Here, we established a novel model to predict the efficacy of immunotherapy regardless of expression.
RNA-seq data from The Cancer Genome Atlas (TCGA)-BLCA project were integrated into fragments per kilobase million (FPKM), and the GTF files were used to annotate and differentiate mRNAs and lncRNAs according to the Ensembl database (http://asia.ensembl.org). The ImmPort portal database (http://www.immport.org) was used to obtain confirmed immune-related genes and ir-lnRNAs by coexpression analysis.
Establishment of DE-ir-lncRNA pairs
The relationship between immune-related genes and all lncRNAs was verified by correlation tests; the highly correlated lncRNAs were considered ir-lncRNAs, with the cutoff value of correlation efficacy being > 0.5 and a P value of < 0.05. Thereafter, the R package “limma” (Bioconductor, USA) was used to detect differentially expressed lncRNAs (DE-lncRNAs), with the thresholds being defined as log fold change (FC) > 2, with a false discovery rate (FDR) < 0.05.
Lasso regression analysis and construction of Cox regression model
For DE-ir-lncRNA pairing, if one of two markers was highly expressed in a sample, the sample was regarded as a highly expressing sample of the two DE-ir-lncRNA markers. DE-ir-lncRNAs were tautologically paired, and a 0 or 1 matrix was constructed as per the following rule: considering that A is equal to lncRNA B plus lncRNA C, A is 1 if the expression of lncRNA B is higher than that of lncRNA C; if not, A is defined as 0. Then, the established matrix was filtered. Pairs were considered unrelated to prognosis as long as the expression value of the lncRNA pair was 0 or 1. DE-lncRNA pairs were deemed to be an applicable match when the expression value was > 20% of the total pairs. The least absolute shrinkage and selection operator (Lasso) regression model  was constructed with a P value of 0.05. The lasso regression was performed for 1000 cycles, and for each cycle, a random stimulation was set up 1000 times. Next, the frequency of each pair in the 1000-time-repeated lasso regression model was recorded, and pairs with frequencies > 100 times were selected for Cox proportional hazard regression analysis as well as the construction of the model. The area under the curve (AUC) of each model was calculated and plotted as a curve. If the curve reached the highest point, indicating the maximum AUC value, the calculation procedure was terminated, and the model was considered the optimal candidate.
We conducted a Kaplan–Meier analysis to validate the accuracy of the risk model using the following R packages : “survival,” “glmnet,” “pbapply,” “survivalROC,” “survminer,” and “heatmap.” In addition, the chi-squared test was used to analyze the relationship between the risk model and clinical characteristics, and the Wilcoxon test was used to evaluate risk score differences among the clinical groups.
Immune infiltration status analysis
We applied novel methods, including TIMER (http://cistrome.org/TIMER/), CIBERSORT, XCELL, QUANTISEQ, MCPcounter, and EPIC, to calculate the immune infiltration status of BLCA. The Wilcoxon signed-rank test was then applied to calculate the differences in infiltrating immune cells between the high- and low-risk groups. Subsequently, the relationship between the risk score values and the immune-infiltrated cells was evaluated using Spearman’s correlation analysis. The significance cutoff was set at P < 0.05, and the R package “ggplot2” was used for this analysis.
Finally, we calculated the half-maximal inhibitory concentration (IC50) of common chemotherapeutic drugs among patients with BCLA in the TCGA-BLCA project. The difference between the high- and low-risk groups was determined using the Wilcoxon test, and results were obtained using the R packages “pRRophetic” and “ggplot2.”
The chi-squared and Fisher’s exact tests were performed to detect the relationship between the risk score and clinical characteristics . The prognostic value of the risk model was assessed by determining the AUC of receiver operating characteristic (ROC) . R environment and Bioconductor packages (version 3.5.5) were used for statistical analysis, and P value < 0.05 was considered statistically significant .
Identification of DE-ir-lncRNAs and construction of two DE-ir-lncRNA pairs
A flow chart of the study is shown in Fig. 1. First, we identified the raw data for BLCA from the TCGA project, which comprised 19 normal and 411 tumor samples. Then, we annotated the transcriptome according to the Ensembl database. Consequently, a total of 1094 ir-lncRNAs were detected (Table S1), among which 109 were identified as DE-ir-lncRNAs, (14 downregulated and 95 upregulated; Fig. 2A). Ultimately, we constructed a 0 or 1 matrix to generate DE-ir-lncRNA pairs. In total, 4896 pairs were constructed, 251 pairs were identified using univariate analysis, and 49 DE-ir-lncRNA pairs were verified by lasso regression model analysis. Then, we established a multi-Cox regression model including 21 DE-lncRNA pairs using the forward method (Fig. 2B).
Establishing a risk assessment model and evaluating the relationship between the model and the prognosis of patients with BLCA
We calculated the AUCs for each ROC curve for the 21 DE-lncRNA pairs (Fig. 3A) and detected the optimal cutoff value, which referred to 1483 using the Akaike information criterion  (AIC) values  (Fig. 3B). Based on the cutoff point, we divided the patients into high- and low-risk groups. To validate the cutoff value, we delineated the 1-, 3-, and 5-year ROC curves, the AUC values of which were over 0.80 (Fig. 3C) and outlined the 5-year ROC curves with other clinical characteristics (Fig. 3D).
Evaluating the relationship between the risk assessment model and clinical characteristics
Based on the cutoff value previously defined, 156 and 244 patients were categorized into the high- and low-risk groups, respectively. The risk score and survival data of each patient are shown in Fig. 4A and B; this result confirms that the clinical prognosis of the low-risk group was superior to that of the high-risk group. Moreover, we observed that patients in the low-risk group had longer survival than those in the high-risk group, according to analysis using the Kaplan–Meier method (P < 0.001) (Fig. 4C). Subsequently, we conducted chi-squared tests to elucidate the relationship between the risk of BLCA and clinical characteristics. The ribbon chart and ladder diagrams established using the Wilcoxon signed-rank test showed that age (Fig. 5B), grade (Fig. 5C), and stage (Fig. 5D) were significantly associated with the risk group (P < 0.001). In addition, age (P < 0.01, hazard ratio (HR) = 1.026, 95% confidence interval (CI) [1.009–1.042]), stage (P < 0.001, HR = 1.564, 95% CI [1.280–1.912]), and risk score (P < 0.001, HR = 1.154, 95% CI [1.126–1.182]) were statistically significant as indicated by univariate Cox regression model analysis (Fig. 5E) and further verified by multivariate Cox regression analysis.
Relationship between the tumor microenvironment and the risk model
After establishing and verifying the risk model, we investigated whether the model was relevant to the tumor immune microenvironment. The high-risk group was more significantly associated with tumor-infiltrating immune cells, such as macrophages, neutrophils, and CD8+ T cells, but negatively associated with myeloid dendritic cells and CD4+ T cells, as verified using the Wilcoxon signed-rank test (Fig. 6A). As ICIs have been used to treat BLCA in clinical practice, we investigated whether the risk model was correlated with ICI-related biomarkers. Overall, high-risk scores were positively correlated with the high expression of discoidin domain receptor tyrosine kinase 2 (DDR2) (P < 0.05, Fig. 6C) and hepatitis A virus cellular receptor 2 (HAVCR2) (P < 0.001, Fig. 6D), whereas lymphocyte activating 3 (LAG3) (P > 0.05, Fig. 6E) and cytotoxic T lymphocyte associated protein 4 (CTLA4) (P > 0.05, Fig. 6A) showed no significant differences.
Relationship between the risk model and clinical chemotherapeutics
In addition to the aforementioned immunotherapy, we identified the potential relationship between the risk model and the efficacy of common chemotherapeutics in treating BLCA. The analysis of the TCGA-BLCA dataset revealed that a high-risk score was associated with a lower IC50 of chemotherapeutics, such as cisplatin (P = 0.00021, Fig. 7A), docetaxel (P < 0.0001, Fig. 7B), and paclitaxel (P < 0.0047, Fig. 7C). In contrast, we found that it was associated with a higher IC50 for metformin (P < 0.001, Fig. 7D) and methotrexate (P < 0.001, Fig. 7E). Collectively, these results demonstrate the predictive value of the proposed DE-lncRNA-based risk model.
It is widely known that RNA expression (including mRNA and lncRNA) have crucial biological functions [22,23,24]. Some of the DE-ir-lncRNAs detected in this study, such as TRPM2-AS, LINC01605, AC104041.1, and UCA1, have been confirmed to play significant roles in BLCA progression. Avgeris et al.  reported that the downregulation of UCA1 was correlated with a higher risk of short-term relapse in BLCA. Tian et al.  reported that TRPM2-AS promoted BLCA by targeting miR-22-3p and regulating the expression of GINS2. Qin et al.  revealed that high LINC01605 expression promoted BLCA progression by upregulating MMP9. Moreover, Lian et al.  established an 8-lncRNA signature, comprising APCDD1L-AS1, FAM225B, LINC00626, LINC00958, LOC100996694, LOC441601, LOC101928111, and ZSWIM8-AS1, as candidate prognostic markers for BLCA. Although various functions of lncRNAs have been proposed [29,30,31,32], single lncRNAs may be biased in predicting the prognosis of patients with BLCA. Furthermore, previous studies [33,34,35,36,37,38] have shown that the combinations of two genetic markers are more accurate than single genes in establishing prognostic cancer models. To date, few studies have confirmed the prognostic value of lncRNA pairs in this setting [39,40,41]. In the present study, we established a prognostic risk model by pairing immune-related genes and constructed a risk model with two lncRNA pairs without adopting their exact expression value. First, we screened the lncRNAs within the TCGA-BLCA dataset, selected the DE-lncRNAs, conducted a coexpression analysis for DE-ir-lncRNAs identification, and validated the obtained DE-ir-lncRNA pairs using a 0 or 1 matrix. Second, we applied a modified lasso penalized regression model, including the procedures of the cross, multiple repetitions of validation, and random stimulations to determine DE-ir-lncRNA pairs. Third, we delineated ROC curves and calculated the AUC values to acquire the optimized model. In addition, we calculated the AIC value of each point on the AUC to detect the best cutoff value to differentiate the high- and low-risk groups among patients with BLCA. Finally, we assessed the relationship between this novel risk model and different clinical parameters.
Preclinical studies have confirmed that increased infiltration of CD4+ or CD8+ immune cells [42,43,44] leads to a better response to ICIs. In the present study, we used various online tools, including CIBERSORT, XCELL, CIBERSORT-ABS, QUANTISEQ, MCPcounter, EPIC, and TIMER, to estimate the tumor-infiltrating cells in patients with BLCA, and analyzed their association with the predicted risk scores. Our results showed that CD4+ T cells, monocytes, macrophages, cancer-associated fibroblasts, and myeloid dendritic cells were enriched in the high-risk group, which may explain why the high-risk group was related to poor prognosis. In addition, correlation analysis demonstrated that the high-risk group was positively correlated with the expression of some immune microenvironmental inhibitory genes, such as HAVCR2 and DDR2, and it had a positive correlation trend with the expression of LAG3.
LINC00665 and some other lncRNAs have been shown to enhance the efficacy of immunotherapy in BLCA [45,46,47]. In addition, Zhang et al.  found that the lncRNA HOTAIR can inhibit 5-fluorouracil sensitivity by mediating MTHFR methylation, and Gu et al. reported that NONHSAT141924 was associated with paclitaxel chemotherapy resistance . Overall, these findings demonstrate that lncRNAs may be related to chemotherapy resistance. Based on this, herein, we explored the relationship between the identified risk group and chemotherapy. Our risk model suggested that the high-risk group was more sensitive to methotrexate and metformin, whereas the low-risk group was more sensitive to cisplatin, docetaxel, and paclitaxel, which was consistent with previous studies [50,51,52].
This study has some limitations. First, the raw data obtained from the TCGA database were relatively insufficient for an initial analysis. Second, external validation was necessary to verify the efficiency of the risk model established in this study. To overcome these limitations, we screened lncRNA pairs using a 0 or 1 matrix, which was optimal in this study. Further studies comprising more clinical samples are underway for further verification of the proposed model. In summary, we defined a novel risk predictive model comprising ir-lncRNAs that does not require the exact expression of the lncRNAs. This may help clinicians identify patients who can benefit from immunotherapy.
This study established a lncRNA pair model with the exact expression to predict the prognosis of patients with bladder cancer, which may have significant value for clinical practice.
Availability of data and materials
All data used in this study are publicly available in an online database.
Akaike information criterion
Area under the curve
Differentially expressed lncRNAs
Half-maximal inhibitory concentration
Immune checkpoint inhibitor
Long non-coding RNAs
Muscle-invasive bladder cancer
Receiver operating characteristic
The Cancer Genome Atlas
Urothelial carcinoma-associated lncRNA
Afonso J, Santos LL, Longatto-Filho A, Baltazar F. Competitive glucose metabolism as a target to boost bladder cancer immunotherapy. Nat Rev Urol. 2020;17(2):77–106. https://doi.org/10.1038/s41585-019-0263-6.
Autio KA, Boni V, Humphrey RW, Naing A. Probody therapeutics: an emerging class of therapies designed to enhance on-target effects with reduced off-tumor toxicity for use in immuno-oncology. Clin Cancer Res. 2020;26(5):984–9. https://doi.org/10.1158/1078-0432.Ccr-19-1457.
Fradet Y, Bellmunt J, Vaughn DJ, Lee JL, Fong L, Vogelzang NJ, et al. Randomized phase III KEYNOTE-045 trial of pembrolizumab versus paclitaxel, docetaxel, or vinflunine in recurrent advanced urothelial cancer: results of >2 years of follow-up. Ann Oncol. 2019;30(6):970–6. https://doi.org/10.1093/annonc/mdz127.
Linxweiler J, Junker K. Extracellular vesicles in urological malignancies: an update. Nat Rev Urol. 2020;17(1):11–27. https://doi.org/10.1038/s41585-019-0261-8.
Snyder A, Nathanson T, Funt SA, Ahuja A, Buros Novik J, Hellmann MD, et al. Contribution of systemic and somatic factors to clinical response and resistance to PD-L1 blockade in urothelial cancer: an exploratory multi-omic analysis. PLoS Med. 2017;14(5):e1002309. https://doi.org/10.1371/journal.pmed.1002309.
Jacqueline Vuky AVB, Bellmunt J, Fang X, Plimack ER, Castellano D, Godwin JL, et al. Long-term outcomes in KEYNOTE-052: phase II study investigating first-line pembrolizumab in cisplatin-ineligible patients with locally advanced or metastatic urothelial cancer. J Clin Oncol. 2020;38(23):2658–66.
Tran L, Xiao JF, Agarwal N, Duex JE, Theodorescu D. Advances in bladder cancer biology and therapy. Nat Rev Cancer. 2021;21(2):104–21. https://doi.org/10.1038/s41568-020-00313-1.
Pan J, Li X, Wu W, Xue M, Hou H, Zhai W, et al. Long non-coding RNA UCA1 promotes cisplatin/gemcitabine resistance through CREB modulating miR-196a-5p in bladder cancer cells. Cancer Lett. 2016;382(1):64–76. https://doi.org/10.1016/j.canlet.2016.08.015.
Chen C, He W, Huang J, Wang B, Li H, Cai Q, et al. LNMAT1 promotes lymphatic metastasis of bladder cancer via CCL2 dependent macrophage recruitment. Nat Commun. 2018;9(1):3826. https://doi.org/10.1038/s41467-018-06152-x.
Chen C, Luo Y, He W, Zhao Y, Kong Y, Liu H, et al. Exosomal long noncoding RNA LNMAT2 promotes lymphatic metastasis in bladder cancer. J Clin Investig. 2020;130(1):404–21. https://doi.org/10.1172/jci130892.
Chen X, Xie R, Gu P, Huang M, Han J, Dong W, et al. Long noncoding RNA LBCS inhibits self-renewal and chemoresistance of bladder cancer stem cells through eEpigenetic silencing of SOX2. Clin Cancer Res. 2019;25(4):1389–403. https://doi.org/10.1158/1078-0432.Ccr-18-1656.
Zhan Y, Chen Z, He S, Gong Y, He A, Li Y, et al. Long non-coding RNA SOX2OT promotes the stemness phenotype of bladder cancer cells by modulating SOX2. Mol Cancer. 2020;19(1):25. https://doi.org/10.1186/s12943-020-1143-7.
Zhan Y, Du L, Wang L, Jiang X, Zhang S, Li J, et al. Expression signatures of exosomal long non-coding RNAs in urine serve as novel non-invasive biomarkers for diagnosis and recurrence prediction of bladder cancer. Mol Cancer. 2018;17(1):142. https://doi.org/10.1186/s12943-018-0893-y.
Zheng R, Du M, Wang X, Xu W, Liang J, Wang W, et al. Exosome-transmitted long non-coding RNA PTENP1 suppresses bladder cancer progression. Mol Cancer. 2018;17(1):143. https://doi.org/10.1186/s12943-018-0880-3.
Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med. 1997;16:385–95.
Robert C, Gentleman VJC, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):80–6.
Ma W, Zhao F, Yu X, Guan S, Suo H, Tao Z, et al. Immune-related lncRNAs as predictors of survival in breast cancer: a prognostic signature. J Transl Med. 2020;18(1):442. https://doi.org/10.1186/s12967-020-02522-6.
Chen Q, Hu L, Huang D, Chen K, Qiu X, Qiu B. Six-lncRNA immune prognostic signature for cervical cancer. Front Genet. 2020;11:533628. https://doi.org/10.3389/fgene.2020.533628.
Lin Y, Pan X, Chen Z, Lin S, Chen S. Identification of an immune-related nine-lncRNA signature predictive of overall survival in colon cancer. Front Genet. 2020;11:318. https://doi.org/10.3389/fgene.2020.00318.
Posada D, Buckley TR. Model selection and model averaging in phylogenetics: advantages of Akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol. 2004;53(5):793–808. https://doi.org/10.1080/10635150490522304.
Vrieze SI. Model selection and psychological theory: a discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). Psychol Methods. 2012;17(2):228–43. https://doi.org/10.1037/a0027127.
Deng S, Wang W, Li X, Zhang P. Common genetic polymorphisms in pre-microRNAs and risk of bladder cancer. World J Surg Oncol. 2015;13:297. https://doi.org/10.1186/s12957-015-0683-6.
Cao W, Ma E, Zhou L, Yuan T, Zhang C. Exploring the FGFR3-related oncogenic mechanism in bladder cancer using bioinformatics strategy. World J Surg Oncol. 2017;15(1):66. https://doi.org/10.1186/s12957-017-1125-4.
Shi HB, Yu JX, Yu JX, Feng Z, Zhang C, Li GY, et al. Diagnostic significance of microRNAs as novel biomarkers for bladder cancer: a meta-analysis of ten articles. World J Surg Oncol. 2017;15(1):147. https://doi.org/10.1186/s12957-017-1201-9.
Avgeris M, Tsilimantou A, Levis PK, Rampias T, Papadimitriou MA, Panoutsopoulou K, et al. Unraveling UCA1 lncRNA prognostic utility in urothelial bladder cancer. Carcinogenesis. 2019;40(8):965–74. https://doi.org/10.1093/carcin/bgz045.
Tian Y, Guan Y, Su Y, Yang T, Yu H. TRPM2-AS promotes bladder cancer by targeting miR-22-3p and regulating GINS2 mRNA expression. Onco Targets Ther. 2021;14:1219–37. https://doi.org/10.2147/OTT.S282151.
Qin Z, Wang Y, Tang J, Zhang L, Li R, Xue J, et al. High LINC01605 expression predicts poor prognosis and promotes tumor progression via up-regulation of MMP9 in bladder cancer. Biosci Rep. 2018;38(5). https://doi.org/10.1042/BSR20180562.
Penghu Lian QW, Zhao Y, Chen C, Sun X, Li H, Deng J, et al. An eight-long non-coding RNA signature as a candidate prognostic biomarker for bladder cancer. Aging. 2019;11:6930–40.
Li L, Huang S, Yao Y, Chen J, Li J, Xiang X, et al. Follistatin-like 1 (FSTL1) is a prognostic biomarker and correlated with immune cell infiltration in gastric cancer. World J Surg Oncol. 2020;18(1):324. https://doi.org/10.1186/s12957-020-02070-9.
Zhou X, Xiao C, Han T, Qiu S, Wang M, Chu J, et al. Prognostic biomarkers related to breast cancer recurrence identified based on logit model analysis. World J Surg Oncol. 2020;18(1):254. https://doi.org/10.1186/s12957-020-02026-z.
Ding Y, Liu N, Chen M, Xu Y, Fang S, Xiang W, et al. Overexpressed pseudogene MT1L associated with tumor immune infiltrates and indicates a worse prognosis in BLCA. World J Surg Oncol. 2021;19(1):133. https://doi.org/10.1186/s12957-021-02231-4.
Hu X, Zhang L, Tian J, Ma J. Long non-coding RNA PART1 predicts a poor prognosis and promotes the malignant progression of pancreatic cancer by sponging miR-122. World J Surg Oncol. 2021;19(1). https://doi.org/10.1186/s12957-021-02232-3.
Hendriks RJ, Dijkstra S, Smit FP, Vandersmissen J, Van de Voorde H, Mulders PFA, et al. Epigenetic markers in circulating cell-free DNA as prognostic markers for survival of castration-resistant prostate cancer patients. Prostate. 2018;78(5):336–42. https://doi.org/10.1002/pros.23477.
Liu X, Wang J, Chen M, Liu S, Yu X, Wen F. Combining data from TCGA and GEO databases and reverse transcription quantitative PCR validation to identify gene prognostic markers in lung cancer. Onco Targets Ther. 2019;12:709–20. https://doi.org/10.2147/OTT.S183944.
Song D, Wang Y, Zhu K, Tian L, Gao Q, Zhou J, et al. DCK is a promising prognostic biomarker and correlated with immune infiltrates in hepatocellular carcinoma. World J Surg Oncol. 2020;18(1):176. https://doi.org/10.1186/s12957-020-01953-1.
Ding J, He X, Cheng X, Cao G, Chen B, Chen S, et al. A 4-gene-based hypoxia signature is associated with tumor immune microenvironment and predicts the prognosis of pancreatic cancer patients. World J Surg Oncol. 2021;19(1):123. https://doi.org/10.1186/s12957-021-02204-7.
Gu X, Zhang Q, Wu X, Fan Y, Qian J. Gene coexpression network approach to develop an immune prognostic model for pancreatic adenocarcinoma. World J Surg Oncol. 2021;19(1):112. https://doi.org/10.1186/s12957-021-02201-w.
Shen X, Yang Z, Feng S, Li Y. Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm. World J Surg Oncol. 2021;19(1):223. https://doi.org/10.1186/s12957-021-02333-z.
Wu Y, Zhang L, He S, Guan B, He A, Yang K, et al. Identification of immune-related LncRNA for predicting prognosis and immunotherapeutic response in bladder cancer. Aging. 2020;12(22):23306–25. https://doi.org/10.18632/aging.104115.
Robertson AG, Kim J, Al-Ahmadie H, Bellmunt J, Guo G, Cherniack AD, et al. Comprehensive molecular characterization of muscle-invasive bladder cancer. Cell. 2017;171(3):540–56 e525. https://doi.org/10.1016/j.cell.2017.09.007.
Zhou M, Zhang Z, Bao S, Hou P, Yan C, Su J, et al. Computational recognition of lncRNA signature of tumor-infiltrating B lymphocytes with potential implications in prognosis and immunotherapy of bladder cancer. Brief Bioinform. 2020. https://doi.org/10.1093/bib/bbaa047.
Pluhara GE, Pennell CA, Olin MR. CD8 T cell-independent immune-mediated mechanisms of antitumor activity +. Crit Rev Immunol. 2015;35(2):153–72.
Ostroumov D, Fekete-Drimusz N, Saborowski M, Kuhnel F, Woller N. CD4 and CD8 T lymphocyte interplay in controlling tumor growth. Cell Mol Life Sci. 2018;75(4):689–713. https://doi.org/10.1007/s00018-017-2686-7.
Tay RE, Richardson EK, Toh HC. Revisiting the role of CD4(+) T cells in cancer immunotherapy-new insights into old paradigms. Cancer Gene Ther. 2021;28(1-2):5–17. https://doi.org/10.1038/s41417-020-0183-x.
Zhang Y, Liu Q, Liao Q. Long noncoding RNA: a dazzling dancer in tumor immune microenvironment. J Exp Clin Cancer Res. 2020;39(1):231. https://doi.org/10.1186/s13046-020-01727-3.
Cong Z, Diao Y, Xu Y, Li X, Jiang Z, Shao C, et al. Long non-coding RNA linc00665 promotes lung adenocarcinoma progression and functions as ceRNA to regulate AKR1B10-ERK signaling by sponging miR-98. Cell Death Dis. 2019;10(2):84. https://doi.org/10.1038/s41419-019-1361-3.
Yang D, Feng W, Zhuang Y, Liu J, Feng Z, Xu T, et al. Long non-coding RNA linc00665 inhibits CDKN1C expression by binding to EZH2 and affects cisplatin sensitivity of NSCLC cells. Mol Ther Nucleic Acids. 2021;23:1053–65. https://doi.org/10.1016/j.omtn.2021.01.013.
Zhang S, Zheng F, Zhang L, Huang Z, Huang X, Pan Z, et al. LncRNA HOTAIR-mediated MTHFR methylation inhibits 5-fluorouracil sensitivity in esophageal cancer cells. J Exp Clin Cancer Res. 2020;39(1):131. https://doi.org/10.1186/s13046-020-01610-1.
Gu M, Zheng W, Zhang M, Dong X, Zhao Y, Wang S, et al. LncRNA NONHSAT141924 promotes paclitaxel chemotherapy resistance through p-CREB/Bcl-2 apoptosis signaling pathway in breast cancer. J Cancer. 2020;11(12):3645–54. https://doi.org/10.7150/jca.39463.
Patel VG, Oh WK, Galsky MD. Treatment of muscle-invasive and advanced bladder cancer in 2020. CA Cancer J Clin. 2020;70(5):404–23. https://doi.org/10.3322/caac.21631.
Mathes J, Rausch S, Todenhöfer T, Stenzl A. Trimodal therapy for muscle-invasive bladder cancer. Expert Rev Anticancer Ther. 2018;18(12):1219–29. https://doi.org/10.1080/14737140.2018.1535314.
Coen JJ, Zhang P, Saylor PJ, Lee CT, Wu CL, Parker W, et al. Bladder preservation with twice-a-day radiation plus fluorouracil/cisplatin or once daily radiation plus gemcitabine for muscle-invasive bladder cancer: NRG/RTOG 0712-A randomized phase II trial. J Clin Oncol. 2019;37(1):44–51. https://doi.org/10.1200/jco.18.00537.
We appreciate the support of the Dean of the Jiaxing Second Hospital (G. Chen).
Financial support was provided by the Jiaxing Science and Technology Bureau (2021AD3017 and 2020AD30084) and the Medical and Health Science and Technology Project of the Zhejiang Province (2021KY354).
Ethics approval and consent to participate
All analyses were based on publicly available online datasets; thus, no ethical approval and patient consent were required.
Consent for publication
The authors declare no potential conflicts of interest concerning the research, authorship, and/or publication of this article.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gao, Z., Wu, D., Zheng, W. et al. Prognostic value of immune-related lncRNA pairs in patients with bladder cancer. World J Surg Onc 19, 304 (2021). https://doi.org/10.1186/s12957-021-02419-8
- Bladder cancer
- Long non-coding RNA
- Risk model