Long noncoding RNAs as potential biomarkers in retinoblastoma: a systematic review and meta-analysis
Cancer Cell International volume 20, Article number: 201 (2020)
Retinoblastoma is the most common malignant rare intraocular tumor of childhood. Long noncoding RNAs (lncRNAs) have been reported participating in its progression, but their significance remains inconclusive. We conducted this systematic review and meta-analysis to explore specific lncRNA biomarker in patients with retinoblastoma.
Materials and methods
Eligible articles were searched from the Pubmed, Web of Science, Embase and the Cochrane library. Hazard ratios (HRs) and odds ratios (ORs) were extracted or calculated to evaluate the relationship between lncRNAs and retinoblastoma. The meta-analysis part was conducted with STATA v.15 software.
A total of 9 articles with 834 retinoblastoma patients are yielded. Heterogeneity among HRs of overall survival (OS) is notably high (I2 = 91.3%, p < 0.001). Subgroup analysis suggests that elevated expression of lncRNA BDNF-AS and MT1JP are favorable factors in OS (pooled HR = 1.89, 95% CI 1.72–2.07, I2 = 0%). Six articles included optic nerve invasion as a clinicopathological outcome and showed a notable correlation (pooled HR = 2.38, 95% CI 1.26–3.50, I2 = 0.0%). We validate our analysis via the public dataset and also sum up the studies of lncRNA BDNF-AS and MT1JP in other cancers.
Differential expression of lncRNAs has been reported in retinoblastoma. Some of them showed potential in retinoblastoma prognosis and progression.
Retinoblastoma is the most common intraocular rare malignant tumor in children . The estimated incidence is approximately 8000 new cases and 3000 deaths in the world annually [2, 3]. Retinoblastoma is not life-threatening as long as being diagnosed early enough and treated appropriately. However, in developing countries, the situation is quite critical, probably due to rather high birth rates and limited medical resources. Approximately 40–70% of patients died in some Asian and African countries . The largest hospital in South Western China used to conduct a retrospective study and found the majority of patients presented at advanced stage . Late diagnosis leads to delayed treatment, which usually means increasing metastasis and decreasing survival rate . Once proliferating into brain through optic nerve, the damage becomes irreparable and fatal. Therefore, effective and early diagnosis is particularly crucial. Exploring novel diagnostic and therapeutic signatures that can be used in retinoblastoma management deserves more study to achieve the aim of eradicating retinoblastoma in the future decade.
Long noncoding RNAs (lncRNAs) are a group of non-coding RNAs that are greater than 200 nucleotides in length. They are widely reported to contribute to tumor progression including proliferation, migration and invasion . Some of them have already been implicated as tumor biomarkers. For example, lncRNA MALAT1 can be applied to predict early and metastatic lung cancers . It is also reported as a cancer suppressor. Professor Ma from famous MD Anderson Cancer Center found it suppressing the metastasis of breast cancer cells . Similarly, emerging studies demonstrate the differential expression of lncRNAs between retinoblastoma and normal tissues, leading us to consider whether they can be potential approaches for early diagnosis and clinical outcomes.
Therefore, we conducted this meta-analysis to evaluate the function of candidate lncRNAs in retinoblastoma in order to find potential prognostic lncRNA markers for further investigation.
Materials and methods
Two investigators (WJL and QDS) performed this systematic literature search independently through Pubmed, Web of Science, Embase and the Cochrane library (from September 10, 2019 to October 10, 2019). The search terms were used with various conjunctions: (“lncRNA” OR “lincRNA” OR “long noncoding RNA” OR “long untranslated RNA” OR “long intergenic noncoding RNA”) AND (“retinoblastoma”). Related citations of every included article were also manually examined. We only included articles in English language and agreement were achieved on each articles by both of the investigators after careful screening and discussion.
Studies from the literature search were included based on the following conditions: (1) LncRNA expression was detected in retinoblastoma tissues; (2) the retinoblastoma tissues were all pathologically diagnosed and classified by ≥ 2 pathologists; (3) the relationship between lncRNA expression and retinoblastoma patient prognosis was investigated, including overall survival and other related clinicopathological parameters; and (4) sufficient data were provided for the hazard ratio (HR) and the 95% confidence interval (95% CI) of survival and clinicopathological parameters. Meanwhile, unsuited studies were excluded if they are: (1) animal studies and reviews; (2) studies based on patient blood and serum; (3) duplicate records; and (4) studies without other sufficient data.
Data extraction and quality assessment
Two investigators extracted requisite information from each eligible study: first authors, publication year, lncRNA type and expression, sample size, sample type, the HR and 95% CI of lncRNA for survival. Engauge Digitizer v.10.11 software and Tierney’s spreadsheet were utilized to obtain the HR when it was not offered directly. HR from multivariate analysis was preferred over that from the univariate analysis as it accounted for confounding factors. Quality assessment of each eligible study was performed based on the acknowledged Newcastle–Ottawa Quality Assessment Scale (NOS) with scores ranging from zero to nine. A study with a NOS score higher than 6 is regarded as high quality.
Datasets of GSE125903 were normalized before calculation. Differential expressing genes (DEGs) were analyzed by R packages of “limma” (http://www.bioconductor.org/packages/release/bioc/html/limma.html). Limma is an open source package on Bioconductor platform. We use Benjamini–Hochberg method to get the adjusted p values, lmfit and eBayes functions to get the different expressed genes. The cutoff of adjusted p values was set to 0.05 and the cutoff of |log2FC| was 1.5. The genes which meet the two conditions above were identified as DEGs. We use ggpubr and ggthemes to draw volcano plots for each datasets, with an x axis of “log2(Fold Change)” and an y axis of “−log10(Adjust p-value)”. The heatmap was done by R package of “pheatmap”.
As different lncRNAs act differently in retinoblastoma, pooling HRs directly would lead to great heterogeneity, we assessed the heterogeneity among studies and conducted the subgroup analysis of HRs based on different lncRNAs. Heterogeneity of the results was estimated by the Q test and I2 statistics. We selected the random pooling model, to be on the safe side. 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 STATA v.15 software.
Characteristics of the included studies
As shown in Fig. 1, a total of 209 articles were identified from online resources. Among them, 9 articles with 834 retinoblastoma patients published between 2015 and 2019 in China were included in our analysis. The relevant lncRNAs include PVT1 , AFAP1-AS1 , HOTAIR , BANCR , H19 , DANCR , LINC00202 , MT1JP  and BDNF-AS . 7 of them used OS as a main outcome and 2 of them used both overall survival (OS) and disease-free survival (DFS) for survival analysis. OS means the percentage of patients still alive at the endpoint while DFS means the percentage of patients without any symptoms of the retinoblastoma. The expression was detected by quantitative polymerase chain reaction (qPCR) in human retinoblastoma tissue specimens. More than half of the studies included applied median value as the cut-off value to stratify patients into high of low expression group. 6 of the articles reported the HRs directly while the remaining used survival curves. Most included articles are of high quality with a rather high NOS score. More details are shown in Table 1.
The relationship between lncRNAs and patient survival
LncRNAs PVT1, AFAP10AS1, HOTAIR, BANCR, H19, DANCR and LINC00202 were up-regulated in retinoblastoma tissues while MT1JP and BDNF-AS were down-regulated. Owing to the rather significant heterogeneity (I2 = 91.3%, p<0.001), a random model was employed to calculate the 95% CI and pooled HR. To decrease the heterogeneity, we performed subgroup analysis. It suggested no significance between overall high expression and poor OS (pooled HR = 1.31, 95% CI 0.43–2.20, I2 = 51.8%; Fig. 2). Individually, high expression of lncRNA AFAP1-AS1 (HR = 3.6, 95% CI 1.33–9.72), BANCR (HR = 2.9, 95% CI 1.05–8.04) and H19 (HR = 2.91, 95% CI 1.02–8.45) are obviously associated with worse OS. In contrast, increased levels of LncRNA BDNF-AS and MT1JP were favorable factors in OS (pooled HR = 1.89, 95% CI 1.72–2.07, I2 = 0%; Fig. 2). Both LncRNA BDNF-AS and MT1JP have been reported decreased in quite a few cancers including gastric cancer, bladder cancer and lung cancer, suggesting their potential role as a rather broad and sensitive tumor suppresser. Two of the articles included provided DFS data and we also found no significance between their overall high expression and poor DFS (pooled HR = 1.98, 95% CI 0.87–3.09, I2 = 0.0%; Fig. 3).
The relationship between lncRNAs and clinicopathological outcomes
The clinicopathological outcomes mentioned by these 9 articles include tumor size, choroidal invasion, gender, laterality, optic nerve invasion and pathologic grade. All lncRNAs had no remarkable relationship with tumor size, choroidal invasion, gender, laterality and pathologic grade (tumor size: p = 0.209, choroidal invasion: p = 0.996, gender: p = 0.844, laterality: p = 0.823 and pathologic grade: p = 0.005). LncRNA HOTAIR was significantly correlated with larger tumor size (HR = 5, 95% CI 2.86–9.09). Besides, it also related with TNM classification (HR = 1.89%, 95% CI 1.06–3.33). 6 articles included applied optic nerve invasion as a clinicopathological outcome and showed a notable correlation (pooled HR = 2.38, 95% CI 1.26–3.50, I2 = 0.0%; Fig. 4). Optic nerve invasion is known as a major risk factor for retinoblastoma mortality . There is a relation between the extent of invasion and metastasis. Once the invasion come to the cut end of the optic nerve, the mortality increases to 50 ~ 80% . However, it remains unsolved how to predict the invasion before it happens or evaluate the involvement during development . Probably these three lncRNAs can provide some new hints. More details are shown in Fig. 5.
Publication bias and sensitivity analysis
Egger’s publication bias plot and Begg’s funnel plot were constructed to explore the potential publication bias in this study. 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. 5). In addition, the sensitivity analysis suggested the robustness of the results because eliminating any study did not change the results significantly (Fig. 6).
Action mechanisms of lncRNAs in retinoblastom
Furthermore, we concentrated on potential targets and pathways of these lncRNAs included in our study (Table 2). Most of them sponge with microRNA to influence retinoblastoma progression.
Validation of the candidate lncRNAs
Right now, we have figured out two potential lncRNAs. For the further validation, we searched the public datasets and found GSE125903 (Fig. 7). It was the first transcriptomic report on pediatric retinoblastoma tumor based on 7 retinoblastoma tissue and 3 control retina tissue . We downloaded its original data and analyzed the 9 lncRNAs discussed in our analysis. The results were showed in a heatmap (Fig. 8). LncRNA MT1JP and BDNF-AS were slightly decreased in retinoblastoma tissues, which was consistent with our analysis. In addition, lncRNA DANCR, PVT1 and H19 also showed significant changes, which deserved more attention.
Retinoblastoma accounts for 3% of all childhood cancers and is one of the most malignant intraocular tumors . Despite of the development of overall health care, the mortality of retinoblastoma still remains relatively high in Asia, especially in China, partly due to the likely delayed diagnosis . Thus, looking for suitable and sensitive screening biomarkers would benefit clinicians a lot in taking actions in time.
In the past few years, various lncRNAs have been reported in the pathogenesis and progression of retinoblastoma. Last year, Pluousiou et al. did a comprehensive review of lncRNAs involved in retinoblastoma , however, right now, no study has described the relation between lncRNAs and patient survival or clinicopathological outcomes systematically.
In our study, we analyzed the role of 9 lncRNAs in retinoblastoma survival and clinical indicators. Seven upregulated lncRNAs showed no overall remarkable relationship with OS. Two downregulated lncRNAs: BDNF-AS and MT1JP predicted favorable OS. We further searched the public dataset to validate our analysis. However, we also found both of them have been reported in quite a few cancers (Table 3), thus their specificity for retinoblastoma needs further verification. High expression of Lnc00202 and BANCR both indicated poorer DFS while taken together, this result is not significant.
Among the 5 main clinicopathological parameters, lncRNAs had significant relation with optic nerve invasion. Optic nerve invasion is one of the most significant risk factor for eye enucleation and mortality and with the increase of involvement, the mortality risk increases, too. Predicting or managing the development of nerve invasion remains a demanding task  while lncRNAs might be somewhat of benefit.
We also summed up the mechanisms of lncRNAs in retinoblastoma, which might be helpful for further comprehensive recognition.
Limitations of our meta-analysis include: (1) All of the studies included are conducted in China. We propose studies outside China to compare and confirm our results. (2) Sample size of most studies included are rather limited and larger researches are still needed. (3) Some of the studies did not provide HRs directly and we collected the necessary information from survival curves, which might be not so accurate. (4) We only collected the targets and pathways of lncRNAs studied in the articles we included. More regulatory mechanisms are available in articles we did not enrolled.
Our study systematically reviewed and analyzed the abnormally expressed lncRNAs in retinoblastoma. Based on all eligible evidence, we found two lncRNAs:BDNF-AS and MT1JP as potential prognostic biomarker for retinoblastoma. Considering our existed limitations and further needs, larger-scale prospective investigations are warranted for validation.
Availability of data and materials
All data are included in this article.
Yahaya JJ, Rugwizangoga B, Mremi A, Munema A. Clinicopathological findings of retinoblastoma: a 10-year experience from a tertiary hospital in Kampala, Uganda. J Ophthalmol. 2019;2019:5829284. https://doi.org/10.1155/2019/5829284(Published 2019 Jun 16).
Dimaras H, Corson TW, Cobrinik D, et al. Retinoblastoma. Nat Rev Dis Primers. 2015;1:15021. https://doi.org/10.1038/nrdp.2015.21(Published 2015 Aug 27).
Kivela T. The epidemiological challenge of the most frequent eye cancer: retinoblastoma, an issue of birth and death. Br J Ophthalmol. 2009;93(9):1129–31. https://doi.org/10.1136/bjo.2008.150292.
Gudiseva HV, Berry JL, Polski A, Tummina SJ, O’Brien JM. Next-generation technologies and strategies for the management of retinoblastoma. Genes. 2019;10(12):1032. https://doi.org/10.3390/genes10121032(Published 2019 Dec 11).
Gao J, Zeng J, Guo B, et al. Clinical presentation and treatment outcome of retinoblastoma in children of South Western China. Medicine. 2016;95(42):e5204. https://doi.org/10.1097/md.0000000000005204.
Dimaras H, Kimani K, Dimba EA, Gronsdahl P, White A, Chan HS. Retinoblastoma. Lancet. 2012;379(9824):1436–46. https://doi.org/10.1016/S0140-6736(11)61137-9(ISSN 0140-6736).
Li J, Li Z, Zheng W, et al. LncRNA-ATB: an indispensable cancer-related long noncoding RNA. Cell Prolif. 2017;50:e12381. https://doi.org/10.1111/cpr.12381.
Bhan A, Soleimani M, Mandal SS. Long noncoding RNA and cancer: a new paradigm. Cancer Res. 2017;77(15):3965–81. https://doi.org/10.1158/0008-5472.can-16-2634.
Kim J, Piao HL, Kim BJ, et al. Long noncoding RNA MALAT1 suppresses breast cancer metastasis. Nat Genet. 2018;50(12):1705–15. https://doi.org/10.1038/s41588-018-0252-3.
Wu XZ, Cui HP, Lv HJ, Feng L. Knockdown of lncRNA PVT1 inhibits retinoblastoma progression by sponging miR-488-3p. Biomed Pharmacother. 2019;112:108627. https://doi.org/10.1016/j.biopha.2019.108627(ISSN 0753-3322).
Hao F, Mou Y, Zhang L, Wang S, Yang Y. LncRNA AFAP1-AS1 is a prognostic biomarker and serves as oncogenic role in retinoblastoma. Biosci Rep. 2018;38(3):BSR20180384. https://doi.org/10.1042/bsr20180384.
Yang G, Fu Y, Lu X, Wang M, Dong H, Li Q. LncRNA HOTAIR/miR-613/c-met axis modulated epithelial-mesenchymal transition of retinoblastoma cells. J Cell Mol Med. 2018;22(10):5083–96. https://doi.org/10.1111/jcmm.13796.
Su S, Gao J, Wang T, Wang J, Li H, Wang Z. Long non-coding RNA BANCR regulates growth and metastasis and is associated with poor prognosis in retinoblastoma. Tumor Biol. 2015;36(9):7205–11. https://doi.org/10.1007/s13277-015-3413-3.
Li L, Chen W, Wang Y, Tang L, Han M. Long non-coding RNA H19 regulates viability and metastasis, and is upregulated in retinoblastoma. Oncol Lett. 2018;15(6):8424–32. https://doi.org/10.3892/ol.2018.8385.
Wang J-X, Yang Y, Li K. Long noncoding RNA DANCR aggravates retinoblastoma through miR-34c and miR-613 by targeting MMP-9. J Cell Physiol. 2018;233:6986–95. https://doi.org/10.1002/jcp.26621.
Yan Guigang, Su Yi, Ma Zhao, Yu Lianzhi, Chen Ning. Long noncoding RNA LINC00202 promotes tumor progression by sponging miR-3619-5p in retinoblastoma. Cell Struct Funct. 2019;44:51–60. https://doi.org/10.1247/csf.18033.
Bi LL, Han F, Zhang XM, Li YY. LncRNA MT1JP acts as a tumor inhibitor via reciprocally regulating Wnt/beta-catenin pathway in retinoblastoma. Eur Rev Med Pharmacol Sci. 2018;22(13):4204–14. https://doi.org/10.26355/eurrev_201807_15414.
Shang W, Yang Y, Zhang J, Wu Q. Long noncoding RNA BDNF-AS is a potential biomarker and regulates cancer development in human retinoblastoma. Biochem Biophys Res Commun. 2018;497(4):1142–8. https://doi.org/10.1016/j.bbrc.2017.01.134(ISSN 0006-291X).
Fabian ID, Onadim Z, Karaa E, et al. The management of retinoblastoma. Oncogene. 2018;37:1551–60. https://doi.org/10.1038/s41388-017-0050-x.
Kim JW. Retinoblastoma. Int Ophthalmol Clin. 2015;55(1):77–96. https://doi.org/10.1097/iio.0000000000000048.
Berry JL, Zolfaghari E, Chen A, Murphree AL, Jubran R, Kim JW. Optic nerve obscuration in retinoblastoma: a risk factor for optic nerve invasion? Ocul Oncol Pathol. 2017;3(4):283–91. https://doi.org/10.1159/000464468.
Rajasekaran S, Nagarajha Selvan LD, Dotts K, et al. Non-coding and coding transcriptional profiles are significantly altered in pediatric retinoblastoma tumors. Front Oncol. 2019;9:221. https://doi.org/10.3389/fonc.2019.00221(Published 2019 Apr 16).
Rao R, Honavar SG. Retinoblastoma. Indian J Pediatr. 2017;84:937–44. https://doi.org/10.1007/s12098-017-2395-0.
Dimaras H, Corson TW. Retinoblastoma, the visible CNS tumor: a review. J Neurosci Res. 2019;97(1):29–44. https://doi.org/10.1002/jnr.24213.
Plousiou M, Vannini I. Non-coding RNAs in retinoblastoma. Front Genet. 2019;10:1155. https://doi.org/10.3389/fgene.2019.01155(Published 2019 Nov 14).
Zhao H, Diao C, Wang X, et al. LncRNA BDNF-AS inhibits proliferation, migration, invasion and EMT in oesophageal cancer cells by targeting miR-214. J Cell Mol Med. 2018;22(8):3729–39. https://doi.org/10.1111/jcmm.13558(published online ahead of print, 2018 Jun 12).
Zhang H, Liu C, Yan T, Wang J, Liang W. Long noncoding RNA BDNF-AS is downregulated in cervical cancer and has anti-cancer functions by negatively associating with BDNF. Arch Biochem Biophys. 2018;646:113.
Zhi H, Lian J. LncRNA BDNF-AS suppresses colorectal cancer cell proliferation and migration by epigenetically repressing GSK-3β expression. Cell Biochem Funct. 2019. https://doi.org/10.1002/cbf.3403.
Li W, Dou Z, We S, Zhu Z, Pan D, Jia Z, Liu H, Wang X, Yu G. Long noncoding RNA BDNF-AS is associated with clinical outcomes and has functional role in human prostate cancer. Biomed Pharmacother. 2018;102:1105–10.
Huang Q, Yang J, He X, Shi S, Xing S. LncRNA BDNF-AS is associated with the malignant status and regulates cell proliferation and apoptosis in osteosarcoma. Biosci Rep. 2018;38(6):BSR20181498. https://doi.org/10.1042/bsr20181498(Published 2018 Nov 15).
Esfandi F, Bouraghi H, Glassy MC, et al. Brain-derived neurotrophic factor downregulation in gastric cancer. J Cell Biochem. 2019;120:17831–7. https://doi.org/10.1002/jcb.29050.
Xu Y, Zhang G, Zou C, Zhang H, Gong Z, Wang W, Zhang W. LncRNA MT1JP suppresses gastric cancer cell proliferation and migration through MT1JP/MiR-214-3p/RUNX3 Axis. Cell Physiol Biochem. 2018;46(6):2445–59. https://doi.org/10.1159/000489651.
Zhang G, Li S, Lu J, et al. LncRNA MT1JP functions as a ceRNA in regulating FBXW7 through competitively binding to miR-92a-3p in gastric cancer. Mol Cancer. 2018;17(1):87. https://doi.org/10.1186/s12943-018-0829-6(Published 2018 May 2).
Zhu D, Zhang X, Lin Y, Liang S, Song Z, Dong C. MT1JP inhibits tumorigenesis and enhances cisplatin sensitivity of breast cancer cells through competitively binding to miR-24-3p. Am J Transl Res. 2019;11(1):245–56 (Published 2019 Jan 15).
Yu H, Wang S, Zhu H, Rao D. LncRNA MT1JP functions as a tumor suppressor via regulating miR-214-3p expression in bladder cancer. J Cell Physiol. 2019. https://doi.org/10.1002/jcp.28274.
Ma J, Yan H, Zhang J, Tan Y, Gu W. Long-chain non-coding RNA (lncRNA) MT1JP suppresses biological activities of lung cancer by regulating miRNA-423-3p/Bim axis. Med Sci Monit. 2019;25:5114–26. https://doi.org/10.12659/msm.914387(Published 2019 Jul 10).
This study was supported by the National Natural Science Foundation of China (Key Program) (81730026, 81425006), Frontier Project of Shanghai Hospital Development Center (SHDC12016105), Science and Technology Commission of Shanghai Municipality(17411953000, 16140900800, 0303N17001).
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.
About this article
Cite this article
Wu, J., Qian, D. & Sun, X. Long noncoding RNAs as potential biomarkers in retinoblastoma: a systematic review and meta-analysis. Cancer Cell Int 20, 201 (2020). https://doi.org/10.1186/s12935-020-01281-0