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.

The brain reacting to COVID-19: analysis of the cerebrospinal fluid proteome, RNA and inflammation.

Patients with COVID-19 can have a variety of neurological symptoms, but the active involvement of central nervous system (CNS) in COVID-19 remains unclear. While routine cerebrospinal fluid (CSF) analyses in patients with neurological manifestations of COVID-19 generally show no or only mild inflammation, more detailed data on inflammatory mediators in the CSF of patients with COVID-19 are scarce. We studied the inflammatory response in paired CSF and serum samples of patients with COVID-19 (n = 38). Patients with herpes simplex virus encephalitis (HSVE, n = 10) and patients with non-inflammatory, non-neurodegenerative neurological diseases (n = 28) served as controls. We used proteomics, enzyme-linked immunoassays, and semiquantitative cytokine arrays to characterize inflammatory proteins. Autoantibody screening was performed with cell-based assays and native tissue staining. RNA sequencing of long-non-coding RNA and circular RNA was done to study the transcriptome. Proteomics on single protein level and subsequent pathway analysis showed similar yet strongly attenuated inflammatory changes in the CSF of COVID-19 patients compared to HSVE patients with, e.g., downregulation of the apolipoproteins and extracellular matrix proteins. Protein upregulation of the complement system, the serpin proteins pathways, and other proteins including glycoproteins alpha-2 and alpha-1 acid. Importantly, calculation of interleukin-6, interleukin-16, and CXCL10 CSF/serum indices suggest that these inflammatory mediators reach the CSF from the systemic circulation, rather than being produced within the CNS. Antibody screening revealed no pathological levels of known neuronal autoantibodies. When stratifying COVID-19 patients into those with and without bacterial superinfection as indicated by elevated procalcitonin levels, inflammatory markers were significantly (p < 0.01) higher in those with bacterial superinfection. RNA sequencing in the CSF revealed 101 linear RNAs comprising messenger RNAs, and two circRNAs being significantly differentially expressed in COVID-19 than in non-neuroinflammatory controls and neurodegenerative patients. Our findings may explain the absence of signs of intrathecal inflammation upon routine CSF testing despite the presence of SARS-CoV2 infection-associated neurological symptoms. The relevance of blood-derived mediators of inflammation in the CSF for neurological COVID-19 and post-COVID-19 symptoms deserves further investigation.

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.

Circular RNAs are abundantly expressed and upregulated during human epidermal stem cell differentiation.

The expression patterns of endogenous circular RNA (circRNA) molecules during epidermal stem cell (EpSC) differentiation have not previously been explored. Here, we show that circRNAs are abundantly expressed in EpSCs and that their expression change dramatically during differentiation in a coordinated manner. Overall, circRNAs are expressed at higher levels in the differentiated cells, and many upregulated circRNAs are derived from developmental genes, including four different circRNAs from DLG1. The observed changes in circRNA expression were largely independent of host gene expression, and circRNAs independently upregulated upon differentiation are more prone to AGO2 binding and have more predicted miRNA binding sites compared to stably expressed circRNAs. In particular, upregulated circRNAs from the HECTD1 and ZNF91 genes have exceptionally high numbers of AGO2 binding sites and predicted miRNA target sites, and circZNF91 contains 24 target sites for miR-23b-3p, which is known to play important roles in keratinocyte differentiation. We also observed that upregulated circRNAs are less likely to be flanked by homologues inverted Alu repeats compared to stably expressed circRNAs. This coincide with DHX9 being upregulated in the differentiated keratinocytes. Finally, none of the circRNAs upregulated upon differentiation were also upregulated upon DNMT3A or DNMT3B knockdown, making it unlikely that epigenetic mechanisms are governing the observed circRNA expression changes. Together, we provide a map of circRNA expression in EpSCs and their differentiated counterparts and shed light on potential function and regulation of differentially expressed circRNAs.

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.

Nanopore-Mediated Sequencing of Circular RNA.

Circular RNAs (circRNAs) constitute a group of RNAs defined by a covalent bond between the 5' and 3' end formed by a unique back-splicing event. Most circRNAs are composed of more than one exon, which are spliced together in a linear fashion. This protocol describes methods to sequence full-length circRNA across the back-splicing junction, allowing unambiguous characterization of circRNA-specific exon-intron structures by long-read sequencing (LRS). Two different sequencing approaches are provided: (1) Global circRNA sequencing (the circNick-LRS strategy) relying on circRNA enrichment from total RNA followed by total circRNA long-read sequencing, and (2) targeted circRNA sequencing (the circPanel-LRS strategy) where a preselected panel of circRNA are sequenced without prior circRNA enrichment. Both methods were originally described in Karim et al. (Rahimi et al., Nat Commun 12: 4825, 2021) where they were applied to characterize the exon-intron structure of >10.000 circRNAs in mouse and human brains.

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.

A Circular RNA Expressed from the FAT3 Locus Regulates Neural Development.

Circular RNAs (circRNAs) are key regulators of cellular processes, are abundant in the nervous system, and have putative regulatory roles during neural differentiation. However, the knowledge about circRNA functions in brain development is limited. Here, using RNA-sequencing, we show that circRNA levels increased substantially over the course of differentiation of human embryonic stem cells into rostral and caudal neural progenitor cells (NPCs), including three of the most abundant circRNAs, ciRS-7, circRMST, and circFAT3. Knockdown of circFAT3 during early neural differentiation resulted in minor transcriptional alterations in bulk RNA analysis. However, single-cell transcriptomics of 30 and 90 days differentiated cerebral organoids deficient in circFAT3 showed a loss of telencephalic radial glial cells and mature cortical neurons, respectively. Furthermore, non-telencephalic NPCs in cerebral organoids showed changes in the expression of genes involved in neural differentiation and migration, including FAT4, ERBB4, UNC5C, and DCC. In vivo depletion of circFat3 in mouse prefrontal cortex using in utero electroporation led to alterations in the positioning of the electroporated cells within the neocortex. Overall, these findings suggest a conserved role for circFAT3 in neural development involving the formation of anterior cell types, neuronal differentiation, or migration.