MicroRNA-335-5p suppresses voltage-gated sodium channel expression and may be a target for seizure control.

There remains an urgent need for new therapies for treatment-resistant epilepsy. Sodium channel blockers are effective for seizure control in common forms of epilepsy, but loss of sodium channel function underlies some genetic forms of epilepsy. Approaches that provide bidirectional control of sodium channel expression are needed. MicroRNAs (miRNA) are small noncoding RNAs which negatively regulate gene expression. Here we show that genome-wide miRNA screening of hippocampal tissue from a rat epilepsy model, mice treated with the antiseizure medicine cannabidiol, and plasma from patients with treatment-resistant epilepsy, converge on a single target-miR-335-5p. Pathway analysis on predicted and validated miR-335-5p targets identified multiple voltage-gated sodium channels (VGSCs). Intracerebroventricular injection of antisense oligonucleotides against miR-335-5p resulted in upregulation of Scn1a, Scn2a, and Scn3a in the mouse brain and an increased action potential rising phase and greater excitability of hippocampal pyramidal neurons in brain slice recordings, consistent with VGSCs as functional targets of miR-335-5p. Blocking miR-335-5p also increased voltage-gated sodium currents and SCN1A, SCN2A, and SCN3A expression in human induced pluripotent stem cell-derived neurons. Inhibition of miR-335-5p increased susceptibility to tonic-clonic seizures in the pentylenetetrazol seizure model, whereas adeno-associated virus 9-mediated overexpression of miR-335-5p reduced seizure severity and improved survival. These studies suggest modulation of miR-335-5p may be a means to regulate VGSCs and affect neuronal excitability and seizures. Changes to miR-335-5p may reflect compensatory mechanisms to control excitability and could provide biomarker or therapeutic strategies for different types of treatment-resistant epilepsy.

A systems approach delivers a functional microRNA catalog and expanded targets for seizure suppression in temporal lobe epilepsy.

Temporal lobe epilepsy is the most common drug-resistant form of epilepsy in adults. The reorganization of neural networks and the gene expression landscape underlying pathophysiologic network behavior in brain structures such as the hippocampus has been suggested to be controlled, in part, by microRNAs. To systematically assess their significance, we sequenced Argonaute-loaded microRNAs to define functionally engaged microRNAs in the hippocampus of three different animal models in two species and at six time points between the initial precipitating insult through to the establishment of chronic epilepsy. We then selected commonly up-regulated microRNAs for a functional in vivo therapeutic screen using oligonucleotide inhibitors. Argonaute sequencing generated 1.44 billion small RNA reads of which up to 82% were microRNAs, with over 400 unique microRNAs detected per model. Approximately half of the detected microRNAs were dysregulated in each epilepsy model. We prioritized commonly up-regulated microRNAs that were fully conserved in humans and designed custom antisense oligonucleotides for these candidate targets. Antiseizure phenotypes were observed upon knockdown of miR-10a-5p, miR-21a-5p, and miR-142a-5p and electrophysiological analyses indicated broad safety of this approach. Combined inhibition of these three microRNAs reduced spontaneous seizures in epileptic mice. Proteomic data, RNA sequencing, and pathway analysis on predicted and validated targets of these microRNAs implicated derepressed TGF-ß signaling as a shared seizure-modifying mechanism. Correspondingly, inhibition of TGF-ß signaling occluded the antiseizure effects of the antagomirs. Together, these results identify shared, dysregulated, and functionally active microRNAs during the pathogenesis of epilepsy which represent therapeutic antiseizure targets.

tsRNAsearch: a pipeline for the identification of tRNA and ncRNA fragments from small RNA-sequencing data.

MOTIVATION: tRNAs were originally considered uni-functional RNA molecules involved in the delivery of amino acids to growing peptide chains on the ribosome. More recently, the liberation of tRNA fragments from tRNAs via specific enzyme cleavage has been characterized. Detection of tRNA fragments in sequencing data is difficult due to tRNA sequence redundancy and the short length of both tRNAs and their fragments. RESULTS: Here, we introduce tsRNAsearch, a Nextflow pipeline for the identification of differentially abundant tRNA fragments and other non-coding RNAs from small RNA-sequencing data. tsRNAsearch is intended for use when comparing two groups of datasets, such as control and treatment groups. tsRNAsearch comparatively searches for tRNAs and ncRNAs with irregular read distribution profiles (a proxy for RNA cleavage) using a combined score made up of four novel methods and a differential expression analysis, and reports the top ranked results in simple PDF and TEXT files. In this study, we used publicly available small RNA-seq data to replicate the identification of tsRNAs from chronic hepatitis-infected liver tissue data. In addition, we applied tsRNAsearch to pancreatic ductal adenocarcinoma (PDAC) and matched healthy pancreatic tissue small RNA-sequencing data. Our results support the identification of miR135b from the original study as a potential biomarker of PDAC and identify other potentially stronger miRNA biomarkers of PDAC. AVAILABILITY AND IMPLEMENTATION: https://github.com/GiantSpaceRobot/tsRNAsearch. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nanopore long-read sequencing of circRNAs.

Circular RNA (circRNA) is a group of highly stable RNA molecules with suggested roles in development and disease. They derive from linear pre-mRNAs when a 5'-splice site splices back to an upstream 3'-splice site in a process termed back-splicing. Most circRNAs are multi-exonic and may contain several thousand nucleotides. The extensive sequence overlap between the linear and circular forms of an RNA means that circRNA identification depends on the detection of back-splice-junction sequence reads that are unique to the circRNA. However, the short-read length obtained using standard next-generation sequencing techniques means that the internal sequence, exon composition and alternative splicing of circRNAs are unknown in many cases. Recently, several labs, including ours, have reported protocols for sequencing of circRNAs using long-read nanopore sequencing and thereby expanded our understanding of circRNA size distribution and internal splicing patterns. Here, we review these protocols and discuss the different approaches taken to study the full length composition of circRNAs.

A microRNA-129-5p/Rbfox crosstalk coordinates homeostatic downscaling of excitatory synapses.

Synaptic downscaling is a homeostatic mechanism that allows neurons to reduce firing rates during chronically elevated network activity. Although synaptic downscaling is important in neural circuit development and epilepsy, the underlying mechanisms are poorly described. We performed small RNA profiling in picrotoxin (PTX)-treated hippocampal neurons, a model of synaptic downscaling. Thereby, we identified eight microRNAs (miRNAs) that were increased in response to PTX, including miR-129-5p, whose inhibition blocked synaptic downscaling in vitro and reduced epileptic seizure severity in vivo Using transcriptome, proteome, and bioinformatic analysis, we identified the calcium pump Atp2b4 and doublecortin (Dcx) as miR-129-5p targets. Restoring Atp2b4 and Dcx expression was sufficient to prevent synaptic downscaling in PTX-treated neurons. Furthermore, we characterized a functional crosstalk between miR-129-5p and the RNA-binding protein (RBP) Rbfox1. In the absence of PTX, Rbfox1 promoted the expression of Atp2b4 and Dcx. Upon PTX treatment, Rbfox1 expression was downregulated by miR-129-5p, thereby allowing the repression of Atp2b4 and Dcx. We therefore identified a novel activity-dependent miRNA/RBP crosstalk during synaptic scaling, with potential implications for neural network homeostasis and epileptogenesis.

SMARTer single cell total RNA sequencing.

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.

Global MicroRNA Profiling in Human Bone Marrow Skeletal-Stromal or Mesenchymal-Stem Cells Identified Candidates for Bone Regeneration.

Bone remodeling and regeneration are highly regulated multistep processes involving posttranscriptional regulation by microRNAs (miRNAs). Here, we performed a global profiling of differentially expressed miRNAs in bone-marrow-derived skeletal cells (BMSCs; also known as stromal or mesenchymal stem cells) during in vitro osteoblast differentiation. We functionally validated the regulatory effects of several miRNAs on osteoblast differentiation and identified 15 miRNAs, most significantly miR-222 and miR-423, as regulators of osteoblastogenesis. In addition, we tested the possible targeting of miRNAs for enhancing bone tissue regeneration. Scaffolds functionalized with miRNA nano-carriers enhanced osteoblastogenesis in 3D culture and retained this ability at least 2 weeks after storage. Additionally, anti-miR-222 enhanced in vivo ectopic bone formation through targeting the cell-cycle inhibitor CDKN1B (cyclin-dependent kinase inhibitor 1B). A number of additional miRNAs exerted additive osteoinductive effects on BMSC differentiation, suggesting that pools of miRNAs delivered locally from an implanted scaffold can provide a promising approach for enhanced bone regeneration.

Comparison of circular RNA prediction tools.

CircRNAs are novel members of the non-coding RNA family. For several decades circRNAs have been known to exist, however only recently the widespread abundance has become appreciated. Annotation of circRNAs depends on sequencing reads spanning the backsplice junction and therefore map as non-linear reads in the genome. Several pipelines have been developed to specifically identify these non-linear reads and consequently predict the landscape of circRNAs based on deep sequencing datasets. Here, we use common RNAseq datasets to scrutinize and compare the output from five different algorithms; circRNA_finder, find_circ, CIRCexplorer, CIRI, and MapSplice and evaluate the levels of bona fide and false positive circRNAs based on RNase R resistance. By this approach, we observe surprisingly dramatic differences between the algorithms specifically regarding the highly expressed circRNAs and the circRNAs derived from proximal splice sites. Collectively, this study emphasizes that circRNA annotation should be handled with care and that several algorithms should ideally be combined to achieve reliable predictions.

miRdentify: high stringency miRNA predictor identifies several novel animal miRNAs.

During recent years, miRNAs have been shown to play important roles in the regulation of gene expression. Accordingly, much effort has been put into the discovery of novel uncharacterized miRNAs in various organisms. miRNAs are structurally defined by a hairpin-loop structure recognized by the two-step processing apparatus, Drosha and Dicer, necessary for the production of mature ∼ 22-nucleotide miRNA guide strands. With the emergence of high-throughput sequencing applications, tools have been developed to identify miRNAs and profile their expression based on sequencing reads. However, as the read depth increases, false-positive predictions increase using established algorithms, underscoring the need for more stringent approaches. Here we describe a transparent pipeline for confident miRNA identification in animals, termed miRdentify. We show that miRdentify confidently discloses more than 400 novel miRNAs in humans, including the first male-specific miRNA, which we successfully validate. Moreover, novel miRNAs are predicted in the mouse, the fruit fly and nematodes, suggesting that the pipeline applies to all animals. The entire software package is available at www.ncrnalab.dk/mirdentify.

Aberrant expression of miR-218 and miR-204 in human mesial temporal lobe epilepsy and hippocampal sclerosis-convergence on axonal guidance.

OBJECTIVE: Mesial temporal lobe epilepsy (MTLE) is one of the most common types of the intractable epilepsies and is most often associated with hippocampal sclerosis (HS), which is characterized by pronounced loss of hippocampal pyramidal neurons. microRNAs (miRNAs) have been shown to be dysregulated in epilepsy and neurodegenerative diseases, and we hypothesized that miRNAs could be involved in the pathogenesis of MTLE and HS. METHODS: miRNA expression was quantified in hippocampal specimens from human patients using miRNA microarray and quantitative real-time polymerase chain reaction RT-PCR, and by RNA-seq on fetal brain specimens from domestic pigs. In situ hybridization was used to show the spatial distribution of miRNAs in the human hippocampus. The potential effect of miRNAs on targets genes was investigated using the dual luciferase reporter gene assay. RESULTS: miRNA expression profiling showed that 25 miRNAs were up-regulated and 5 were down-regulated in hippocampus biopsies of MTLE/HS patients compared to controls. We showed that miR-204 and miR-218 were significantly down-regulated in MTLE and HS, and both were expressed in neurons in all subfields of normal hippocampus. Moreover, miR-204 and miR-218 showed strong changes in expression during fetal development of the hippocampus in pigs, and we identified four target genes, involved in axonal guidance and synaptic plasticity, ROBO1, GRM1, SLC1A2, and GNAI2, as bona fide targets of miR-218. GRM1 was also shown to be a direct target of miR-204. SIGNIFICANCE: miR-204 and miR-218 are developmentally regulated in the hippocampus and may contribute to the molecular mechanisms underlying the pathogenesis of MTLE and HS.

A distant cis acting intronic element induces site-selective RNA editing.

Transcripts have been found to be site selectively edited from adenosine-to-inosine (A-to-I) in the mammalian brain, mostly in genes involved in neurotransmission. While A-to-I editing occurs at double-stranded structures, other structural requirements are largely unknown. We have investigated the requirements for editing at the I/M site in the Gabra-3 transcript of the GABA(A) receptor. We identify an evolutionarily conserved intronic duplex, 150 nt downstream of the exonic hairpin where the I/M site resides, which is required for its editing. This is the first time a distant RNA structure has been shown to be important for A-to-I editing. We demonstrate that the element also can induce editing in related but normally not edited RNA sequences. In human, thousands of genes are edited in duplexes formed by inverted repeats in non-coding regions. It is likely that numerous such duplexes can induce editing of coding regions throughout the transcriptome.

Integrative network analysis of miRNA-mRNA expression profiles during epileptogenesis in rats reveals therapeutic targets after emergence of first spontaneous seizure.

Epileptogenesis is the process by which a normal brain becomes hyperexcitable and capable of generating spontaneous recurrent seizures. The extensive dysregulation of gene expression associated with epileptogenesis is shaped, in part, by microRNAs (miRNAs) - short, non-coding RNAs that negatively regulate protein levels. Functional miRNA-mediated regulation can, however, be difficult to elucidate due to the complexity of miRNA-mRNA interactions. Here, we integrated miRNA and mRNA expression profiles sampled over multiple time-points during and after epileptogenesis in rats, and applied bi-clustering and Bayesian modelling to construct temporal miRNA-mRNA-mRNA interaction networks. Network analysis and enrichment of network inference with sequence- and human disease-specific information identified key regulatory miRNAs with the strongest influence on the mRNA landscape, and miRNA-mRNA interactions closely associated with epileptogenesis and subsequent epilepsy. Our findings underscore the complexity of miRNA-mRNA regulation, can be used to prioritise miRNA targets in specific systems, and offer insights into key regulatory processes in epileptogenesis with therapeutic potential for further investigation.

Differential microRNA editing may drive target pathway switching in human temporal lobe epilepsy.

MicroRNAs have emerged as important regulators of the gene expression landscape in temporal lobe epilepsy. The mechanisms that control microRNA levels and influence target choice remain, however, poorly understood. RNA editing is a post-transcriptional mechanism mediated by the adenosine acting on RNA (ADAR) family of proteins that introduces base modification that diversifies the gene expression landscape. RNA editing has been studied for the mRNA landscape but the extent to which microRNA editing occurs in human temporal lobe epilepsy is unknown. Here, we used small RNA-sequencing data to characterize the identity and extent of microRNA editing in human temporal lobe epilepsy brain samples. This detected low-to-high editing in over 40 of the identified microRNAs. Among microRNA exhibiting the highest editing was miR-376a-3p, which was edited in the seed region and this was predicted to significantly change the target pool. The edited form was expressed at lower levels in human temporal lobe epilepsy samples. We modelled the shift in editing levels of miR-376a-3p in human-induced pluripotent stem cell-derived neurons. Reducing levels of the edited form of miR-376a-3p using antisense oligonucleotides resulted in extensive gene expression changes, including upregulation of mitochondrial and metabolism-associated pathways. Together, these results show that differential editing of microRNAs may re-direct targeting and result in altered functions relevant to the pathophysiology of temporal lobe epilepsy and perhaps other disorders of neuronal hyperexcitability.

Unlocking the diagnostic power of plasma extracellular vesicle miR-200 family in pancreatic ductal adenocarcinoma.

BACKGROUND: Distinguishing benign from malignant pancreaticobiliary disease is challenging because of the absence of reliable biomarkers. Circulating extracellular vesicles (EVs) have emerged as functional mediators between cells. Their cargos, including microRNAs (miRNAs), are increasingly acknowledged as an important source of potential biomarkers. This multicentric, prospective study aimed to establish a diagnostic plasma EV-derived miRNA signature to discriminate pancreatic ductal adenocarcinoma (PDAC) from benign pancreaticobiliary disease. METHODS: Plasma EVs were isolated using size exclusion chromatography (SEC) and characterised using nanoparticle tracking analysis, electron microscopy and Western blotting. EV-RNAs underwent small RNA sequencing to discover differentially expressed markers for PDAC (n = 10 benign vs. 10 PDAC). Candidate EV-miRNAs were then validated in a cohort of 61 patients (n = 31 benign vs. 30 PDAC) by RT-qPCR. Logistic regression and optimal thresholds (Youden Index) were used to develop an EV-miR-200 family model to detect cancer. This model was tested in an independent cohort of 95 patients (n = 30 benign, 33 PDAC, and 32 cholangiocarcinoma). RESULTS: Small RNA sequencing and RT-qPCR showed that EV-miR-200 family members were significantly overexpressed in PDAC vs. benign disease. Combined expression of the EV-miR-200 family showed an AUC of 0.823. In an independent validation cohort, application of this model showed a sensitivity, specificity and AUC of 100%, 88%, and 0.97, respectively, for diagnosing PDAC. CONCLUSIONS: This is the first study to validate plasma EV-miR-200 members as a clinically-useful diagnostic biomarker for PDAC. Further validation in larger cohorts and clinical trials is essential. These findings also suggest the potential utility in monitoring response and/or recurrence.

ciRS-7 and miR-7 regulate ischemia-induced neuronal death via glutamatergic signaling.

Brain functionality relies on finely tuned regulation of gene expression by networks of non-coding RNAs (ncRNAs) such as the one composed by the circular RNA ciRS-7 (also known as CDR1as), the microRNA miR-7, and the long ncRNA Cyrano. We describe ischemia-induced alterations in the ncRNA network both in vitro and in vivo and in transgenic mice lacking ciRS-7 or miR-7. Our data show that cortical neurons downregulate ciRS-7 and Cyrano and upregulate miR-7 expression during ischemia. Mice lacking ciRS-7 exhibit reduced lesion size and motor impairment, while the absence of miR-7 alone results in increased ischemia-induced neuronal death. Moreover, miR-7 levels in pyramidal excitatory neurons regulate neurite morphology and glutamatergic signaling, suggesting a potential molecular link to the in vivo phenotype. Our data reveal the role of ciRS-7 and miR-7 in modulating ischemic stroke outcome, shedding light on the pathophysiological function of intracellular ncRNA networks in the brain.

Knockdown of Circular RNAs Using LNA-Modified Antisense Oligonucleotides.

Circular RNAs (circRNAs) constitute an abundant class of covalently closed noncoding RNA molecules that are formed by backsplicing from eukaryotic protein-coding genes. Recent studies have shown that circRNAs can act as microRNA or protein decoys, as well as transcriptional regulators. However, the functions of most circRNAs are still poorly understood. Because circRNA sequences overlap with their linear parent transcripts, depleting specific circRNAs without affecting host gene expression remains a challenge. In this study, we assessed the utility of LNA-modified antisense oligonucleotides (ASOs) to knock down circRNAs for loss-of-function studies. We found that, while most RNase H-dependent gapmer ASOs mediate effective knockdown of their target circRNAs, some gapmers reduce the levels of the linear parent transcript. The circRNA targeting specificity can be enhanced using design-optimized gapmer ASOs, which display potent and specific circRNA knockdown with a minimal effect on the host genes. In summary, our results demonstrate that LNA-modified ASOs complementary to backsplice-junction sequences mediate robust knockdown of circRNAs in vitro and, thus, represent a useful tool to explore the biological roles of circRNAs in loss-of-function studies in cultured cells and animal models.

Brain cell-specific origin of circulating microRNA biomarkers in experimental temporal lobe epilepsy.

The diagnosis of epilepsy is complex and challenging and would benefit from the availability of molecular biomarkers, ideally measurable in a biofluid such as blood. Experimental and human epilepsy are associated with altered brain and blood levels of various microRNAs (miRNAs). Evidence is lacking, however, as to whether any of the circulating pool of miRNAs originates from the brain. To explore the link between circulating miRNAs and the pathophysiology of epilepsy, we first sequenced argonaute 2 (Ago2)-bound miRNAs in plasma samples collected from mice subject to status epilepticus induced by intraamygdala microinjection of kainic acid. This identified time-dependent changes in plasma levels of miRNAs with known neuronal and microglial-cell origins. To explore whether the circulating miRNAs had originated from the brain, we generated mice expressing FLAG-Ago2 in neurons or microglia using tamoxifen-inducible Thy1 or Cx3cr1 promoters, respectively. FLAG immunoprecipitates from the plasma of these mice after seizures contained miRNAs, including let-7i-5p and miR-19b-3p. Taken together, these studies confirm that a portion of the circulating pool of miRNAs in experimental epilepsy originates from the brain, increasing support for miRNAs as mechanistic biomarkers of epilepsy.

Spatio-temporal regulation of ADAR editing during development in porcine neural tissues.

Editing by ADAR enzymes is essential for mammalian life. Still, knowledge of the spatio-temporal editing patterns in mammals is limited. By use of 454 amplicon sequencing we examined the editing status of 12 regionally extracted mRNAs from porcine developing brain encompassing a total of 64 putative ADAR editing sites. In total 24 brain tissues, dissected from up to five regions from embryonic gestation day 23, 42, 60, 80, 100 and 115, were examined for editing.   Generally, editing increased during embryonic development concomitantly with an increase in ADAR2 mRNA level. Notably, the Gria2 (GluR-B) Q/R site, reported to be ~100% edited in previous studies, is only 54% edited at embryonic day 23. Transcripts with multiple editing sites in close proximity to each other exhibit coupled editing and an extraordinary incidence of long-range coupling of editing events more than 32 kb apart is observed for the kainate glutamate receptor 2 transcript, Grik2. Our study reveals complex spatio-temporal ADAR editing patterns of coordinated editing events that may play important roles in the development of the mammalian brain.

Expression of Circ_Satb1 Is Decreased in Mesial Temporal Lobe Epilepsy and Regulates Dendritic Spine Morphology.

Mesial temporal lobe epilepsy (mTLE) is a chronic disease characterized by recurrent seizures that originate in the temporal lobes of the brain. Anti-epileptic drugs (AEDs) are the standard treatment for managing seizures in mTLE patients, but are frequently ineffective. Resective surgery is an option for some patients, but does not guarantee a postoperative seizure-free period. Therefore, further insight is needed into the pathogenesis of mTLE to enable the design of new therapeutic strategies. Circular RNAs (circRNAs) have been identified as important regulators of neuronal function and have been implicated in epilepsy. However, the mechanisms through which circRNAs contribute to epileptogenesis remain unknown. Here, we determine the circRNA transcriptome of the hippocampus and cortex of mTLE patients by using RNA-seq. We report 333 differentially expressed (DE) circRNAs between healthy individuals and mTLE patients, of which 23 circRNAs displayed significant adjusted p-values following multiple testing correction. Interestingly, hippocampal expression of circ_Satb1, a circRNA derived from special AT-rich sequence binding protein 1 (SATB1), is decreased in both mTLE patients and in experimental epilepsy. Our work shows that circ_Satb1 displays dynamic patterns of neuronal expression in vitro and in vivo. Further, circ_Satb1-specific knockdown using CRISPR/CasRx approaches in hippocampal cultures leads to defects in dendritic spine morphology, a cellular hallmark of mTLE. Overall, our results identify a novel epilepsy-associated circRNA with disease-specific expression and previously unidentified cellular effects that are relevant for epileptogenesis.

Nanopore sequencing of brain-derived full-length circRNAs reveals circRNA-specific exon usage, intron retention and microexons.

Circular RNA (circRNA) is a class of covalently joined non-coding RNAs with functional roles in a wide variety of cellular processes. Their composition shows extensive overlap with exons found in linear mRNAs making it difficult to delineate their composition using short-read RNA sequencing, particularly for long and multi-exonic circRNAs. Here, we use long-read nanopore sequencing of nicked circRNAs (circNick-LRS) and characterize a total of 18,266 and 39,623 circRNAs in human and mouse brain, respectively. We further develop an approach for targeted long-read sequencing of a panel of circRNAs (circPanel-LRS), eliminating the need for prior circRNA enrichment and find >30 circRNA isoforms on average per targeted locus. Our data show that circRNAs exhibit a large number of splicing events such as novel exons, intron retention and microexons that preferentially occur in circRNAs. We propose that altered exon usage in circRNAs may reflect resistance to nonsense-mediated decay in the absence of translation.