Table of Contents Author Guidelines Submit a Manuscript
BioMed Research International
Volume 2018 (2018), Article ID 8518563, 9 pages
Research Article

Profiling and Bioinformatic Analyses Indicate Differential circRNA and miRNA/isomiR Expression and Interactions

1Department of Bioinformatics, School of Geographic and Biologic Information, Nanjing University of Posts and Telecommunications, Nanjing 210023, China
2Department of Cardiothoracic Surgery, The Third Affiliated Hospital of Soochow University, Changzhou 213003, China
3Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing 211166, China

Correspondence should be addressed to Yang Zhao and Qianyun Wang

Received 18 October 2017; Revised 28 December 2017; Accepted 1 January 2018; Published 25 February 2018

Academic Editor: Ernesto Picardi

Copyright © 2018 Li Guo 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.


As a novel class of noncoding RNAs, circular RNAs (circRNAs) have been reported to play a role in various biological processes. Some circRNAs may serve as microRNA (miRNA) sponges, regulating transcription or splicing. Herein, we investigated the expression profiles and interactions of miRNAs/isomiRs and circRNAs in male patients with esophageal cancer. We found that some miRNA genes generated two deregulated miRNA products (miR-#-5p and miR-#-3p), and these products were consistently abnormally expressed. Some circRNAs were predicted to be miRNA sponges for specific miRNAs. Some of these typically showed opposing expression patterns in cancer tissues: one upregulated and the other downregulated. Although fewer miRNAs were predicted to interact with circRNAs, the number of predicted interactions would be substantially increased if detailed isomiRs were involved. High sequence similarity across multiple isomiRs suggested that they might interact with circRNAs, similar to the interaction of homologous miRNAs with circRNAs. At the isomiR level, due to the characteristics of the sequences and expression patterns involved, the cross-talk between different ncRNAs is complicated despite simplification of the isomiRs involved through clustering. We expect that our results may provide methods for further study of the cross-talk among ncRNAs and elucidate their biological roles in human diseases.

1. Introduction

Esophageal cancer is a common cancer with high cancer mortality [1]. The impact of this disease is most severe in southern Africa, eastern Africa, and eastern Asia [2]. In China, the incidence of esophageal cancer in some areas is high [3]. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), may play crucial roles in the regulation of the occurrence and development of many cancers. Indeed, some ncRNAs have been proposed as biomarkers for the diagnosis of many cancers [4, 5].

Among the ncRNAs, miRNAs have been shown to play crucial regulatory roles in a variety of biological processes, particularly those related to pathological and physiological pathways [6, 7]. An increasing number of studies have shown that miRNA loci can generate miRNA variants (isomiRs) [815]. It has further been shown that these isomiRs are functional molecules that have expression and sequence heterogeneities [12, 1625]. Some small RNAs have been shown to have important regulatory roles in occurrence and development of esophageal cancer. For example, miR-330-3p regulates PDCD4 and might have an oncogenic role in esophageal squamous cell carcinoma [26], circuit of miR-31, SOX4, and EZH2 increases tumor progression in invasive esophageal carcinomas [27], and miR-10b may target KLF4 and thus promote the migration and invasion of tumorous cells [28]. These findings suggest that the small RNAs play important roles in esophageal cancer. However, few studies have investigated isomiRs, despite the fact that the canonical miRNA sequence is only one specific variant of multiple possible isomiRs [14, 29].

Another class of recently discovered animal RNAs, circRNAs, forms covalently closed continuous loops lacking 5′ and 3′ ends; these are also potential regulatory molecules [30]. Similar to miRNAs, circRNAs are thought to be crucial to cancer development and have thus been widely studied: a series of dysregulated circRNAs have been identified as important to the development of radiation resistance [31], circ_002059 is significantly downregulated in gastric carcinoma and may be a stable and novel biomarker [32], and circRNA_001569 may target miR-145 and thus contribute to colorectal cancer cell proliferation and invasion [33].

An understanding of miRNA-circRNA interactions is important for understanding the roles of ncRNAs in general, especially with respect to cross-talk among ncRNAs. Some circRNAs have been shown to have important regulatory roles in miRNA activity via sponging miRNAs. For example, the circRNA CDR1 acts as a miR-7 super-sponge [30, 3436]. Similar to miRNAs, the expression profiles of circRNAs are specific to different cell types and tissues [30, 37], indicating potential regulatory functions and roles in cellular processes. Some studies have suggested that circRNAs may be potential biomarkers for disease diagnosis [32, 38], based on their potential biological roles and on their stable features, such as closed RNA loops.

Although miRNAs and circRNAs have been studied in many cancers, few have been investigated at the isomiR level, despite the fact that multiple isomiRs are the detailed products from miRNA locus. Herein, we aim to understand the expression profiles of miRNA/isomiR and circRNA in esophageal cancer using deep sequencing data, focusing on their interactions and potential roles in cancer development. We also propose an approach for the analysis of miRNA/isomiR-circRNA cross-talk.

2. Materials and Methods

2.1. Patient Samples

A total of 5 paired tissue samples (tumor tissue and matched adjacent nontumor normal tissue) were taken from 5 male patients (age: ) with squamous esophageal cancer (in all cases, the distance between the tumor and nontumor normal tissue was >5 cm, Table S1). The 5 patients underwent surgical resection without prior radiotherapy or chemotherapy, in the Department of Cardiothoracic Surgery at the Third Affiliated Hospital of Soochow University, Changzhou, China. All samples from patients were HPV-negative and HPV has been suggested to play an important role in the development of esophageal cancer [39]. Informed consent was obtained from each subject. This study was approved by the ethics committee of Soochow University.

2.2. High-Throughput Sequencing and Analysis

Total RNA was extracted from each sample using Trizol (Invitrogen, USA). Purified RNA was quantitatively and qualitatively analyzed with NanoDrop (Thermo Company, USA). CircRNA and small RNA sequencing libraries were constructed with HiSeq 2000 sequencing platform following the manufacturer’s instructions.

CIRCexplorer software (version 1.1.7) [40] was used to identify circRNAs based on the fusion gene. After trimming adaptors, the small sequence fragments were mapped onto the human genome using Bowtie [41]. miRNAs/isomiRs were identified based on the canonical miRNA and pre-miRNA sequences annotated in miRBase version 21.0 ( [42]. IsomiRs were identified according to the location of the canonical miRNA in pre-miRNA with 4 nt [43], and further analyses were performed from the miRNA and isomiR levels, respectively. Differentially expressed circRNA and miRNA/isomiR profiles were obtained with edgeR, particularly those that were highly expressed. In the small RNA analysis, the canonical miRNA sequence and isomiRs from the specific miRNA locus were analyzed, respectively.

2.3. miRNA/isomiR-circRNA Interaction Analysis

The potential binding sites of abnormally expressed miRNAs were predicted using MiRanda [44] based on deregulated circRNAs. The potential relationships of different RNAs were constructed using Cytoscape version 3.4.0 [45]. It is possible to predict specific interactions at the isomiR levels, as well as based on the canonical miRNA sequence, or via clustering analysis of multiple isomiRs based on their close seed sequences (close sequence relationships can be obtained based on miRNA locus or independent miRNA locus) (Figure 1). Potential miRNA-circRNA interactions were also predicted and obtained from starBase v2.0 [46].

Figure 1: Flowchart for miRNA/isomiR-circRNA interaction. (a) Flowchart. Cancer and adjacent-normal samples were obtained from 5 male patients, and small RNA expression profiles and circRNA expression profiles were obtained via sequencing. (b) The detailed integrative analysis of small RNAs and circRNAs from the miRNA and isomiR levels.
2.4. Functional Analysis of Differentially Expressed miRNAs and isomiRs

For the top differentially expressed miRNAs (the canonical miRNA and isomiRs), we performed functional analysis using CapitalBio Molecule Annotation System V4.0 (MAS, Target mRNAs of the miRNAs were either identified from experimentally validated annotated in miRTarBase [47] or predicted with TargetScan [48]. Based on the same predictive methods, the target mRNAs of the isomiRs were predicted after clustering multiple isomiRs using seed sequences. We also further discussed the potential relationships among miRNAs/isomiRs, mRNAs, and circRNAs based on their functional relationships.

2.5. Statistical Analysis

Differentially expressed circRNAs and miRNA/isomiRs were obtained using edgeR based on paired t-test results with FDR correction: that is, when value 0.05 and the , we considered the RNA a deregulated species. Relative expressions were described using mean ± SD () based on normalized sequence counts. Specifically, the relative expression of the isomiR was further estimated based on percentage of multiple isomiRs in the specific miRNA locus, and those abundantly expressed isomiR species were specifically concerned in relevant analysis. Distributions of different small RNAs were estimated using Venny [49].

3. Results

3.1. Overview of circRNA and miRNA/isomiR Expression Profiles Showed Abundant ncRNA Expression

Although some circRNAs were detected in both of the paired tumor and normal samples, some were expressed differentially in the tumor samples as compared to the normal samples based on the detected circRNAs from sequencing data (Figures S1 and S2). Among the small ncRNAs with abundant detection, many miRNA/isomiRs showed significant abnormal expression (Figure S3), and the isoform profiles were consistent with the relative expression levels of the miRNA locus.

We identified several pre-miRNAs capable of generating two mature miRNAs from both 5p and 3p arms. Interestingly, in deregulated miRNAs, based on canonical miRNAs, we found that some pairs of miR-#-5p and miR-#-3p showed consistent downregulation, including miR-204-5p (, ) and miR-204-3p (, ). This pattern was not detected in upregulated miRNA species.

3.2. Deregulated Small RNAs Contributed to Basic Biological Progress Based on Functional Analysis

Further functional analysis suggested that deregulated small RNAs might contribute to some basic biological progress via the negative regulation of target mRNAs (Figure 2). For example, these small RNAs might be involved in processes such as single-multicellular organism process, system development, cell development, and cell differentiation. Although the role of miRNAs in these pathways was indirect, their flexible regulatory contributions might be important to many relevant biological processes.

Figure 2: Functional enrichment analysis based on deregulated miRNAs. (a) Enriched pathways based on GO (blue column: some of the top biological processes; green column: some of the top cellular components; red column: some of the top molecular functions). (b) Enriched pathways based on KEGG.
3.3. Predicted miRNA-circRNA Interaction Based on the Canonical miRNA

Based on known mature miRNAs and predicted circRNAs, we identified miRNA-circRNA pairs that potentially interact. Of these, some miRNAs were homologous with high sequence similarity (some nucleotides might vary, but these varied nucleotides were not located in the functional seed sequence). Homologous miRNAs might have been generated by homologous miRNA genes or by multiple pre-miRNAs. A total of 93 predicted circRNAs and 1,399 miRNAs might interact. Of these, 51% of them were oppositely expressed: one was upregulated as the other was downregulated. Due to their many-to-many interactions, many RNAs were repeatedly estimated. We found that 511 miRNAs and 59 circRNAs were downregulated. The miRNA-circRNA network based on potential relationships indicated that some circRNAs and miRNAs had a star-like distribution, suggesting potential crucial network positions (Figure 3).

Figure 3: Predicted miRNA-circRNA interaction network. Square indicates circRNA, and circle indicates miRNA. Red indicates significantly upregulated RNAs, green indicates significantly downregulated RNAs, and light blue indicates no significant change in expression.

A series of miRNA-circRNA interactions was obtained from starBase [46]. Based on the collected miRNA-circRNA interactions and the deregulated RNA molecules, we collected several abnormally expressed miRNAs, including 11 downregulated miRNAs (miR-124-3p, miR-129-5p, miR-135a-5p, miR-153-3p, miR-204-5p, miR-208a-3p, miR-211-5p, miR-218-5p, miR-488-3p, miR-490-3p, and miR-504-5p) and 1 upregulated miRNA (miR-373-3p). Most of the deregulated miRNAs with potential circRNA interactions were downregulated. Several circRNAs were detected in male patients with esophageal squamous cancer, mainly including circ-AGTPBP1, circ-ANKRD17, circ-GMPS, circ-PAPD4, and circ-SP140L. Interestingly, except for miR-373-3p:circ-SP140L (both upregulated), other miRNA:circRNA pairs had opposing deregulated expressions: in the cancer samples, miRNAs were downregulated and their paired circRNAs were upregulated (Table 1). Among the 97 deregulated circRNAs identified in starBase, 21 have been suggested to be competing endogenous RNAs (ceRNAs) of miRNAs. The paired miRNAs of these deregulated circRNAs were screened, focusing on the deregulated miRNAs without statistical difference (Table S2). A total of 19 pairs of potential circRNA:miRNA interactions were detected, and 14 of them (73.68%) were found to have opposing expression. Each member of these pairs might be involved in repeated circRNA and miRNA. For example, circ-SP140L was found to interact with 7 different miRNAs derived from several miRNA gene families.

Table 1: Selected deregulated miRNAs and their interacted deregulated circRNAs.

We also found that target gene of miRNA may generate circRNA as ceRNA of the miRNA. For example, circ-GMPS and miR-217 might interact, and GMPS was also a potential target of miR-217 according to TargetScan [50].

3.4. isomiR-circRNA Interaction at the isomiR Levels

Based on the miRNA-circRNA pairs identified, we further investigated potential isomiR-circRNA interaction at the isomiR levels. For the dominantly expressed miRNA gene, both miR-#-5p and miR-#-3p might generate multiple isomiRs with various expressions (Figure 4(a)). Due to the expression characteristics of isomiRs, most isomiRs sequences are similar to the canonical miRNA sequence (Figure 4(b)), despite differences in expression and involvement with additional nontemplate nucleotides. Based on the potential binding characteristics of miRNA and circRNA, all of the relevant isomiRs from a given miRNA locus can also bind to the circRNA (Figure 4(c)). Thus, at the isomiR level, it is possible to detect many isomiR:circRNA pairs based on close sequence relationships. In order to survey these interactions, multiple isomiRs could first be clustered based on miRNA loci or different miRNA loci (some homologous miRNA loci generate similar isomiR sequences) (Figure 1(b)). Clustered specific miRNA sequences with enriched expressions and sequences could be used to predict miRNA-circRNA interactions. Thus, many types of small RNAs may interact with a specific circRNA, and these relationships are complex.

Figure 4: An example of isomiR expression and interaction with circRNA. (a-b) Hsa-mir-204 has two mature products, but miR-204-5p is the dominant miRNA generated. From each miRNA locus, multiple isomiR can be detected with various expression patterns and sequences. However, some of these isomiRs are more dominant and others are rare. (c) Sequence distribution patterns and interactions with circRNA. Multiple isomiRs may interact with the same circRNA, as long as the sequence in the binding region is homologous.

4. Discussion

Bioinformatic analysis indicates that 276 miRNAs bind to circRNAs according to starBase [46], and 2,603 mature miRNAs have been reported and validated according to miRBase (version 21.0) [42]. Many miRNA gene families are involved in these interactions; this may be because homologous miRNAs always have higher similarity. Only 10.60% of all identified miRNAs potentially interact with circRNAs. Most miRNAs are not involved in miRNA-circRNA interactions, while miRNA-mRNA interactions are quite popular. miRNA-circRNA interactions are not specific, and many small RNAs may interact with different molecules. These interactions may be important to regulatory processes, particularly those regulated by the small RNA. Some circRNAs that have been identified as miRNA sponges binding miRNA and reducing miRNA activity include ciRS-7/CDR1 [30, 34], circ-Sry [34], circ-ITCH [51, 52], circ-Foxo3 [51], circ-HRCR [53], circ-HIPK3 [54], and hsa_circ_001569 [33]. The ability of circRNA to act as miRNA sponges may also be important for the regulatory function of small RNAs. Herein, we aim to understand the two noncoding RNAs expression profiles and their interactions in male patients with esophagus cancer, focusing on the interactions at the isomiR level.

Based on sequence data, we collected several miRNA-circRNA pairs. Based on widespread miRNA gene families and gene clusters, many miRNA sequences are highly similar any may thus have similar expression patterns and functional relationships. Even though few miRNA-circRNA interactions were predicted when compared to the large numbers of RNA species, the details of the corresponding relationships vary substantially. However, similar to miRNA-mRNA and miRNA-miRNA interactions, the actual functional binding of the two ncRNAs in vitro remains unclear. Indeed, competitive relationships may also exist between multiple isomiRs. Of these isomiRs, many are repeatedly estimated due to high sequence similarity. miRNA-circRNA pairs are not as common as miRNA-mRNA pairs, but the role of potential miRNA sponges should be not ignored. For example, ciRS-7 abrogates miR-7 [55], circRNA_010567 suppress miR-141 [56], and circRNA_100290 may act as a sponge for the miR-29 gene family [57].

The miRNA-circRNA network that we constructed demonstrates the complex interaction patterns between the two ncRNAs (Figure 2). This network would increase in complexity complex if we included the detailed isomiRs. A given miRNA locus may generate multiple isomiRs with difference sequence and expression patterns. These isomiRs mainly differ at the 3′ ends. Thus, these isomiRs may potentially interact with circRNAs as well as their canonical miRNA sequence and isomiRs with high sequence similarity and similar expression relationships that are similar to homologous miRNAs [29]. Therefore, at the isomiR level, more complex interactions with circRNAs are possible, and miRNA sponges may absorb many types of small RNAs with similar sequences. The interesting phenomenon seems to enlarge interaction profiles as well as current pattern: a circRNA may be miRNA sponge for homologous miRNAs in a miRNA gene family, while multiple isomiRs can be considered a special larger gene family with high sequence similarity. At the isomiR level, the isomiR-circRNA interaction may be quite common and the network may be quite complex. Indeed, except for the interaction between isomiRs and circRNAs, interactions between miRNAs or isomiRs are also quite common, such as coordinated or restricted interactions, and these horizontal and vertical interactions increase RNA cross-talk. Moreover, based on the deregulated miRNAs and circRNAs we identified, the paired circRNAs may be derived from target gene of miRNA, although this seems rare compared to the number of miRNA-cirRNA interactions. CircRNA may suppress miRNA to increase the expression level of the miRNA target [34], simultaneously generating ceRNAs to control miRNAs to ensure gene expression.

In conclusion, based on the small RNA and circRNA expression profiles in male patients with esophagus cancer, some deregulated RNA species may show potential functional relationship, and their expression patterns are always opposed. Although fewer miRNAs are predicted to interact with circRNAs, the interaction between two ncRNAs is more common if detailed analysis is performed at the isomiR level. Similar to homologous miRNAs, isomiRs increase analysis and cross-talk complexity.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Li Guo and Liang Zheng equally contributed to this work.


This work was supported by the National Natural Science Foundation of China (nos. 81373102 and 61771251), the Key Project of Social Development in Jiangsu Province (no. BE2016773), the National Natural Science Foundation of Jiangsu (no. BK20171443), and Top-Notch Academic Programs Project of Jiangsu Higher Education Institutions (TAPP: PPZY2015A067). The work was also supported by the Qing-lan Project of Jiangsu Province and the Excellent Young Teacher Project of Nanjing Medical University, Changzhou High-Level Medical Talents Training Project (no. 2016CZBJ043), and the Applied Basic Research Program of Changzhou (CJ20159030) and Sponsored by NUPTSF (nos. NY215068 and NY217100) and the Priority Academic Program Development of Jiangsu Higher Education Institution (PAPD).

Supplementary Materials

Figure S1: clustering analysis showed expression distribution of circRNAs across different samples. Figure S2: Venn’s distributions of expressed circRNAs between paired normal and tumor samples. Figure S3: clustering analysis and volcano distribution of miRNAs. Table S1: patients characteristics of 5 male patients with esophageal squamous cancer. Table S2: selected deregulated circRNAs and their interacted deregulated miRNAs. (Supplementary Materials)


  1. P. C. Enzinger and R. J. Mayer, “Esophageal Cancer,” The New England Journal of Medicine, vol. 349, no. 23, pp. 2241–2252, 2003. View at Publisher · View at Google Scholar · View at Scopus
  2. P. M. Forde and R. J. Kelly, “Chemotherapeutic and targeted strategies for locally advanced and metastatic esophageal cancer,” Journal of Thoracic Oncology, vol. 8, no. 6, pp. 673–684, 2013. View at Publisher · View at Google Scholar · View at Scopus
  3. G. D. Tran, X.-D. Sun, C. C. Abnet et al., “Prospective study of risk factors for esophageal and gastric cancers in the Linxian General Population Trial cohort in China,” International Journal of Cancer, vol. 113, no. 3, pp. 456–463, 2005. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Boro, D. Bauer, W. Born, and B. Fuchs, “Plasma levels of miRNA-155 as a powerful diagnostic marker for dedifferentiated liposarcoma,” American Journal of Cancer Research, vol. 6, no. 2, pp. 544–552, 2016. View at Google Scholar · View at Scopus
  5. F. R. Kulcheski, A. P. Christoff, and R. Margis, “Circular RNAs are miRNA sponges and can be used as a new class of biomarker,” Journal of Biotechnology, vol. 238, pp. 42–51, 2016. View at Publisher · View at Google Scholar
  6. G. A. Calin and C. M. Croce, “MicroRNA signatures in human cancers,” Nature Reviews Cancer, vol. 6, no. 11, pp. 857–866, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. A. Esquela-Kerscher and F. J. Slack, “Oncomirs—microRNAs with a role in cancer,” Nature Reviews Cancer, vol. 6, no. 4, pp. 259–269, 2006. View at Publisher · View at Google Scholar · View at Scopus
  8. J. Li, Z. Yang, B. Yu, J. Liu, and X. Chen, “Methylation protects miRNAs and siRNAs from a 3-end uridylation activity in Arabidopsis,” Current Biology, vol. 15, no. 16, pp. 1501–1507, 2005. View at Publisher · View at Google Scholar · View at Scopus
  9. P. Landgraf, M. Rusu, R. Sheridan et al., “A mammalian microRNA expression atlas based on small RNA library sequencing,” Cell, vol. 129, no. 7, pp. 1401–1414, 2007. View at Publisher · View at Google Scholar · View at Scopus
  10. R. D. Morin, M. D. O'Connor, M. Griffith et al., “Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells (Genome Research (2008) 18, (610-621)),” Genome Research, vol. 19, no. 5, p. 958, 2009. View at Google Scholar · View at Scopus
  11. J. G. Reid, A. K. Nagaraja, F. C. Lynn et al., “Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5-seed/cleavage/anchor regions and stabilize predicted mmu-let-7a:mRNA duplexes,” Genome Research, vol. 18, no. 10, pp. 1571–1581, 2008. View at Publisher · View at Google Scholar · View at Scopus
  12. S. L. Fernandez-Valverde, R. J. Taft, and J. S. Mattick, “Dynamic isomiR regulation in Drosophila development,” RNA, vol. 16, no. 10, pp. 1881–1888, 2010. View at Publisher · View at Google Scholar · View at Scopus
  13. L. W. Lee, S. Zhang, A. Etheridge et al., “Complexity of the microRNA repertoire revealed by next-generation sequencing,” RNA, vol. 16, no. 11, pp. 2170–2180, 2010. View at Publisher · View at Google Scholar · View at Scopus
  14. C. T. Neilsen, G. J. Goodall, and C. P. Bracken, “IsomiRs—the overlooked repertoire in the dynamic microRNAome,” Trends in Genetics, vol. 28, no. 11, pp. 544–549, 2012. View at Publisher · View at Google Scholar · View at Scopus
  15. L. Guo and F. Chen, “A challenge for miRNA: Multiple isomiRs in miRNAomics,” Gene, vol. 544, no. 1, pp. 1–7, 2014. View at Publisher · View at Google Scholar · View at Scopus
  16. S. Lu, Y. H. Sun, and V. L. Chiang, “Adenylation of plant miRNAs,” Nucleic Acids Research, vol. 37, no. 6, pp. 1878–1885, 2009. View at Publisher · View at Google Scholar · View at Scopus
  17. A. M. Burroughs, Y. Ando, M. J. L. De Hoon et al., “A comprehensive survey of 3 animal miRNA modification events and a possible role for 3 adenylation in modulating miRNA targeting effectiveness,” Genome Research, vol. 20, no. 10, pp. 1398–1410, 2010. View at Publisher · View at Google Scholar · View at Scopus
  18. S. Li, Y. Liao, M. Ho, K. Tsai, C. Lai, and W. Lin, “miRNA arm selection and isomiR distribution in gastric cancer,” BMC Genomics, vol. 13, supplement 1, p. S13, 2012. View at Publisher · View at Google Scholar · View at Scopus
  19. F. Llorens, M. Bañez-Coronel, L. Pantano et al., “A highly expressed miR-101 isomiR is a functional silencing small RNA,” BMC Genomics, vol. 14, no. 1, article no. 104, 2013. View at Publisher · View at Google Scholar · View at Scopus
  20. F. Llorens, M. Hummel, L. Pantano et al., “Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor,” BMC Genomics, vol. 14, no. 1, article 371, 2013. View at Publisher · View at Google Scholar · View at Scopus
  21. M. K. McGahon, J. M. Yarham, A. Daly et al., “Distinctive Profile of IsomiR Expression and Novel MicroRNAs in Rat Heart Left Ventricle,” PLoS ONE, vol. 8, no. 6, Article ID e65809, 2013. View at Publisher · View at Google Scholar · View at Scopus
  22. S. Babapoor, E. Fleming, R. Wu, and S. S. Dadras, “A novel miR-451a isomiR, associated with amelanotypic phenotype, acts as a tumor suppressor in melanoma by retarding cell migration and invasion,” PLoS ONE, vol. 9, no. 9, p. e107502, 2014. View at Publisher · View at Google Scholar · View at Scopus
  23. A. Schamberger and T. I. Orbán, “3 IsomiR species and DNA contamination influence reliable quantification of microRNAs by stem-loop quantitative PCR,” PLoS ONE, vol. 9, no. 8, Article ID e106315, 2014. View at Publisher · View at Google Scholar · View at Scopus
  24. G. C. Tan, E. Chan, A. Molnar et al., “5 isomiR variation is of functional and evolutionary importance,” Nucleic Acids Research, vol. 42, no. 14, pp. 9424–9435, 2014. View at Publisher · View at Google Scholar · View at Scopus
  25. M. Karali, M. Persico, M. Mutarelli et al., “High-resolution analysis of the human retina miRNome reveals isomiR variations and novel microRNAs,” Nucleic Acids Research, vol. 44, no. 4, pp. 1525–1540, 2016. View at Publisher · View at Google Scholar · View at Scopus
  26. H. Meng, K. Wang, X. Chen et al., “MicroRNA-330-3p functions as an oncogene in human esophageal cancer by targeting programmed cell death 4,” American Journal of Cancer Research, vol. 5, no. 3, pp. 1062–1075, 2015. View at Google Scholar
  27. R. B. Koumangoye, T. Andl, K. J. Taubenslag et al., “SOX4 interacts with EZH2 and HDAC3 to suppress microRNA-31 in invasive esophageal cancer cells,” Molecular Cancer, vol. 14, no. 1, article no. 24, 2015. View at Publisher · View at Google Scholar · View at Scopus
  28. Y. Tian, A. Luo, and Y. Cai, “MicroRNA-10b promotes migration and invasion through KLF4 in human esophageal cancer cell lines,” The Journal of Biological Chemistry, vol. 285, no. 11, pp. 7986–7994, 2010. View at Publisher · View at Google Scholar · View at Scopus
  29. L. Guo and T. Liang, “MicroRNAs and their variants in an RNA world: implications for complex interactions and diverse roles in an RNA regulatory network,” Briefings in Bioinformatics, p. bbw124. View at Publisher · View at Google Scholar
  30. S. Memczak, M. Jens, A. Elefsinioti et al., “Circular RNAs are a large class of animal RNAs with regulatory potency,” Nature, vol. 495, no. 7441, pp. 333–338, 2013. View at Publisher · View at Google Scholar · View at Scopus
  31. H. Su, F. Lin, X. Deng et al., “Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells,” Journal of Translational Medicine, vol. 14, no. 1, article no. 225, 2016. View at Publisher · View at Google Scholar · View at Scopus
  32. P. Li, S. Chen, H. Chen et al., “Using circular RNA as a novel type of biomarker in the screening of gastric cancer,” Clinica Chimica Acta, vol. 444, pp. 132–136, 2015. View at Publisher · View at Google Scholar · View at Scopus
  33. H. Xie, X. Ren, S. Xin et al., “Emerging roles of circRNA_001569 targeting miR-145 in the proliferation and invasion of colorectal cancer,” Oncotarget , vol. 7, no. 18, pp. 26680–26691, 2016. View at Publisher · View at Google Scholar · View at Scopus
  34. T. B. Hansen, T. I. Jensen, B. H. Clausen et al., “Natural RNA circles function as efficient microRNA sponges,” Nature, vol. 495, no. 7441, pp. 384–388, 2013. View at Publisher · View at Google Scholar · View at Scopus
  35. H. Xu, S. Guo, W. Li, and P. Yu, “The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells,” Scientific Reports, vol. 5, Article ID 12453, 2015. View at Publisher · View at Google Scholar · View at Scopus
  36. L. Yu, X. Gong, L. Sun, Q. Zhou, B. Lu, and L. Zhu, “The circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 expression,” PLoS ONE, vol. 11, no. 7, Article ID e0158347, 2016. View at Publisher · View at Google Scholar · View at Scopus
  37. J. Salzman, R. E. Chen, M. N. Olsen, P. L. Wang, and P. O. Brown, “Cell-type specific features of circular RNA expression,” PLoS Genetics, vol. 9, no. 9, Article ID e1003777, 2013. View at Publisher · View at Google Scholar · View at Scopus
  38. S. J. Conn, K. A. Pillman, and G. J. Goodall, “Identification of nuclear RNA-binding proteins as critical factors in circular RNA (circRNA) biogenesis,” Febs Journal, vol. 281, pp. 33-34, 2014. View at Google Scholar
  39. F. Chang, S. Syrjänen, Q. Shen, H. Ji, and K. Syrjänen, “Human papillomavirus (HPV) DNA in esophageal precancer lesions and squamous cell carcinomas from China,” International Journal of Cancer, vol. 45, no. 1, pp. 21–25, 1990. View at Publisher · View at Google Scholar · View at Scopus
  40. X.-O. Zhang, H.-B. Wang, Y. Zhang, X. Lu, L.-L. Chen, and L. Yang, “Complementary sequence-mediated exon circularization,” Cell, vol. 159, no. 1, pp. 134–147, 2014. View at Publisher · View at Google Scholar · View at Scopus
  41. B. Langmead and S. L. Salzberg, “Fast gapped-read alignment with Bowtie 2,” Nature Methods, vol. 9, no. 4, pp. 357–359, 2012. View at Publisher · View at Google Scholar · View at Scopus
  42. A. Kozomara and S. Griffiths-Jones, “miRBase: annotating high confidence microRNAs using deep sequencing data,” Nucleic Acids Research, vol. 42, pp. D68–D73, 2013. View at Publisher · View at Google Scholar
  43. L. Guo, Q. Yang, and J. Lu, “A comprehensive survey of miRNA repertoire and 3 addition events in the placentas of patients with pre-eclampsia from high-throughput sequencing,” PLoS ONE, vol. 6, no. 6, Article ID e21072, 2011. View at Publisher · View at Google Scholar · View at Scopus
  44. D. Betel, A. Koppal, P. Agius, C. Sander, and C. Leslie, “Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites,” Genome Biology, vol. 11, no. 8, article R90, 2010. View at Publisher · View at Google Scholar · View at Scopus
  45. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software Environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at Publisher · View at Google Scholar · View at Scopus
  46. J.-H. Li, S. Liu, H. Zhou, L.-H. Qu, and J.-H. Yang, “StarBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein—RNA interaction networks from large-scale CLIP-Seq data,” Nucleic Acids Research, vol. 42, no. 1, pp. D92–D97, 2014. View at Publisher · View at Google Scholar · View at Scopus
  47. S.-D. Hsu, F.-M. Lin, W.-Y. Wu et al., “MiRTarBase: a database curates experimentally validated microRNA-target interactions,” Nucleic Acids Research, vol. 39, supplement 1, pp. D163–D169, 2011. View at Publisher · View at Google Scholar · View at Scopus
  48. B. P. Lewis, I.-H. Shih, M. W. Jones-Rhoades, D. P. Bartel, and C. B. Burge, “Prediction of mammalian microRNA targets,” Cell, vol. 115, no. 7, pp. 787–798, 2003. View at Publisher · View at Google Scholar · View at Scopus
  49. J. C. Oliveros, VENNY. An interactive tool for comparing lists with Venn Diagrams., Available,
  50. V. Agarwal, G. W. Bell, J.-W. Nam, and D. P. Bartel, “Predicting effective microRNA target sites in mammalian mRNAs,” eLife, vol. 4, no. 2015, Article ID e05005, 2015. View at Publisher · View at Google Scholar · View at Scopus
  51. G. Huang, H. Zhu, Y. Shi, W. Wu, H. Cai, and X. Chen, “Cir-ITCH plays an inhibitory role in colorectal cancer by regulating the Wnt/β-Catenin Pathway,” PLoS ONE, vol. 10, no. 6, article e0131225, 2015. View at Publisher · View at Google Scholar · View at Scopus
  52. F. Li, L. Zhang, W. Li et al., “Circular RNA ITCH has inhibitory effect on ESCC by suppressing the Wnt/β-catenin pathway,” Oncotarget , vol. 6, no. 8, pp. 6001–6013, 2015. View at Publisher · View at Google Scholar · View at Scopus
  53. K. Wang, B. Long, F. Liu et al., “A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223,” European Heart Journal, vol. 37, no. 33, pp. 2602–2611, 2016. View at Publisher · View at Google Scholar
  54. Q. Zheng, C. Bao, W. Guo et al., “Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs,” Nature Communications, vol. 7, Article ID 11215, 2016. View at Publisher · View at Google Scholar
  55. H. Pan, T. Li, Y. Jiang et al., “Overexpression of Circular RNA ciRS-7 Abrogates the Tumor Suppressive Effect of miR-7 on Gastric Cancer via PTEN/PI3K/AKT Signaling Pathway,” Journal of Cellular Biochemistry, 2017. View at Publisher · View at Google Scholar · View at Scopus
  56. B. Zhou and J.-W. Yu, “A novel identified circular RNA, circRNA_010567, promotes myocardial fibrosis via suppressing miR-141 by targeting TGF-β1,” Biochemical and Biophysical Research Communications, vol. 487, no. 4, pp. 769–775, 2017. View at Publisher · View at Google Scholar · View at Scopus
  57. 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. View at Publisher · View at Google Scholar · View at Scopus