Guoming Su
Qili He2†- 1Department of Pharmacy and Laboratory, Sichuan Nursing Vocational College, Chengdu, China
- 2Institute of Toxicological Detection, Sichuan Center for Disease Control and Prevention, Chengdu, China
- 3Guangdong Provincial Key Laboratory of Medical Molecular Diagnostics, Guangdong Medical University, Dongguan, China
Background: Increasing evidence shows that dysregulated expression of long non-coding RNAs (lncRNAs) can serve as diagnostic or prognostic markers in bladder cancer. The aim of this study was to evaluate the clinical values of dysregulated lncRNAs in bladder cancer.
Methods: Eligible studies were systematically searched in PubMed, Embase, and Web of Science databases from inception to December 2017. Odds ratios (OR) were calculated to investigate the correlation between lncRNAs and clinicopathological parameters. Pooled hazard ratios (HR) and 95% confidence interval (CI) were calculated to explore the prognostic value of lncRNAs in bladder cancer. Pooled diagnostic parameters were also calculated to estimate the performance of lncRNAs in diagnosing bladder cancer. All statistical analyses were performed by using STATA 13.1 program.
Results: A total of 37 relevant studies were included to the present systematic review according to the inclusion and exclusion criteria, including 26 on clinicopathological parameters, 19 on prognosis, and 7 on diagnosis. For clinicopathological parameters, MALAT1 expression was significantly associated with lymph node metastasis (OR = 2.731; 95% CI: 1.409–5.292; p = 0.003), and high-level expression of XIST was related to larger tumor size (OR = 2.473; 95% CI: 1.159–5.276; p = 0.019) and higher TNM stage (OR = 0.400; 95% CI, 0.184–0.868; p = 0.020). For the prognostic values, the most significant association was observed between increased expressions of SPRY4-IT1 and poor overall survival (OS) (HR = 3.716; 95% CI: 2.084–6.719; p < 0.001); high MALAT1 expression was significantly associated with poor OS (HR = 1.611; 95% CI: 1.076–2.412; p = 0.020). For the diagnostic values, UCA1 expression profile achieved a combined AUC of 0.92, with sensitivity of 0.84 and specificity of 0.89 in distinguishing patients with bladder cancer from non-cancerous controls.
Conclusions: In summary, systematic review elaborated that abnormal lncRNAs expression can serve as potential markers for prognostic evaluation in bladder cancer patients. In addition, the diagnostic meta-analysis concluded that abnormally expressed UCA1 can function as potential diagnostic markers for bladder cancer.
Introduction
Bladder cancer ranks as the ninth most frequently-diagnosed cancer worldwide and it is estimated that nearly 500,000 cases are diagnosed annually worldwide (Antoni et al., 2017). Despite improvements in current clinical treatment such as surgery, radiation therapy, and chemotherapy, 50–70% of patients are relapsed within the next 5 years (Terracciano et al., 2017). Therefore, it is urgent to find novel markers for diagnosis at early stage and identify effective therapeutic targets for improving the survival rate of patients with bladder cancer.
Long non-coding RNAs (lncRNAs) are generally defined as RNA transcripts longer than 200 nucleotides that lack an open reading frame. Recently, increasing evidences show that lncRNAs play important roles in various cancers, which influence all the “hallmarks of cancer” (Gutschner and Diederichs, 2012). It is reported that lncRNAs are involved in various cell biological processes, such as tumor cell proliferation, apoptosis, invasion, and metastasis (Hansji et al., 2014; Terracciano et al., 2017). So the aberrant expression patterns of lncRNAs are correlated with cancer diagnosis and prognosis and serve as predictors of patient outcomes. For example, LncRNA H19 expression was up-regulated and closely related to TNM cancer stages in patients with gastric cancer, which can serve as a potential non-invasive diagnostic biomarker in gastric cancer (Hashad et al., 2016). Sun et al. (2016) indicated that the lncRNA antisense non-coding RNA in the INK4 locus (ANRIL) was up-regulated in colorectal cancer tissues, which was associated with the survival rate of patients with colorectal cancer. LncRNA-activated by TGF-β (lncRNA-ATB) was significantly up-regulated in hepatocellular carcinoma metastases and associated with poor prognosis (Yuan et al., 2014).
Up to now, it was reported that some lncRNAs was aberrantly expressed in bladder cancer, such as HULC (Wang J. et al., 2017), MALAT1 (Li et al., 2017), and SNHG16 (Cao et al., 2017). Hu R. G. et al. (2017) found that the lncRNA cancer susceptibility candidate 8 (CASC8) was significantly down-regulated in bladder cancers and associated with the advanced stage of bladder cancer patients, overexpression of which remarkably suppressed the bladder cancer cell proliferation. Hepatocellular carcinoma up-regulated long non-coding RNA (HULC) promoted bladder cancer cells proliferation and inhibited apoptosis via regulation of ZIC2 and PI3K/AKT signaling pathway (Wang J. et al., 2017). LncRNA urothelial cancer-associated 1 (UCA1) promoted bladder cancer cell migration and invasion via hsa-miR-145/ZEB1/2/FSCN1 pathway (Xue et al., 2016). Recently several studies have investigated the prognostic and diagnostic value of lncRNAs in bladder cancer. However, most studies examining the clinical values of aberrantly expressed lncRNAs was limited by the small sample size or a single lncRNA. Therefore, we conducted a systematic review and meta-analysis to evaluate the clinicopathological, prognostic, and diagnostic roles of multiple lncRNAs expression in patients with bladder cancer.
Methods
Publication Search
The present systematic review was performed according to the PRISMA Statement (Moher et al., 2010) (see Table S1 PRISMA Checklist) and Cochrane Collaboration guidelines (http://handbook.cochrane.org/). We searched the Pubmed, Embase and Web of Science to identify relevant studies until December 21, 2017. The search strategies were based on combinations of the following key words: (“long non-coding RNA,” “lncRNA,” “lincRNA,” “long ncRNA,” “long intergenic non-coding RNA”) AND (“bladder”) AND (“cancer,” “carcinoma,” “neoplasm,” “tumor,” “tumors,” “tumor,” “tumors,” “malignancy,” “metastasis”). In addition, the references of eligible studies and relevant systematic reviews were checked for other eligible studies. We provided the detailed search strategies and results in the Table S2.
Selection Criteria
The included studies met the following criteria: (1) patients in the study were diagnosed with bladder cancer; (2) studies investigated the association between lncRNAs and bladder cancer; (3) sample size was no less than 40 cases; (4) for clinicopathological studies, the correlation between lncRNAs and clinicopathological parameters of patients with bladder cancer was performed, and the expression level of lncRNAs was divided into high or low groups; (5) for prognostic studies, the correlation between lncRNAs and survival was performed and the primary endpoints as overall survival (OS), disease free survival (DFS), cancer-specific survival (CSS) or recurrence-free survival (RFS) were clearly defined, then Kaplan–Meier survival curves or sufficient original data was provided to extract hazard ratio (HR) with 95% confidence interval (CI); (6) for diagnostic studies, diagnostic accuracy of lncRNAs for bladder cancer was performed, and sufficient data was provided for constructing the diagnostic two-by-two tables.
The exclusion criteria were: (1) overlapping or duplicate data; (2) lack of essential information; (3) letter, review article, case report, and conference abstract; and (4) non-English papers and non-human studies.
Data Quality Assessment and Extraction
Data were extracted by two authors independently from included studies using a predefined data extraction form. Then another author verified them and any discrepancies were resolved by consensus. The following information were collected: (1) basic information: first author's name, publication year, country, study design, patient population, lncRNAs, expression, sample size, tumor type, detected sample, detection method, and cutoff value; (2) clinicopathological parameters: gender, age, tumor size, tumor number, histological grade, TNM stage, tumor stage T, and lymph node metastasis; (3) prognostic information: follow-up months, outcome of survival analysis, and HR with corresponding 95% CI; (4) diagnostic information: sensitivity, specificity, area under the curve (AUC), sample sizes for diagnostic analysis, and data for two-by-two tables [true positive (TP), false positive (FP), true negative (TN), and false negative (FN)]. For studies that showed only Kaplan-Meier survival curve, HR with their 95% CI was calculated by using Engauge Digitizer version 4.1 and Tierney's method (Tierney et al., 2007).
We assessed the methodological quality of prognostic studies using the Newcastle-Ottawa Scale (NOS) tool that was extracted and modified from previous studies (Gao et al., 2017). The NOS scores ranged from 0 to 8, and a study with the higher scores indicated better methodological quality. Moreover, the Quality Assessment of Diagnostic Accuracy Studies (QUADAS) list was used to systematically assess the quality of all the included diagnostic studies (Whiting et al., 2003). Fourteen items from the QUADAS list were applied to each article, with an answer of “yes,” “no,” or “unclear.” The answer “yes” obtained a score of 1, whereas “no” or “unclear” gained a score of 0, and the full score was 14. If a cumulative score is higher than 8, the study will be considered as low risk of bias.
Statistical Analysis
Heterogeneity among the included studies was assessed by using the Cochrane's Q-test and I2 statistics. If heterogeneity (p < 0.05 or I2 > 50%) was statistically significant among studies, the random-effect model was chosen for the meta-analysis; otherwise, the fixed-effect model was used. Odds ratio (OR) with 95% CI was used to evaluate association between lncRNAs expression and clinicopathological parameters. Pooled HR with 95% CI was calculated to summarize the effect between lncRNAs and survival in patients with bladder cancer. For the diagnostic meta-analysis, correlated diagnostic accuracy indexes were computed as follows: sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR), diagnostic odds ratio (DOR), summary receiver operating characteristic (SROC) curve, and AUC. Publication bias was detected using Deeks' regression test of asymmetry (Deeks et al., 2005). All statistical analyses were performed by using STATA 13.1 program (Stata Corpotion, College Station, Texas, USA).
Results
Study Selection
As shown in the flow diagram (Figure 1), we identified 737 records in the electronic databases, including Pubmed, Embase, and Web of Science. Firstly, 283 duplicate records were excluded using EndNote X8. With the inclusion and exclusion criteria, 396 records were excluded by reviewing titles and abstracts. Subsequently, the 58 remaining full-text articles were assessed. Among 58 articles, 21 were excluded from the quantitative synthesis for the reasons depicted in Figure 1. No additional studies were identified by a manual search of the references of the original studies. Finally, the remaining 37 articles were eligible for the systematic review (Wang et al., 2006; He et al., 2013, 2016a,b; Fan et al., 2014; Li et al., 2014, 2017; Srivastava et al., 2014; Yan et al., 2014; Chen et al., 2015, 2017; Eissa et al., 2015a,b; Milowich et al., 2015; Zhao F. J. et al., 2015; Zhao X. L. et al., 2015; Duan et al., 2016; Iliev et al., 2016; Zhan et al., 2016a,b, 2017a,b; Zhang et al., 2016, 2017; Cao et al., 2017; Cui et al., 2017; Dudek et al., 2017; Du et al., 2017; Hu Y. Y. et al., 2017; Liao et al., 2017; Lin et al., 2017; Tolkach et al., 2017; Wu et al., 2017; Xiong et al., 2017; Yang et al., 2017; Ye et al., 2017; Zhuang et al., 2017), including 26 studies for clinicopathological parameters, 19 studies for prognosis, and 7 studies for diagnosis.
Figure 1. Flow diagram of study selection process. lncRNAs, long non-coding RNAs; OS, overall survival; DFS, disease free survival; CSS, cancer-specific survival; RFS, recurrence-free survival.
Clinicopathological Parameters
Table 1 summarized the main characteristics of studies on the association between lncRNAs and clinicopathological parameters. All the selected studies on clinicopathological parameters were from China, with 20/26 (76.9%) being published between 2016 and 2017. The systematic review of clinicopathological parameters was performed in 1,896 patients with bladder cancer, including urinary bladder cancer, urothelial carcinoma of the bladder, muscle-invasive bladder cancer, bladder transitional cell carcinomas, and non-muscle invasive bladder cancer. Twenty-four lncRNAs were described in the 26 studies involved in clinicopathological parameters. The expression of MALAT1 (Fan et al., 2014; Li et al., 2017), ASAP1-IT1 (Yang et al., 2017), SPRY4-IT1 (Zhao X. L. et al., 2015), lncRNA-n336928 (Chen et al., 2015), linc-UBC1 (He et al., 2013), SUMO1P3 (Zhan et al., 2016b), HNF1A-AS1 (Zhan et al., 2017b), CCEPR (Zhan et al., 2017a), linc00346 (Ye et al., 2017), XIST (Hu Y. Y. et al., 2017; Xiong et al., 2017), ZEB2-AS1 (Wu et al., 2017), ZEB1-AS1 (Lin et al., 2017), PVT1 (Cui et al., 2017), ABHD11-AS1 (Chen et al., 2017), SNHG16 (Cao et al., 2017), lncRNA-UNMIBC (Zhang et al., 2016), PANDAR (Zhan et al., 2016a), AATBC (Zhao F. J. et al., 2015), and HOTAIR (Yan et al., 2014) were up-regulated in bladder cancer patients, while the expression of lncRNA-LET (Zhuang et al., 2017), lncRNA-LOWEG (Liao et al., 2017), NBAT1 (Du et al., 2017), BANCR (He et al., 2016b), and MIR31HG (He et al., 2016a) were down-regulated. Only one study reported that down-regulated lncRNA-LOWEG were significantly associated with gender of patients (Liao et al., 2017). The results of these studies indicated that 24 lncRNAs were not significantly correlated with age of patients and tumor number. Two studies claimed that up-regulated ZEB2-AS1 (Wu et al., 2017) and XIST (Hu Y. Y. et al., 2017) were significantly related to tumor size. Dysregulated SPRY4-IT1 (Zhao X. L. et al., 2015), lncRNA-n336928 (Chen et al., 2015), SUMO1P3 (Zhan et al., 2016b), HNF1A-AS1 (Zhan et al., 2017b), CCEPR (Zhan et al., 2017a), ZEB1-AS1 (Lin et al., 2017), NBAT1 (Du et al., 2017), PVT1 (Cui et al., 2017), ABHD11-AS1 (Chen et al., 2017), PANDAR (Zhan et al., 2016a), AATBC (Zhao F. J. et al., 2015), and HOTAIR (Yan et al., 2014) were significantly associated with histological grade. Dysregulated ASAP1-IT1 (Yang et al., 2017), lncRNA-LET (Zhuang et al., 2017), XIST (Hu Y. Y. et al., 2017; Xiong et al., 2017), ABHD11-AS1 (Chen et al., 2017), SNHG16 (Cao et al., 2017), BANCR (He et al., 2016b), and MIR31HG (He et al., 2016a) were significantly associated with TNM stage. Furthermore, SPRY4-IT1 (Zhao X. L. et al., 2015), MALAT1 (Li et al., 2017), linc-UBC1 (He et al., 2013), lncRNA-LET (Zhuang et al., 2017), ZEB2-AS1 (Wu et al., 2017), XIST (Hu Y. Y. et al., 2017), PVT1 (Cui et al., 2017), and SNHG16 (Cao et al., 2017) were significantly associated with lymph node metastasis status in patients with bladder cancer.
Table 1. Summary of the comparison for the p-values of the association between lncRNAs and clinicopathological parameters.
Two lncRNAs (MALAT1 and XIST) were investigated in two studies, respectively. Therefore, we performed a meta-analysis to evaluate the association between MALAT1 and XIST and clinicopathological parameters. For the MALAT1, we combined two studies with a total of three groups according to different clinicopathological parameters (Figure 2). Heterogeneity was observed in two groups (Gender, I2 = 81.0%, p = 0.022; Tumor stage T, I2 = 50.8%, p = 0.154); therefore, a random effect model was used for the quantitative pooling. The results revealed that high MALAT1 expression was significantly associated with lymph node metastasis (OR = 2.731; 95% CI: 1.409–5.292; p = 0.003). However, expression of MALAT1 was not significantly associated with gender of patients (OR = 1.748; 95% CI: 0.440–6.951; p = 0.428) and tumor stage T (OR = 0.501; 95% CI: 0.225–1.120; p = 0.092). For the XIST, we combined two studies with a total of three groups according to different clinicopathological parameters (Figure 3). After analysis using fixed effect model, our results revealed that high expression of XIST was significantly associated with larger tumor size (OR = 2.473; 95% CI: 1.159–5.276; p = 0.019). In addition, high expression of XIST was related to higher TNM stage (OR = 0.400; 95% CI: 0.184–0.868; p = 0.020). However, expression of XIST was not significantly associated with tumor number (OR = 0.859; 95% CI: 0.413–1.783; p = 0.682). Publication bias was not assessed because only two studies investigated the same lncRNA MALAT1 or XIST that were pooled into the meta-analysis.
Figure 2. Forest plots of odds ratios (OR) for the association between MALAT1 expression and clinicopathological features in bladder cancer patients.
Figure 3. Forest plots of odds ratios (OR) for the association between XIST expression and clinicopathological features in bladder cancer patients.
Prognosis
Nineteen studies on prognosis were eligible for the final analysis, with 13/19 (68.4%) being published between 2016 and 2017. Most of the eligible studies were conducted in Chinese populations (84.2%), followed by German (5.3%), Czech (5.3%), and Netherlander (5.3%). Additionally, NOS scores indicated that 17 (89.5%) of the 19 eligible studies were not <7 (Table S3). Summary of lncRNAs used as prognostic biomarkers of bladder cancer was presented in Table 2. The systematic review of prognosis was performed in 1,604 patients with bladder cancer, including urinary bladder cancer, urothelial carcinoma of the bladder, muscle-invasive bladder cancer, bladder transitional cell carcinomas, and non-muscle invasive bladder cancer. Fifteen studies containing 17 lncRNAs were available to evaluate the relationship between abnormally expressed lncRNAs and OS of bladder cancer patients. HRs and their corresponding 95% CI were produced from the eligible studies. An observed HR <1 implied that the patients with increased lncRNAs expression had a better survival. Conversely, an observed HR > 1 implied that the patients with increased lncRNAs expression had a worse survival. Increased expressions of MALAT1 (Fan et al., 2014; Li et al., 2017), ASAP1-IT1 (Yang et al., 2017), SPRY4-IT1 (Zhao X. L. et al., 2015), TUG1 (Iliev et al., 2016), lncRNA-n336928 (Chen et al., 2015), GHET1 (Li et al., 2014), linc-UBC1 (He et al., 2013), YRNA5 (Tolkach et al., 2017), XIST (Hu Y. Y. et al., 2017), PVT1 (Cui et al., 2017), and SNHG16 (Cao et al., 2017) were significantly correlated with poor prognosis in OS, along with decreased expressions of CAT1297 (Dudek et al., 2017), lncRNA-LET, YRNA1 (Tolkach et al., 2017), YRNA3 (Tolkach et al., 2017), YRNA4 (Tolkach et al., 2017), and NBAT1 (Du et al., 2017). It was worth noting that the most significant association was observed between increased expressions of SPRY4-IT1 and poor OS (HR = 3.716; 95% CI: 2.084–6.719; p < 0.001). One study containing 4 lncRNAs were available to evaluate the relationship between abnormally expressed lncRNAs and CSS of bladder cancer patients. Increased expressions of YRNA4 (Tolkach et al., 2017) and YRNA5 (Tolkach et al., 2017) were significantly correlated with poor prognosis in CSS, along with decreased expressions of YRNA1 (Tolkach et al., 2017) and YRNA3 (Tolkach et al., 2017). Three studies containing 3 lncRNAs were available to evaluate the relationship between abnormally expressed lncRNAs and RFS of bladder cancer patients. Increased expressions of lncRNA-UNMIBC (Zhang et al., 2016) and HOTAIR (Yan et al., 2014) were significantly correlated with poor prognosis in RFS, along with decreased expression of MEG3 (Duan et al., 2016). It was also worth noting that the most significant association was observed between increased expressions of HOTAIR and poor RFS (HR = 4.712; 95% CI: 2.894–8.714; p < 0.001). Only one study was available to evaluate the relationship between abnormally expressed lncRNAs and DFS of bladder cancer patients. Decreased expression of GAS5 (Zhang et al., 2017) was significantly correlated with poor prognosis in DFS.
Table 2. Summary of lncRNAs used as prognostic biomarkers of bladder cancer.
Two studies investigated the relationship between the expression of MALAT1 and OS in a total number of 215 bladder cancer patients. Therefore, we carried out a meta-analysis on the association between abnormally expressed MALAT1 and the OS of bladder cancer patients. There was not statistically significant heterogeneity among studies (I2 = 29.9%, p = 0.232). After analysis using fixed effect model, our result suggested that high MALAT1 expression was significantly correlated with poor prognosis in OS (HR = 1.611; 95% CI: 1.076–2.412; p = 0.020) (Figure 4). Publication bias was not assessed because only two studies investigated the same lncRNA MALAT1 that were pooled into the meta-analysis.
Figure 4. Forest plots of studies evaluating hazard ratios of up-regulated MALAT1 expression and the overall survival (OS) of bladder cancer patients.
Diagnosis
In the diagnosis category, summary of lncRNAs used as diagnostic biomarkers of bladder cancer are presented in Table 3. Eight lncRNAs were described in the seven studies providing complete diagnostic data. The expression of SNHG16 (Duan et al., 2016), MALAT1 (Duan et al., 2016), and UCA1 (Wang et al., 2006) were up-regulated in bladder cancer tissue, while the expression of YRNA1 (Tolkach et al., 2017), YRNA3 (Tolkach et al., 2017), YRNA4 (Tolkach et al., 2017), YRNA5 (Tolkach et al., 2017), and MEG3 (Duan et al., 2016) were down-regulated. Diagnostic accuracy differed greatly between different lncRNAs tests. We found that UCA1 test had the highest sensitivity (91.5%) and specificity (96.5%) in the study conducted by Eissa et al. (Eissa et al., 2015a). Quality assessment of diagnostic studies included in the systematic review by the QUADAS tool was presented in Table S4, which showed that QUADAS scores of prognostic studies ranged from 9 to 12, indicating the high quality of included studies.
Table 3. Summary of lncRNAs used as diagnostic biomarkers of bladder cancer.
Since five studies investigated diagnostic value of UCA1 for bladder cancer, a meta-analysis was performed based on UCA1. There was statistically significant heterogeneity in pooled sensitivity (I2 = 79.41%, p < 0.01) and pooled specificity (I2 = 90.43%, p < 0.01), then, a random-effect model was chosen for the generation of pooled indexes. The pooled sensitivity of UCA1 was 0.84 (95% CI: 0.76–0.89), specificity was 0.89 (95% CI: 0.78–0.95), PLR was 7.81 (95% CI: 3.45–17.68), NLR was 0.18 (95% CI: 0.11–0.30), and DOR was 39.65 (95% CI: 10.40–151.12). Forest plots of pooled sensitivity and specificity for UCA1 was presented in Figure 5. The AUC of SROC curve based on summary sensitivity and specificity were 0.92 (95% CI: 0.89–0.94) (Figure 6), indicating a moderate accuracy for the diagnostic test. To obtain the post-test probability, we performed a simulation of an environment that had a prevalence of 20% for bladder cancer, with base on the included studies. Incorporating this evidence in a Fagan's nomogram (Figure 7), it appeared that the positive post-test probability was 66% and the negative post-test probability 4%. Deeks' funnel plot asymmetry test was performed to check publication bias in this meta-analysis. The result indicated that no significant bias was found (t = −0.59; p = 0.598). The shape of the funnel plot was presented in Figure 8, without any evidence of obvious asymmetry. Therefore, no obvious publication bias existed in the meta-analysis of diagnostic studies.
Figure 5. Forest plot of the sensitivity (Left) and specificity (Right) of UCA1 for the diagnosis of bladder cancer.
Figure 6. The summary receiver operator characteristic (SROC) curve based on UCA1. SECS, sensitivity; SPEC, specificity; AUC, area under the curve.
Figure 7. Fagan's nomogram for the calculation of post-test probability.
Figure 8. Funnel graph for the assessment of potential publication bias of the diagnostic studies.
Discussion
Increasing evidences have indicated that abnormally expressed lncRNAs were correlated with clinical outcomes for patients with bladder cancer in recent years. Multiple lncRNAs were highlighted as potential diagnostic and prognostic biomarkers for bladder cancer and shown to be potential new targets for cancer drugs. Eissa et al. (2015b) found that there was a significant difference between bilharzial benign and malignant cases regarding urinary lncRNA-UCA1 expression, and suggested that UCA1-nanoassay was a valid test for direct detection of urine UCA1 for bladder cancer detection. Fan et al. (2014) revealed that MALAT1 level was higher in primary tumors that subsequently metastasized than those in non-metastatic bladder cancer, and suggested that MALAT1 inhibition may represent a promising therapeutic option for suppressing bladder cancer progression. While Li et al. (2017) demonstrated that high tumor stage, positive lymph nodes, and high MALAT1 expression were independent prognostic indicators for OS of bladder cancer patients, and suggested that high MALAT1 expression could be considered as a potential therapeutic target of bladder cancer. However, most studies examining the clinical values of aberrantly expressed lncRNAs was limited by relatively small sample size or single lncRNA, which may result in inconsistent biological conclusions. Therefore, we performed a comprehensive systematic review and meta-analysis to systematically evaluate the clinical values of various lncRNAs in bladder cancer.
In the present systematic review, we investigated the relationship between multiple lncRNAs and clinicopathological parameters of bladder cancer. Most of the included studies suggested that multiple lncRNAs might be used as potential biomarkers of histological grade, TNM stage, tumor stage T, and lymph node metastasis. Many lncRNAs were identified in multiple different studies but only 2 (MALAT1 and XIST) were found to be studied in more than one study. For the MALAT1, Li et al. (2017) reported that high expression of MALAT1 was closely associated with higher probability of lymph node metastasis; but Fan et al. (2014) did not find statistically significant association between increased MALAT1 expression and lymph node metastasis. Pooled result showed that high expression of MALAT1 was significantly associated with lymph node metastasis, suggesting that MALAT1 can serve as a valuable biomarker for predicting lymph node metastasis status. For the XIST, high expression of XIST was related to larger tumor size and higher TNM stage after we pooled OR, suggesting that XIST can serve as a valuable biomarker for predicting tumor size and TNM stage in patients with bladder cancer. However, we did not provide enough information about bladder cancer patients. Some studies reported that patients were treated by surgery only and did not receive radiotherapy, chemotherapy, or other therapy before surgery; others did not mention treatment. Clinicopathological data such as gender, age, histological grade, TNM stage, or lymph nodal status were obtained at the same time. TNM stage was reported only in some studies and ranged from 0 to IV. LncRNAs expression was detected using by quantitative real-time PCR in tissue sample which were obtained and immediately frozen at the same time. Then study evaluated the association between lncRNAs and clinicopathological parameters. So whether the patients can be cured was not reported in the most studies included.
We explored the prognostic role of multiple lncRNAs in bladder cancer. For OS, we found that the increased expressions of 11 lncRNAs were related to poor prognosis in bladder cancer, while the decreased expressions of 6 lncRNAs were related to poor prognosis. Among them, SPRY4-IT1 (Zhao X. L. et al., 2015) exhibited the highest HR of 3.72, while NBAT1 (Du et al., 2017) exhibited the lowest HR of 0.41. Meta-analysis of different lcRNAs with prognosis was not performed, because we thought that meta-analysis of random lcRNAs doesn't make sense on a scientific level. However, MALAT1 was investigated in two studies (Fan et al., 2014; Li et al., 2017), so we carried out a meta-analysis on the association between abnormally expressed MALAT1 and the OS of bladder cancer patients. Our result suggested that high expression of MALAT1 was significantly correlated with poor prognosis in OS among patients with bladder cancer. A previous meta-analysis conducted by Tian and Xu (2015) reported that MALAT1 expression was an independent prognostic marker for OS in patients with cancer using univariate and multivariate analyses, those findings in consist with our results. Therefore, high MALAT1 expression can serve as an independent prognostic factor for OS of bladder cancer patients and can be considered as a potential therapeutic target of bladder cancer.
Diagnostic accuracy of multiple lncRNAs tests was explored in the present systematic review. Different lncRNAs tests differed in their sensitivity and specificity. Five studies investigated diagnostic value of lncRNA UCA1 for bladder cancer, but diagnostic accuracy differed greatly among these studies. Meta-analysis is a method of summarizing discrepant data on the accuracy of diagnostic tests. So a meta-analysis was performed based on UCA1. As a result, the overall pooled sensitivity and specificity of UCA1 for bladder cancer were 0.84 and 0.89, respectively, along with an AUC value of 0.92, suggesting that the diagnostic accuracy of UCA1 was moderate. However, there was statistically significant heterogeneity in pooled sensitivity and pooled specificity. We supposed that the patient populations, sample size, and the different cut-off value might be the potential source. But data were not enough to illustrate main source of heterogeneity. Thus, a meta-regression analysis is urgently needed to assess how study-specific attributes caused heterogeneity. The DOR shows the correlation between diagnostic efficiency and the disease, which has better discriminatory test performance with an extremely higher value (Glas et al., 2003). In current study, the DOR value was calculated to be 39.65, suggesting a moderate diagnostic accuracy of UCA1 for bladder cancer diagnosis. In addition, a lower NLR and a higher PLR values show a better diagnostic performance. In this study, the pooled PLR and NLR values were calculated to be 7.81 and 0.18 for UCA1, respectively, also suggesting a moderate diagnostic accuracy. These data suggested that UCA1 expression test showed a moderate diagnostic accuracy for bladder cancer. Therefore, UCA1 can be considered as a potential biomarker to assist in the diagnosis of bladder cancer. The diagnostic significance of UCA1 expression in bladder cancer has been investigated by meta-analysis in recent studies (Wang Z. et al., 2017; Zhen et al., 2017). Compared with a previous meta-analysis (Wang Z. et al., 2017), we found that there was similar sensitivity, but there was slightly higher specificity, DOR, and AUC in our results. The different in results may be due to the fact that non-English papers were excluded from our meta-analysis.
There are several limitations in our systematic review. Firstly, most of the included studies were conducted in Chinese populations, so the conclusion of this study might not be extended to all populations. Secondly, mutiple lncRNAs were used to evaluate the clinicopathological parameters, prognosis, and diagnosis of bladder cancer, so we may overestimate the clinical values of single lncRNA. Thirdly, HR and its 95% CI were extracted from Kaplan–Meier survival curves in seven studies, which may be less reliable than those directly acquired from survival data. Finally, although our searches were extensive and were not limited by language, language bias should not be completely avoided because of all included studies written in English.
In conclusion, systematic review elaborated that abnormal lncRNAs expression can serve as potential markers for prognostic evaluation in bladder cancer patients. More importantly, the diagnostic meta-analysis concluded that abnormally expressed UCA1 can function as potential diagnostic markers for bladder cancer. However, most lncRNAs were investigated in a single study, whose results were not enough to illustrate the clinical value of lncRNAs completely. So larger-size and higher-quality studies need to be conducted to validate the clinical value of single lncRNA in patients with bladder cancer.
Author Contributions
GS and QH designed this study. GS, QH, and JW participated in study selection and data extraction. GS, QH, and JW performed statistical analysis. GS and QH wrote and reviewed the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer VC and handling Editor declared their shared affiliation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00652/full#supplementary-material
References
Antoni, S., Ferlay, J., Soerjomataram, I., Znaor, A., Jemal, A., and Bray, F. (2017). Bladder cancer incidence and mortality: a global overview and recent trends. Eur. Urol. 71, 96–108. doi: 10.1016/j.eururo.2016.06.010
Cao, X., Xu, J., and Yue, D. (2017). LncRNA-SNHG16 predicts poor prognosis and promotes tumor proliferation through epigenetically silencing p21 in bladder cancer. Cancer Gene Ther. 25, 10–17. doi: 10.1038/s41417-017-0006-x
Chen, M., Li, J., Zhuang, C., and Cai, Z. (2017). Increased lncRNA ABHD11-AS1 represses the malignant phenotypes of bladder cancer. Oncotarget 8, 28176–28186. doi: 10.18632/oncotarget.14945
Chen, T., Xie, W., Xie, L., Sun, Y., Zhang, Y., Shen, Z., et al. (2015). Expression of long noncoding RNA lncRNA-n336928 is correlated with tumor stage and grade and overall survival in bladder cancer. Biochem. Biophys. Res. Commun. 468, 666–670. doi: 10.1016/j.bbrc.2015.11.013
Cui, Y., Liu, L. F., Chen, J. B., Li, C., Cheng, X., and Zu, X. B. (2017). Increased expression of long non-coding RNA PVT1 correlates with clinical progression and poor prognosis in bladder cancer. Int. J. Clin. Exp. Pathol. 10, 3265–3271.
Deeks, J. J., Macaskill, P., and Irwig, L. (2005). The performance of tests of publication bias and other sample size effects in systematic reviews of diagnostic test accuracy was assessed. J. Clin. Epidemiol. 58, 882–893. doi: 10.1016/j.jclinepi.2005.01.016
Du, D., Zhang, Z., Yin, Y. F., Min, M. Y., Du, J. H., Fu, J., et al. (2017). Decreased expression of long noncoding RNA NBAT1 indicates a poor prognosis and promotes cell proliferation and EMT in bladder cancer. Int. J. Clin. Exp. Med. 10, 9214–9221.
Duan, W., Du, L., Jiang, X., Wang, R., Yan, S., Xie, Y., et al. (2016). Identification of a serum circulating IncRNA panel for the diagnosis and recurrence prediction of bladder cancer. Oncotarget 7, 78850–78858. doi: 10.18632/oncotarget.12880
Dudek, A. M., Boer, S. J., Boon, N., Witjes, J. A., Kiemeney, L., and Verhaegh, G. W. (2017). Identification of long non-coding RNAs that stimulate cell survival in bladder cancer. Oncotarget 8, 34442–34452. doi: 10.18632/oncotarget.16284
Eissa, S., Matboli, M., Essawy, N. O., and Kotb, Y. M. (2015a). Integrative functional genetic-epigenetic approach for selecting genes as urine biomarkers for bladder cancer diagnosis. Tumour Biol. 36, 9545–9552. doi: 10.1007/s13277-015-3722-6
Eissa, S., Matboli, M., Essawy, N. O., Shehta, M., and Kotb, Y. M. (2015b). Rapid detection of urinary long non-coding RNA urothelial carcinoma associated one using a PCR-free nanoparticle-based assay. Biomarkers 20, 212–217. doi: 10.3109/1354750X.2015.1062918
Fan, Y., Shen, B., Tan, M., Mu, X., Qin, Y., Zhang, F., et al. (2014). TGF-beta-induced upregulation of malat1 promotes bladder cancer metastasis by associating with suz12. Clin. Cancer Res. 20, 1531–1541. doi: 10.1158/1078-0432.CCR-13-1455
Gao, Q., Xie, H., Zhan, H., Li, J., Liu, Y., and Huang, W. (2017). Prognostic values of long noncoding RNA GAS5 in various carcinomas: an updated systematic review and meta-analysis. Front. Physiol. 8:814. doi: 10.3389/fphys.2017.00814
Glas, A. S., Lijmer, J. G., Prins, M. H., Bonsel, G. J., and Bossuyt, P. M. (2003). The diagnostic odds ratio: a single indicator of test performance. J. Clin. Epidemiol. 56, 1129–1135. doi: 10.1016/S0895-4356(03)00177-X
Gutschner, T., and Diederichs, S. (2012). The hallmarks of cancer: a long non-coding RNA point of view. RNA Biol. 9, 703–719. doi: 10.4161/rna.20481
Hansji, H., Leung, E. Y., Baguley, B. C., Finlay, G. J., and Askarian-Amiri, M. E. (2014). Keeping abreast with long non-coding RNAs in mammary gland development and breast cancer. Front. Genet. 5:379. doi: 10.3389/fgene.2014.00379
Hashad, D., Elbanna, A., Ibrahim, A., and Khedr, G. (2016). Evaluation of the role of circulating long non-coding RNA H19 as a promising novel biomarker in plasma of patients with gastric cancer. J. Clin. Lab. Anal. 30, 1100–1105. doi: 10.1002/jcla.21987
He, A., Chen, Z., Mei, H., and Liu, Y. (2016a). Decreased expression of LncRNA MIR31HG in human bladder cancer. Cancer Biomarkers 17, 231–236. doi: 10.3233/CBM-160635
He, A., Liu, Y., Chen, Z., Li, J., Chen, M., Liu, L., et al. (2016b). Over-expression of long noncoding RNA BANCR inhibits malignant phenotypes of human bladder cancer. J. Exp. Clin. Cancer Res. 35:125. doi: 10.1186/s13046-016-0397-9
He, W., Cai, Q., Sun, F., Zhong, G., Wang, P., Liu, H., et al. (2013). linc-UBC1 physically associates with polycomb repressive complex 2 (PRC2) and acts as a negative prognostic factor for lymph node metastasis and survival in bladder cancer. Biochim. Biophys. Acta 1832, 1528–1537. doi: 10.1016/j.bbadis.2013.05.010
Hu, R., Zhong, P., Xiong, L., and Duan, L. (2017). Long noncoding RNA cancer susceptibility candidate 8 suppresses the proliferation of bladder cancer cells via regulating glycolysis. DNA Cell Biol. 36, 767–774. doi: 10.1089/dna.2017.3785
Hu, Y., Deng, C., Zhang, H., Zhang, J., Peng, B., and Hu, C. (2017). Long non-coding RNA XIST promotes cell growth and metastasis through regulating miR-139-5p mediated Wnt/beta-catenin signaling pathway in bladder cancer. Oncotarget 8, 94554–94568. doi: 10.18632/oncotarget.21791
Iliev, R., Kleinova, R., Juracek, J., Dolezel, J., Ozanova, Z., Fedorko, M., et al. (2016). Overexpression of long non-coding RNA TUG1 predicts poor prognosis and promotes cancer cell proliferation and migration in high-grade muscle-invasive bladder cancer. Tumour Biol. 37, 13385–13390. doi: 10.1007/s13277-016-5177-9
Li, C., Cui, Y., Liu, L. F., Ren, W. B., Li, Q. Q., Zhou, X., et al. (2017). High expression of long noncoding RNA MALAT1 indicates a poor prognosis and promotes clinical progression and metastasis in bladder cancer. Clin. Genitourin. Cancer 15, 570–576. doi: 10.1016/j.clgc.2017.05.001
Li, L. J., Zhu, J. L., Bao, W. S., Chen, D. K., Huang, W. W., and Weng, Z. L. (2014). Long noncoding RNA GHET1 promotes the development of bladder cancer. Int. J. Clin. Exp. Pathol. 7, 7196–7205.
Liao, X. H., He, A. B., Zhong, J. H., Wang, F., Mei, H. B., Wu, J. T., et al. (2017). Over-expression of long noncoding RNA LOWEG inhibits cell migration in human bladder cancer. Int. J. Clin. Exp. Pathol. 10, 7249–7255.
Lin, J., Zhan, Y., Liu, Y., Chen, Z., Liang, J., Li, W., et al. (2017). Increased expression of ZEB1-AS1 correlates with higher histopathological grade and promotes tumorigenesis in bladder cancer. Oncotarget 8, 24202–24212. doi: 10.18632/oncotarget.15527
Milowich, D., Le Mercier, M., De Neve, N., Sandras, F., Roumeguere, T., Decaestecker, C., et al. (2015). Diagnostic value of the UCA1 test for bladder cancer detection: a clinical study. Springerplus 4:349. doi: 10.1186/s40064-015-1092-6
Moher, D., Liberati, A., Tetzlaff, J., and Altman, D. G. (2010). Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Int. J. Surg. 8, 336–341. doi: 10.1016/j.ijsu.2010.02.007
Srivastava, A. K., Singh, P. K., Rath, S. K., Dalela, D., Goel, M. M., and Bhatt, M. L. (2014). Appraisal of diagnostic ability of UCA1 as a biomarker of carcinoma of the urinary bladder. Tumour Biol. 35, 11435–11442. doi: 10.1007/s13277-014-2474-z
Sun, Y., Zheng, Z. P., Li, H., Zhang, H. Q., and Ma, F. Q. (2016). ANRIL is associated with the survival rate of patients with colorectal cancer, and affects cell migration and invasion in vitro. Mol. Med. Rep. 14, 1714–1720. doi: 10.3892/mmr.2016.5409
Terracciano, D., Ferro, M., Terreri, S., Lucarelli, G., D'Elia, C., Musi, G., et al. (2017). Urinary long noncoding RNAs in nonmuscle-invasive bladder cancer: new architects in cancer prognostic biomarkers. Transl. Res. 184, 108–117. doi: 10.1016/j.trsl.2017.03.005
Tian, X., and Xu, G. (2015). Clinical value of lncRNA MALAT1 as a prognostic marker in human cancer: systematic review and meta-analysis. BMJ Open 5:e008653. doi: 10.1136/bmjopen-2015-008653
Tierney, J. F., Stewart, L. A., Ghersi, D., Burdett, S., and Sydes, M. R. (2007). Practical methods for incorporating summary time-to-event data into meta-analysis. Trials 8:16. doi: 10.1186/1745-6215-8-16
Tolkach, Y., Stahl, A. F., Niehoff, E. M., Zhao, C. M., Kristiansen, G., Muller, S. C., et al. (2017). YRNA expression predicts survival in bladder cancer patients. BMC Cancer 17:749. doi: 10.1186/s12885-017-3746-y
Wang, J., Ma, W., and Liu, Y. (2017). Long non-coding RNA HULC promotes bladder cancer cells proliferation but inhibits apoptosis via regulation of ZIC2 and PI3K/AKT signaling pathway. Cancer Biomark 20, 425–434. doi: 10.3233/cbm-170188
Wang, X. S., Zhang, Z., Wang, H. C., Cai, J. L., Xu, Q. W., Li, M. Q., et al. (2006). Rapid identification of UCA1 as a very sensitive and specific unique marker for human bladder carcinoma. Clin. Cancer Res. 12, 4851–4858. doi: 10.1158/1078-0432.CCR-06-0134
Wang, Z., Wang, X., Zhang, D., Yu, Y., Cai, L., and Zhang, C. (2017). Long non-coding RNA urothelial carcinoma-associated 1 as a tumor biomarker for the diagnosis of urinary bladder cancer. Tumour Biol. 39:1010428317709990. doi: 10.1177/1010428317709990
Whiting, P., Rutjes, A. W., Reitsma, J. B., Bossuyt, P. M., and Kleijnen, J. (2003). The development of QUADAS: a tool for the quality assessment of studies of diagnostic accuracy included in systematic reviews. BMC Med. Res. Methodol. 3:25. doi: 10.1186/1471-2288-3-25
Wu, X., Yan, T., Wang, Z., Wu, X., Cao, G., and Zhang, C. (2017). LncRNA ZEB2-AS1 promotes bladder cancer cell proliferation and inhibits apoptosis by regulating miR-27b. Biomed. Pharmacother. 96, 299–304. doi: 10.1016/j.biopha.2017.08.060
Xiong, Y., Wang, L., Li, Y., Chen, M. F., He, W., and Qi, L. (2017). The long non-coding RNA XIST interacted with MiR-124 to modulate bladder cancer growth, invasion and migration by targeting Androgen Receptor (AR). Cell. Physiol. Biochem. 43, 405–418. doi: 10.1159/000480419
Xue, M., Pang, H., Li, X., Li, H., Pan, J., and Chen, W. (2016). Long non-coding RNA urothelial cancer-associated 1 promotes bladder cancer cell migration and invasion by way of the hsa-miR-145-ZEB1/2-FSCN1 pathway. Cancer Sci. 107, 18–27. doi: 10.1111/cas.12844
Yan, T. H., Lu, S. W., Huang, Y. Q., Que, G. B., Chen, J. H., Chen, Y. P., et al. (2014). Upregulation of the long noncoding RNA HOTAIR predicts recurrence in stage Ta/T1 bladder cancer. Tumour Biol. 35, 10249–10257. doi: 10.1007/s13277-014-2344-8
Yang, L., Xue, Y., Liu, J., Zhuang, J., Shen, L., Shen, B., et al. (2017). Long noncoding RNA ASAP1-IT1 promotes cancer stemness and predicts a poor prognosis in patients with bladder cancer. Neoplasma 64, 847–855. doi: 10.4149/neo_2017_606
Ye, T., Ding, W., Wang, N., Huang, H., Pan, Y., and Wei, A. (2017). Long noncoding RNA linc00346 promotes the malignant phenotypes of bladder cancer. Biochem. Biophys. Res. Commun. 491, 79–84. doi: 10.1016/j.bbrc.2017.07.045
Yuan, J. H., Yang, F., Wang, F., Ma, J. Z., Guo, Y. J., Tao, Q. F., et al. (2014). A long noncoding RNA activated by TGF-beta promotes the invasion-metastasis cascade in hepatocellular carcinoma. Cancer Cell 25, 666–681. doi: 10.1016/j.ccr.2014.03.010
Zhan, Y., Lin, J., Liu, Y., Chen, M., Chen, X., Zhuang, C., et al. (2016a). Up-regulation of long non-coding RNA PANDAR is associated with poor prognosis and promotes tumorigenesis in bladder cancer. J. Exp. Clin. Cancer Res. 35:83. doi: 10.1186/s13046-016-0354-7
Zhan, Y., Liu, Y., Wang, C., Lin, J., Chen, M., Chen, X., et al. (2016b). Increased expression of SUMO1P3 predicts poor prognosis and promotes tumor growth and metastasis in bladder cancer. Oncotarget 7, 16038–16048. doi: 10.18632/oncotarget.6946
Zhan, Y., Li, Y., Guan, B., Chen, X., Chen, Z., He, A., et al. (2017a). Increased expression of long non-coding RNA CCEPR is associated with poor prognosis and promotes tumorigenesis in urothelial bladder carcinoma. Oncotarget 8, 44326–44334. doi: 10.18632/oncotarget.17872
Zhan, Y. H., Li, Y. F., Guan, B., Wang, Z. C., Peng, D., Chen, Z. C., et al. (2017b). Long non-coding RNA HNF1A-AS1 promotes proliferation and suppresses apoptosis of bladder cancer cells through up-regulating Bcl-2. Oncotarget 8, 76656–76665. doi: 10.18632/oncotarget.20795
Zhang, H., Guo, Y., Song, Y., and Shang, C. (2017). Long noncoding RNA GAS5 inhibits malignant proliferation and chemotherapy resistance to doxorubicin in bladder transitional cell carcinoma. Cancer Chemother. Pharmacol. 79, 49–55. doi: 10.1007/s00280-016-3194-4
Zhang, S., Zhong, G., He, W., Yu, H., Huang, J., and Lin, T. (2016). lncRNA up-regulated in nonmuscle invasive bladder cancer facilitates tumor growth and acts as a negative prognostic factor of recurrence. J. Urol. 196, 1270–1278. doi: 10.1016/j.juro.2016.05.107
Zhao, F., Lin, T., He, W., Han, J., Zhu, D., Hu, K., et al. (2015). Knockdown of a novel lincRNA AATBC suppresses proliferation and induces apoptosis in bladder cancer. Oncotarget 6, 1064–1078. doi: 10.18632/oncotarget.2833
Zhao, X. L., Zhao, Z. H., Xu, W. C., Hou, J. Q., and Du, X. Y. (2015). Increased expression of SPRY4-IT1 predicts poor prognosis and promotes tumor growth and metastasis in bladder cancer. Int. J. Clin. Exp. Pathol. 8, 1954–1960.
Zhen, S., Hua, L., Liu, Y. H., Sun, X. M., Jiang, M. M., Chen, W., et al. (2017). Inhibition of long non-coding RNA UCA1 by CRISPR/Cas9 attenuated malignant phenotypes of bladder cancer. Oncotarget 8, 9634–9646. doi: 10.18632/oncotarget.14176
Keywords: bladder cancer, long non-coding RNAs, clinicopathological parameters, prognosis, diagnosis, systematic review
Citation: Su G, He Q and Wang J (2018) Clinical Values of Long Non-coding RNAs in Bladder Cancer: A Systematic Review. Front. Physiol. 9:652. doi: 10.3389/fphys.2018.00652
Received: 21 February 2018; Accepted: 14 May 2018;
Published: 30 May 2018.
Edited by:
Marina Pajic, Garvan Institute of Medical Research, AustraliaReviewed by:
Venessa T. Chin, Garvan Institute of Medical Research, AustraliaAdnan Nagrial, University of Sydney, Australia
Copyright © 2018 Su, He and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guoming Su, suguoming@yahoo.com
†These authors have contributed equally to this work.