Research Article | Open Access
Zhaojun Wang, Haifeng Li, Fajun Li, Xin Su, Junhang Zhang, "Bioinformatics-Based Identification of a circRNA-miRNA-mRNA Axis in Esophageal Squamous Cell Carcinomas", Journal of Oncology, vol. 2020, Article ID 8813800, 9 pages, 2020. https://doi.org/10.1155/2020/8813800
Bioinformatics-Based Identification of a circRNA-miRNA-mRNA Axis in Esophageal Squamous Cell Carcinomas
Background. Esophageal squamous cell carcinoma (ESCC) has a poor prognosis due to the lack of early disease symptoms. Using bioinformatics tools, this study aimed to discover differentially expressed nonprotein-coding RNAs and genes with potential prognostic relevance in ESCC. Methods. Two microRNAs (miRNAs) and one circular RNA (circRNA) microarray datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differential expression of miRNAs (DEMs) and circRNAs (DECs) was, respectively, identified in ESCC tissue and compared to adjacent healthy tissue. Further analysis was performed using the miRNA microarray datasets, where miRTarBase was used to predict which messenger RNAs (mRNAs) was present. This was followed by protein-protein interaction (PPI) network, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) analyses. Moreover, cytoHubba and UALCAN were used to predict the important nodes and perform patient survival analysis, respectively. The miRNA-associated circRNAs were predicted using the ENCORI website. The interaction between DECs and the predicted circRNAs was also determined. A circRNA-miRNA-mRNA axis was constructed. Results. Associated with RAP1B and circ_0052867, two miRNAs (miR-133b and miR-139-5p) were identified as being differentially expressed and downregulated across the two datasets. Finally, the circ_0052867/miR-139-5p/RAP1B regulatory axis was established. Conclusion. This study provides support for the possible mechanisms of disease progression in ESCC.
Globally, esophageal carcinoma is the eighth most common type of cancer, and the sixth most dominant cause of cancer death [1, 2]. Esophageal squamous cell carcinoma (ESCC) is the major pathology (60–70%) of esophageal cancer (ESCA) . Despite great advances in surgical and radiation therapy, the five-year survival rate of ESCC patients is still under 25%. This is because obvious clinical symptoms of early-stage ESCA are lacking, and many patients are therefore already in an advanced disease stage at the time of diagnosis . A better prognosis is related to the early diagnosis of disease . While molecular markers have been regarded as prognostic markers of ESCC in recent years [5, 6], these biomarkers are not always applicable for early diagnosis or prognosis prediction. Therefore, there is an urgent need to identify new biomarkers to improve the diagnosis and prognosis of ESCC. Biomarkers are stably expressed in plasma, serum, and other body fluids and include circRNAs, miRNAs, and genes .
The occurrence and progression of ESCC are related to genetic factors such as genomic amplifications, mutations, insertions, and deletions, as well as tumor epigenetics. The latter includes noncoding RNA, DNA methylation, and histone acetylation [7, 8]. MiRNAs are endogenous, single-stranded, noncoding RNA polymers formed from 20–22 nucleotides and encoded by a single nuclear DNA. They are mainly involved in transcriptional and posttranscriptional regulation of gene expression and bind to complementary sequences of target mRNA. Decreased expression of tumor suppressor miRNA-associated genes and overexpression of oncogenic miRNAs may lead to tumor formation . Recently, miRNAs have been accepted as good prognostic or diagnostic biomarkers in cancer research . Several studies have assessed the possible role of miRNA in identifying the potential for hepatocellular carcinoma metastasis and tumor recurrence , and it has also been shown that miRNAs can act on different pathways associated with chemotherapeutic drug sensitivity or resistance in several tumors [12, 13]. Sun et al. found that the level of serum miRNA-1290 could be a useful diagnostic and prognostic marker in ESCC patients . Similarly, miR-455-3p was found to play an antioncogenic role by downregulating FAM83F and could be an independent clinical biomarker for predicting the outcome of ESCC patients . One study indicated that decreased miR-145 expression levels could be a useful prognostic marker and could be used to predict overall survival in ESCC patients . The ciRS-7 molecule was shown to be a novel prognostic biomarker and therapeutic target for ESCC .
This study aimed to identify differentially expressed miRNAs (DEMs) using GEO datasets, predict which mRNAs were associated to the DEMs, and perform functional enrichment and protein-protein interaction (PPI) network analysis to detect potential disease-associated target genes. The key ESCC-associated genes were identified by survival analysis. The intersection between predicted miRNAs, circRNAs, and differentially expressed circRNAs (DECs) was determined against GEO data, and a potential circRNA-miRNA-mRNA network for ESCC was subsequently established.
2. Materials and Methods
2.1. Source of the Microarray Data
The public GEO database (http://www.ncbi.nlm.nih.gov/geo) is a functional genomic database from which three ESCC microarray expression datasets (GSE97051, GSE59973, and GSE131969) were downloaded. The circRNA dataset (GSE131969) included data from three ESCC tissue and three adjacent noncancerous tissue samples. The samples were tested against GPL19978 Agilent-069978 Arraystar Human CircRNA microarray V1 data. The miRNA dataset (GSE97051) included data from seven ESCC tissue and seven nonpathogenic esophageal tissue samples that were tested against GPL21572 Agilent-067406 Human CBC lncRNA + mRNA microarray V4.0 data. The miRNA dataset (GSE59973) included three ESCC tissue and three nonpathogenic esophageal tissue samples that were tested against data from the GPL16770 Agilent-031181 Unrestricted_Human_miRNA_V16.0_microarray system.
2.2. Identification of DECs and DEMs
Differential analysis was conducted for GSE97051, GSE59973, and GSE131969 to obtain DEMs and DECs between ESCC and healthy control tissue samples by using GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/) . The interaction across miRNA microarray datasets was performed through Venny (http://bioinformatics.psb.ugent.be/webtools/Venn/). Using ENCORI (http://starbase.sysu.edu.cn/index.php), which is an open-source platform for studying RNA interactomes , the expression and overall survival of DEMs was examined. The criterion for DECs was an adjusted value <0.05 and , while the criterion for DEMs was an adjusted value <0.05 and . Volcano plots of the DEMs were created using Sangerbox (http://sangerbox.com/AllTools?tool_id=9699135). The miRTarBase was used to predict potentially significant mRNAs (http://mirtarbase.cuhk.edu.cn/php/index.php).
2.3. Functional Enrichment Analysis of Differentially Expressed Genes
GO and KEGG were used to identify the gene functions of individual genomic products and associated pathway information. GO evaluated the cellular component (CC), biological process (BP), and molecular function (MF) of each element [20, 21]. Annotation, visualization, and integrated discovery were performed using the online DAVID database tool (https://david.ncifcrf.gov/). A value <0.05 was considered to be significant.
2.4. Construction of Protein-Protein Interaction (PPI) Networks and Validation of the Hub Genes
The STRING database was applied for PPI network construction. Using Cytoscape via the CytoHubba app, the top 50 candidates were considered to be hub genes. To determine the expression of the 50 hub genes in ESCC tissue and to understand the effect(s) of hub genes on prognosis, the expression and survival analyses of hub genes were performed using UALCAN (http://ualcan.path.uab.edu/index.html). The miRNA-mRNA subnetwork of interest was determined. was considered to indicate statistical significance.
2.5. Validation of the circRNAs and Construction of the Competitive Endogenous RNA Network
StarBase 3.0 was introduced to predict linkage between the miRNAs and circRNAs. Using Venny against the GEO dataset, DECs of interest were predicted. The intersection between circRNAs and these DECs indicated the target circRNAs in ESCC tissue. Finally, a circRNA-miRNA-mRNA network was constructed, and the data were visualized using Cytoscape software 3.7.1.
3.1. Two Potential miRNAs Selected for ESCC
To find the potential miRNA biomarkers in ESCC, two miRNA datasets (GSE97051 and GSE59973) were selected (Table 1). As per the selection criteria, and as shown in Figure 1(d), two downregulated DEMs (miR-133b and miR-139-5p) were isolated. Next, the expression and overall survival of two downregulated DEMs in ESCA data were examined using ENCORI. Both miR-133b () and miR-139-5p () showed significantly low expression in cancer samples, but no effect on overall patient survival was seen in ESCA (Figure 2).
Note. T: tumor; N: normal.
3.2. Prediction and Analysis of mRNAs Binding to the Potential ESCC-Associated miRNAs
The miRTarBase database was used to identify which mRNAs might bind to the two predicted miRNAs. In total, 87 mRNAs potentially bind to miR-133b, and 105 mRNAs potentially bind to miR-139-5p. In order to identify key nodes in the PPI network, the top 50 hub genes were selected using CytoHubba based on the STRING database (Figure 3(c) and Table S1). For better visualization, the miRNA-mRNA connecting network was constructed using Cytoscape (Figure 3(d)).
3.3. Enrichment Analysis of the Target Genes
KEGG enrichment and GO analyses were performed for the target genes of the two DEMs (Figure 2). In the BP category, the “extrinsic apoptotic signaling pathway in absence of ligand,” “negative regulation of anoikis,” and “negative regulation of apoptotic process” were enriched. While “protein binding” was enriched according to the MF category, enrichment of the “nucleoplasm,” “nucleus,” and “cytosol” was shown in the CC category (Figure 3(b)). The KEGG pathway analysis revealed enrichment of the “focal adhesion,” “estrogen signaling,” “sphingolipid signaling,” and “PI3K-Akt signaling” pathways (Figure 3(a)).
3.4. Identification of mRNA Influences Patient Survival
To analyze the 50 mRNAs selected by the PPI network, survival analysis was performed using UALCAN. It was found that only RAP1B expression was significantly increased in ESCA tumor samples (). Additionally, the higher the expression of RAP1B, the worse the patient’s survival (). The miR-139-5p/RAP1B subnetwork was therefore determined.
3.5. Validation of the circRNAs and Construction of the Competitive Endogenous RNA Network
ENCORI was utilized to identify the 3766 potential circRNAs targeted by miR-139-5p. The intersection of these circRNAs and the predicted DECs found circ_0052867 to be of interest (Figure 4(c)). Finally, the circ_0052867/miR-139-5p/RAP1B network was constructed (Figure 4(d)) and was concluded to be a circRNA-miRNA-mRNA regulatory axis.
Due to the absence of obvious clinical symptoms in the early stages of ESCC and the lack of sensitive early diagnostic detection methods, most patients with ESCC are diagnosed at an advanced disease stage. Consequently, the best treatment opportunity is lost. Stein et al. reported that the five-year survival rate for early ESCA was 62.9% following surgery . Nevertheless, the validity and sensitivity in detecting traditional tumor markers are insufficient for early ESCC diagnoses . There is therefore an urgent need for improved biomarkers for early detection of ESCC.
Playing an important role in pathological and physiological processes, noncoding RNA comprises approximately 98% of the human genome . Many studies have demonstrated that competitive endogenous (ceRNAs) are able to act as sponges for miRNAs , while miRNAs have been shown to play an important role in tumorigenesis and its prognosis . Li et al. found that the survival time of high miRNA-506 expression was significantly shorter than that of low miRNA-506 expression . Downregulated miRNA-718 expression may also serve as a potential diagnostic marker for ESCC .
The circRNAs are newly found noncoding RNAs that participate in cancer development  and are currently being investigated as potential cancer biomarkers [30, 31]. Leading to poor prognosis in ESCC, upregulated circRNA_100876 expression can accelerate cell proliferation and metastasis , while the upregulation of circRNA_100873 is associated with increased lymphatic metastasis in ESCC . A further study reported that circGSK3 overexpression enhances multiple metastasis and invasion of the tumor . The underlying mechanisms of circRNAs and miRNAs in ESCC patient survival remain unknown.
Real-time quantitative PCR (q-PCR) is a classical method to detect abnormal noncoding RNA expression in ESCC tissue. Next generation sequencing techniques and available microarray data provide researchers with a more complete noncoding RNA expression profile between ESCC tumor and nonpathogenic tissue samples. For example, Yang et al. detected 39 dysregulated miRNAs (28 downregulated and 11 upregulated) in tumor/nonpathogenic ESCC tissue . Serum miRNA expression profiles were detected in 52 ESCC patients, and three miRNAs were identified as biomarkers for the diagnosis of ESCC .
In the present study, the intersection between two GEO databases confirmed miR-133b and miR-139-5p to be of interest, with an miRNA-mRNA-circRNA interaction being predicted. According to ceRNA theory, the expressions of miRNAs should be negatively correlated with the expressions of its targeted mRNAs and circRNAs. In this study, two downregulated DEMs were predicted to be associated with upregulated mRNA and circRNA targets. The PPI network was used to identify hub genes, and survival analysis was performed in order to find the key genes. Subsequently, RAP1B and circ_0052867 were found to have an impact on the survival of ESCC patients. Consequently, we constructed the circ_0052867/miR-139-5p/RAP1B ceRNA regulatory network, which may be a potential biomarker for ESCC patient survival. To our knowledge, this is the first study to report the expression of circ_0052867 in ESCC. While it has been found that the expression level of miR-139-5P in 11 cases of esophageal cancer is significantly higher than that of adjacent paracarcinoma tissue . Yang et al. reported that hsa-miR-139-5p was downregulated in ESCC tumor tissue compared to nontumor tissue samples . Jia et al. found that the expression of RAP1B was upregulated in ESCC tissue, and RAP1B promoted the growth, migration, and metastasis of the ESCC cells . However, Zhang et al. concluded that miR-518b may play an anticancer role in the development of ESCC by targeting RAP1B . Nevertheless, the exact mechanism remains to be discovered.
There were still some limitations to this study. Firstly, the findings of this study require further experimental verification. Secondly, the sample size is too small, and there are only a limited number of samples in each microarray dataset. Lastly, as mRNAs that influence patient survival and miRNAs that alter expression in ESCC tissue were identified in a database which includes ESCA tumor samples, there may be bias towards ESCC patients in this analysis.
The data used to support this study are available in TCGA database and GEO database.
Conflicts of Interest
The authors declare no conflicts of interest in this work.
Zhaojun Wang wrote the manuscript. Junhang Zhang conceived the study and provided guidance. All other authors contributed to the analysis and reviewed the manuscript. Zhaojun Wang and Haifeng Li contributed equally to this work.
The authors thank International Science Editing (http://www.internationalscienceediting.com) for editing this manuscript.
Table S1: 50 hub genes selected using CytoHubba ranked by MCC method. (Supplementary Materials)
- W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
- R. Siegel, J. Ma, Z. Zou, and A. Jemal, “Cancer statistics, 2014,” CA: A Cancer Journal for Clinicians, vol. 64, no. 1, pp. 9–29, 2014.
- M. Maghsudlu, E. Farashahi Yazd, and T. Amiriani, “Increased expression of MiR-27a and MiR-24-2 in esophageal squamous cell carcinoma,” Journal of Gastrointestinal Cancer, vol. 51, no. 1, pp. 227–233, 2020.
- A. Pennathur, M. K. Gibson, B. A. Jobe, and J. D. Luketich, “Oesophageal carcinoma,” The Lancet, vol. 381, no. 9864, pp. 400–412, 2013.
- W. Cao, W. Wu, F. Shi et al., “Integrated analysis of long noncoding RNA and coding RNA expression in esophageal squamous cell carcinoma,” International Journal of Genomics, vol. 2013, Article ID 480534, 10 pages, 2013.
- S. Mehta, A. Shelling, A. Muthukaruppan et al., “Predictive and prognostic molecular markers for cancer medicine,” Therapeutic Advances in Medical Oncology, vol. 2, no. 2, pp. 125–148, 2010.
- C.-Y. Li, W.-W. Zhang, J.-L. Xiang, X.-H. Wang, J. Li, and J.-L. Wang, “Identification of microRNAs as novel biomarkers for esophageal squamous cell carcinoma: a study based on the cancer genome atlas (TCGA) and bioinformatics,” Chinese Medical Journal, vol. 132, no. 18, pp. 2213–2222, 2019.
- C. Kandoth, M. D. Mclellan, F. Vandin et al., “Mutational landscape and significance across 12 major cancer types,” Nature, vol. 502, no. 7471, pp. 333–339, 2013.
- O. Brunetti, A. Russo, A. Scarpa et al., “MicroRNA in pancreatic adenocarcinoma: predictive/prognostic biomarkers or therapeutic targets?” Oncotarget, vol. 6, no. 27, pp. 23323–23341, 2015.
- F. R. Yang, H. J. Li, T. T. Li et al., “Prognostic value of microrna-15a in human cancers: a meta-analysis and bioinformatics,” BioMed Research International, vol. 2019, Article ID 2063823, 12 pages, 2019.
- A. Gnoni, D. Santini, M. Scartozzi et al., “Hepatocellular carcinoma treatment over sorafenib: epigenetics, microRNAs and microenvironment. Is there a light at the end of the tunnel?” Expert Opinion on Therapeutic Targets, vol. 19, no. 12, pp. 1623–1635, 2015.
- K. Danza, N. Silvestris, G. Simone et al., “Role of miR-27a, miR-181a and miR-20b in gastric cancer hypoxia-induced chemoresistance,” Cancer Biology & Therapy, vol. 17, no. 4, pp. 400–406, 2016.
- R. Pinto, S. Strippoli, S. De Summa et al., “MicroRNA expression in BRAF-mutated and wild-type metastatic melanoma and its correlation with response duration to BRAF inhibitors,” Expert Opinion on Therapeutic Targets, vol. 19, no. 8, pp. 1027–1035, 2015.
- H. Sun, L. Wang, Q. Zhao, and J. Dai, “Diagnostic and prognostic value of serum miRNA-1290 in human esophageal squamous cell carcinoma,” Cancer Biomarkers, vol. 25, no. 4, pp. 381–387, 2019.
- H. Yang, Y. N. Wei, J. Zhou et al., “MiR-455-3p acts as a prognostic marker and inhibits the proliferation and invasion of esophageal squamous cell carcinoma by targeting FAM83F,” European Review for Medical and Pharmacological Sciences, vol. 21, no. 14, pp. 3200–3206, 2017.
- W. Jin, W. Luo, W. Fang et al., “miR-145 expression level in tissue predicts prognosis of patients with esophageal squamous cell carcinoma,” Pathology—Research and Practice, vol. 215, no. 6, Article ID 152401, 2019.
- R. C. Li, S. Ke, F. K. Meng et al., “CiRS-7 promotes growth and metastasis of esophageal squamous cell carcinoma via regulation of miR-7/HOXB13,” Cell Death & Disease, vol. 9, no. 8, Article ID 838, 2018.
- Y. Qin, L. Li, F. Wang et al., “Knockdown of mir-135b sensitizes colorectal cancer cells to oxaliplatin-induced apoptosis through increase of FOXO1,” Cellular Physiology and Biochemistry, vol. 48, no. 4, pp. 1628–1637, 2018.
- W. Wang, W. Lou, B. Ding et al., “A novel mRNA-miRNA-lncRNA competing endogenous RNA triple sub-network associated with prognosis of pancreatic cancer,” Aging, vol. 11, no. 9, pp. 2610–2627, 2019.
- M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene ontology: tool for the unification of biology. The gene ontology consortium,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
- J. Du, Z. Yuan, Z. Ma, J. Song, X. Xie, and Y. Chen, “KEGG-PATH: kyoto encyclopedia of genes and genomes-based pathway analysis using a path analysis model,” Molecular BioSystems, vol. 10, no. 9, pp. 2441–2447, 2014.
- H. J. Stein, M. Feith, B. L. Bruecher, J. Naehrig, M. Sarbia, and J. R. Siewert, “Early esophageal cancer: pattern of lymphatic spread and prognostic factors for long-term survival after surgical resection,” Annals of Surgery, vol. 242, pp. 566–573, 2005.
- H. Zhang, H. Li, Q. Ma, F.-Y. Yang, and T.-Y. Diao, “Predicting malignant transformation of esophageal squamous cell lesions by combined biomarkers in an endoscopic screening program,” World Journal of Gastroenterology, vol. 22, no. 39, pp. 8770–8778, 2016.
- Y. Zhang, Y. Xu, L. Feng et al., “Comprehensive characterization of lncRNA-mRNA related ceRNA network across 12 major cancers,” Oncotarget, vol. 7, no. 39, pp. 64148–64167, 2016.
- W.-H. Xue, Z.-R. Fan, L.-F. Li et al., “Construction of an oesophageal cancer-specific ceRNA network based on miRNA, lncRNA, and mRNA expression data,” World Journal of Gastroenterology, vol. 24, no. 1, pp. 23–34, 2018.
- S. Zhu, M.-L. Si, H. Wu, and Y.-Y. Mo, “MicroRNA-21 targets the tumor suppressor gene tropomyosin 1 (TPM1),” Journal of Biological Chemistry, vol. 282, no. 19, pp. 14328–14336, 2007.
- S.-P. Li, H.-X. Su, D. Zhao, and Q.-L. Guan, “Plasma miRNA-506 as a prognostic biomarker for esophageal squamous cell carcinoma,” Medical Science Monitor, vol. 22, pp. 2195–2201, 2016.
- L. Sun, S. Dong, C. Dong et al., “Predictive value of plasma miRNA-718 for esophageal squamous cell carcinoma,” Cancer Biomarkers, vol. 16, no. 2, pp. 265–273, 2016.
- X. Li, L. Yang, and L.-L. Chen, “The biogenesis, functions, and challenges of circular RNAs,” Molecular Cell, vol. 71, no. 3, pp. 428–442, 2018.
- L. Chen, S. Zhang, J. Wu et al., “circRNA_100290 plays a role in oral cancer by functioning as a sponge of the miR-29 family,” Oncogene, vol. 36, no. 32, pp. 4551–4561, 2017.
- W. W. Du, W. Yang, X. Li et al., “A circular RNA circ-DNMT1 enhances breast cancer progression by activating autophagy,” Oncogene, vol. 37, no. 44, pp. 5829–5842, 2018.
- S. Cao, G. Chen, L. Yan, L. Li, and X. Huang, “Contribution of dysregulated circRNA_100876 to proliferation and metastasis of esophageal squamous cell carcinoma,” OncoTargets and Therapy, vol. 11, pp. 7385–7394, 2018.
- B. Zheng, Z. Wu, S. Xue et al., “hsa_circRNA_100873 upregulation is associated with increased lymphatic metastasis of esophageal squamous cell carcinoma,” Oncology letters, vol. 18, no. 6, pp. 6836–6844, 2019.
- X. Hu, D. Wu, X. He et al., “circGSK3β promotes metastasis in esophageal squamous cell carcinoma by augmenting β-catenin signaling,” Molecular Cancer, vol. 18, no. 1, p. 160, 2019.
- H. Yang, H. Su, N. Hu et al., “Integrated analysis of genome-wide miRNAs and targeted gene expression in esophageal squamous cell carcinoma (ESCC) and relation to prognosis,” BMC Cancer, vol. 20, no. 1, p. 388, 2020.
- D. Zheng, Y. Ding, Q. Ma et al., “Identification of serum MicroRNAs as novel biomarkers in esophageal squamous cell carcinoma using feature selection algorithms,” Frontiers in Oncology, vol. 8, p. 674, 2018.
- W. Jiao, J. Zhang, Y. Wei et al., “MiR-139-5p regulates VEGFR and downstream signaling pathways to inhibit the development of esophageal cancer,” Digestive and Liver Disease, vol. 51, no. 1, pp. 149–156, 2019.
- M. Yang, R. Liu, J. Sheng et al., “Differential expression profiles of microRNAs as potential biomarkers for the early diagnosis of esophageal squamous cell carcinoma,” Oncology Reports, vol. 29, no. 1, pp. 169–176, 2013.
- Z. Jia, Y. Yang, Z. Dengyan et al., “RAP1B, a DVL2 binding protein, activates Wnt/beta-catenin signaling in esophageal squamous cell carcinoma,” Gene, vol. 611, pp. 15–20, 2017.
- M. Zhang, S. Zhou, L. Zhang et al., “miR-518b is down-regulated, and involved in cell proliferation and invasion by targeting Rap1b in esophageal squamous cell carcinoma,” FEBS Letters, vol. 586, no. 19, pp. 3508–3521, 2012.
Copyright © 2020 Zhaojun Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.