==== Front Noncoding RNA ResNoncoding RNA ResNon-coding RNA Research2468-0540KeAi Publishing S2468-0540(17)30047-110.1016/j.ncrna.2018.02.002ArticleCircular RNAs in the cardiovascular system Gomes Clarissa P.C. a1Salgado-Somoza Antonio a1Creemers Esther E. bDieterich Christoph cLustrek Mitja dDevaux Yvan yvan.devaux@lih.lua∗Cardiolinc™ network 2a Cardiovascular Research Unit, Luxembourg Institute of Health, Luxembourg, Luxembourgb Experimental Cardiology, Academic Medical Center, Amsterdam, The Netherlandsc German Center for Cardiovascular Research, University Hospital Heidelberg, Heidelberg, Germanyd Department of Intelligent Systems, Jozef Stefan Institute, Ljubljana, Slovenia∗ Corresponding author. Luxembourg Institute of Health, Cardiovascular Research Unit, 84, Val Fleuri, L-1526, Luxembourg, Luxembourg. yvan.devaux@lih.lu1 These authors contributed equally to the work. 2 www.cardiolinc.org. 25 2 2018 3 2018 25 2 2018 3 1 1 11 31 10 2017 16 1 2018 22 2 2018 © 2018 The Authors2018This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).Until recently considered as rare, circular RNAs (circRNAs) are emerging as important regulators of gene expression. They are ubiquitously expressed and represent a novel branch of the family of non-coding RNAs. Recent investigations showed that circRNAs are regulated in the cardiovascular system and participate in its physiological and pathological development. In this review article, we will provide an overview of the role of circRNAs in cardiovascular health and disease. After a description of the biogenesis of circRNAs, we will summarize what is known of the expression, regulation and function of circRNAs in the cardiovascular system. We will then address some technical aspects of circRNAs research, discussing how artificial intelligence may aid in circRNAs research. Finally, the potential of circRNAs as biomarkers of cardiovascular disease will be addressed and directions for future research will be proposed. Keywords Non-coding RNAsCircular RNAsCardiovascular systemCardiovascular diseaseArtificial intelligenceBiomarkerAbbreviations circRNAs, circular RNAsCRISPR, clustered regularly interspaced short palindromic repeatsCV, cardiovascularDCM, dilated cardiomyopathyEMT, epithelial-mesenchymal transitionlncRNAs, long non-coding RNAsmiRNAs, microRNAsncRNAs, non-coding RNAsRNA-seq, RNA sequencingRPAD, RNase R treatment followed by polyadenylation and poly(A)+ RNA depletionRT-qPCR, reverse transcription quantitative polymerase chain reaction ==== Body 1 Introduction February 2001 marked the completion of the sequencing of the human genome [1,2]. Eleven years later, researchers from the Encode Project Consortium revealed that, while 80% of our genome is transcribed, only 2% of expressed RNAs encode proteins [3]. This landmark discovery implied that the vast majority of our genes are transcribed into non-protein coding RNAs, more simply named non-coding RNAs (ncRNAs). Multiple categories of ncRNAs have been characterized and classified. The most widely used classification is based on the size of ncRNAs, the shorter ones (<200 nucleotides) being classified as small ncRNAs, while those exceeding 200 nucleotides are grouped as long non-coding RNAs (lncRNAs) (Fig. 1). Inside the lncRNAs family, both linear and circular forms of RNAs can be found. Although circular RNAs (circRNAs) were previously thought to be rare as compared to linear messenger RNAs, recent work showed that circRNAs can be relatively abundant and be the predominant transcripts of hundreds of genes [4]. These circular molecules are ubiquitously expressed in human and more than 30,000 different circRNAs have been identified so far [5]. They are regulated in several disease states and contribute to disease development and progression. The role of circRNAs in the cardiovascular system in only emerging.Fig. 1 The non-coding RNA family. Non-coding RNAs are classified depending on their size with an arbitrary cut-off of 200 nucleotides. Abbreviations: ncRNAs: noncoding RNAs; lncRNA: long noncoding RNA; eRNA: enhancer RNA; circRNA: circular RNA; rRNA: ribosomal RNA; snRNA: small nuclear RNA; scaRNA: small cajal body-specific RNA; snoRNA: small nucleolar RNA; tiRNA: tRNA-derived stress-induced small RNA; tRNA: transfer RNA; RNAi: RNA interference; siRNA: small interfering RNA; piRNA: piwi-interacting RNA; miRNA: microRNA; rasiRNA: repeat associated small interfering RNA. Fig. 1 Cardiovascular (CV) disease is still the main cause of death and disabilities worldwide [6,7]. While research has made enormous progress in terms of healthcare and treatment of CV disease, the use of RNAs for its treatment and diagnosis is facing technical challenges that slow down their clinical application. Recent developments revealed the potential of microRNAs (miRNAs) for the management of CV disease patients [8,9]. LncRNAs have also emerged as potential diagnostic and therapeutic targets of CV disease [10,11]. Yet, the knowledge of the role of circRNAs in the diseased heart and their value as CV disease markers and therapeutic targets are still in their infancy. This review article is part of a Special Issue on “Non-coding RNAs in Cardiovascular Health and Disease” and will focus on circRNAs in the CV system. After a description of the mechanisms leading to circRNA biogenesis, we will gather the knowledge of the expression, regulation and function of circRNAs in the CV system. We will then address some methodological aspects of circRNAs research and we will provide some considerations of the potential of artificial intelligence to aid in circRNAs research. The value of circRNAs as biomarkers of CV disease will be delineated and we will finally propose directions for future research. 2 Biogenesis of circRNAs CircRNAs are produced by the spliceosome machinery by direct back-splicing of exons or introns within pre-mRNA (Fig. 2) [12] or by exon skipping [13]. As a result of a direct back-splicing event, an exon is not associated with an adjacent downstream exon, as seen in linear splicing, but with an upstream exon or intron. This gives rise to circular RNA molecules harboring exons or introns that are out-of-order from the genomic context [14,15]. Exon skipping involves a lariat that undergoes internal splicing, removing the intron and generating a circRNA. Both proposed mechanisms involve the spliceosome machinery, but the mechanism by which the spliceosome selects only certain exons to circularize is largely unknown. It is however known that many circRNAs in mammals are processed from exons with long flanking introns, containing reverse complementary sequences [16,17]. These sequences can be inverted repeats, such as very abundant Alu elements (∼300 base pairs repeats), or non-repetitive complementary sequences [16]. The common feature is that the intronic sequences base-pair to one another, thereby bringing the splice sites into close proximity of each other, allowing the spliceosome to perform a back-splicing reaction. An alternative mechanism to bridge flanking introns is through interaction with RNA binding proteins that bind cis-regulatory regions in these introns.Fig. 2 Biogenesis of circRNAs. The spliceosome machinery, which normally catalyses linear splicing of pre-mRNA, can also perform a back-splicing reaction between two exons, resulting in the formation of a circRNA. Back-splicing uses the same canonical splicing machinery and canonical splice sites as needed for linear splicing. Mechanistically, back-splicing requires that the donor and acceptor site of the back-spliced exons are brought in close proximity to each other, which can be accomplished by direct RNA base-pairing of reverse complementary sequences in the introns flanking the back-spliced exons, or by the interaction of RNA binding proteins (RBPs) that dock on these flanking introns. Fig. 2 Among the hundreds of RNA binding proteins characterized so far [18], the first one identified as a regulator of circRNA production was Muscleblind (MBL), a zinc finger protein which is known to regulate alternative splicing in muscle [14] and for which dysregulation can lead to pathologies, particularly myotonic dystrophies (reviewed in Ref. [19]). Interestingly, MBL regulates circRNA production from its own pre-mRNA. Numerous MBL-binding sites are present in the introns flanking the back-spliced junction of circMBL, which turned out to be essential for the biogenesis of circMBL. This raises the possibility of a sophisticated control mechanism for MBL protein levels. When MBL protein levels are high, it decreases the production of its own mRNA by promoting circMBL production. The circRNA can subsequently sponge out the excess MBL protein by binding to it [14] (more details about this function in the next section). The second RNA binding protein shown to regulate circRNA production by pre-mRNA binding was Quaking (QKI) [20], a splicing factor that has been implicated in diseases including ataxia, schizophrenia, and glioblastoma [21]. The insertion of synthetic QKI-binding sites into introns was found to be sufficient to produce circRNAs. Mechanistically it was shown that QKI binds to two flanking introns and subsequently dimerizes to bridge the circularizing exons. During epithelial-mesenchymal transition (EMT) it was shown that about one-third of the circRNAs were regulated by QKI [20]. Whether QKI is also required for circRNA formation in the CV system is currently unknown. Among the post-transcriptional modifications that increase molecular diversity and can lead to alternative splicing and changes in gene expression levels, called RNA editing, the prevailing modification in higher eukaryotes is the single base substitution of adenosine to inosine mediated by the adenosine deaminase acting on RNA (ADAR) enzyme family [22,23]. A recent study revealed that the RNA-editing enzyme ADAR1 functions as an antagonist of circRNA production [24,25]. Circularized exons are normally between long introns, often termed bracketing introns. It was found that introns bracketing circRNAs are enriched for RNA editing events and that knockdown of ADAR1 significantly enhanced circRNA expression in HEK293 cells. Repression of RNA circularization is probably not a general function of ADAR1, as only a few circRNAs were affected by ADAR1 knockdown [25]. RNA Binding Motif protein 20 (RBM20) is yet another splicing factor that has recently been implicated in circRNA production, specifically from the titin gene (TTN) [26]. 3 Function of circRNAs in the CV system In the following lines, we will summarize the current knowledge of the roles of circRNAs in CV pathophysiology, focusing on their mechanisms of action and interaction with partner molecules. A schematic representation of the different functions of circRNAs is provided in Fig. 3.Fig. 3 Functions of circRNAs. a) In the presence of internal ribosomal enter sites and a suitable open reading frame, circRNAs could be translated into proteins. b) Generation of circRNAs by back-splicing of the precursor RNA can lead to a regulation of the linear/circular ratio of the hosting gene. c) CircRNAs can sponge microRNAs to act as a cytoplasmic reservoir of microRNAs. The scavenging of microRNAs removes the repression of target RNAs leading to an increase in their translation (if mRNAs) or activity (lncRNAs). d) CircRNAs can act as a scaffold for cytoplasmic proteins, retain certain transcription factors in the cytoplasm or serve as a vehicle for the transport of these molecules. e) CircRNAs could undergo degradation after being targeted by microRNAs. f) CircRNAs could function as a scaffold for transcription factors leading them to specific locations of the genome or directly interact with the DNA to regulate transcription. Fig. 3 3.1 circRNAs: coding or non-coding? First, we would like to discuss the “non-coding RNA” feature of circRNAs. While circRNAs are widely considered as lacking protein-coding potential, recent reports questioned this concept. Translation of circRNAs (Fig. 3a) was suggested by ribosome footprinting data in several animal models (from flies to murine ones), which revealed that a number of circRNAs were consistently found to be associated with translating ribosomes [27]. In another study in murine and human myoblasts, cZNF609 (also known as myocardial infarction-associated circRNA – MICRA [28]) was reported to contain an open reading frame allowing for protein production after overexpression [29]. However, the protein-coding potential of circRNAs in naïve conditions may differ from artificial up-regulation and thus their capability to encode proteins is still a matter of debate. 3.2 Expression level As discussed in section 2, splicing events can determine the expression levels of circRNAs (Fig. 3b). Thus, the balance between circular and linear isoforms originated from the same hosting gene is modified by general splicing factors and the regulation of back splicing mechanisms [14]. 3.3 Sponge for miRNAs More attention has been brought to the ability of circRNAs to bind miRNAs (Fig. 3c), ∼22 nucleotide ncRNAs that regulate transcription [8,9]. The circRNA antisense of the cerebellar degeneration-related protein 1 (Cdr1as, also known as ciRS-7) harbours more than 60 potential binding sites for miR-7 [30,31]. The sequestration of miR-7 by Cdr1as, a process termed “sponge effect”, relates to insulin production and secretion by beta cells [32]. Furthermore, Cdr1as is elevated after myocardial infarction in mice and its effect on apoptosis seems to be regulated via miR-7 [33]. Also in mice, the heart-related circRNA (Hrcr) scavenges the pro-hypertrophic miR-223 and thus blocks the isoproterenol-induced cardiomyocyte hypertrophy [34]. Although miR-223 is up-regulated in human failing hearts [34], whether it is sequestered by HRCR in human is still unknown. In serum samples from patients with carotid artery dissection, the ratio between the circRNA circR-284 and miR-221 is altered due to an increase of circR-284 and a decrease of miR-221 [35]. In thoracic aortic dissection patients, an increase of matrix metalloproteinase 9 expression is observed and bioinformatics tools predicted that hsa_circ_101283 could contribute to it by binding three miRNAs (mir-320a, miR-583-5p and miR-138-5p) [36]. In human aortic smooth muscle cells, circACTA2 regulates alpha-actin expression by interacting with miR-548f-5p [37]. In mouse cardiac fibroblasts, circRNA_000203 increases collagen production by binding miR-26b-5p [38]. The mitochondrial fission and apoptosis-related circRNA (MFACR) sponges miR-652-3p, upregulating its target, MTP18, consequently regulating mitochondrial fission and cardiomyocyte apoptosis in mice [39]. Finally, cZNF609 appears to regulate vascular endothelial dysfunction by sponging miR-615-5p [40]. The therapeutic potential of these findings remains to be demonstrated. 3.4 Scaffold CircRNAs could also behave as scaffolds, either drive other elements in the cell towards their targets or block them to prevent their action (Fig. 3d). As an example of the latter, circ-Foxo3 sequesters proteins like Id1, E2f1, Fak, and Hif1a in the cytoplasm of cardiomyocytes and cardiac fibroblasts, hence avoiding their translocation to the nucleus and contributing to the ageing of those cells [41]. In a mouse cell line, circ-Foxo3 showed an affinity for cell cycle-related proteins such as Cdk2 and p21, and to a lesser extent Cdk6, p16 and p27, also suggesting a role of this circRNA in cell proliferation [42]. 3.5 Other putative functions of circRNAs The miRNA-sponge function of circRNAs described before is typically presumed for cytoplasmic circRNAs. Of note, although cZNF292 is cytoplasmic, a miRNA-sponge effect could not be demonstrated using high-throughput sequencing of RNA isolated by cross-linking and immunoprecipitation [43]. In addition to the sequestration of miRNAs by circRNAs, it can be foreseen that miRNAs could bind circRNAs to induce their degradation (Fig. 3e), similar to the down-regulation of linear RNAs. However, this needs further investigation. Some RNA molecules are able to target genomic DNA in a sequence-specific manner. This characteristic is the basis of the gene editing mechanism named Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), where a small guide RNA places the CRISPR associated protein (Cas9) at a specific genomic location [44]. Together with their scaffold ability, one could speculate that circRNAs can act as transcription regulators, either alone or in combination with other RNA or protein partners (Fig. 3f). Therefore, circRNAs could target specific DNA sequences and modulate their transcription, as it was previously described for lncRNAs [10,45]. 4 Expression and regulation of circRNAs in the CV system Although circRNAs have been described at least as early as 1976 as viroids [46] and then in the 1990s in eukaryotes [47,48], interest in their function and application has been recently renewed [30,[49], [50], [51]]. RNA sequencing data of non-polyadenylated transcripts allowed the realization that circRNA expression is widespread [16,24,30,[50], [51], [52]]. Studies revealed that at least 14% of expressed genes in human fibroblasts yield over 25,000 circRNAs [51]. CircRNAs are abundant and conserved in eukaryotes, have been detected in different developmental stages and cell types, and their expression can also be highly cell-type dependent [30,50,51,53]. These observations together with their specific changes in pathological states (reviewed in Ref. [54]) suggest their expression is a regulated process [4]. Investigations of circRNA expression during human development revealed dynamic changes [55,56], particularly a trend in upregulation [56]. Interestingly, the most expressed gene during foetal development was the circular isoform of NCX1 in the heart, a calcium transport gene essential for cardiac development. In this study, several other circRNAs were more expressed than their linear mRNAs [56]. This observation is consistent with other reports that indicate that the circRNA expression is not necessarily correlated with the expression of its linear mRNA from which it derives, varying according to the cell type [25,51,53,56]. Two contributions towards a better understanding of the potential role of circRNAs in the CV system were the characterization of these molecules' expression in both developing and adult hearts [57,58]. One study detected over 9,000 candidate circRNAs for each species using RNA sequencing (RNA-seq) of heart tissue from humans (failing and non-failing), mice (sham or after transverse aortic constriction) and rats (neonatal and adult) [57]. The other study identified also by RNA-seq over 15,000 circRNAs in humans and over 3,000 in mice. The group analysed tissue from failing and non-failing hearts from both species and from human embryonic stem cell-derived cardiomyocytes across different stages of differentiation [58]. The most expressed circRNAs, which were more abundant than their linear counterparts, were produced from critical cardiac genes, such as TTN, ryanodine receptor 2 (RYR2 – encoding a protein found primarily in cardiac muscle and involved in calcium release from the cytoplasm), and dystrophin (DMD – encoding a component of the dystrophin-glycoprotein complex which links the cytoskeleton to the extracellular matrix and which genetic alterations cause Duchenne muscular dystrophy) [58]. The titin gene encodes titin (TTN), a giant protein involved in muscle contraction. RBM20, involved in circRNA biogenesis, is expressed in the heart and skeletal muscle, where it regulates alternative splicing of multiple genes, including TTN [59]. Mutations in the RBM20 gene lead to a severe and early onset dilated cardiomyopathy (DCM), and aberrant splicing of TTN is believed to be the main reason why RBM20 mutation carriers develop DCM [[59], [60], [61]]. However, a recent study showed that RBM20 is also required for the production of a large set of circRNAs that arise from the TTN gene [26]. The authors profiled circRNAs of human left ventricle samples and control from hypertrophic or dilated cardiomyopathy patients, finding 826 back-splice junctions in common between the groups, where 80 were from the titin transcript. The RBM20-binding sites were found to be enriched in the introns flanking the titin back-splice junctions. These RBM20-dependent TTN circRNAs specifically arise from a region within TTN transcript: the so-called I-band, which is known to undergo extensive alternative splicing. Interestingly, in the hearts of Rbm20 knockout mice, and also in human RBM20 mutation carriers, circRNAs from TTN's I-band were completely absent, while at the same time, the corresponding exons were abundantly included in the linear TTN transcript. This implies a mechanism in which exons that are spliced out of the TTN pre-mRNA can serve as a substrate for the production of circRNAs [26]. In conclusion, this study shows that TTN circRNA production is intimately connected to alternative splicing of TTN, and this is in line with other studies demonstrating that the more an exon is circularized, the less it is represented in the linearly processed mRNA [14,62]. A microarray study contributed to the identification of 63 differentially expressed circRNAs from 1,163 circRNA candidates in murine tissues from normal and failing hearts due to myocardial infarction [63]. Twelve of these differentially expressed circRNAs were validated by quantitative PCR, including Cdr1as. The authors observed an upregulation in failing compared to normal hearts, a result consistent with another report showing that Cdr1as was overexpressed in mice hearts with increased infarct size and cardiomyocytes under hypoxia. Additionally, Cdr1as acted as a sponge to miR-7a in myocardial cells, suppressing its activity [33]. Few other circRNAs have been reported correlated with cardiovascular disease. In an analysis of plasma from myocardial infarction patients, the circRNA circRNA_081881 was down-regulated more than ten-fold compared to plasma from healthy volunteers and was predicted to bind miR-548, which targets the heart-protective gene PPARγ [64]. The circRNA produced from the forkhead family of transcription factors Foxo3, circ-Foxo3, was shown to promote cardiac senescence in hearts of aged patients and mice, where it was highly expressed. Silencing of the circRNA inhibited senescence of mouse embryonic fibroblasts inhibited senescence suppressing cardiomyopathy symptoms [41]. As diabetes mellitus is a prime risk factor for cardiovascular disease and cardiac fibrosis is one of the main pathologies of diabetic cardiomyopathies, studies investigating the role of circRNAs in such conditions are emerging. A report with diabetic mice hearts and Angiotensin II-induced cardiac fibrosis investigated the modulation of circRNA_000203 on fibrosis-associated genes expression. This circRNA was found to be upregulated compared to controls and to interact with the anti-fibrotic miR-26b-5p, enhancing the proliferation of cardiac fibroblast [65]. High levels of another circRNA_010567 were found in diabetic mice myocardium and cardiac fibroblasts, negatively regulating miR-141, and thus upregulating its target, TGF-β1. These three molecules play a role in regulating the pathogenesis of myocardial fibrosis in the studied model [66]. In vascular endothelial cells, computational analysis of RNA-seq data resulted in the identification of over 7,000 circRNAs, among which cZNF292 was shown to be highly expressed and up-regulated under hypoxic stress [43]. Importantly, cZNF292 was not a mere byproduct of RNA splicing but was functionally involved in angiogenic sprouting [43]. 5 Methodological aspects of circRNAs research Most circRNAs escaped identification for decades because they are not easily separated from other RNAs by size (as miRNAs can be), they do not have polyadenylation or free 3′ or 5′ ends, and they have exons arranged in a different order (back-splice), thus being considered as artefacts [12]. However, technological advances in RNA analyses and bioinformatic tools enabled circRNAs to be appreciated as an important type of ncRNAs and potential regulators of gene expression. The events that give rise to circRNAs specify certain properties of these molecules, such as the presence of a non-linear splice junction, the absence of a polyA tail and the resistance to exonuclease treatment, which distinguish them from linear mRNAs. These characteristics can be used to detect novel circRNAs. 5.1 Detection of back-splicing events Back-splicing junctions can be readily detected from ribosomal RNA (rRNA) depleted whole transcriptome sequencing data as opposed to conventional polyA + sequencing data (Fig. 4). The detection of circRNAs is a multi-stage bioinformatics process. Briefly, short sequencing reads, which are indicative of back-splicing events and, hence, do not align co-linearly to their genomic locus of origin, are identified in rRNA-depleted whole transcriptome sequencing data. Depending on tissue and condition, this will lead to a list of potential back-splicing junction candidates, which may represent hundreds to thousands of potential circRNAs. The quality of these predictions is typically assessed by their absence of polyA sequencing data (i.e. false positives) and by their relative enrichment in RNase R treated over mock-treated rRNA-depleted whole transcriptome sequencing data (Fig. 4) [67]. An alternative RNase R-based approach, which is generally applicable for all species where no rRNA removal strategy exists, has been proposed and termed RPAD (RNase R treatment followed by polyadenylation and poly(A)+ RNA depletion) [68]. Briefly, total RNA is digested by RNase R and remaining linear RNAs receive a polyA tail in a second step, which is then used to deplete them from the RNA pool. This method is more efficient to eliminate linear RNAs than treatment with RNase R alone, facilitating the detection of novel circRNAs [68].Fig. 4 Different sequencing strategies to detect back-splicing junctions and to recover the internal structure of circular RNAs. Fig. 4 A number of software tools exist for circRNA recognition through the detection of back-splicing junctions (Table 1). A comparison of the benefits and drawbacks of some of these tools has been recently reported [69]. However, publicly available downstream analysis tools for the reconstruction of circRNA species and additional in-depth analyses are rare.Table 1 Prediction tools for the discovery of circRNAs from RNA sequencing (RNA-seq) datasets. Table 1Tool name Website Description Aligner Reference Acfs https://github.com/arthuryxt/acfs Identifies back-splicing junctions from RNA-seq data BWA-MEM [84] circExplorer http://yanglab.github.io/CIRCexplorer/ Identifies junction reads from back-spliced exons and intron lariats TopHat Fusion [16] circRNA_finder https://github.com/bioxfu/circRNAFinder Pipeline to find circRNAs from RNA-seq data STAR [85] CIRI2 https://sourceforge.net/projects/ciri/files/CIRI2/ Detects back-splicing junction reads using efficient maximum likelihood estimation based on multiple seed matching BWA-MEM [86] DCC https://github.com/dieterich-lab/ Detects back-splicing junctions from RNA-seq data and applies filters to integrate data across replicate sets STAR [67] find_circ https://github.com/marvin-jens/find_circ Detects back-spliced sequencing reads from RNA-seq data bowtie2 [30] KNIFE https://github.com/lindaszabo/KNIFE Statistically-based splicing detection for circular and linear isoforms from RNA-Seq data bowtie and bowtie2 [56] MapeSplice http://www.netlab.uky.edu/p/bioinfo/MapSplice2 Maps RNA-seq data to reference genome to discover splice junctions bowtie [87] NCLScan https://github.com/TreesLab/NCLscan Pipeline to identify intra and intergenic non-colinear transcripts (including circRNA) from paired-end RNA-seq data BWA [88] PTESFinder https://sourceforge.net/projects/ptesfinder-v1/ Pipeline for identifying post-transcriptional exon shuffling events from high-throughput RNA-seq data Bowtie [89] Sailfish-cir https://github.com/zerodel/sailfish-cir Estimates the relative abundance of circRNAs from RNA-seq data RapMap [90] Segemehl http://www.bioinf.uni-leipzig.de/Software/segemehl/ Maps short sequencer reads to reference genomes, detecting mismatches, insertions and deletions Segemehl [91] SUPeR-seq https://www.illumina.com/science/sequencing-method-explorer/kits-and-arrays/super-seq.html Method to sequence both polyadenylated and non-polyadenylated RNAs from individual cells bowtie2 [52] UROBORUS https://github.com/WGLab/UROBORUS Pipeline to identify circRNA from total RNA-seq data without RNase R treatment TopHat Bowtie [92] For example, Circtools (https://github.com/dieterich-lab/circtools), a modular, Python/R-based framework for circRNA-related tools that unifies several essential functions in a one-step software solution, aims to address the issue of circRNA detection and quantification. Circtools' base modules cover detection of circRNAs [67], a beta-binomial model for testing relative abundance changes of circular vs. linear RNA species from the same host gene (CircTest), and to fully reconstruct the exon-intron structure of circRNAs from long sequencing reads (FUCHS) [70]. This framework was used to identify circRNAs and host gene dynamics in cardiac development and disease [55]. 5.2 Alternative splicing analysis Oftentimes, several alternative circRNAs emerge from the same gene locus. Similar to alternative linear splicing, circRNAs could share the same back-splicing site, yet differ in their internal exon composition, that is, by alternative exon skipping. Moreover, alternative back-splicing may occur within the same host gene by using alternative splice site acceptors or donors. Two software packages are specifically designed for resolving this question, namely FUCHS [70] and CIRI-AS [71]. 5.3 Assessment of circRNAs by RT-qPCR Often used to validate RNA-seq data, reverse transcription quantitative polymerase chain reaction (RT-qPCR) has been adapted to specifically assess the expression of circRNAs. Primer design is a necessary first step in validating and quantitating circRNAs. Having considered alternative circRNA splicing, one needs to keep in mind that sequences flanking a particular back-splicing junction may represent one or more circRNAs. A comparison to the expression of the host gene could become a complex endeavour too if multiple isoforms exist. Herein, Circtools supports the experimentalist by its primer design module. 6 Artificial intelligence for circRNA research Artificial intelligence methods, particularly machine learning, are increasingly used to analyse biological data. The large amount of experimental data generated by high-throughput methods makes this both possible and necessary, since machine learning requires a lot of training data, but is also often the most effective approach for finding interesting hidden relationships. So far there have been few applications of artificial intelligence to circRNAs since they have not been researched as intensely as other molecules such as linear RNAs or proteins. However, we expect new utilisations of artificial intelligence methods to arise following the continuous improvement of the knowledge related to the role of circRNAs in biology and disease. Both existing and potential applications of artificial intelligence to circRNA are discussed in this chapter. 6.1 CircRNA recognition The first problem that requires sophisticated computational methods is circRNA recognition as discussed in section 5. This is usually tackled with heuristic methods gathered in Table 1. These tools take RNA-seq data as input and use two main approaches to recognize circRNAs: one that is based on detecting back-splices in RNA-seq data, and another that requires genomic information (explained in Ref. [12]). With increasing numbers of experimentally verified circRNAs, it is becoming feasible to develop machine-learning methods to recognize them. A method called PredcircRNA, as well as the web tool WebCircRNA that implements a version of this method distinguishes circRNAs from other lncRNAs [72,73]. The study describing it started with a dataset of 14,084 circRNAs and 19,722 linear lncRNAs. From these, 188 features were extracted: graph describing RNA structure and sequence, conservation score, nucleotide composition, and other features such as single nucleotide polymorphism density at miRNA binding sites. Each group of features was transformed via a Gaussian kernel, which was fused using lp-norm multiple kernel learning to build a model to recognize circRNAs. The evaluation showed an accuracy of around 80% (depending on the dataset). For a general overview of machine learning from RNA-seq data (of which PredcircRNA is an example), refer to [74]. 6.2 Binding site recognition A possible application of artificial intelligence to circRNAs is to recognize miRNAs or protein binding sites [75]. Here again, to our knowledge, no applications of artificial intelligence to this problem has been reported but computational methods have been used to recognize RNA-binding proteins and binding sites in other contexts [76]. These methods largely rely on machine learning from protein sequence and structure. It is so far still unclear to what degree they apply to circRNA and whether they can be used to recognize binding to circular as opposed to linear RNAs. Unlike most work on RNA binding, which deals with proteins and their binding sites, one framework detects binding sites on RNA. The novelty of this approach is that it uses not only the primary RNA sequence and secondary structure but also tertiary structure [77]. In addition, instead of relying on traditional machine learning algorithms, the framework uses deep learning, an approach that builds deep artificial neural networks and has recently proven superior to traditional machine learning for many problems. It seems likely that such an approach could be used to recognize binding sites on circRNAs as well. 7 Biomarker value of circRNAs Any newly discovered type of RNA molecules has been evaluated for its biomarker potential. Bioinformatics and computational biology approaches based on the concepts of systems biology allowed the discovery of messenger RNAs [78], miRNAs [8,9,79] and lncRNAs [10,11,45] as biomarkers of CV disease. The blood transcriptome has been revealed as a rich reservoir of potential biomarkers of CV diseases [80]. Yet, the biomarker value of circRNAs for CV disease has only recently been suggested. A study in patients with acute myocardial infarction identified the circRNA MICRA which showed an association with the development of heart failure in the few months following the event [28]. This association was validated in two independent patient cohorts. Importantly, MICRA provided a significant incremental value to predict heart failure development over existing risk factors and biomarkers [28,81]. Even if, so far, no other circRNA has been reported to have a diagnostic or a prognostic value for CV disease, it is expected that the resistance of circRNAs to exonuclease degradation confers them with an interesting biomarker potential that needs to be further explored. 8 Considerations for future work The knowledge of the role of circRNAs in CV health and disease is still scarce and more research is needed to reach a better knowledge of their potential as therapeutic targets and biomarkers of CV disease. The following considerations may be pondered when designing future research projects.• Discovery of circRNAs. Different prediction algorithms available for free download from several databases have been developed for the discovery of circRNAs from RNA-seq datasets (Table 1). Use of RNase R-digested total RNA as input for RNA-seq allows targeting circRNAs for increased specificity and sensitivity of detection. • Technical aspects of circRNAs detection. Since circRNAs are produced as an alternative of linear RNA transcripts, it is important to ensure that what is thought of as a circRNA is effectively in a circular form. This verification can be achieved using RNase R digestion and divergent primers or probes spanning the junction site in RT-qPCR experiments. • Biogenesis of circRNAs. A thorough characterization of the mechanisms responsible for circRNAs biogenesis will help to understand their regulation in disease states and may lead to the identification of novel therapeutic targets. • Role of RNA binding proteins in circRNA biogenesis. Since there are hundreds of RNA binding proteins known, an important challenge for the future will be to identify those proteins specifically involved in circRNA biogenesis. • Function of circRNAs. Several databases are available to search for a potential function of circRNAs (Table 2). One of the most interesting, since it links circRNAs to disease states, might be the Circ2Traits database.Table 2 Databases to search for a potential function of circRNAs. Table 2Database name Website Description Reference circBase http://circbase.org/ Merged and unified public datasets of circRNAs [93] circ2Traits http://gyanxet-beta.com/circdb/ Disease-circRNA associations according to a network of predicted interactions between disease-associated miRNAs and protein coding, long non-coding and circRNA genes, as well as the disease-and traits-associated variations in circRNA loci [94] CircNet http://circnet.mbc.nctu.edu.tw/ Tissue-specific circRNA expression profiles and circRNA-miRNA-gene regulatory networks [95] starBase http://starbase.sysu.edu.cn/mirCircRNA.php Predicted miRNA-circRNA interactions overlapped with CLIP-Seq data [96] CLIP: cross-linking immunoprecipitation. • Mechanisms of action. Apart from their presumed miRNA sponge effect, the mechanisms of action of circRNAs are only partly understood and deserve further molecular studies. • Protein-coding potential. Further investigations are required to clarify the potential of circRNAs to encode proteins, as recently suggested in several reports [27,29,82,83]. Some bioinformatics tools have been developed to evaluate the protein-coding potential of circRNAs (Table 3), which may be helpful for future research.Table 3 Bioinformatic tools to evaluate the protein-coding potential of circRNAs. Table 3Tool name Website Aim Reference CircPro http://bis.zju.edu.cn/CircPro/ Predicts the protein-coding potential of circRNAs and discovers junction reads from Ribo-Seq data [97] circRNADB http://202.195.183.4:8000/circrnadb/circRNADb.php A comprehensive database for human circRNAs with protein-coding annotations (IRES, ORF, protein domains) [5] CPAT http://lilab.research.bcm.edu/cpat/ Imputes the probability of an RNA to encode a peptide [98] IRESite http://iresite.org/ Aligns validated IRES to circRNA sequences [99] ORF Finder https://www.ncbi.nlm.nih.gov/orffinder/ Searches for open reading frames in DNA sequences Pfam http://pfam.xfam.org/ Collects protein families with sequence alignments [100] PhyloCSF https://github.com/mlin/PhyloCSF/wiki/ Determines whether a nucleotide sequence represents a protein-coding region [101] CAPT: coding-potential assessment tool; IRES: internal ribosomal entry site; ORF: open reading frame. • Artificial intelligence. As research on circRNAs advance, we envision that new solutions based on artificial intelligence methods will help to identify circRNAs, their interacting molecules, their biological function, and their role in disease. • Biomarker value of circRNAs. The potential of circRNAs to be used as biomarkers is highly related to their relative stability in body fluids due to resistance to exonucleases. However, as any novel biomarker, circRNAs will have to satisfy stringent criteria to be clinically useful: reliable detection, automatable detection, robust predictive value (specificity and sensitivity), incremental value to existing markers, extensive validation in independent and large patient cohorts, and cost-effectiveness. Acknowledgements C. P. C. Gomes is funded by the Eurostars E! 9686 MIPROG project. A. S. S. is funded by the Ministry of Higher Education and Research and the National Research Fund of Luxembourg (grant # C14/BM/8225223). Y. Devaux is supported by the National Research Fund (grant # INTER/EUROSTARS/15/10282117) and the Ministry of Higher Education and Research of Luxembourg. ==== Refs References 1 Lander E.S. Linton L.M. Birren B. Nusbaum C. Zody M.C. Baldwin J. Devon K. Dewar K. Doyle M. FitzHugh W. Funke R. Gage D. Harris K. Heaford A. Howland J. Kann L. Lehoczky J. LeVine R. McEwan P. McKernan K. Meldrim J. Mesirov J.P. Miranda C. Morris W. Naylor J. Raymond C. Rosetti M. Santos R. Sheridan A. Sougnez C. Stange-Thomann N. Stojanovic N. Subramanian A. Wyman D. Rogers J. Sulston J. Ainscough R. Beck S. Bentley D. Burton J. Clee C. Carter N. Coulson A. Deadman R. Deloukas P. Dunham A. Dunham I. Durbin R. French L. Grafham D. Gregory S. Hubbard T. Humphray S. Hunt A. Jones M. Lloyd C. McMurray A. Matthews L. Mercer S. Milne S. Mullikin J.C. Mungall A. Plumb R. Ross M. Shownkeen R. Sims S. Waterston R.H. Wilson R.K. Hillier L.W. McPherson J.D. Marra M.A. Mardis E.R. Fulton L.A. Chinwalla A.T. Pepin K.H. Gish W.R. Chissoe S.L. Wendl M.C. Delehaunty K.D. Miner T.L. Delehaunty A. Kramer J.B. Cook L.L. Fulton R.S. Johnson D.L. Minx P.J. Clifton S.W. Hawkins T. Branscomb E. Predki P. Richardson P. Wenning S. Slezak T. Doggett N. Cheng J.F. Olsen A. Lucas S. Elkin C. Uberbacher E. Frazier M. Gibbs R.A. Muzny D.M. Scherer S.E. Bouck J.B. Sodergren E.J. Worley K.C. Rives C.M. Gorrell J.H. Metzker M.L. Naylor S.L. Kucherlapati R.S. Nelson D.L. Weinstock G.M. Sakaki Y. Fujiyama A. Hattori M. Yada T. Toyoda A. Itoh T. Kawagoe C. Watanabe H. Totoki Y. Taylor T. Weissenbach J. Heilig R. Saurin W. Artiguenave F. Brottier P. Bruls T. Pelletier E. Robert C. Wincker P. Smith D.R. Doucette-Stamm L. Rubenfield M. Weinstock K. Lee H.M. Dubois J. Rosenthal A. Platzer M. Nyakatura G. Taudien S. Rump A. Yang H. Yu J. Wang J. Huang G. Gu J. Hood L. Rowen L. Madan A. Qin S. Davis R.W. Federspiel N.A. Abola A.P. Proctor M.J. Myers R.M. Schmutz J. Dickson M. Grimwood J. Cox D.R. Olson M.V. Kaul R. Raymond C. Shimizu N. Kawasaki K. Minoshima S. Evans G.A. Athanasiou M. Schultz R. Roe B.A. Chen F. Pan H. Ramser J. Lehrach H. Reinhardt R. McCombie W.R. de la Bastide M. Dedhia N. Blocker H. Hornischer K. Nordsiek G. Agarwala R. Aravind L. Bailey J.A. Bateman A. Batzoglou S. Birney E. Bork P. Brown D.G. Burge C.B. Cerutti L. Chen H.C. Church D. Clamp M. Copley R.R. Doerks T. Eddy S.R. Eichler E.E. Furey T.S. Galagan J. Gilbert J.G. Harmon C. Hayashizaki Y. Haussler D. Hermjakob H. Hokamp K. Jang W. Johnson L.S. Jones T.A. Kasif S. Kaspryzk A. Kennedy S. Kent W.J. Kitts P. Koonin E.V. Korf I. Kulp D. Lancet D. Lowe T.M. McLysaght A. Mikkelsen T. Moran J.V. Mulder N. Pollara V.J. Ponting C.P. Schuler G. Schultz J. Slater G. Smit A.F. Stupka E. Szustakowski J. Thierry-Mieg D. Thierry-Mieg J. Wagner L. Wallis J. Wheeler R. Williams A. Wolf Y.I. Wolfe K.H. Yang S.P. Yeh R.F. Collins F. Guyer M.S. Peterson J. Felsenfeld A. Wetterstrand K.A. Patrinos A. Morgan M.J. de Jong P. Catanese J.J. Osoegawa K. Shizuya H. Choi S. Chen Y.J. Initial sequencing and analysis of the human genome Nature 409 6822 2001 860 921 11237011 2 Venter J.C. Adams M.D. Myers E.W. Li P.W. Mural R.J. Sutton G.G. Smith H.O. Yandell M. Evans C.A. Holt R.A. Gocayne J.D. Amanatides P. Ballew R.M. Huson D.H. Wortman J.R. Zhang Q. Kodira C.D. Zheng X.H. Chen L. Skupski M. Subramanian G. Thomas P.D. Zhang J. Gabor Miklos G.L. Nelson C. Broder S. Clark A.G. Nadeau J. McKusick V.A. Zinder N. Levine A.J. Roberts R.J. Simon M. Slayman C. Hunkapiller M. Bolanos R. Delcher A. Dew I. Fasulo D. Flanigan M. Florea L. Halpern A. Hannenhalli S. Kravitz S. Levy S. Mobarry C. Reinert K. Remington K. Abu-Threideh J. Beasley E. Biddick K. Bonazzi V. Brandon R. Cargill M. Chandramouliswaran I. Charlab R. Chaturvedi K. Deng Z. Di Francesco V. Dunn P. Eilbeck K. Evangelista C. Gabrielian A.E. Gan W. Ge W. Gong F. Gu Z. Guan P. Heiman T.J. Higgins M.E. Ji R.R. Ke Z. Ketchum K.A. Lai Z. Lei Y. Li Z. Li J. Liang Y. Lin X. Lu F. Merkulov G.V. Milshina N. Moore H.M. Naik A.K. Narayan V.A. Neelam B. Nusskern D. Rusch D.B. Salzberg S. Shao W. Shue B. Sun J. Wang Z. Wang A. Wang X. Wang J. Wei M. Wides R. Xiao C. Yan C. Yao A. Ye J. Zhan M. Zhang W. Zhang H. Zhao Q. Zheng L. Zhong F. Zhong W. Zhu S. Zhao S. Gilbert D. Baumhueter S. Spier G. Carter C. Cravchik A. Woodage T. Ali F. An H. Awe A. Baldwin D. Baden H. Barnstead M. Barrow I. Beeson K. Busam D. Carver A. Center A. Cheng M.L. Curry L. Danaher S. Davenport L. Desilets R. Dietz S. Dodson K. Doup L. Ferriera S. Garg N. Gluecksmann A. Hart B. Haynes J. Haynes C. Heiner C. Hladun S. Hostin D. Houck J. Howland T. Ibegwam C. Johnson J. Kalush F. Kline L. Koduru S. Love A. Mann F. May D. McCawley S. McIntosh T. McMullen I. Moy M. Moy L. Murphy B. Nelson K. Pfannkoch C. Pratts E. Puri V. Qureshi H. Reardon M. Rodriguez R. Rogers Y.H. Romblad D. Ruhfel B. Scott R. Sitter C. Smallwood M. Stewart E. Strong R. Suh E. Thomas R. Tint N.N. Tse S. Vech C. Wang G. Wetter J. Williams S. Williams M. Windsor S. Winn-Deen E. Wolfe K. Zaveri J. Zaveri K. Abril J.F. Guigo R. Campbell M.J. Sjolander K.V. Karlak B. Kejariwal A. Mi H. Lazareva B. Hatton T. Narechania A. Diemer K. Muruganujan A. Guo N. Sato S. Bafna V. Istrail S. Lippert R. Schwartz R. Walenz B. Yooseph S. Allen D. Basu A. Baxendale J. Blick L. Caminha M. Carnes-Stine J. Caulk P. Chiang Y.H. Coyne M. Dahlke C. Mays A. Dombroski M. Donnelly M. Ely D. Esparham S. Fosler C. Gire H. Glanowski S. Glasser K. Glodek A. Gorokhov M. Graham K. Gropman B. Harris M. Heil J. Henderson S. Hoover J. Jennings D. Jordan C. Jordan J. Kasha J. Kagan L. Kraft C. Levitsky A. Lewis M. Liu X. Lopez J. Ma D. Majoros W. McDaniel J. Murphy S. Newman M. Nguyen T. Nguyen N. Nodell M. Pan S. Peck J. Peterson M. Rowe W. Sanders R. Scott J. Simpson M. Smith T. Sprague A. Stockwell T. Turner R. Venter E. Wang M. Wen M. Wu D. Wu M. Xia A. Zandieh A. Zhu X. The sequence of the human genome Science 291 5507 2001 1304 1351 11181995 3 E.P. Consortium An integrated encyclopedia of DNA elements in the human genome Nature 489 7414 2012 57 74 22955616 4 Barrett S.P. Salzman J. Circular RNAs: analysis, expression and potential functions Development 143 11 2016 1838 1847 27246710 5 Chen X. Han P. Zhou T. Guo X. Song X. Li Y. circRNADb: a comprehensive database for human circular RNAs with protein-coding annotations Sci. Rep. 6 2016 34985 27725737 6 Benjamin E.J. Blaha M.J. Chiuve S.E. Cushman M. Das S.R. Deo R. de Ferranti S.D. Floyd J. Fornage M. Gillespie C. Isasi C.R. Jimenez M.C. Jordan L.C. Judd S.E. Lackland D. Lichtman J.H. Lisabeth L. Liu S. Longenecker C.T. Mackey R.H. Matsushita K. Mozaffarian D. Mussolino M.E. Nasir K. Neumar R.W. Palaniappan L. Pandey D.K. Thiagarajan R.R. Reeves M.J. Ritchey M. Rodriguez C.J. Roth G.A. Rosamond W.D. Sasson C. Towfighi A. Tsao C.W. Turner M.B. Virani S.S. Voeks J.H. Willey J.Z. Wilkins J.T. Wu J.H. Alger H.M. Wong S.S. Muntner P. C. American Heart Association Statistics, S. Stroke Statistics, Heart Disease and Stroke Statistics-2017 Update A report from the american heart association Circulation 135 10 2017 e146 e603 28122885 7 Townsend N. Wilson L. Bhatnagar P. Wickramasinghe K. Rayner M. Nichols M. Cardiovascular disease in Europe: epidemiological update 2016 Eur. Heart J. 37 42 2016 3232 3245 27523477 8 Beermann J. Piccoli M.T. Viereck J. Thum T. Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches Physiol. Rev. 96 4 2016 1297 1325 27535639 9 Goretti E. Wagner D.R. Devaux Y. miRNAs as biomarkers of myocardial infarction: a step forward towards personalized medicine? Trends Mol. Med. 20 12 2014 716 725 25457620 10 Gomes C.P.C. Spencer H. Ford K.L. Michel L.Y.M. Baker A.H. Emanueli C. Balligand J.L. Devaux Y. Cardiolinc n. The function and therapeutic potential of long non-coding RNAs in cardiovascular development and disease Mol. Ther. Nucleic Acids 8 2017 494 507 28918050 11 Greco S. Salgado Somoza A. Devaux Y. Martelli F. Long noncoding RNAs and cardiac disease Antioxidants Redox Signal. August 2017 (ahead of print). 10.1089/ars.2017.7126 12 Jeck W.R. Sharpless N.E. Detecting and characterizing circular RNAs Nat. Biotechnol. 32 5 2014 453 461 24811520 13 Petkovic S. Muller S. RNA circularization strategies in vivo and in vitro Nucleic Acids Res. 43 4 2015 2454 2465 25662225 14 Ashwal-Fluss R. Meyer M. Pamudurti N.R. Ivanov A. Bartok O. Hanan M. Evantal N. Memczak S. Rajewsky N. Kadener S. circRNA biogenesis competes with pre-mRNA splicing Mol. Cell 56 1 2014 55 66 25242144 15 Starke S. Jost I. Rossbach O. Schneider T. Schreiner S. Hung L.H. Bindereif A. Exon circularization requires canonical splice signals Cell Rep. 10 1 2015 103 111 25543144 16 Zhang X.O. Wang H.B. Zhang Y. Lu X. Chen L.L. Yang L. Complementary sequence-mediated exon circularization Cell 159 1 2014 134 147 25242744 17 Liang D. Wilusz J.E. Short intronic repeat sequences facilitate circular RNA production Gene Dev. 28 20 2014 2233 2247 25281217 18 Castello A. Fischer B. Eichelbaum K. Horos R. Beckmann B.M. Strein C. Davey N.E. Humphreys D.T. Preiss T. Steinmetz L.M. Krijgsveld J. Hentze M.W. Insights into RNA biology from an atlas of mammalian mRNA-binding proteins Cell 149 6 2012 1393 1406 22658674 19 Konieczny P. Stepniak-Konieczna E. Sobczak K. MBNL proteins and their target RNAs, interaction and splicing regulation Nucleic Acids Res. 42 17 2014 10873 10887 25183524 20 Conn S.J. Pillman K.A. Toubia J. Conn V.M. Salmanidis M. Phillips C.A. Roslan S. Schreiber A.W. Gregory P.A. Goodall G.J. The RNA binding protein quaking regulates formation of circRNAs Cell 160 6 2015 1125 1134 25768908 21 Chenard C.A. Richard S. New implications for the QUAKING RNA binding protein in human disease J. Neurosci. Res. 86 2 2008 233 242 17787018 22 Gott J.M. Emeson R.B. Functions and mechanisms of RNA editing Annu. Rev. Genet. 34 2000 499 531 11092837 23 Nishikura K. Functions and regulation of RNA editing by ADAR deaminases Annu. Rev. Biochem. 79 2010 321 349 20192758 24 Ivanov A. Memczak S. Wyler E. Torti F. Porath H.T. Orejuela M.R. Piechotta M. Levanon E.Y. Landthaler M. Dieterich C. Rajewsky N. Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals Cell Rep. 10 2 2015 170 177 25558066 25 Rybak-Wolf A. Stottmeister C. Glazar P. Jens M. Pino N. Giusti S. Hanan M. Behm M. Bartok O. Ashwal-Fluss R. Herzog M. Schreyer L. Papavasileiou P. Ivanov A. Ohman M. Refojo D. Kadener S. Rajewsky N. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed Mol. Cell 58 5 2015 870 885 25921068 26 Khan M.A. Reckman Y.J. Aufiero S. van den Hoogenhof M.M. van der Made I. Beqqali A. Koolbergen D.R. Rasmussen T.B. van der Velden J. Creemers E.E. Pinto Y.M. RBM20 regulates circular RNA production from the titin gene Circ. Res. 119 9 2016 996 1003 27531932 27 Pamudurti N.R. Bartok O. Jens M. Ashwal-Fluss R. Stottmeister C. Ruhe L. Hanan M. Wyler E. Perez-Hernandez D. Ramberger E. Shenzis S. Samson M. Dittmar G. Landthaler M. Chekulaeva M. Rajewsky N. Kadener S. Translation of CircRNAs Mol. Cell 66 1 2017 9 21 e7 28344080 28 Vausort M. Salgado-Somoza A. Zhang L. Leszek P. Scholz M. Teren A. Burkhardt R. Thiery J. Wagner D.R. Devaux Y. Myocardial infarction-associated circular RNA predicting left ventricular dysfunction J. Am. Coll. Cardiol. 68 11 2016 1247 1248 27609688 29 Legnini I. Di Timoteo G. Rossi F. Morlando M. Briganti F. Sthandier O. Fatica A. Santini T. Andronache A. Wade M. Laneve P. Rajewsky N. Bozzoni I. Circ-ZNF609 is a circular RNA that can Be translated and functions in myogenesis Mol. Cell 66 1 2017 22 37 e9 28344082 30 Memczak S. Jens M. Elefsinioti A. Torti F. Krueger J. Rybak A. Maier L. Mackowiak S.D. Gregersen L.H. Munschauer M. Loewer A. Ziebold U. Landthaler M. Kocks C. le Noble F. Rajewsky N. Circular RNAs are a large class of animal RNAs with regulatory potency Nature 495 7441 2013 333 338 23446348 31 Hansen T.B. Jensen T.I. Clausen B.H. Bramsen J.B. Finsen B. Damgaard C.K. Kjems J. Natural RNA circles function as efficient microRNA sponges Nature 495 7441 2013 384 388 23446346 32 Xu H. Guo S. Li W. Yu P. The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells Sci. Rep. 5 2015 12453 26211738 33 Geng H.H. Li R. Su Y.M. Xiao J. Pan M. Cai X.X. Ji X.P. The circular RNA Cdr1as promotes myocardial infarction by mediating the regulation of miR-7a on its target genes expression PLoS One 11 3 2016 e0151753 26998750 34 Wang K. Long B. Liu F. Wang J.X. Liu C.Y. Zhao B. Zhou L.Y. Sun T. Wang M. Yu T. Gong Y. Liu J. Dong Y.H. Li N. Li P.F. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223 Eur. Heart J. 37 33 2016 2602 2611 26802132 35 Bazan H.A. Hatfield S.A. Brug A. Brooks A.J. Lightell D.J. Jr. Woods T.C. Carotid plaque rupture is accompanied by an increase in the ratio of serum circR-284 to miR-221 levels Circ. Cardiovasc. Genet. 10 4 2017 36 Zou M. Huang C. Li X. He X. Chen Y. Liao W. Liao Y. Sun J. Liu Z. Zhong L. Bin J. Circular RNA expression profile and potential function of hsa_circRNA_101238 in human thoracic aortic dissection Oncotarget 8 47 2017 37 Sun Y. Yang Z. Zheng B. Zhang X.H. Zhang M.L. Zhao X.S. Zhao H.Y. Suzuki T. Wen J.K. A novel regulatory mechanism of smooth muscle alpha-actin expression by NRG-1/circACTA2/miR-548f-5p axis Circ. Res. 121 6 2017 628 635 28698179 38 Tang W.H. Hazen S.L. The gut microbiome and its role in cardiovascular diseases Circulation 135 11 2017 1008 1010 28289004 39 Wang K. Gan T.Y. Li N. Liu C.Y. Zhou L.Y. Gao J.N. Chen C. Yan K.W. Ponnusamy M. Zhang Y.H. Li P.F. Circular RNA mediates cardiomyocyte death via miRNA-dependent upregulation of MTP18 expression Cell Death Differ. 24 6 2017 1111 1120 28498369 40 Liu C. Yao M.D. Li C.P. Shan K. Yang H. Wang J.J. Liu B. Li X.M. Yao J. Jiang Q. Yan B. Silencing of circular RNA-ZNF609 ameliorates vascular endothelial dysfunction Theranostics 7 11 2017 2863 2877 28824721 41 Du W.W. Yang W. Chen Y. Wu Z.K. Foster F.S. Yang Z. Li X. Yang B.B. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses Eur. Heart J. 38 18 2017 1402 1412 26873092 42 Du W.W. Yang W. Liu E. Yang Z. Dhaliwal P. Yang B.B. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2 Nucleic Acids Res. 44 6 2016 2846 2858 26861625 43 Boeckel J.N. Jae N. Heumuller A.W. Chen W. Boon R.A. Stellos K. Zeiher A.M. John D. Uchida S. Dimmeler S. Identification and characterization of hypoxia-regulated endothelial circular RNA Circ. Res. 117 10 2015 884 890 26377962 44 Doetschman T. Georgieva T. Gene editing with CRISPR/cas9 RNA-directed nuclease Circ. Res. 120 5 2017 876 894 28254804 45 Devaux Y. Zangrando J. Schroen B. Creemers E.E. Pedrazzini T. Chang C.P. Dorn G.W. 2nd Thum T. Heymans S. Cardiolinc n. Long noncoding RNAs in cardiac development and ageing Nat. Rev. Cardiol. 12 7 2015 415 425 25855606 46 Sanger H.L. Klotz G. Riesner D. Gross H.J. Kleinschmidt A.K. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures Proc. Natl. Acad. Sci. U.S.A. 73 11 1976 3852 3856 1069269 47 Cocquerelle C. Mascrez B. Hetuin D. Bailleul B. Mis-splicing yields circular RNA molecules Faseb. J. : Official Publ. Fed. Fed. Am. Soc. Exp. Biol. 7 1 1993 155 160 48 Nigro J.M. Cho K.R. Fearon E.R. Kern S.E. Ruppert J.M. Oliner J.D. Kinzler K.W. Vogelstein B. Scrambled exons Cell 64 3 1991 607 613 1991322 49 Hansen T.B. Wiklund E.D. Bramsen J.B. Villadsen S.B. Statham A.L. Clark S.J. Kjems J. miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA EMBO J. 30 21 2011 4414 4422 21964070 50 Salzman J. Gawad C. Wang P.L. Lacayo N. Brown P.O. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types PLoS One 7 2 2012 e30733 22319583 51 Jeck W.R. Sorrentino J.A. Wang K. Slevin M.K. Burd C.E. Liu J. Marzluff W.F. Sharpless N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats Rna 19 2 2013 141 157 23249747 52 Fan X. Zhang X. Wu X. Guo H. Hu Y. Tang F. Huang Y. Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos Genome Biol. 16 2015 148 26201400 53 Salzman J. Chen R.E. Olsen M.N. Wang P.L. Brown P.O. Cell-type specific features of circular RNA expression PLoS Genet. 9 9 2013 e1003777 24039610 54 Greene J. Baird A.M. Brady L. Lim M. Gray S.G. McDermott R. Finn S.P. Circular RNAs: biogenesis, function and role in human diseases Front Mol Biosci 4 2017 38 28634583 55 Siede D. Rapti K. Gorska A.A. Katus H.A. Altmuller J. Boeckel J.N. Meder B. Maack C. Volkers M. Muller O.J. Backs J. Dieterich C. Identification of circular RNAs with host gene-independent expression in human model systems for cardiac differentiation and disease J. Mol. Cell. Cardiol. 109 2017 48 56 28676412 56 Szabo L. Morey R. Palpant N.J. Wang P.L. Afari N. Jiang C. Parast M.M. Murry C.E. Laurent L.C. Salzman J. Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development Genome Biol. 16 2015 126 26076956 57 Werfel S. Nothjunge S. Schwarzmayr T. Strom T.M. Meitinger T. Engelhardt S. Characterization of circular RNAs in human, mouse and rat hearts J. Mol. Cell. Cardiol. 98 2016 103 107 27476877 58 Tan W.L. Lim B.T. Anene-Nzelu C.G. Ackers-Johnson M. Dashi A. See K. Tiang Z. Lee D.P. Chua W.W. Luu T.D. Li P.Y. Richards A.M. Foo R.S. A landscape of circular RNA expression in the human heart Cardiovasc. Res. 113 3 2017 298 309 28082450 59 Guo W. Schafer S. Greaser M.L. Radke M.H. Liss M. Govindarajan T. Maatz H. Schulz H. Li S. Parrish A.M. Dauksaite V. Vakeel P. Klaassen S. Gerull B. Thierfelder L. Regitz-Zagrosek V. Hacker T.A. Saupe K.W. Dec G.W. Ellinor P.T. MacRae C.A. Spallek B. Fischer R. Perrot A. Ozcelik C. Saar K. Hubner N. Gotthardt M. RBM20, a gene for hereditary cardiomyopathy, regulates titin splicing Nat. Med. 18 5 2012 766 773 22466703 60 Haas J. Frese K.S. Peil B. Kloos W. Keller A. Nietsch R. Feng Z. Muller S. Kayvanpour E. Vogel B. Sedaghat-Hamedani F. Lim W.K. Zhao X. Fradkin D. Kohler D. Fischer S. Franke J. Marquart S. Barb I. Li D.T. Amr A. Ehlermann P. Mereles D. Weis T. Hassel S. Kremer A. King V. Wirsz E. Isnard R. Komajda M. Serio A. Grasso M. Syrris P. Wicks E. Plagnol V. Lopes L. Gadgaard T. Eiskjaer H. Jorgensen M. Garcia-Giustiniani D. Ortiz-Genga M. Crespo-Leiro M.G. Deprez R.H. Christiaans I. van Rijsingen I.A. Wilde A.A. Waldenstrom A. Bolognesi M. Bellazzi R. Morner S. Bermejo J.L. Monserrat L. Villard E. Mogensen J. Pinto Y.M. Charron P. Elliott P. Arbustini E. Katus H.A. Meder B. Atlas of the clinical genetics of human dilated cardiomyopathy Eur. Heart J. 36 18 2015 1123–35a 61 Maatz H. Jens M. Liss M. Schafer S. Heinig M. Kirchner M. Adami E. Rintisch C. Dauksaite V. Radke M.H. Selbach M. Barton P.J. Cook S.A. Rajewsky N. Gotthardt M. Landthaler M. Hubner N. RNA-binding protein RBM20 represses splicing to orchestrate cardiac pre-mRNA processing J. Clin. Invest. 124 8 2014 3419 3430 24960161 62 Kelly S. Greenman C. Cook P.R. Papantonis A. Exon skipping is correlated with exon circularization J. Mol. Biol. 427 15 2015 2414 2417 25728652 63 Wu H.J. Zhang C.Y. Zhang S. Chang M. Wang H.Y. Microarray expression profile of circular RNAs in heart tissue of mice with myocardial infarction-induced heart failure Cell. Physiol. Biochem. Int. J. Exp. Cell. Physiol. Biochem. Pharmacol. 39 1 2016 205 216 64 Deng Y.-Y. Zhang W. She J. Zhang L. Chen T. Zhou J. Yuan Z. GW27-e1167 Circular RNA Related to PPARg Function as ceRNA of microRNA in human acute myocardial infarction J. Am. Coll. Cardiol. 68 16 Supplement 2016 C51 C52 65 Tang C.M. Zhang M. Huang L. Hu Z.Q. Zhu J.N. Xiao Z. Zhang Z. Lin Q.X. Zheng X.L. Yang M. Wu S.L. Cheng J.D. Shan Z.X. CircRNA_000203 enhances the expression of fibrosis-associated genes by derepressing targets of miR-26b-5p, Col1a2 and CTGF, in cardiac fibroblasts Sci. Rep. 7 2017 40342 28079129 66 Zhou B. Yu J.W. A novel identified circular RNA, circRNA_010567, promotes myocardial fibrosis via suppressing miR-141 by targeting TGF-beta1 Biochem. Biophys. Res. Commun. 487 4 2017 769 775 28412345 67 Cheng J. Metge F. Dieterich C. Specific identification and quantification of circular RNAs from sequencing data Bioinformatics 32 7 2016 1094 1096 26556385 68 Panda A.C. De S. Grammatikakis I. Munk R. Yang X. Piao Y. Dudekula D.B. Abdelmohsen K. Gorospe M. High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs Nucleic Acids Res. 45 12 2017 e116 28444238 69 Zeng X. Lin W. Guo M. Zou Q. A comprehensive overview and evaluation of circular RNA detection tools PLoS Comput. Biol. 13 6 2017 e1005420 28594838 70 Metge F. Czaja-Hasse L.F. Reinhardt R. Dieterich C. FUCHS-towards full circular RNA characterization using RNAseq PeerJ 5 2017 e2934 28265491 71 Gao Y. Wang J. Zheng Y. Zhang J. Chen S. Zhao F. Comprehensive identification of internal structure and alternative splicing events in circular RNAs Nat. Commun. 7 2016 12060 27350239 72 Pan X. Xiong K. PredcircRNA: computational classification of circular RNA from other long non-coding RNA using hybrid features Mol. Biosyst. 11 8 2015 2219 2226 26028480 73 C.f.n.-c.R.i.T.a. Health, WebCircRNA: Assessing the Circular RNA Potential of Coding and Noncoding RNA. 74 Jabeen A. Ahmad N. Raza K. Machine learning-based state-of-the-art methods for the classification of RNA-seq data Dey N. Ashour A.S. Borra S. Classification in BioApps 2018 Springer International Publishing 75 Dudekula D.B. Panda A.C. Grammatikakis I. De S. Abdelmohsen K. Gorospe M. CircInteractome: a web tool for exploring circular RNAs and their interacting proteins and microRNAs RNA Biol. 13 1 2016 34 42 26669964 76 Si J. Cui J. Cheng J. Wu R. Computational prediction of RNA-binding proteins and binding sites Int. J. Mol. Sci. 16 11 2015 25952 77 Zhang S. Zhou J. Hu H. Gong H. Chen L. Cheng C. Zeng J. A deep learning framework for modeling structural features of RNA-binding protein targets Nucleic Acids Res. 44 4 2016 e32 78 Azuaje F.J. Dewey F.E. Brutsaert D.L. Devaux Y. Ashley E.A. Wagner D.R. Systems-based approaches to cardiovascular biomarker discovery Circ. Cardiovasc. Genet. 5 3 2012 360 367 22715280 79 Devaux Y. Vausort M. McCann G.P. Zangrando J. Kelly D. Razvi N. Zhang L. Ng L.L. Wagner D.R. Squire I.B. MicroRNA-150: a novel marker of left ventricular remodeling after acute myocardial infarction, circulation Cardiovasc. Genet. 6 3 2013 290 298 80 Devaux Y. Creemers E.E. Boon R.A. Werfel S. Thum T. Engelhardt S. Dimmeler S. Squire I. n. Cardiolinc, Circular RNAs in heart failure Eur. J. Heart Fail. 19 6 2017 701 709 28345158 81 Salgado-Somoza A. Zhang L. Vausort M. Devaux Y. The circular RNA MICRA for risk stratification after myocardial infarction Int. J. Cardiol. Heart Vasc. 17 2017 33 36 29159270 82 Abe N. Matsumoto K. Nishihara M. Nakano Y. Shibata A. Maruyama H. Shuto S. Matsuda A. Yoshida M. Ito Y. Abe H. Rolling circle translation of circular RNA in living human cells Sci. Rep. 5 2015 16435 26553571 83 Yang Y. Fan X. Mao M. Song X. Wu P. Zhang Y. Jin Y. Yang Y. Chen L.L. Wang Y. Wong C.C. Xiao X. Wang Z. Extensive translation of circular RNAs driven by N6-methyladenosine Cell Res. 27 5 2017 626 641 28281539 84 You X. Conrad T.O. Acfs: accurate circRNA identification and quantification from RNA-Seq data Sci. Rep. 6 2016 38820 27929140 85 Westholm J.O. Miura P. Olson S. Shenker S. Joseph B. Sanfilippo P. Celniker S.E. Graveley B.R. Lai E.C. Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation Cell Rep. 9 5 2014 1966 1980 25544350 86 Gao Y. Zhang J. Zhao F. Circular RNA identification based on multiple seed matching Briefings Bioinf. bbx014 2017 87 Wang K. Singh D. Zeng Z. Coleman S.J. Huang Y. Savich G.L. He X. Mieczkowski P. Grimm S.A. Perou C.M. MacLeod J.N. Chiang D.Y. Prins J.F. Liu J. MapSplice: accurate mapping of RNA-seq reads for splice junction discovery Nucleic Acids Res. 38 18 2010 e178 88 Chuang T.J. Wu C.S. Chen C.Y. Hung L.Y. Chiang T.W. Yang M.Y. NCLscan: accurate identification of non-co-linear transcripts (fusion, trans-splicing and circular RNA) with a good balance between sensitivity and precision Nucleic Acids Res. 44 3 2016 e29 89 Izuogu O.G. Alhasan A.A. Alafghani H.M. Santibanez-Koref M. Elliott D.J. Jackson M.S. PTESFinder: a computational method to identify post-transcriptional exon shuffling (PTES) events BMC Bioinf. 17 2016 31 90 Li M. Xie X. Zhou J. Sheng M. Yin X. Ko E.A. Zhou T. Gu W. Quantifying circular RNA expression from RNA-seq data using model-based framework Bioinformatics 33 14 2017 2131 2139 28334396 91 Hoffmann S. Otto C. Doose G. Tanzer A. Langenberger D. Christ S. Kunz M. Holdt L.M. Teupser D. Hackermuller J. Stadler P.F. A multi-split mapping algorithm for circular RNA, splicing, trans-splicing and fusion detection Genome Biol. 15 2 2014 R34 92 Song X. Zhang N. Han P. Moon B.S. Lai R.K. Wang K. Lu W. Circular RNA profile in gliomas revealed by identification tool UROBORUS Nucleic Acids Res. 44 9 2016 e87 93 Glazar P. Papavasileiou P. Rajewsky N. circBase: a database for circular RNAs Rna 20 11 2014 1666 1670 25234927 94 Ghosal S. Das S. Sen R. Basak P. Chakrabarti J. Circ2Traits: a comprehensive database for circular RNA potentially associated with disease and traits Front. Genet. 4 2013 283 24339831 95 Liu Y.C. Li J.R. Sun C.H. Andrews E. Chao R.F. Lin F.M. Weng S.L. Hsu S.D. Huang C.C. Cheng C. Liu C.C. Huang H.D. CircNet: a database of circular RNAs derived from transcriptome sequencing data Nucleic Acids Res. 44 D1 2016 D209 D215 26450965 96 Li J.H. Liu S. Zhou H. Qu L.H. Yang J.H. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data Nucleic Acids Res. 42 Database issue 2014 D92 D97 24297251 97 Meng X. Chen Q. Zhang P. Chen M. CircPro: an integrated tool for the identification of circRNAs with protein-coding potential Bioinformatics 33 20 2017 3314 3316 29028266 98 Wang L. Park H.J. Dasari S. Wang S. Kocher J.P. Li W. CPAT: coding-Potential Assessment Tool using an alignment-free logistic regression model Nucleic Acids Res. 41 6 2013 e74 99 Mokrejs M. Masek T. Vopalensky V. Hlubucek P. Delbos P. Pospisek M. IRESite–a tool for the examination of viral and cellular internal ribosome entry sites Nucleic Acids Res. 38 Database issue 2010 D131 D136 19917642 100 Finn R.D. Coggill P. Eberhardt R.Y. Eddy S.R. Mistry J. Mitchell A.L. Potter S.C. Punta M. Qureshi M. Sangrador-Vegas A. Salazar G.A. Tate J. Bateman A. The Pfam protein families database: towards a more sustainable future Nucleic Acids Res. 44 D1 2016 D279 D285 26673716 101 Lin M.F. Jungreis I. Kellis M. PhyloCSF: a comparative genomics method to distinguish protein coding and non-coding regions Bioinformatics 27 13 2011 i275 i282 21685081