- Open Access
Long noncoding RNAs as potential biomarkers and therapeutic targets in gallbladder cancer: a systematic review and meta-analysis
Cancer Cell International volume 19, Article number: 169 (2019)
Mounting evidence has shown that long noncoding RNAs (lncRNAs) can play a substantial role in gallbladder cancer (GBC) development as tumor promotors or suppressors, and their abnormal expression is relevant to GBC patient outcomes. We completed this systematic review and meta-analysis to explore the clinical significance and mechanisms of lncRNAs in GBC.
We conducted a comprehensive literature search and selected eligible records according to the inclusion and exclusion criteria. Hazard ratios (HRs) and odds ratios (ORs) were extracted or calculated to estimate the relationships of high lncRNA expression with GBC patient survival and clinical outcomes.
Eighteen studies were identified as eligible for this systematic review and meta-analysis. Heterogeneity among HRs of overall survival (OS) was notably high (I2 = 86.2%, p < 0.001). Subgroup analysis suggested that overexpression of lncRNAs in a group that is upregulated in GBC showed a significant association with poor OS (HR = 2.454, 95% CI 2.004–3.004, I2 = 0%). Conversely, overexpression of lncRNAs in a downregulated group was markedly related to good OS (HR = 0.371, 95% CI 0.267–0.517, I2 = 0%). High expression levels of lncRNA AFAP1-AS1, MALAT1 and ROR were positively correlated with tumor size. Expression of lncRNA LET, LINC00152 and HEGBC exhibited a positive correlation with high T status. LncRNA LINC00152, HEGBC, MALAT1 and ROR showed a marked correlation with positive lymph node metastasis (LNM), while lncRNA GCASPC, MEG3, LET and UCA1 had the opposite effect. High expression levels of lncRNA HEGBC, PAGBC, PVT1 and UCA1 predicted high tumor node metastasis (TNM) stages, while lncRNA LET, GCASPC and MEG3 indicated low TNM stages. We also summarized the mechanisms of lncRNAs in GBC.
Aberrant expression of several lncRNAs was indicative of the prognosis of GBC patients, and lncRNAs showed promise as biomarkers and therapeutic targets for GBC.
Gallbladder cancer (GBC) is the most common malignant extrahepatic biliary duct tumor with high malignancy . Statistics show an estimated 12,360 new gallbladder cancer and other biliary cancer cases and 3960 related deaths in the United States in 2019 . Only when the tumor grows up to a large size do symptoms appear, such as pain, jaundice and fever . Because of the asymptomatic characteristics in the early stage, most patients are diagnosed late and miss the best opportunity for surgical resection, which is the only possible curative treatment for GBC. The 5-year survival rates are 85.9% and 56.1% for T1 (tumor was localized at the mucosal or lamina propria) and T2 (tumor invaded the subserosal layer) tumors, but only 19.2% and 14.1% for T3 and T4 tumors (tumor invasion was over the serosa) . Thus, exploring novel biomarkers and therapeutic targets for GBC is a top priority.
Long noncoding RNAs (lncRNAs) play vital roles in cell processes, such as cell proliferation, differentiation, DNA damage response, chromosomal imprinting, etc. via transcriptional, posttranscriptional or epigenetic regulation of gene expression, rather than functioning as coding RNAs [5, 6]. Abnormal lncRNA expression is tightly correlated with tumorigenesis and progression, and lncRNAs act as not only tumor promotors but also suppressors. For example, Di Huang et al.  found that NF-κB-interacting lncRNA (NKILA) increased T cell sensitivity to activation-induced cell death (AICD) by inhibiting NF-κB, which promoted immune evasion in breast cancer. A p53-regulated lncRNA (loc285194), a tumor suppressor, could inhibit colon cancer cell growth and proliferation by repressing miR-211 .
Wu et al.  demonstrated that in gallbladder cancer, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) promoted GBC cell proliferation and metastasis through the activation of the ERK/MAPK pathway. In the following years, substantial numbers of lncRNAs related to GBC development and patient outcomes were discovered, such as LET , GCASPC , ANRIL , etc. The results of these studies have not been systematically reviewed and analyzed. Thus, we thoroughly collected all relevant studies and summarized the results to explore the potential prognostic value of lncRNAs in GBC.
Materials and methods
A systematic literature search was conducted by two authors (Yuan Zhong and Xiaochao Wu) through PubMed, Web of Science, Embase and the Cochrane Library (up to March 21, 2019). The following search terms were used with Boolean conjunctions: (“lncRNA” OR “lincRNA” OR “long noncoding RNA” OR “long untranslated RNA” OR “long intergenic noncoding RNA”) AND (“gallbladder cancer” OR “gallbladder carcinoma” OR “carcinoma of gallbladder” OR “gallbladder tumor” OR “cancer of gallbladder” OR “gallbladder neoplasm”).
Studies from the literature search were included based on the following criteria: (1) LncRNA expression was detected in gallbladder cancer tissues; (2) the diagnosis of gallbladder cancer was confirmed by pathology and histology; (3) the relationship between lncRNA expression and gallbladder cancer patient prognosis was investigated, including survival and clinicopathological parameters; and (4) sufficient data were provided for the hazard ratio (HR) and the 95% confidence interval (95% CI) of survival or the odds ratio (OR) and the 95% CI of clinicopathological parameters. Meanwhile, unsuited studies were excluded by the following criteria: (1) animal studies, reviews, editorials and letters; (2) duplicate records; and (3) studies without sufficient data.
Data extraction and quality assessment
Two investigators extracted requisite information from each eligible study: first authors, publication year, country, lncRNA type, sample size (high/low), sample type, lncRNA detection method, the HR and 95% CI of lncRNA for survival and patient number for histological grade, TNM stage, LNM, etc. Engauge Digitizer v.10.11 software and Tierney’s spreadsheet were utilized to obtain the HR when it was not reported directly . The HR from multivariate analysis was preferred over that from the univariate analysis because the multivariate analysis considered the influence of confounding factors. The quality assessment of each eligible study was performed in accordance with the Newcastle–Ottawa Quality Assessment Scale (NOS) .
As different lncRNAs could act inversely in GBC, pooling HRs and ORs directly would lead to great heterogeneity. Thus, we assessed the heterogeneity among studies and conducted the subgroup analysis of HRs based on the role of the lncRNA in gallbladder cancer. The heterogeneity of the results was estimated by the Q test and I2 statistics. When I2 ≥ 50%, the random pooling model was selected. In contrast, the fixed pooling model was used. HR > 1 suggested a significant association of lncRNA overexpression with poor survival, and HR < 1 indicated that high lncRNA expression predicted long survival . All these calculations were completed with the assistance of STATA v.15 software.
Characteristics of the included studies
As shown in Fig. 1, a total of 18 articles published between 2015 and 2019 with 1001 gallbladder patients were included in our analysis. The 15 relevant lncRNAs were as follows: LET , GCASPC , LINC00152 [16, 17], MEG3 [12, 18], ANRIL , MALAT1 , MINCR , AFAP1-AS1 , UCA1 , PAGBC , Loc344887 , ROR , HEGBC , H19 [27, 28] and PVT1 . The 12 studies included one study of 702 gallbladder cancer patients that used OS to estimate patient survival and 3 studies of 217 patients that used both OS and disease-free survival (DFS). All lncRNA expression levels were detected by quantitative real-time polymerase chain reaction (qRT-PCR). All studies were conducted in China and all the samples were tissues. Furthermore, 4 studies reported HRs directly. More details are shown in Table 1.
The relationship between lncRNAs and patient survival
LncRNAs LET, GCASPC and MEG3 were downregulated in gallbladder cancer tissues and acted as cancer suppressors, while the remaining 12 lncRNAs (LINC00152, ANRIL, MALAT1, MINCR, AFAP1-AS1, UCA1, PAGBC, Loc344887, ROR, HEGBC, H19, PVT1), were upregulated in cancer tissues and promoted gallbladder cancer progression. Heterogeneity among HRs of OS was markedly high (I2 = 86.2%, p < 0.001; Fig. 2). The subgroup analysis suggested that high expression levels of lncRNAs in the upregulation subgroup were significantly related to poor OS (pooled HR = 2.454, 95% CI 2.004–3.004, I2 = 0%). In contrast, increased levels of LET, MEG3 (Fig. 3) and GCASPC were favorable factors in OS (pooled HR = 0.371, 95% CI 0.267–0.517, I2 = 0%; Fig. 2). We also found that high expression levels of LET and GCASPC were markedly associated with DFS (pooled HR = 0.426, 95% CI 0.302–0.601; Fig. 4).
The relationship between lncRNAs and clinicopathological outcomes
The results indicated that increased expression of AFAP1-AS1, MALAT1 and ROR was significantly correlated with large tumor size, and overexpression of LET, LINC00152 and HEGBC indicated advanced T status. LET and MEG3 showed a positive association with better histological grade. Moreover, LINC00152, HEGBC, MALAT1 and ROR exhibited a notable correlation with positive LNM (LINC00152: OR = 4.500, 95% CI 1.166–17.373; HEGBC: OR = 30.750, 95% CI 10.266–92.103; MALAT1: OR = 11.000, 95% CI 1.998–60.572; ROR: OR = 15.889, 95% CI 2.652–95.208). In contrast, GCASPC, MEG3, LET and UCA1 were favorable factors for LNM (GCASPC: OR = 0.297, 95% CI 0.103–0.860; MEG3: OR = 0.201, 95% CI 0.057–0.713; LET: OR = 0.360, 95% CI 0.162–0.802; UCA1: OR = 0.242, 95% CI 0.062–0.942). Four upregulated lncRNAs (HEGBC, PAGBC, PVT1, UCA1) were unfavorable factors for TNM stage, while three lncRNAs (LET, GCASPC, MEG3) showed an opposite trend. All lncRNAs had no remarkable relationship with age or sex (age: p = 0.146 and sex: p = 0.848). More details are shown in Fig. 5.
Both LINC00152 and H19 were explored by two studies on their association with clinicopathological features. As presented in Fig. 6, no obvious relevance of H19 to clinical parameters was found. LINC00152 overexpression was remarkably related to advanced T status (OR = 5.045, 95% CI 1.769–14.384) and positive LNM (OR = 4.639, 95% CI 1.737–12.390) but not TNM stage (p = 0.862).
Publication bias and sensitivity analysis
Egger’s publication bias plot and Begg’s funnel plot were constructed to explore the potential publication bias. Two plots did not display apparent asymmetry (Egger’s test: p = 0.749 and Begg’s test: p = 0.553), which showed no significant publication bias (Fig. 7). In addition, the sensitivity analysis suggested the robustness of the results because eliminating any study did not change the results significantly (Fig. 8).
Action mechanisms of lncRNAs in gallbladder cancer
In addition, we concentrated on potential targets and pathways of the included lncRNAs in gallbladder cancer, as presented in Table 2. MEG3 and UCA1 could both target Enhancer of Zeste Homolog 2 (EZH2), despite their differential roles in gallbladder cancer. Furthermore, MINCR could stimulate EZH2 expression via targeting miR-26a-5p. MEG3 and ANRIL acted inversely on p53 and CyclinD1, which was consistent with their roles in GBC.
Gallbladder cancer is characterized by a difficult early diagnosis, strong invasion and poor prognosis. Although complete resection can eliminate gallbladder tumors confined to the mucosa, the majority of patients are diagnosed at a late stage, which means that tumors have disseminated to nearby organs or lymph nodes or distant body sites. Patients with advanced GBC lose their opportunity for curative radical cholecystectomy, while radiotherapy or chemotherapy tend to be limited and palliative because of their low sensitivity. Thus, investigations for novel biomarkers and therapeutic targets for GBC make sense practically. Noncoding RNAs including lncRNA, microRNA and circRNA play vital roles in GBC tumorigenesis and progression [30, 31]. Increased studies have reported the relationship between aberrant lncRNA expression and GBC, but the results have not yet been systematically evaluated or analyzed. Therefore, we aimed to summarize and analyze the results to obtain more accurate conclusions.
In this systematic review and meta-analysis, we explored the role of lncRNAs in patient survival and clinicopathological parameters. The heterogeneity among HRs of OS was notably high, possibly because of differential expression of lncRNAs in GBC tissues according to nearly non-existent heterogeneity in both subgroups. Overexpression of GBC promoters (MALAT1, MINCR, ROR, AFAP1-AS1, PAGBC, LINC00152, UCA1, HEGBC and PVT1) predicted poor OS, while high expression of GBC suppressors (LET, GCASPC and MEG3) indicated favorable OS. In addition, the expression levels of LET and GCASPC were positively related to long DFS.
With regard to clinicopathological outcomes, high expression of AFAP1-AS1, MALAT1 and ROR corresponded to a large tumor size. High expression of LET, LINC00152 and HEGBC predicted high T status. LET and MEG3 acted as indicators for low histological grade, while other lncRNAs showed no significance. LINC00152, HEGBC, MALAT1 and ROR functioned as markers for positive LNM, while GCASPC, MEG3, LET and UCA1 indicated negative LNM. In consideration of the tumor-promoting effects of UCA1, more studies with a large sample size need to be conducted to verify the conclusion. HEGBC, PAGBC, PVT1 and UCA1 predicted high TNM stage, and LET, GCASPC and MEG3 predicted low TNM stage. H19 showed no association with any parameter according to two studies with a small sample size. The overexpression of LINC00152 predicted high T status and positive LNM.
We also focused on the mechanisms of lncRNAs in GBC, which might be helpful for further studies. LINC00152 could promote GBC development via both the PI3K/AKT pathway and the LINC00152/miR-138/HIF-1α pathway. MEG3 could suppress GBC proliferation and metastasis through targeting both EZH2 and Large Tumor Suppressor 2 (LATS2). PAGBC plays an oncogenic role in GBC by targeting miR-133b and miR-511 and activating AKT/mTOR pathway.
Limitations need to be mentioned in this systematic review and meta-analysis. First, because all the included studies were from China, the conclusions might not be applicable to other populations. Second, estimates that exaggerate the magnitude of the effect might arise from the small sample size of most studies. Third, 11 HRs of OS were calculated from the K–M curves rather than being reported directly.
Our study first systematically reviewed and estimated the association of abnormal lncRNA expression with GBC patient survival and clinical outcomes. We confirmed several lncRNAs, whether GBC promoters or GBC suppressors, as potential biomarkers and therapeutic targets for GBC. In consideration of the limitations of this study, more large-scale and high-quality studies on various ethnic populations are needed to form better conclusions of the value of lncRNAs in GBC.
Availability of data and materials
All data are included in this article.
long noncoding RNA
- 95% CI:
95% confidence interval
quantitative real-time polymerase chain reaction
metastasis-associated lung adenocarcinoma transcript 1
lymph node metastasis
the Newcastle–Ottawa Quality Assessment Scale
maternally expressed gene 3
antisense non-coding RNA at the INK4 locus
actin filament-associated protein 1 antisense RNA 1
plasmacytoma variant translocation 1
high expressed in gallbladder cancer
low expression in tumor
gallbladder cancer-associated suppressor of pyruvate carboxylase lncRNA
urothelial carcinoma associated 1
MYC-induced long noncoding RNA
regulator of reprogramming
prognosis-associated gallbladder cancer
large tumor suppressor 2
v-AKT murine thymoma viral oncogene homologue 2
Ebata T, Ercolani G, Alvaro D, Ribero D, Di Tommaso L, Valle JW. Current status on cholangiocarcinoma and gallbladder cancer. Liver Cancer. 2016;6(1):59–65.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69(1):7–34.
Baiu I, Visser B. Gallbladder cancer. JAMA. 2018;320(12):1294.
Shindoh J, de Aretxabala X, Aloia TA, Roa JC, Roa I, Zimmitti G, Javle M, Conrad C, Maru DM, Aoki T, et al. Tumor location is a strong predictor of tumor progression and survival in T2 gallbladder cancer: an international multicenter study. Ann Surg. 2015;261(4):733–9.
Khorkova O, Hsiao J, Wahlestedt C. Basic biology and therapeutic implications of lncRNA. Adv Drug Deliv Rev. 2015;87:15–24.
Schmitt AM, Garcia JT, Hung T, Flynn RA, Shen Y, Qu K, Payumo AY, Peres-da-Silva A, Broz DK, Baum R, et al. An inducible long noncoding RNA amplifies DNA damage signaling. Nat Genet. 2016;48(11):1370–6.
Huang D, Chen J, Yang L, Ouyang Q, Li J, Lao L, Zhao J, Liu J, Lu Y, Xing Y, et al. NKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death. Nat Immunol. 2018;19(10):1112–25.
Liu Q, Huang J, Zhou N, Zhang Z, Zhang A, Lu Z, Wu F, Mo YY. LncRNA loc285194 is a p53-regulated tumor suppressor. Nucleic Acids Res. 2013;41(9):4976–87.
Wu XS, Wang XA, Wu WG, Hu YP, Li ML, Ding Q, Weng H, Shu YJ, Liu TY, Jiang L, et al. MALAT1 promotes the proliferation and metastasis of gallbladder cancer cells by activating the ERK/MAPK pathway. Cancer Biol Ther. 2014;15(6):806–14.
Ma MZ, Kong X, Weng MZ, Zhang MD, Qin YY, Gong W, Zhang WJ, Quan ZW. Long non-coding RNA-LET is a positive prognostic factor and exhibits tumor-suppressive activity in gallbladder cancer. Mol Carcinog. 2015;54(11):1397–406.
Ma MZ, Zhang Y, Weng MZ, Wang SH, Hu Y, Hou ZY, Qin YY, Gong W, Zhang YJ, Kong X, et al. Long noncoding RNA GCASPC, a target of miR-17-3p, negatively regulates pyruvate carboxylase-dependent cell proliferation in gallbladder cancer. Cancer Res. 2016;76(18):5361–71.
Liu B, Shen ED, Liao MM, Hu YB, Wu K, Yang P, Zhou L, Chen WD. Expression and mechanisms of long non-coding RNA genes MEG3 and ANRIL in gallbladder cancer. Tumour Biol. 2016;37(7):9875–86.
Tierney JF, Stewart LA, Ghersi D, Burdett S, Sydes MR. Practical methods for incorporating summary time-to-event data into meta-analysis. Trials. 2007;8:16.
Lo CK, Mertz D, Loeb M. Newcastle–Ottawa scale: comparing reviewers’ to authors’ assessments. BMC Med Res Methodol. 2014;14:45.
Zhou Q, Liu J, Quan J, Liu W, Tan H, Li W. lncRNAs as potential molecular biomarkers for the clinicopathology and prognosis of glioma: a systematic review and meta-analysis. Gene. 2018;668:77–86.
Cai Q, Wang ZQ, Wang SH, Li C, Zhu ZG, Quan ZW, Zhang WJ. Upregulation of long non-coding RNA LINC00152 by SP1 contributes to gallbladder cancer cell growth an d tumor metastasis via PI3K/AKT pathway. Am J Transl Res. 2016;8(10):4068–81.
Cai Q, Wang Z, Wang S, Weng M, Zhou D, Li C, Wang J, Chen E, Quan Z. Long non-coding RNA LINC00152 promotes gallbladder cancer metastasis and epithelial–mesenchymal transition by regulating HIF-1α via miR-138. Open Biol. 2017. https://doi.org/10.1098/rsob.160247.
Jin L, Cai Q, Wang S, Wang S, Mondal T, Wang J, Quan Z. Long noncoding RNA MEG3 regulates LATS2 by promoting the ubiquitination of EZH2 and inhibits proliferation and invasion in gallbladder cancer. Cell Death Dis. 2018;9(10):1017.
Wang SH, Zhang WJ, Wu XC, Zhang MD, Weng MZ, Zhou D, Wang JD, Quan ZW. Long non-coding RNA Malat1 promotes gallbladder cancer development by acting as a molecular sponge to regulate miR-206. Oncotarget. 2016;7(25):37857–67.
Wang SH, Yang Y, Wu XC, Zhang MD, Weng MZ, Zhou D, Wang JD, Quan ZW. Long non-coding RNA MINCR promotes gallbladder cancer progression through stimulating EZH2 expression. Cancer Lett. 2016;380(1):122–33.
Ma F, Wang SH, Cai Q, Zhang MD, Yang Y, Ding J. Overexpression of LncRNA AFAP1-AS1 predicts poor prognosis and promotes cells proliferation and invasion in gallbladder cancer. Biomed Pharmacother. 2016;84:1249–55.
Cai Q, Jin L, Wang S, Zhou D, Wang J, Tang Z, Quan Z. Long non-coding RNA UCA1 promotes gallbladder cancer progression by epigenetically repressing p21 and E-cadherin expression. Oncotarget. 2017;8(29):47957–68.
Wu XS, Wang F, Li HF, Hu YP, Jiang L, Zhang F, Li ML, Wang XA, Jin YP, Zhang YJ, et al. LncRNA-PAGBC acts as a microRNA sponge and promotes gallbladder tumorigenesis. EMBO Rep. 2017;18(10):1837–53.
Wu XC, Wang SH, Ou HH, Zhu B, Zhu Y, Zhang Q, Yang Y, Li H. The NmrA-like family domain containing 1 pseudogene Loc344887 is amplified in gallbladder cancer and promotes epithelial–mesenchymal transition. Chem Biol Drug Des. 2017;90(3):456–63.
Wang SH, Zhang MD, Wu XC, Weng MZ, Zhou D, Quan ZW. Overexpression of LncRNA-ROR predicts a poor outcome in gallbladder cancer patients and promotes the tumor cells proliferation, migration, and invasion. Tumour Biol. 2016;37(9):12867–75.
Yang L, Gao Q, Wu X, Feng F, Xu K. Long noncoding RNA HEGBC promotes tumorigenesis and metastasis of gallbladder cancer via forming a positive feedback loop with IL-11/STAT3 signaling pathway. J Exp Clin Cancer Res. 2018;37(1):186.
Wang SH, Wu XC, Zhang MD, Weng MZ, Zhou D, Quan ZW. Upregulation of H19 indicates a poor prognosis in gallbladder carcinoma and promotes epithelial–mesenchymal transition. Am J Cancer Res. 2016;6(1):15–26.
Wang SH, Wu XC, Zhang MD, Weng MZ, Zhou D, Quan ZW. Long noncoding RNA H19 contributes to gallbladder cancer cell proliferation by modulated miR-194-5p targeting AKT2. Tumour Biol. 2016;37(7):9721–30.
Chen J, Yu Y, Li H, Hu Q, Chen X, He Y, Xue C, Ren F, Ren Z, Li J, et al. Long non-coding RNA PVT1 promotes tumor progression by regulating the miR-143/HK2 axis in gallbladder cancer. Mol Cancer. 2019;18(1):33.
Kai D, Yannian L, Yitian C, Dinghao G, Xin Z, Wu J. Circular RNA HIPK3 promotes gallbladder cancer cell growth by sponging microRNA-124. Biochem Biophys Res Commun. 2018;503(2):863–9.
Jin YP, Hu YP, Wu XS, Wu YS, Ye YY, Li HF, Liu YC, Jiang L, Liu FT, Zhang YJ, et al. miR-143-3p targeting of ITGA6 suppresses tumour growth and angiogenesis by downregulating PLGF expression via the PI3K/AKT pathway in gallbladder carcinoma. Cell Death Dis. 2018;9(2):182.
Wang SH, Zhang WJ, Wu XC, Weng MZ, Zhang MD, Cai Q, Zhou D, Wang JD, Quan ZW. The lncRNA MALAT1 functions as a competing endogenous RNA to regulate MCL-1 expression by sponging miR-363-3p in gallbladder cancer. J Cell Mol Med. 2016;20(12):2299–308.
Sun KK, Hu PP, Xu F. Prognostic significance of long non-coding RNA MALAT1 for predicting the recurrence and metastasis of gallbladder cancer and its effect on cell proliferation, migration, invasion, and apoptosis. J Cell Biochem. 2018;119(4):3099–110.
This study was supported by the Project of Standard Diagnosis and Treatment of Key Disease of Jiangsu Province (BE2015722).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.