Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Author Interviews
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Journal of Human Immunity
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Journal of Human Immunity
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Author Interviews
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow LSA on Bluesky
  • Follow lsa Template on Twitter
Resource
Transparent Process
Open Access

Sequence and expression levels of circular RNAs in progenitor cell types during mouse corticogenesis

View ORCID ProfileMartina Dori, Leila Haj Abdullah Alieh, View ORCID ProfileDaniel Cavalli, Simone Massalini, View ORCID ProfileMathias Lesche, View ORCID ProfileAndreas Dahl, View ORCID ProfileFederico Calegari  Correspondence email
Martina Dori
1CRTD-Center for Regenerative Therapies Dresden, School of Medicine, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Martina Dori
Leila Haj Abdullah Alieh
1CRTD-Center for Regenerative Therapies Dresden, School of Medicine, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel Cavalli
1CRTD-Center for Regenerative Therapies Dresden, School of Medicine, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel Cavalli
Simone Massalini
1CRTD-Center for Regenerative Therapies Dresden, School of Medicine, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mathias Lesche
2DRESDEN-concept Genome Center c/o Center for Molecular and Cellular Bioengineering, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mathias Lesche
Andreas Dahl
2DRESDEN-concept Genome Center c/o Center for Molecular and Cellular Bioengineering, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Andreas Dahl
Federico Calegari
1CRTD-Center for Regenerative Therapies Dresden, School of Medicine, Technische Universität Dresden, Dresden, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Federico Calegari
  • For correspondence: federico.calegari@tu-dresden.de
Published 29 March 2019. DOI: 10.26508/lsa.201900354
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

Circular (circ) RNAs have recently emerged as a novel class of transcripts whose identification and function remain elusive. Among many tissues and species, the mammalian brain is the organ in which circRNAs are more abundant and first evidence of their functional significance started to emerge. Yet, even within this well-studied organ, annotation of circRNAs remains fragmentary, their sequence is unknown, and their expression in specific cell types was never investigated. Overcoming these limitations, here we provide the first comprehensive identification of circRNAs and assessment of their expression patterns in proliferating neural stem cells, neurogenic progenitors, and newborn neurons of the developing mouse cortex. Extending the current knowledge about the diversity of this class of transcripts by the identification of nearly 4,000 new circRNAs, our study is the first to provide the full sequence information and expression patterns of circRNAs in cell types representing the lineage of neurogenic commitment. We further exploited our data by evaluating the coding potential, evolutionary conservation, and biogenesis of circRNAs that we found to arise from a specific subclass of linear mRNAs. Our study provides the arising field of circRNA biology with a powerful new resource to address the complexity and potential biological significance of this new class of transcripts.

Introduction

In the last few decades, the field of RNA biology has witnessed impressive developments. Fuelled by new sequencing technologies, these included the comprehensive annotation of micro- and long noncoding (lnc) RNAs in various organisms and tissues, the characterization of RNA modifications and the new field of epitranscriptomic and the discovery of an entirely new class of noncoding RNAs: circular (circ) RNAs (Kosik, 2013).

CircRNAs are transcripts whose 3′ and 5′ ends are covalently linked in a nonlinear manner resulting in a so-called backsplice junction (Lasda & Parker, 2014; Vicens & Westhof, 2014). The lack of a 3′ poly(A) tail and 5′ capping provides this class of RNAs resistance to exonuclease activity and, thus, an average longer half-live as compared with linear RNAs (Suzuki et al, 2006; Vincent & Deutscher, 2006; Jeck et al, 2013). Transcripts with these characteristics have long been known, but until recently, circRNAs were primarily found in viruses (Kos et al, 1986), and although some reports indicated their origin also from eukaryotic genomes (Nigro et al, 1991; Capel et al, 1993; Zaphiropoulos, 1996), these were still considered a rarity or a byproduct of splicing with no specific function. This view was completely changed very recently after the identification of thousands of circRNAs, including some with regulatory functions during brain development (Salzman et al, 2012; Hansen et al, 2013; Jeck et al, 2013; Memczak et al, 2013; Piwecka et al, 2017).

Despite their abundance, predicting circRNAs remains burdensome and typically relies on bioinformatic tools identifying sequences across backsplice junctions from RNA sequencing data obtained upon depletion of ribosomal RNA (Szabo & Salzman, 2016). Although this has resulted in the prediction of thousands of potential circRNAs in cell lines or whole organs of many species (Salzman et al, 2012; Jeck et al, 2013; Memczak et al, 2013; Westholm et al, 2014; Filippenkov et al, 2015; Jakobi et al, 2016; Kristensen et al, 2018), this approach based on ribosomal RNA depletion has the critical limitation that reads that do not map on a backsplice junction cannot be assigned to a circular, as opposed to a linear, transcript, resulting in the exclusion of the overwhelming majority of the sequencing data. In turn, this makes it impossible to reliably reconstruct neither the full sequence nor the expression level of large pools of circRNAs. Overcoming these limitations, the use of the exonuclease RNase R triggers the digestion of linear RNAs, thereby allowing the isolation of circRNAs (Suzuki et al, 2006). To date, this strategy was applied to a few cell lines or whole organs (Jeck et al, 2013; Ashwal-Fluss et al, 2014; Bachmayr-Heyda et al, 2015; Jakobi et al, 2016; Yang et al, 2017), but the expression patterns of circRNAs in specific cell types in physiological conditions is not known, and a comprehensive reconstruction of their sequence is yet to be reported.

In addition to a poor classification, barely a handful of circRNAs have been suggested to be functionally relevant. For example, and beside their use as biomarkers in various diseases from cancer to diabetes (Bahn et al, 2015; Memczak et al, 2015; Abu & Jamal, 2016; Kulcheski et al, 2016), a study concluded that at least some circRNA, such as circ-ZNF609, may retain coding potential (Legnini et al, 2017). In addition, Drosophila’s circMbl was found to interact with MBL protein resulting from the linear form of the same transcript (Ashwal-Fluss et al, 2014). Finally, and perhaps as the only circRNA-mediated molecular mechanism underlying a specific cellular effect, Cdr1as/ciRS-7 was shown to act as a sponge for miR-7 during development (Hansen et al, 2013; Memczak et al, 2013), and its depletion altered synaptic transmission in adulthood (Piwecka et al, 2017).

From these and other studies, it emerged that the mammalian brain is the organ most enriched in circRNAs (Rybak-Wolf et al, 2015; Veno et al, 2015). Hence, given the scarce knowledge about circRNA function and lack of studies reporting their sequence and expression in specific cell types, we here decided to exploit a double reporter mouse line previously characterized by our group and allowing the isolation of proliferating neural stem cells, neurogenic progenitors, and newborn neurons based on the combinatorial expression of RFP and GFP, respectively (Aprea et al, 2013).

Specifically, during cortical development, proliferative progenitors (PPs) progressively switch from divisions that expand their population to divisions that generate more committed, differentiative progenitors (DPs), which in turn are consumed to generate newborn neurons (N) (Lui et al, 2011; Taverna et al, 2014). Hence, to identify the three coexisting subpopulations of PP, DP, and N during mouse corticogenesis, our group has generated a double reporter mouse line expressing (i) RFP under the control of the Btg2 promoter and identifying the switch of PP to DP and (ii) GFP under the control of the Tubb3 promoter as a marker of newborn N (Aprea et al, 2013). Validating this approach, the combinatorial expression of the two reporters allowed the isolation of PP (RFP−/GFP−) from DP (RFP+/GFP−) and N (GFP+, irrespective of RFP) and identification and validation of transcription factors, lncRNAs, and epigenetic modifications functionally involved in neurogenic commitment (Aprea et al, 2013, 2015; Artegiani et al, 2015; Noack et al, 2019). Hence, we here decided to further exploit this mouse line and provide the arising field of circRNA biology with the first resource identifying circRNAs sequence and expression patterns in specific cell types of the developing mammalian cortex.

Results

Comprehensive identification of cell type–specific circRNAs of the mouse cortex

To identify cell type–specific circRNAs, we FAC-sorted PP, DP, and N, each in three biological replicates, from the mouse cortex at embryonic day (E) 14.5 as previously described (Aprea et al, 2013). Total RNA was then treated with RNase R to degrade linear RNAs, which we found to be very efficient in reducing the levels of linear transcript down to undetectable levels even among those with the highest expression levels and most stable predicted secondary structure (Fig S1A). This was then followed by 150-bp single-end, strand-specific, high-throughput sequencing. The reads were then aligned to the reference mouse genome (mm9) and unmapped reads used to identify the circularizing, backsplice junctions predictive of putative circRNAs (Fig S1B).

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Related to Fig 1.

(A) Validation of RNase R treatment on highly stable expressed transcript. Graphs represent the cycle threshold values (Ct; top) or calculated FC (bottom) upon quantitative real-time PCR of total RNA from E14.5 mouse cortex, with or without prior RNase R treatment (as indicated). For linear transcripts, regular primers were designed, whereas divergent primers were used for circEzh2. N = 3, n = 3; bars = SDs (*P < 0.05; **P < 0.01; ***P < 0.001; t test). (B) Sequencing quality was assessed by comparing the overall statistics (left) with the only other study using RNase R treatment before sequencing (Jeck et al, 2013). Total number of reads for each sample is reported (right). (C) DNA gel electrophoresis upon PCR amplification of 10 circRNAs predicted to include introns. Divergent primers were designed to include also the intronic sequence. Full-length gel is available in Fig S2. (D) Expressed circRNAs were identified as genic when overlapping in the sense strand of genes and the start/end base of the circRNA coincided with the start/end base of an annotated exon (top). If genic, but with start/end base falling within and intron or exon, we grouped the circular with the antisense and intergenic ones and generically termed them as other (bottom, top to bottom, respectively). (E) Primer designing strategy. Divergent primers were designed for specific detection of the circular-specific backsplice junction (top), whereas convergent, regular primer pairs were used to identify linear, non-circularized, mRNA (bottom).

As a first step to assess the sequence of circRNAs and their expression in cortical cell types, we collected their genomic coordinates putting together all replicates of PP, DP, and N and obtaining an initial set of 6,033 putative circRNAs (Fig 1A). Then, according to their genomic locations, we selected genic transcripts whose start and end sites coincided with the annotated start and end site of an exon and separated them from the remaining transcripts that included genic transcripts with ends not coinciding with exons, antisense or intergenic ones (Fig S1D). Within the former, we separated sequences belonging to introns from exons, whereas the latter were considered as a single exon. Finally, we calculated the relative reads per kilo base per million (RPKM) value of each intron and exon (Fig 1A and see the Materials and Methods section).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. circRNA prediction, sequence determination, and validation.

(A) Strategy for the identification of circRNAs and determination of their sequence and expression in cortical cell types. (B) DNA gel electrophoresis upon PCR amplification of 10 circRNAs identified as in (A). Full gels are available in Fig S2. Note the detection of bands with a molecular weight consistent with the expected exonic, but not intronic, regions. (C) Graph representing the cycle threshold (Ct) values upon RT-PCR of total RNA from the E14.5 mouse cortex, with or without a prior treatment with RNase R (as indicated) and using divergent primers for 30 predicted circRNAs together with linear and circular controls (red dashed line). Dashed black line indicates the Ct threshold of detection. N = 2, n = 3; bars = SDs (*P < 0.05; **P < 0.01; ***P < 0.001; t test). (D, E) Distribution of CiCo’s entries overlapping circRNAs deposited in circBase and (E) Venn diagram representing the proportion of circBase- and CiCo-specific circRNAs as well as common ones.

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Full-length gel with PCR products.

Bold circRNA IDs represent the original panel for each circRNA, whereas small-sized IDs are either for circRNA for which a second PCR was necessary to obtain enough material for cloning or unspecific product (determined after sequencing of the cloned PCR product). (A) Gel shown in Fig 1B. (B) Gel shown in Fig S1C.

Given the need to establish an unbiased minimum threshold of RPKM to define “expression,” we next selected 10 predicted genic circRNAs from our dataset, including six that were not described by previous studies. We then cloned and sequenced these circRNAs from RNase R–treated lysates from the E14.5 mouse lateral cortex and found that in all cases, the predicted exon(s) (2–15 for each circRNA; 45 in total) were included in their sequence, whereas not a single intron (35 in total) could be detected (Fig 1B). Hence, we chose the highest RPKM value previously calculated among these predicted, but not detected, introns as a minimum threshold to define expression (RPKM > 3.5). As a result, this allowed us to redefine among the original list of 6,033 putative circRNAs the 5,073 fulfilling our criteria for expression (Supplemental Data 1). Furthermore, we selected form this list six genic circRNAs for which introns were predicted as part of their sequence and designed divergent primers to validate their existence. Again, while confirming the presence of exons in these circRNAs, we were unable to confirm intronic sequences (Fig S1C), leading us to conclude that the genic circRNAs in our list are primarily exonic and that reads mapping on intronic locations may derive from lariats. As a result, we decided to exclude intronic sequences from any subsequent analyses (these sequences were, however, left in the list of expressed features provided in Supplemental Data 1).

Supplemental Data 1.

Expressed feature. Tab-delimited file containing exon/intron coordinates, the custom circID assigned and the type of feature expressed (exon or intron) as well as the average RPKM for each cell population. CiCo (Circular RNA from lateral Cortex). The first six columns are in BED format. In addition, the overlapped gene name (or "Intergenic"), the relative ENSEMBL gene ID and the number of exon included (in case of genic circRNAs; for antisense and for genic with unusual start/end point no exon is indicated) are reported. In the last two columns are recorded the circBase ID (if there is a correspondence) and the minimum level of overlap between the circRNAs. miRNA seed prediction. Tab-delimited file containing the total number of seeds predicted, the relative miRNA, the circRNA ID on which the seed was predicted and the length of the circRNA.[LSA-2019-00354_Supplemental_DataS1.xlsx]

Next, as a validation of our approach, we selected 30 circRNAs among this list of 5,073, including nine belonging to the bottom 30% and four to the bottom 10% in expression levels. Subsequent real time PCRs (RT-PCRs) were performed from the E14.5 mouse cortex, with or without RNase R treatment, using either divergent primers spanning over the backsplice junctions (Fig S1E) or convergent primers for two linear transcripts used as internal negative control of the enzymatic treatment (GAPDH and Ezh2_linear). This confirmed the presence of the vast majority (80%) of the selected circRNAs (Fig 1C). Importantly, failure to detect some circRNAs was independent from their predicted expression levels pointing out more a suboptimal choice of primers than false positives in our analysis. From here, we then reconstructed the comprehensive list, sequence, and expression levels of the 5,073 circRNAs detected in PP, DP, and N and representing the neurogenic lineage during corticogenesis that we named CiCo for CircRNAs of the mouse Cortex (Fig 1A and Supplemental Data 1).

We next compared CiCo with the most complete resource of circRNAs available: circBase (Glazar et al, 2014). In particular, because of the stage- and tissue-specific expression of circRNAs, we selected from circBase murine circRNAs predicted from whole embryos and brains. Considering total or partial overlaps, down to the single nucleotide, as known circRNAs we found that our dataset extended this list by nearly three-fold adding to the ∼1,900 transcripts of circBase ∼3,700 new circRNAs of CiCo (Fig 1D) for a total of ∼5,000 unique circRNAs (Fig 1E). As previously shown in the case of lncRNAs (Aprea et al, 2013, 2015), this high rate of novel circRNAs found in our study relative to previous reports highlights the power of cell type–specific analyses, allowing the identification of transcripts enriched in specific cell populations and being diluted out when considering bulk tissues or whole organs. In turn, this suggests a new layer of complexity in circRNA biology as these appeared to be not only stage- and tissue specific but also within the same developmental stage and tissue, to retain highly specific expression in individual cell types, which we sought to dissect next.

General features and differential expression of CiCo

Previous studies reported that ∼80% of the predicted circRNAs overlap in the sense strand with genes (Jeck et al, 2013; Memczak et al, 2013; Salzman et al, 2012). In our dataset, 97% of the expressed sequences overlapped genes with the remaining 3%, including either antisense or intergenic circRNAs (Fig 2A). This seemed in contrast with the similar proportion of sense and antisense transcripts among lncRNAs (Aprea & Calegari, 2015), which can be explained by the overall lower expression of this class of transcripts biasing against the detection of circRNAs derived from them. Regarding length distribution and exon density, we found that most circRNAs (90%) were less than 1 kb long, primarily 250–500 bp, and including on average 2–3 exons (Fig 2B). As expected, a linear correlation was found between the length of circRNAs and the number of exons that they included, with no specific bias in their distribution across or within chromosomes (data not shown).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. General features of CiCo.

(A) Number and proportion of CiCo and their genomic features showing they almost entirely (97%) derive from genic regions. (B) Length (left) and exon number (for genic transcripts, right) distribution of circRNAs. (C) Differential expression (50% i.e.: FC ≥1.5 or ≤0.67; no FDR being applied) of CiCo expressed in PP (grey), DP (red), and N (green). Abundance and proportion of circRNAs detected in each cell type are indicated and represented proportionally to the area of circles and pattern of expression.

Next, we analysed CiCo’s expression profiles of PP, DP, and N. In particular, we focused on the differentially expressed circRNAs considering a 50% threshold (i.e., a fold change [FC] by ≥1.5 or ≤0.67 for up- or down-regulation, respectively) (Fig 2C). Whereas a substantial proportion (42%) of circRNAs showed no significant change in expression among cell types, subdividing differentially expressed circRNAs into the possible patterns of up- or down-regulation during the neurogenic lineage revealed a distribution that was remarkably similar to that found in linear transcripts (Aprea et al, 2013) (Fig 2C). In addition, when analysing the biological functions and processes of the parental genes enriched in each group, we found that circRNAs up-regulated during neurogenesis were primarily associated with synaptogenesis and neuronal development. Conversely, down-regulated circRNAs were associated to cell cycle and regulation of transcription (data not shown).

Among transcripts showing specific expression patterns, an intriguing group of circRNAs emerged that was transiently up- or down-regulated specifically in DP compared with both PP and N (65 and 183, respectively) (Fig 2C). As previously shown by our group (Aprea et al, 2013), these patterns of expression account for a small, underrepresented proportion of differentially expressed transcripts (1–2%). Yet, at least among mRNAs and lncRNAs, many from this small subset of genes revealed to play key roles in corticogenesis (Aprea et al, 2013; Artegiani et al, 2015), raising the possibility that this also applies to circRNAs here identified for the first time.

Properties of circRNAs: miRNA-sponging potential, translation, and evolutionary conservation

Our novel assessment of the sequence of circRNAs gave us the possibility to gain new insights into their putative function(s). Since sponging of miRNAs by Cdr1as was the first suggested role for a circRNA (Hansen et al, 2013; Memczak et al, 2013), we searched within CiCo for miRNA seed sequences. When assessing the highest possible number of seed sequences present in any given circRNA that would potentially target any given miRNA, we found that the highest number was reached by Cdr1as (CiCo_mm9_circ_006933) for which 99 seeds specific for miR-7b and, among others, a single seed for miR-671 were predicted as previously described (Memczak et al, 2013). Despite this extreme example, only a handful of circRNAs displayed a number of seeds for a single miRNA that exceeded five (Fig 3A). Moreover, since this analysis did not account for the potential of circRNAs to sponge more miRNAs independently, we expanded our search to include all possible miRNA binding to a single circRNA with no regard to the number of seeds. With this approach, we identified in CiCo_mm9_circ_006344 the circRNA with the highest number of potentially binding miRNAs for a total of nearly 1,400 different miRNAs of which four (miR-1187, 466i, 466k, and 669c) with a high number of seeds (46, 43, 37, and 30 respectively; Fig 3A and Supplemental Data 1). However, such a high promiscuity in miRNA-binding sites and limitations in the bioinformatic tools to predict sponging properties makes it difficult to infer the biological significance of this circRNA.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Seed sequences, coding potential, and evolutionary conservation of circRNAs.

(A) Distribution of the maximum number of seed sequences for any miRNA found across CiCo. Note that the overwhelming majority of circRNAs only displayed few (1–4) seeds. (B) Total ORFs (left) and proportion of those with a CAI value above expected (right) found in CiCo or random control datasets (***P < 0.001; chi-squared test). (C) Whiskers–box plots showing the distribution of phyloP scores for CiCo sequences (orange) and relative shuffled controls (grey) (**P < 0.01; ***P < 0.001; t test).

In addition to miRNA sponging, a recent study has suggested that at least one circRNA, circ-ZNF609, is translated into a peptide (Legnini et al, 2017). Hence, we searched for possible open reading frames (ORFs) among all circRNAs and found that a remarkably high proportion (nearly 4,000; i.e., >70%) contained at least one putative ORFs with ≥150 nt in length (Fig 3B, left). To increase our confidence in this result, we assessed the Codon Adaptation Index (CAI) measuring the relative codon usage as a function of gene expression (Sharp & Li, 1987). To this end, we used CAIcal (Puigbo et al, 2008) to compute the index of each predicted ORF within circRNAs and referred this to the expected value calculated from a pool of 500 shuffled sequences reflecting the whole ORF population (eCAI = 0.744) that was used as a significance threshold. We found that 72% of all predicted ORFs and belonging to ∼50% of all circRNAs had an index higher than expected (Fig 3B, right), raising the possibility that these are potentially translated. However, 97% of CiCo’s sequences derive from annotated genes, implying that this abundance of ORFs might simply reflect their sharing of some of, if perhaps not all, the features of coding genes. To assess this, we repeated our analyses (ORF prediction and CAI calculation), but this time considering random sequences generated not only from the entire genome but also from the coding transcriptome as two independent negative controls. Although we could find a comparable number of putative ORFs from CiCo and our two negative control datasets (5,966, 3,239, and 4,372, respectively), only 52% of all random genomic ORFs and 66% of the transcriptomic ones were found to have a CAI value above the eCAI thresholds. Remarkably, when comparing the fraction of significant and nonsignificant ORFs within CiCo with the two control datasets, we found highly significant differences with both (chi-squared test, two-tailed P-value < 0.001) (Fig 3B, right) suggesting that the coding potential of CiCo is not only an inherited feature resulting from their origin from coding genes but might potentially account for their function.

We next assessed the evolutionary conservation of CiCo. To this end, we considered not only the sequence of the mature circRNAs but also their flanking regions (200 bp up- and down-stream) as these might be important for their biogenesis (Jeck et al, 2013; Ashwal-Fluss et al, 2014). Again to obtain an appropriate reference as negative control, we generated for each circRNA and flanking regions a random sequence of the same size and reflecting the features of such circRNA. Specifically, we used any random genomic region for intergenic circRNAs or transcriptome-specific and intronic sequences for genic circRNAs and their flanking regions, respectively. We next computed the average conservation score of all these regions, making use of the phyloP score (a per-base value obtained by the alignment of murine genome against 30 other vertebrates, UCSC Genome Browser [Casper et al, 2018]) and compared it with the one calculated for the respective random sequences. We found that genic circRNAs were more conserved than their reference random sequences and that despite an average lower conservation score, their flanking regions were also significantly more conserved (Fig 3C; top). As this again can be influenced by circRNA origin from expressed genes, we next assessed the conservation of the subgroup of intergenic circRNAs that, by definition, are not overlapping any annotated gene. Strikingly, both the circRNAs themselves as well as their 200-bp up- and down-stream regions were also significantly more conserved than their random counterparts (Fig 3C; bottom). Hence, similarly to their coding potential, the conservation of circRNAs seems not only an inherited feature resulting from their sharing of sequences with expressed genes but may actually reflect their function.

CircRNA biogenesis is independent from alternative splicing

Finally, because nearly all (97%) CiCo transcripts are overlapping genes, we asked if the abundance of circRNAs correlated with that of their respective linear mRNA, which could reveal the mechanisms underlying their biogenesis. In this context, we sought possible mechanisms by which biogenesis of circRNAs may be controlled. We speculated that this may occur at the level of transcription by the synthesis of different pre-RNAs that ultimately mature into circular, rather than linear, transcripts. Alternatively, a circular or linear RNA may result by splicing of a single pre-RNA precursor. Although not mutually exclusive, the two mechanisms may occur either un-specifically without any regulation at the level of individual transcripts or cell types or, alternatively, only apply to a specific subclass of genes. Ultimately, distinguishing between these possibilities may help understanding the biological significance of circRNAs.

To this end, we took advantage of the previous poly(A) transcriptome assessment of PP, DP, and N reported by our group (Aprea et al, 2013, 2015) and compared the expression levels of each circRNA and its corresponding linear counterpart in each cell type. We found a strong positive correlation between the expression of mRNAs and circRNAs (Pearson’s score ∼0.8 in all cell populations) (Fig 4A), implying that an increase in a circRNA did not result in a decrease in its linear form. This in turn led us to reject the hypothesis of regulation at the level of transcription of different pre-RNAs. Next, if synthesis of a circular versus linear RNA should be regulated by splicing from a common transcript, then an increase in such circRNA should inversely correlate with the usage of the exon(s) in common with the linear transcript. To assess this, we performed differential exon usage analyses of PP, DP, and N (to be described elsewhere) and assessed whether exon(s) used in common by both a circular and linear transcript would display an inverse FC from one cell population to the following (e.g., from PP to DP and from DP to N). Surprisingly, barely 0.5% of all circRNA–exon pairs showed an opposite FC pattern with regard to their linear counterparts, whereas the vast majority shared the same pattern of either up- or down-regulation (Fig 4B). In turn, this questioned a mechanistic link between mRNA splicing per se and circRNA biogenesis. To address this, we next selected all genes resulting in the expression of at least two linear isoforms (4,146 genes out of 19,802 expressed) resulting from alternative splicing as one specific mechanism of splicing and investigated what proportion of circular RNA was produced among these genes. We found that only 49% of circRNAs were generated from alternatively spliced genes and that among these, only ∼0.2% shared the very same exon(s) in all three cell populations with the differentially spliced linear transcript (Fig 4C).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. CiCo relation to linear transcripts.

(A) Scatterplot log2(RPKM) of CiCo and mRNA counterparts in cell types (colours as in Fig 2C). R2 for each regression line is indicated. (B) Number and proportion of circRNAs having either an opposite or the same FC pattern as the exon of the linear counterpart (50% i.e.,: FC ≥1.5 or ≤0.67; no FDR being applied; colours as in Fig 2C). (C) Number of expressed circRNAs that result from differentially spliced genes during the neurogenic lineage.

Taken together, our data indicate that the overall expression of circRNAs correlates with that of their linear RNAs but is not a byproduct limited to transcripts and exons undergoing alternative splicing and, hence, can be regulated independently by mechanisms that are still unknown.

Discussion

Here, we provided the arising field of circRNA biology with the first resource describing their sequence and cell-specific expression during mammalian cortical development. This has allowed us to increase by threefold the number of circRNAs currently known, pointing out their full diversity and specificity in cell populations representing the lineage to neurogenic commitment. More importantly, our work has allowed us to reveal novel features of these elusive transcripts that are important to infer their significance and that could not be deduced based on the previous knowledge of backsplice junctions alone.

We found that genic circRNAs are primarily composed by exons of genes encoding for mRNAs or lncRNAs and that a strong correlation existed in the levels of expression of circular and linear pairs of transcripts. Our data were inconsistent with both the synthesis of circular-specific, pre-RNAs during transcription as well as the involvement of alternative splicing as a mechanism underlying their biogenesis. Although more studies are needed to address these aspects, certain features emerging from our study seemed to support the notion that circRNAs may be a generic byproduct of splicing, whereas others highlighted their specificity both in terms of biogenesis and putative biological significance. To start with, CiCo transcripts were found to derive only from a specific group of linear transcripts and consistently included certain exons but not others. In addition, both the coding potential and evolutionary conservation of CiCo revealed to be much higher than expected by chance even after accounting for their origin from exonic regions.

Although it is clear that a significant proportion of circRNAs might nonetheless lack a function despite these features, our study provides the field with a new resource for the identification of biologically relevant circRNAs, fuelling future studies on their molecular regulation and role.

Materials and Methods

Animal care, cell sorting, and RNA extraction

Animal experiments were approved by the Landesdirektion Sachsen (24-9168.11-1/41 and TVV 39/2015) and carried out in accordance with the relevant guidelines and regulation. Pregnant Btg2RFP/Tubb3GFP double heterozygous mice were anaesthetised with isoflurane and euthanized through cervical dislocation. Brains from E14.5 embryos were collected and the lateral cortex isolated after removal of the meninges and ganglionic eminences. The Neural Tissue Dissociation kit with papain (Miltenyi Biotech) was used according to the manufacturer’s protocol to obtain a cell suspension for FAC-sorting or RNA extraction. The cells were resuspended in ice-cold PBS and supplied with 10 μl of 7-AAD to assess cell viability (BD Pharmigen) and sorted with a BD Aria III FACS as previously described (Aprea et al, 2013). Sorted cells were collected in PBS, centrifuged (600 g, 5 min), and RNA extracted using Quick RNA Mini Prep (Zymo Research) according to the manufacturer’s protocol.

CircRNA sequencing, annotation, and validation

Total RNA was denatured for 3 min at 70°C, then treated with RNase R (Epicentre) for 1 h at 40°C, and finally supplied with DNase I (Invitrogen) for 15 min at room temperature. The reaction was then cleaned with RNA Clean & Concentrator (Zymo Research) and cDNA libraries prepared using the NEB Next Ultra Directional RNA Library Prep kit without mRNA enrichment. Effectiveness of RNase R treatment was assessed by quantitative real time PCR of six transcripts with the highest free energy per nucleotide in their most stable predicted secondary structure according to RNAfold (-p --noLP --temp = 40; primers in Table S1C). Samples were sequenced on an Illumina HiSeq 2500 with a read length of 150 bp and resulting reads were aligned using gsnap (Wu & Nacu, 2010) with mm9 and ENSEMBL Genes (v67 [Flicek et al, 2013]) as the reference genome and trancriptome, respectively. Besides adapter removal, reads were not trimmed, and overall quality was evaluated by the percentages of mapped and unmapped reads that resulted in similar values to the only comparable study reported to date (Jeck et al, 2013) (Fig S1B). To predict putative circRNA, unmapped reads were retrieved and then analysed using the “find_circ” pipeline (version 1) with default parameters (Memczak et al, 2013). No filter was applied on the number of reads identifying the circularising junction. For the genomic location, predicted circRNAs were overlapped with mouse ENSEMBL Genes (v67 [Flicek et al, 2013]) and considered them genic when their start/end base coincided with the start/end of annotated exon(s). circRNA overlapping genes but with other start/end points were grouped together with intergenic and antisense ones (Fig 2A, termed as “other”). Overlap with circBase was carried out using bedtools (Quinlan & Hall, 2010) and intersecting predicted circRNAs with the .bed file relative to mouse (Glazar et al, 2014). Different levels of minimal reciprocal overlap were set within bedtools options to take into account also single-nucleotide overlap between data. For validation, 1 µg of total RNA from un-sorted lateral cortex was treated as described previously for sequencing, with one sample digested either with RNase R or water and then converted into cDNA using 200 U of reverse transcriptase and 50 ng of random hexamers according to the SuperScript III (Invitrogen) kit. The resulting cDNA was diluted either 1:10 for RNase R or 1:100 in case of water treatment and 1 μl used as a PCR template with divergent primers (designed with Primer3 ([Untergasser et al, 2012]) and spanning over the circularising junction. All primer pairs were tested with the In-Silico PCR tool from UCSC to minimize byproducts (Kent et al, 2002) (primers in Table S1A).

Table S1 Primer List. (A) List of divergent primers for all the tested circRNAs in Fig 1C. (B) Primers used for cloning and sequencing in Fig 1B. (C) Primers for linear transcript used in RNase R efficacy test, Fig S1A. (D) List of divergent primers for intron validation, Fig S1C.

Sequence prediction, differential expression, and conservation

First, exonic and intronic sequences of previously annotated genic circRNAs were separated. Next, we retrieved the coordinates of these exon(s) from ENSEMBL Genes (v67 [Flicek et al, 2013]), whereas intron coordinates were obtained through command line starting from the exonic ones. circRNAs annotated as intergenic, antisense, or genic with unusual start/end points were considered as a single exon. Exonic and intronic coordinates were kept separate, and featureCounts (Liao et al, 2014) was run using them as reference files to obtain a per-feature read count on which RPKM values were calculated as follows: RPKM=read countR/L, where R=library size×10−6 and L=circ length×10−3. RPKM were averaged across biological replicates. Threshold for expression was set as the highest RPKM of the predicted, but not detected, introns within 10 validated circRNAs, with additional six circRNA used to validate the absence of introns from which we redefined the set of expressed features (primers in Table S1B and D). For the differential expression analysis, we run again featureCounts (Liao et al, 2014) using the set of expressed features as reference and specifying a meta-feature count to automatically sum the reads from different features of the same circRNA. The resulting table was analysed with DESeq2 (Love et al, 2014), considering FCs by ≥1.5 or ≤0.67 (no FDR applied). Clustering analysis was used to evaluate enriched terms using DAVID (Huang da et al, 2009). Finally, circRNA conservation was assessed by computing the phyloP score for the expressed sequence as well as for 200 bp up- and down-stream of each circRNA. To compare our sequences, we generated a shuffled version with the same sizes using bedtools (shuffleBed) and accounting for their genomic location. In particular, for genic circRNAs, the shuffled sequences were obtained from a reference file that included exons and/or introns as appropriate and derived from ENSEMBL (v67). In the case of intergenic circRNAs (both sequence and flanking regions), the shuffled sequences were chosen from the entire genome. BEDOPS (Neph et al, 2012) was used to sort regions per genomic location and to split into a per-base coordinate file. BEDOPS was again used to retrieve the per-base phyloP score from the UCSC table containing the conservation score for mouse versus 30 vertebrate genomes (phyloP30wayAll). For shuffled sequences relative to genic circRNAs, we built a custom phyloP score table in which the values relative to all exons were included. We finally reconstructed the score per sequence by averaging the phyloP values for each circRNA. These steps were performed for all files (original and shuffled sequences).

Prediction of miRNA seeds, ORFs, and exon usage correlation

Seed prediction was performed using miRanda (John et al, 2004) with default parameters and with input the FASTA file for the circRNA sequences obtained with bedtools (getfasta, with the expressed feature coordinates file as input) and miRBase-downloaded mouse mature miRNA sequences (v22) (Kozomara & Griffiths-Jones, 2014). The output generated was then parsed to a more manageable form, and the seed category was added according to Bartel (Bartel, 2009). To count the number of seed for each miRNA present on a circRNA, we first subset the parsed miRanda output by seed category (removed everything that did not have a recognized seed type), alignment score (≥150), and by free energy (∆G ≤ −19) (Bartel, 2009), and then we counted and reported the number of unique miRNA-circRNA combination and sorted by number of occurrences (top to bottom). We predicted the presence of ORF using the standalone version of National Center for Biotechnology Information’s ORFfinder (Sayers et al, 2011) using the FASTA sequences of the expressed circRNA as input. We restricted the search by requiring a minimum ORF length of 150 nt and setting the start codon as ATG only; we also ignored nested ORFs and searched only on the strand of the sequence itself. We then supplied the resulting sequences to the CAIcal web server (Puigbo et al, 2008) together with the codon usage for mouse (Nakamura et al, 2000). Expected CAI value (eCAI) was calculated several times using the Markov method with 95% of confidence and the final value obtained by averaging all the calculated ones. As controls, we generated two shuffled dataset using as reference either the entire genome or the transcriptome only. For this two control datasets, the same ORF prediction and CAI calculation were performed. Differences of our dataset with the two control ones were assessed through chi-squared test with one degree of freedom. Exon usage data were obtained by poly(A) enrichment and PE sequencing of our three cell populations (to be reported elsewhere). For the full mRNA, we selected the RPKM values for the linear corresponding to the circRNAs. For linear versus circRNA expression, we computed the average RPKM for each cell population and compared the log2(RPKM), plotting their relative abundance and computing an overall Pearson’s correlation score. In case of linear exon versus circRNA, the corresponding read counts were analysed with DESeq2 to compute the FC of the exon shared by circRNAs and linear transcript.

Data availability

All custom scripts can be obtained upon request. Sequencing data generated during the present study are available at Gene Expression Omnibus repository (GSE117009).

Acknowledgements

We thank the facilities at Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) and Center for Regenerative Therapies Dresden (CRTD) for maintenance of mouse lines, sequencing, and FACS. We are grateful to the members of Silvio Bicciato’s Lab (University of Modena) and Dr. Michael Hiller (MPI-CBG) for suggestions and discussion. We also thank Prof. Katja Nowick and Maria Beatriz Walter Costa for insights into RNA structure prediction. This work was funded by the CRTD, the School of Medicine of the TU-Dresden, the Italian Epigenomics Flagship Project (Epigen), the Deutsche Forschungsgemeinschaft (DFG) grant DFG CA 893/9-1, and a fellowship from the Dresden International Graduate School for Biomedicine and Bioengineering (DIGS-BB) awarded to M Dori.

Author Contributions

  • M Dori: conceptualization, data curation, formal analysis, validation, investigation, methodology, and writing—original draft, review, and editing.

  • L Haj Abdullah Alieh: validation and investigation.

  • D Cavalli: validation and investigation.

  • S Massalini: validation and investigation.

  • M Lesche: data curation.

  • A Dahl: data curation.

  • F Calegari: conceptualization, supervision, funding acquisition, and writing—original draft, review, and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received February 21, 2019.
  • Revision received March 22, 2019.
  • Accepted March 22, 2019.
  • © 2019 Dori et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Abu N,
    2. Jamal R
    (2016) Circular RNAs as promising biomarkers: A mini-review. Front Physiol 7: 355. doi:10.3389/fphys.2016.00355
    OpenUrlCrossRef
  2. ↵
    1. Aprea J,
    2. Calegari F
    (2015) Long non-coding RNAs in corticogenesis: Deciphering the non-coding code of the brain. EMBO J 34: 2865–2884. doi:10.15252/embj.201592655
    OpenUrlCrossRef
  3. ↵
    1. Aprea J,
    2. Lesche M,
    3. Massalini S,
    4. Prenninger S,
    5. Alexopoulou D,
    6. Dahl A,
    7. Hiller M,
    8. Calegari F
    (2015) Identification and expression patterns of novel long non-coding RNAs in neural progenitors of the developing mammalian cortex. Neurogenesis (Austin) 2: e995524. doi:10.1080/23262133.2014.995524
    OpenUrlCrossRef
  4. ↵
    1. Aprea J,
    2. Prenninger S,
    3. Dori M,
    4. Ghosh T,
    5. Monasor LS,
    6. Wessendorf E,
    7. Zocher S,
    8. Massalini S,
    9. Alexopoulou D,
    10. Lesche M, et al.
    (2013) Transcriptome sequencing during mouse brain development identifies long non-coding RNAs functionally involved in neurogenic commitment. EMBO J 32: 3145–3160. doi:10.1038/emboj.2013.245
    OpenUrlCrossRef
  5. ↵
    1. Artegiani B,
    2. de Jesus Domingues AM,
    3. Bragado Alonso S,
    4. Brandl E,
    5. Massalini S,
    6. Dahl A,
    7. Calegari F
    (2015) Tox: A multifunctional transcription factor and novel regulator of mammalian corticogenesis. EMBO J 34: 896–910. doi:10.15252/embj.201490061
    OpenUrlCrossRef
  6. ↵
    1. Ashwal-Fluss R,
    2. Meyer M,
    3. Pamudurti NR,
    4. Ivanov A,
    5. Bartok O,
    6. Hanan M,
    7. Evantal N,
    8. Memczak S,
    9. Rajewsky N,
    10. Kadener S
    (2014) circRNA biogenesis competes with pre-mRNA splicing. Mol Cell 56: 55–66. doi:10.1016/j.molcel.2014.08.019
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bachmayr-Heyda A,
    2. Reiner AT,
    3. Auer K,
    4. Sukhbaatar N,
    5. Aust S,
    6. Bachleitner-Hofmann T,
    7. Mesteri I,
    8. Grunt TW,
    9. Zeillinger R,
    10. Pils D
    (2015) Correlation of circular RNA abundance with proliferation: Exemplified with colorectal and ovarian cancer, idiopathic lung fibrosis, and normal human tissues. Sci Rep 5: 8057. doi:10.1038/srep08057
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bahn JH,
    2. Zhang Q,
    3. Li F,
    4. Chan TM,
    5. Lin X,
    6. Kim Y,
    7. Wong DT,
    8. Xiao X
    (2015) The landscape of microRNA, Piwi-interacting RNA, and circular RNA in human saliva. Clin Chem 61: 221–230. doi:10.1373/clinchem.2014.230433
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Bartel DP
    (2009) MicroRNAs: Target recognition and regulatory functions. Cell 136: 215–233. doi:10.1016/j.cell.2009.01.002
    OpenUrlCrossRefPubMed
  10. ↵
    1. Capel B,
    2. Swain A,
    3. Nicolis S,
    4. Hacker A,
    5. Walter M,
    6. Koopman P,
    7. Goodfellow P,
    8. Lovell-Badge R
    (1993) Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 73: 1019–1030. doi:10.1016/0092-8674(93)90279-y
    OpenUrlCrossRefPubMed
  11. ↵
    1. Casper J,
    2. Zweig AS,
    3. Villarreal C,
    4. Tyner C,
    5. Speir ML,
    6. Rosenbloom KR,
    7. Raney BJ,
    8. Lee CM,
    9. Lee BT,
    10. Karolchik D, et al.
    (2018) The UCSC Genome Browser database: 2018 update. Nucleic Acids Res 46: D762–D769. doi:10.1093/nar/gkx1020
    OpenUrlCrossRefPubMed
  12. ↵
    1. Filippenkov IB,
    2. Sudarkina OY,
    3. Limborska SA,
    4. Dergunova LV
    (2015) Circular RNA of the human sphingomyelin synthase 1 gene: Multiple splice variants, evolutionary conservatism and expression in different tissues. RNA Biol 12: 1030–1042. doi:10.1080/15476286.2015.1076611
    OpenUrlCrossRef
  13. ↵
    1. Flicek P,
    2. Ahmed I,
    3. Amode MR,
    4. Barrell D,
    5. Beal K,
    6. Brent S,
    7. Carvalho-Silva D,
    8. Clapham P,
    9. Coates G,
    10. Fairley S, et al.
    (2013) Ensembl 2013. Nucleic Acids Res 41: D48–D55. doi:10.1093/nar/gks1236
    OpenUrlCrossRefPubMed
  14. ↵
    1. Glazar P,
    2. Papavasileiou P,
    3. Rajewsky N
    (2014) circBase: A database for circular RNAs. RNA 20: 1666–1670. doi:10.1261/rna.043687.113
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Hansen TB,
    2. Jensen TI,
    3. Clausen BH,
    4. Bramsen JB,
    5. Finsen B,
    6. Damgaard CK,
    7. Kjems J
    (2013) Natural RNA circles function as efficient microRNA sponges. Nature 495: 384–388. doi:10.1038/nature11993
    OpenUrlCrossRefPubMed
  16. ↵
    1. Huang da W,
    2. Sherman BT,
    3. Lempicki RA
    (2009) Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37: 1–13. doi:10.1093/nar/gkn923
    OpenUrlCrossRefPubMed
  17. ↵
    1. Jakobi T,
    2. Czaja-Hasse LF,
    3. Reinhardt R,
    4. Dieterich C
    (2016) Profiling and validation of the circular RNA repertoire in adult murine hearts. Genomics Proteomics Bioinformatics 14: 216–223. doi:10.1016/j.gpb.2016.02.003
    OpenUrlCrossRefPubMed
  18. ↵
    1. Jeck WR,
    2. Sorrentino JA,
    3. Wang K,
    4. Slevin MK,
    5. Burd CE,
    6. Liu J,
    7. Marzluff WF,
    8. Sharpless NE
    (2013) Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 19: 141–157. doi:10.1261/rna.035667.112
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. John B,
    2. Enright AJ,
    3. Aravin A,
    4. Tuschl T,
    5. Sander C,
    6. Marks DS
    (2004) Human MicroRNA targets. PLoS Biol 2: e363. doi:10.1371/journal.pbio.0020363
    OpenUrlCrossRefPubMed
  20. ↵
    1. Kent WJ,
    2. Sugnet CW,
    3. Furey TS,
    4. Roskin KM,
    5. Pringle TH,
    6. Zahler AM,
    7. Haussler D
    (2002) The human genome browser at UCSC. Genome Res 12: 996–1006. doi:10.1101/gr.229102
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Kos A,
    2. Dijkema R,
    3. Arnberg AC,
    4. van der Meide PH,
    5. Schellekens H
    (1986) The hepatitis delta (delta) virus possesses a circular RNA. Nature 323: 558–560. doi:10.1038/323558a0
    OpenUrlCrossRefPubMed
  22. ↵
    1. Kosik KS
    (2013) Molecular biology: Circles reshape the RNA world. Nature 495: 322–324. doi:10.1038/nature11956
    OpenUrlCrossRefPubMed
  23. ↵
    1. Kozomara A,
    2. Griffiths-Jones S
    (2014) miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res 42: D68–D73. doi:10.1093/nar/gkt1181
    OpenUrlCrossRefPubMed
  24. ↵
    1. Kristensen LS,
    2. Okholm TLH,
    3. Veno MT,
    4. Kjems J
    (2018) Circular RNAs are abundantly expressed and upregulated during human epidermal stem cell differentiation. RNA Biol 15: 280–291. doi:10.1080/15476286.2017.1409931
    OpenUrlCrossRef
  25. ↵
    1. Kulcheski FR,
    2. Christoff AP,
    3. Margis R
    (2016) Circular RNAs are miRNA sponges and can be used as a new class of biomarker. J Biotechnol 238: 42–51. doi:10.1016/j.jbiotec.2016.09.011
    OpenUrlCrossRef
  26. ↵
    1. Lasda E,
    2. Parker R
    (2014) Circular RNAs: Diversity of form and function. RNA 20: 1829–1842. doi:10.1261/rna.047126.114
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Legnini I,
    2. Di Timoteo G,
    3. Rossi F,
    4. Morlando M,
    5. Briganti F,
    6. Sthandier O,
    7. Fatica A,
    8. Santini T,
    9. Andronache A,
    10. Wade M, et al.
    (2017) Circ-ZNF609 is a circular RNA that can Be translated and functions in myogenesis. Mol Cell 66: 22–37 e29. doi:10.1016/j.molcel.2017.02.017
    OpenUrlCrossRefPubMed
  28. ↵
    1. Liao Y,
    2. Smyth GK,
    3. Shi W
    (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30: 923–930. doi:10.1093/bioinformatics/btt656
    OpenUrlCrossRefPubMed
  29. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  30. ↵
    1. Lui JH,
    2. Hansen DV,
    3. Kriegstein AR
    (2011) Development and evolution of the human neocortex. Cell 146: 18–36. doi:10.1016/j.cell.2011.06.030
    OpenUrlCrossRefPubMed
  31. ↵
    1. Memczak S,
    2. Jens M,
    3. Elefsinioti A,
    4. Torti F,
    5. Krueger J,
    6. Rybak A,
    7. Maier L,
    8. Mackowiak SD,
    9. Gregersen LH,
    10. Munschauer M, et al.
    (2013) Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 495: 333–338. doi:10.1038/nature11928
    OpenUrlCrossRefPubMed
  32. ↵
    1. Memczak S,
    2. Papavasileiou P,
    3. Peters O,
    4. Rajewsky N
    (2015) Identification and characterization of circular RNAs as a new class of putative biomarkers in human blood. PLoS One 10: e0141214. doi:10.1371/journal.pone.0141214
    OpenUrlCrossRefPubMed
  33. ↵
    1. Nakamura Y,
    2. Gojobori T,
    3. Ikemura T
    (2000) Codon usage tabulated from international DNA sequence databases: Status for the year 2000. Nucleic Acids Res 28: 292. doi:10.1093/nar/28.1.292
    OpenUrlCrossRefPubMed
  34. ↵
    1. Neph S,
    2. Kuehn MS,
    3. Reynolds AP,
    4. Haugen E,
    5. Thurman RE,
    6. Johnson AK,
    7. Rynes E,
    8. Maurano MT,
    9. Vierstra J,
    10. Thomas S, et al.
    (2012) BEDOPS: High-performance genomic feature operations. Bioinformatics 28: 1919–1920. doi:10.1093/bioinformatics/bts277
    OpenUrlCrossRefPubMed
  35. ↵
    1. Nigro JM,
    2. Cho KR,
    3. Fearon ER,
    4. Kern SE,
    5. Ruppert JM,
    6. Oliner JD,
    7. Kinzler KW,
    8. Vogelstein B
    (1991) Scrambled exons. Cell 64: 607–613. doi:10.1016/0092-8674(91)90244-s
    OpenUrlCrossRefPubMed
  36. ↵
    1. Noack F,
    2. Pataskar A,
    3. Schneider M,
    4. Buchholz F,
    5. Tiwari VK,
    6. Calegari F
    (2019) Assessment and site-specific manipulation of DNA (hydroxy-)methylation during mouse corticogenesis. Life Sci Alliance 2: e201900331. doi:10.26508/lsa.201900331
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Piwecka M,
    2. Glazar P,
    3. Hernandez-Miranda LR,
    4. Memczak S,
    5. Wolf SA,
    6. Rybak-Wolf A,
    7. Filipchyk A,
    8. Klironomos F,
    9. Cerda Jara CA,
    10. Fenske P, et al.
    (2017) Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science 357: eaam8526. doi:10.1126/science.aam8526
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Puigbo P,
    2. Bravo IG,
    3. Garcia-Vallve S
    (2008) CAIcal: A combined set of tools to assess codon usage adaptation. Biol Direct 3: 38. doi:10.1186/1745-6150-3-38
    OpenUrlCrossRefPubMed
  39. ↵
    1. Quinlan AR,
    2. Hall IM
    (2010) BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 26: 841–842. doi:10.1093/bioinformatics/btq033
    OpenUrlCrossRefPubMed
  40. ↵
    1. Rybak-Wolf A,
    2. Stottmeister C,
    3. Glazar P,
    4. Jens M,
    5. Pino N,
    6. Giusti S,
    7. Hanan M,
    8. Behm M,
    9. Bartok O,
    10. Ashwal-Fluss R, et al.
    (2015) Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell 58: 870–885. doi:10.1016/j.molcel.2015.03.027
    OpenUrlCrossRefPubMed
  41. ↵
    1. Salzman J,
    2. Gawad C,
    3. Wang PL,
    4. Lacayo N,
    5. Brown PO
    (2012) Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS One 7: e30733. doi:10.1371/journal.pone.0030733
    OpenUrlCrossRefPubMed
  42. ↵
    1. Sayers EW,
    2. Barrett T,
    3. Benson DA,
    4. Bolton E,
    5. Bryant SH,
    6. Canese K,
    7. Chetvernin V,
    8. Church DM,
    9. DiCuccio M,
    10. Federhen S, et al.
    (2011) Database resources of the national center for biotechnology information. Nucleic Acids Res 39: D38–D51. doi:10.1093/nar/gkq1172
    OpenUrlCrossRefPubMed
  43. ↵
    1. Sharp PM,
    2. Li WH
    (1987) The Codon Adaptation Index: A measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res 15: 1281–1295. doi:10.1093/nar/15.3.1281
    OpenUrlCrossRefPubMed
  44. ↵
    1. Suzuki H,
    2. Zuo Y,
    3. Wang J,
    4. Zhang MQ,
    5. Malhotra A,
    6. Mayeda A
    (2006) Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing. Nucleic Acids Res 34: e63. doi:10.1093/nar/gkl151
    OpenUrlCrossRefPubMed
  45. ↵
    1. Szabo L,
    2. Salzman J
    (2016) Detecting circular RNAs: Bioinformatic and experimental challenges. Nat Rev Genet 17: 679–692. doi:10.1038/nrg.2016.114
    OpenUrlCrossRef
  46. ↵
    1. Taverna E,
    2. Gotz M,
    3. Huttner WB
    (2014) The cell biology of neurogenesis: Toward an understanding of the development and evolution of the neocortex. Annu Rev Cell Dev Biol 30: 465–502. doi:10.1146/annurev-cellbio-101011-155801
    OpenUrlCrossRefPubMed
  47. ↵
    1. Untergasser A,
    2. Cutcutache I,
    3. Koressaar T,
    4. Ye J,
    5. Faircloth BC,
    6. Remm M,
    7. Rozen SG
    (2012) Primer3: New capabilities and interfaces. Nucleic Acids Res 40: e115. doi:10.1093/nar/gks596
    OpenUrlCrossRefPubMed
  48. ↵
    1. Veno MT,
    2. Hansen TB,
    3. Veno ST,
    4. Clausen BH,
    5. Grebing M,
    6. Finsen B,
    7. Holm IE,
    8. Kjems J
    (2015) Spatio-temporal regulation of circular RNA expression during porcine embryonic brain development. Genome Biol 16: 245. doi:10.1186/s13059-015-0801-3
    OpenUrlCrossRefPubMed
  49. ↵
    1. Vicens Q,
    2. Westhof E
    (2014) Biogenesis of circular RNAs. Cell 159: 13–14. doi:10.1016/j.cell.2014.09.005
    OpenUrlCrossRefPubMed
  50. ↵
    1. Vincent HA,
    2. Deutscher MP
    (2006) Substrate recognition and catalysis by the exoribonuclease RNase R. J Biol Chem 281: 29769–29775. doi:10.1074/jbc.m606744200
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Westholm JO,
    2. Miura P,
    3. Olson S,
    4. Shenker S,
    5. Joseph B,
    6. Sanfilippo P,
    7. Celniker SE,
    8. Graveley BR,
    9. Lai EC
    (2014) Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation. Cell Rep 9: 1966–1980. doi:10.1016/j.celrep.2014.10.062
    OpenUrlCrossRefPubMed
  52. ↵
    1. Wu TD,
    2. Nacu S
    (2010) Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics 26: 873–881. doi:10.1093/bioinformatics/btq057
    OpenUrlCrossRefPubMed
  53. ↵
    1. Yang X,
    2. Yuan W,
    3. Tao J,
    4. Li P,
    5. Yang C,
    6. Deng X,
    7. Zhang X,
    8. Tang J,
    9. Han J,
    10. Wang J, et al.
    (2017) Identification of circular RNA signature in bladder cancer. J Cancer 8: 3456–3463. doi:10.7150/jca.19345
    OpenUrlCrossRef
  54. ↵
    1. Zaphiropoulos PG
    (1996) Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: Correlation with exon skipping. Proc Natl Acad Sci U S A 93: 6536–6541. doi:10.1073/pnas.93.13.6536
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Download PDF
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Sequence and expression levels of circular RNAs in progenitor cell types during mouse corticogenesis
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Circular RNA assessment in the developing mouse cortex
Martina Dori, Leila Haj Abdullah Alieh, Daniel Cavalli, Simone Massalini, Mathias Lesche, Andreas Dahl, Federico Calegari
Life Science Alliance Mar 2019, 2 (2) e201900354; DOI: 10.26508/lsa.201900354

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Circular RNA assessment in the developing mouse cortex
Martina Dori, Leila Haj Abdullah Alieh, Daniel Cavalli, Simone Massalini, Mathias Lesche, Andreas Dahl, Federico Calegari
Life Science Alliance Mar 2019, 2 (2) e201900354; DOI: 10.26508/lsa.201900354
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
Issue Cover

In this Issue

Volume 2, No. 2
April 2019
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Development
  • Molecular Biology
  • Neuroscience

Related Articles

  • No related articles found.

Cited By...

  • MicroRNA profiling of mouse cortical progenitors and neurons reveals miR-486-5p as a regulator of neurogenesis
  • microRNA profiling of mouse cortical progenitors and neurons reveals miR-486-5p as a novel regulator of neurogenesis
  • Google Scholar

More in this TOC Section

  • Comparison of mitochondrial imaging in aging C. elegans
  • RNA profile by single cell analysis of severe dengue in mice
  • Pangenome-based disease genomics
Show more Resource

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Bluesky
  • X/Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2025 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.