Abstract
To investigate the diagnostic value of urothelial carcinoma–associated 1 as a urine biomarker in urinary bladder cancer patients by performing a comprehensive meta-analysis. A comprehensive literature search was conducted by the databases PubMed, Embase, Cochrane Library, China Knowledge Resource Integrated, and Web of Science. The quality of eligible studies was scored with the Quality Assessment of Diagnostic Accuracy Studies. The bivariate meta-analysis model was used to pool the sensitivity, specificity, likelihood ratio, and diagnostic odds ratio. Receiver operating characteristic curves and hierarchical summary receiver operating characteristic models were employed to check the overall test performance in this meta-analysis. Seven publications involving 678 patients and 563 controls were included in this meta-analysis. The pooled sensitivity was 0.84 (95% confidence interval: 0.80–0.88), specificity was 0.87 (95% confidence interval: 0.75–0.94), positive likelihood ratio was 6.5 (95% confidence interval: 3.10–13.62), negative likelihood ratio was 0.18 (95% confidence interval: 0.13–0.25), and diagnostic odds ratio was 36 (95% confidence interval: 13–99). The area under the summary receiver operating characteristic curve was 0.89 (95% confidence interval: 0.86–0.91). Our results indicated that urothelial carcinoma–associated 1 was a potential diagnostic biomarker with good specificity and sensitivity in urinary bladder cancer. Further prospective studies with larger cohorts are necessary to evaluate the diagnostic accuracy of urothelial carcinoma–associated 1 for urinary bladder cancer.
Introduction
Urinary bladder cancer (UBC) is one of the most common malignancies with the 4th and 10th leading cause of cancer-related deaths in males and females, respectively. 1 In 2015, 74,000 new UBC cases in the United States with a male/female sex ratio of 2.5:1 were reported. 2 The histologic subtype of UBC is quite diverse, such as urothelial carcinoma, squamous cell carcinoma, adenocarcinoma, and small cell carcinoma. The most common subtype of UBC is transitional cell cancer or more properly urothelial cell carcinoma. 3 However, the outcome of UBC patients is not favorable. The majority of patients with UBC could be diagnosed at early stages, but up to 50% cases recur after treatment, and 15%–40% cases progress to muscle-invasive diseases. 4 Screening cystoscopy, voided urinary cytology (VUC), and random bladder biopsies remain the most commonly used methods for initial diagnosis of UBC. 5 As their complex preparatory procedures, high cost or unreliable result of single approach; for example, the sensitivity of cytology was 0.51. 6 Thus, an urgent requirement is needed to identify an alternative urine biomarker to the current diagnostic system.
Long non-coding RNAs (lncRNAs) are longer than 200 nucleotides encoding no protein, which are involved in a variety of biological functions and are dysregulated in various cancers. 7 The aberrant expression of lncRNAs is associated with a broad range of cellular processes, such as cell growth, survival, migration, invasion, and differentiation.8,9 In recent years, increasing studies have highlighted the role of lncRNAs in tumorigenesis and suggested that this group of genes might be used as biomarkers in cancer. 10
The human urothelial carcinoma–associated 1 (UCA1), a 2314-bp lncRNA located on human chromosome 19, 11 was upregulated in many cancers, such as hepatocellular cancer, colorectal cancer, gastric cancer, esophageal squamous cell carcinoma, and lung cancer.12–16 Especially, it was first identified in bladder cancer and thought to be involved in embryogenesis and UBC progression. 17 So, we analyzed transcriptome profiling in human bladder cancer tissues and normal tissue from The Cancer Genome Atlas (TCGA), aiming to validate UCA1 differentially expressed in UBC.
We wonder whether UCA1 differentially expressed in urine sediments so that it could be used as a novel urine biomarker for the detection of UBC. As we found that recent studies have assessed the diagnostic value of UCA1 in UBC, the inconsistencies or discrepancies about diagnostic accuracy exist.6,17–22 For example, Wang et al. showed that the specificity was 91.8% in the diagnosis of bladder cancer, but Milowich et al. showed the specificity was 64.6%. Therefore, we performed this meta-analysis to evaluate the diagnostic efficiency of UCA1 in patients with UBC.
Materials and methods
Identification of differentially expressed genes associated with bladder cancer
For identification of differentially expressed genes, we used edgeR algorithm 23 based on transcriptome profiling downloaded from TCGA. The Benjamini–Hochberg (BH) multiple testing method was used to adjust the p value. The thresholds were the false discovery rate (FDR) < 0.05 and |log2FC| > 2.
Publication search
PubMed, Embase, Cochrane Library, China Knowledge Resource Integrated (CNKI), and Web of Science databases were searched (up to 14 August 2016) for relevant articles that estimated the diagnostic value of UCA1 in UBC. The search terms for the literature retrieval were “UCA1” or “urothelial cancer associated 1” or “CUDR” or “UCAT1” or “LINC00178” or “onco-lncRNA-36” or “NCRNA00178” and “Bladder Neoplasm” or “Bladder Tumor” or “Bladder Cancer” or “Bladder carcinoma” or “Malignant Tumor of Urinary Bladder.” Reference lists of relevant articles were also reviewed to identify potential eligible studies.
Inclusion and exclusion criteria
All the eligible studies included in this meta-analysis had to meet the following criteria: (1) UBC was diagnosed by histopathological confirmation, (2) the expression of UCA1 in urine sediments was tested, (3) healthy individuals or patients with benign disease were included in the control group, and (4) sufficient data in the studies were enough to construct a two-by-two table. Publications were excluded based on any of the following criteria: (1) unqualified data, (2) duplicate publications, and (3) animal research, letters, editorials, expert opinions, case reports, or laboratory articles. The retrieved literature were independently assessed by two reviewers according to the preestablished criteria, and discrepancies were solved through discussions and consensus.
Data extraction and quality assessment
Data were retrieved from the eligible studies independently by three investigators, according to the inclusion and exclusion criteria. Disagreements were discussed with the corresponding author in conference. The characteristics of data extracted for this systematic review included first author, publication year, country of publication, cancer type, sample type, test method, sample size, reference test, cutoff value, and data for two-by-two tables (true positive (TP), false positive (FP), true negative (TN), and false negative (FN)).
The Quality Assessment of Diagnostic Accuracy Studies (QUADAS) list was used to systematically assess the quality of the articles included in this diagnostic meta-analysis. 24 Specifically, 14 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 with a total score of 14.
Statistical analysis
Data analysis was performed using Stata 13.0 software. 25 For this diagnostic meta-analysis, we retrieved TP, FP, TN, and FN as bivariate data directly or through recalculation on the basis of relative data from each eligible study. 26 Subsequently, we calculated the sensitivity, specificity, diagnostic odds ratio (DOR), positive likelihood ratio (PLR), and negative likelihood ratio (NLR) and constructed their forest plots. Then, we generated summary receiver operator characteristic (SROC) curves analysis and calculated the area under the ROC curves (AUC) to evaluate the associations of lncRNA UCA1 and UBC. 27 In case of difference from different diagnostic threshold effect in the traditional SROC, we plotted study-specific sensitivity and specificity by applying the hierarchical summary ROC model (HSROC). 28
To detect cutoff threshold effects, the Spearman correlation coefficient was employed to assess the relationship between sensitivity and specificity. We used Cochran’s Q and I2 statistics to evaluate the heterogeneity across this study, I2 more than 50% indicated the existence of significant heterogeneity. 29 Random-effects model was employed when heterogeneity existed; otherwise, the fixed-effects model was applied.
Results
UCA1 expression levels in human bladder cancer tissue
An analysis was performed in the whole-genome gene expression profiles including 410 tumor tissues and 19 normal tissues from TCGA using edgeR algorithm. When comparing cancer tissues with normal tissues in the data set, we found 1156 differentially expressed genes being assigned for bladder cancer (Figure 1(a)) including UCA1. The cutoff value of the adjusted p values was set to 0.05 and the minimum fold change (FC) to 2. We identified that the expression level of UCA1 in bladder cancer tissues is higher than in normal bladder tissues (Figure 1(b)).

(a) Hierarchical clustering of differentially expressed genes: rows represent genes and columns represent samples. For a gene, red represents higher expression level, blue represents lower expression level, and white represents medial expression level for all samples. For a sample, red means tumor tissues and blue means tumor adjacent tissues from patients with UBC. (b) Relative expression of UCA1 in the tumor tissue compared to normal tissues.
Literature search
A total of 5758 studies were initially retrieved. After excluding duplicates, 2011 articles were excluded. After screening the title and abstract carefully, 3703 articles were excluded for unrelated to UCA1, bladder cancer, reviews, letters, or not original articles. After further inspection of the full articles, 37 articles were excluded for lacking available diagnostic outcome or sufficient clinical data in their records. Finally, 7 articles were considered eligible for this meta-analysis upon the given criteria. The whole procedure of this study selection is shown in Figure 2.

Flow diagram of study selection process based on inclusion and exclusion criteria.
Study characteristics and quality assessment
Seven publications involving 1241 cases and 563 controls were analyzed.6,17–22 In these selected studies, all of the individuals had been histopathologically confirmed. And, the 563 controls included healthy volunteers and patients with benign disease that had never been diagnosed with malignancy. The studies enrolled in this systemic review were conducted in China, Belgium, Egypt, and India. These included publications published between 2006 and 2015 and investigated the diagnostic value of UCA1 as a urine biomarker for patients with UBC. Furthermore, all the studies used real-time quantitative polymerase chain reaction (RT-qPCR) to assess the expression of UCA1 in urine samples. The main characteristics of the eligible studies are presented in Table 1. We used the QUADAS assessment tool to assess the quality of the seven articles with 14 items. The analysis scores are shown in Table 2 with the total point 14.
Main characteristics of the eligible studies.
NA: not available; TCC: transitional cell carcinoma; RT-qPCR: real-time quantitative polymerase chain reaction; ROC: receiver operating characteristic curve; TP: true positive; FP: false positive; TN: true negative; and FN: false negative.
QUADAS assessment for the eligible studies.
QUADAS: Quality Assessment of Diagnostic Accuracy Studies.
N for “no,” Y for “yes,” U for “unclear.”
Data analysis
The forest plot showed the sensitivity ranged from 0.79 to 0.91 (pooled: 0.84, 95% confidence interval (CI): 0.80–0.88), while specificity ranged from 0.60 to 0.97 (pooled: 0.87, 95% CI: 0.75–0.94). The heterogeneity existed in pooled sensitivity and specificity in the seven studies (I2 = 39.45% and I2 = 91.55%, respectively), suggesting significant heterogeneity in sensitivity and specificity (Figure 3).

Forest plots of sensitivity and specificity for UCA1 test accuracy in the diagnosis of bladder cancer.
Since compared with the specificity and sensitivity, likelihood ratios (LRs) are considered to be more comprehensive and clinically meaningful. We calculated the PLR and NLR for the diagnostic performance of UCA1 (Figure 4(a)). The pooled PLR was 6.5 (95% CI: 3.10–13.62), indicating that the UBC patients have more than a six-fold probability to express UCA1 in comparison with healthy individuals. The pooled NLR was 0.18 (95% CI: 0.13–0.25), which means that the probability of individuals with UBC was 18% when the UCA1 test is negative. PLR > 10 or NLR < 0.1 suggests high diagnostic accuracy (Figure 4(b)). It was also noted that DOR was 36 (95% CI: 13–99), indicating that UCA1 can be used as a good biomarker for UBC diagnosis.

(a) Forest plots of positive and negative likelihood ratios for UCA1 test accuracy in the diagnosis of bladder cancer and (b) pooled LNR and LPR of the diagnostic studies of UCA1 test.
The estimates of sensitivity and specificity of all the included studies were also shown in a SROC curve, together with the summary ROC point (pooled sensitivity against pooled specificity). The value of AUC was near to 1, which was considered statistically significant. The area under the summary ROC curve (AUC) was 0.89 (95% CI: 0.86–0.91). The 95% prediction contour and 95% confidence contour were also plotted (Figure 5(a)).

(a) Summary receiver operating characteristic curves for UCA1 in the diagnosis of bladder cancer and (b) hierarchical summary receiver operating characteristics (HSROC) curve for UCA1 in the diagnosis of bladder cancer.
The HSROC curve of these selected studies was consistent with the results from the bivariate model. The summary operating point estimate of sensitivity and specificity was 0.84 (95% CI: 0.80–0.88) and 0.87 (95% CI: 0.75–0.94), respectively. The value of β was 2.48, and the p value was 0.013. This result showed that the HSROC was asymmetrical. The value of γ, which helps distinguish UBC patients from healthy individuals, was 4.32 (95% CI: 2.59–6.05). The 95% prediction region and 95% confidence region were also plotted. This result indicates that UCA1 is a relatively accurate diagnostic urine marker for UBC (Figure 5(b)).
To validate the stability of outcomes in this meta-analysis, sensitivity analysis showed no significant change after omitting any of the included studies (Table 3).
The influence of each trial for the outcome of the meta-analysis.
DOR: diagnostic odds ratio; CI, confidence interval.
Threshold effect and heterogeneity
The threshold effect was the major cause of heterogeneity, which leads to differences in sensitivity and specificity. In this meta-analysis, the Spearman correlation coefficient of logarithm sensitivity and 1-specificity was −0.46 with p = 0.30 (p > 0.05), suggesting that the heterogeneity was not result of the threshold effect. The I2 of the heterogeneity test was 87%, indicating significant heterogeneity. Because there were only seven articles included, meta-regression analysis and subgroup analysis cannot be used to search sources of heterogeneity.
Discussion
Recent studies have revealed several novel diagnostic biomarkers for UBC, including UCA1 as a urine biomarker. UCA1 is an lncRNA, 10 which is localized to nuclear speckles. The UCA1 transcript is a non-coding RNA with more than 1400 nt, which expressed from chromosome 19p13.12. 30 For the reason that UCA1 expression sequences were conserved across several species, it had been indicated to play a regulatory role in the molecular biology of cancer and other human diseases. 10 Now, with the widespread application of new technologies, a growing number of studies have proven that UCA1 is upregulated in several cancers. Abnormal UCA1 expression was strongly correlated to clinicpathologic characteristics, including lymph node metastasis, chemoresistance, and overall survival.31,32 Recently, UCA1 has been found to be required for mitotic progression and controls cell cycle progression in human cells.30,33 In bladder cancer, UCA1 level was upregulated compared with adjacent normal tissues. Yang et al. 34 found that UCA1 regulates cell cycle progression through cAMP response element–binding protein (CREB) via PI3K-AKT–dependent signaling pathways. Zhen et al. 35 demonstrated that UCA1 has an oncogenic role in promoting bladder cancer progression. Wang et al. 10 showed that UCA1 expression increases invasion and drug resistance of bladder cancer cells.
The UCA1 has attracted the researchers’ attention not only because of the application as biomarkers for cancer prognosis and predictors for metastasis 33 but also because of the potential use as a valuable diagnostic urine biomarker. In early studies, the associations between UBC and several biomarkers such as nuclear matrix protein 22 (NMP22), 36 microRNA-96 (miR-96), 37 soluble fragments of cytokeratin 19 (CYFRA 21-1) 38 were identified. However, the diagnostic value of UBC was unsatisfactory (the specificity or specificity was low). The studies included in our meta-analysis used RT-qPCR to test urine UCA1 for detecting UBC, which was a non-invasive way with good sensitivity and specificity. In combination with cytology, it may improve the detection of UBC with a non-invasive diagnostic test and may increase the interval between cystoscope and replace cystoscopies altogether. 21 Wang et al. 17 showed that a UCA1 assay was highly specific (91.8%, 78 of 85) and very sensitive (80.9%, 76 of 94) in the diagnosis of bladder cancer. However, Milowich et al. 19 indicated that the efficiency of the UCA1 test for detecting primary and recurring bladder cancer was unsatisfactory, the sensitivity was 83.9%, and the specificity was 64.6%. Thus, the present meta-analysis was performed to provide an integrated and up-to-date evaluation of the diagnostic and clinical values of UCA1 as a urinary marker for UBC.
In this meta-analysis, the combined sensitivity and specificity were 0.84 (95% CI: 0.80–0.88) and 0.87 (95% CI: 0.75–0.94), respectively, which implied good sensitivity and specificity. The DOR value is defined as the ratio of TP odds to FP odds reflecting the extent of the association between the diagnostic results and the disease. The DOR value of UCA1 was 36, suggesting that UCA1 can serve as a promising diagnostic biomarker for UBC. In addition, LR combines the stability of sensitivity and specificity to provide an omnibus index of test performance and is more clinically valuable than DOR. 39 PLR greater than 10 or NLR less than 0.1 generate large and often conclusive changes from pre-test to post-test probability. Due to the threshold effect among studies, the SROC curve appears to show a better method to assess the summary diagnostic accuracy of the discrimination between cases and controls rather than pooled sensitivity, specificity, PLR, or DOR. In this study, the AUC was 0.89, suggesting a relatively high level of diagnostic accuracy.
Heterogeneity is a potential problem that can influence the incorporation effect and the interpretation of the meta-analysis results. 40 Although we have set strict inclusion and exclusion criteria to gain eligible studies, heterogeneity still exists due to the existence of potential confounding factors. A prime cause of heterogeneity is the threshold effect, which arises because some studies did not supply the cutoff values. So, Spearman rank correlation test was employed to analyze the threshold effect, and the Spearman correlation coefficient was −0.46 with p = 0.30 (p > 0.05). This result demonstrated that threshold effect was not a certain source of heterogeneity. As there were only seven articles included, meta-regression analysis and subgroup analysis cannot be used. Therefore, the other factors such as ethnicity and test method were not evaluated as sources of heterogeneity.
Zhen et al. 35 also conducted an impressive meta-analysis on the diagnostic value of UCA1 for UBC. They reported a pooled sensitivity of 83% and a pooled specificity of 83%. Their results gave support to our conclusion. In comparison with their included data, our data of included original studies were only from screening groups (data of follow-up were not included), which aim to present UCA1 as a diagnostic test in a screening population (for detecting primary BC). Moreover, in this present meta-analysis, HSROC curve, sensitivity analysis, and Spearman rank correlation were performed.
Several limitations in the present meta-analysis should be emphasized. First, only seven studies were included in this meta-analysis; Deeks’ funnel plots cannot be used to assess the presence of publication bias in this meta-analysis, it has weakened the reliability of this meta-analysis’ results to some extent. So, larger-size and better design studies are needed to confirm our results. Second, not all of the studies reported the cutoff values of UCA1. Further comprehensive studies are needed to solve this problem. Third, only publications in English or Chinese were brought into, researches in other languages might be missed. Fourth, the blind detection and blind judgment methods can minimize the tendencies of diagnostic trials; however, all the studies did not report using blind detection.
Conclusion
Our comprehensive analysis validates that UCA1 may act as a clinical potential urine biomarker in UBC diagnosis with a high AUC value. Furthermore, combining UCA1 with other conventional tests may improve the detection of UBC.
Footnotes
Acknowledgements
C.Z. conceived and designed the experiments. Z.W. and X.W. performed the experiments, analyzed the data, and drafted the paper. D.Z., Y.Y., and L.C. contributed to the reagents/materials/analysis tools.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethical standard
As the study was a systemic review of published data and meta-analysis of the pooled data, we did not apply for the approval of Institutional Review Board.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by the National Natural Science Foundation of China (grant numbers 81171996 and 81272289), the Wujieping Medical Foundation (grant number 320.6750.13252), and Heilongjiang Province Health and Family Planning Commission Science Foundation (grant number 2014-287).
