Best practice standards for circular RNA research.

Circular RNAs (circRNAs) are formed in all domains of life and via different mechanisms. There has been an explosion in the number of circRNA papers in recent years; however, as a relatively young field, circRNA biology has an urgent need for common experimental standards for isolating, analyzing, expressing and depleting circRNAs. Here we propose a set of guidelines for circRNA studies based on the authors' experience. This Perspective will specifically address the major class of circRNAs in Eukarya that are generated by a spliceosome-catalyzed back-splicing event. We hope that the implementation of best practice principles for circRNA research will help move the field forward and allow a better functional understanding of this fascinating group of RNAs.

The biogenesis, biology and characterization of circular RNAs.

Circular RNAs (circRNAs) are covalently closed, endogenous biomolecules in eukaryotes with tissue-specific and cell-specific expression patterns, whose biogenesis is regulated by specific cis-acting elements and trans-acting factors. Some circRNAs are abundant and evolutionarily conserved, and many circRNAs exert important biological functions by acting as microRNA or protein inhibitors ('sponges'), by regulating protein function or by being translated themselves. Furthermore, circRNAs have been implicated in diseases such as diabetes mellitus, neurological disorders, cardiovascular diseases and cancer. Although the circular nature of these transcripts makes their detection, quantification and functional characterization challenging, recent advances in high-throughput RNA sequencing and circRNA-specific computational tools have driven the development of state-of-the-art approaches for their identification, and novel approaches to functional characterization are emerging.

circRAB3IP modulates cell proliferation by reorganizing gene expression and mRNA processing in a paracrine manner.

Circular RNAs are an endogenous long-lived and abundant noncoding species. Despite their prevalence, only a few circRNAs have been dissected mechanistically to date. Here, we cataloged nascent RNA-enriched circRNAs from primary human cells and functionally assigned a role to circRAB3IP in sustaining cellular homeostasis. We combined "omics" and functional experiments to show how circRAB3IP depletion deregulates hundreds of genes, suppresses cell cycle progression, and induces senescence-associated gene expression changes. Conversely, excess circRAB3IP delivered to endothelial cells via extracellular vesicles suffices for accelerating their division. We attribute these effects to an interplay between circRAB3IP and the general splicing factor SF3B1, which can affect transcript variant expression levels of cell cycle-related genes. Together, our findings link the maintenance of cell homeostasis to the presence of a single circRNA.

Alzheimer's disease-induced phagocytic microglia express a specific profile of coding and non-coding RNAs.

INTRODUCTION: Alzheimer's disease (AD) is a neurodegenerative disease and the main cause of dementia in the elderly. AD pathology is characterized by accumulation of microglia around the beta-amyloid (Aß) plaques which assumes disease-specific transcriptional signatures, as for the disease-associated microglia (DAM). However, the regulators of microglial phagocytosis are still unknown. METHODS: We isolated Aß-laden microglia from the brain of 5xFAD mice for RNA sequencing to characterize the transcriptional signature in phagocytic microglia and to identify the key non-coding RNAs capable of regulating microglial phagocytosis. Through spatial sequencing, we show the transcriptional changes of microglia in the AD mouse brain in relation to Aß proximity. RESULTS: Finally, we show that phagocytic messenger RNAs are regulated by miR-7a-5p, miR-29a-3p and miR-146a-5p microRNAs and segregate the DAM population into phagocytic and non-phagocytic states. DISCUSSION: Our study pinpoints key regulators of microglial Aß clearing capacity suggesting new targets for future therapeutic approaches.

RNA-Seq profiling of leukocytes reveals a sex-dependent global circular RNA upregulation in multiple sclerosis and 6 candidate biomarkers.

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system, with higher prevalence in women, that leads to neurological disability. The disease course and clinical phenotype are highly variable, and therefore, biomarkers for the diagnosis, classification, monitoring of the disease and treatment assessment are needed. Studies have shown a dysregulation in the coding and non-coding RNAs and proposed some as biomarkers. However, still none of them have reached the clinical practice. Recently, circular RNAs (circRNAs) have emerged as new players in the transcriptome that hold a great potential as biomarkers in several diseases. Leukocytes from 30 MS patients and 20 healthy controls (HCs) were RNA-sequenced to study the linear and circular transcriptome. Differential expression analysis was performed by DESeq, and circRNA candidates were studied in a second cohort (70 MS and 46 HC) by RT-qPCR and in paired samples drawn during the relapse and remission phases (20 patients). Among the differentially expressed circRNAs, 96.1% are upregulated in patients compared with controls, but similar circRNA profiles are found between MS types. The same upregulation trend was observed in females but not in males or in the linear transcriptome. The upregulation of 6 circRNAs was validated, and a change in their expression was found between relapse and remission. The 6 circRNAs showed a good performance to discriminate patients from HC with a combined area under the curve of 0.852. There is global, specific and sex-dependent increase of circRNA expression in MS, and 6 circRNAs are proposed as potential biomarkers.

Natural RNA circles function as efficient microRNA sponges.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that act by direct base pairing to target sites within untranslated regions of messenger RNAs. Recently, miRNA activity has been shown to be affected by the presence of miRNA sponge transcripts, the so-called competing endogenous RNA in humans and target mimicry in plants. We previously identified a highly expressed circular RNA (circRNA) in human and mouse brain. Here we show that this circRNA acts as a miR-7 sponge; we term this circular transcript ciRS-7 (circular RNA sponge for miR-7). ciRS-7 contains more than 70 selectively conserved miRNA target sites, and it is highly and widely associated with Argonaute (AGO) proteins in a miR-7-dependent manner. Although the circRNA is completely resistant to miRNA-mediated target destabilization, it strongly suppresses miR-7 activity, resulting in increased levels of miR-7 targets. In the mouse brain, we observe overlapping co-expression of ciRS-7 and miR-7, particularly in neocortical and hippocampal neurons, suggesting a high degree of endogenous interaction. We further show that the testis-specific circRNA, sex-determining region Y (Sry), serves as a miR-138 sponge, suggesting that miRNA sponge effects achieved by circRNA formation are a general phenomenon. This study serves as the first, to our knowledge, functional analysis of a naturally expressed circRNA.

The agotrons: Gene regulators or Argonaute protectors?

Over the last decades, it has become evident that highly complex networks of regulators govern post-transcriptional regulation of gene expression. A novel class of Argonaute (Ago)-associated RNA molecules, the agotrons, was recently shown to function in a Drosha- and Dicer-independent manner, hence bypassing the maturation steps required for canonical microRNA (miRNA) biogenesis. Agotrons are found in most mammals and associate with Ago as ∼100 nucleotide (nt) long RNA species. Here, we speculate on the functional and biological relevance of agotrons: (i) agotrons could serve as non-promiscuous miRNA-like regulators with reduced off-targeting or (ii) agotrons could encompass other putative functions, such as protecting Ago proteins from taking up aberrant short RNAs or by rescuing and stabilizing otherwise unloaded Ago-proteins from degradation. Collectively, agotrons have emerged as a novel class of interesting non-coding RNA molecules, but their full functional potential and biological impact still remain to be disclosed.

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.

Circular RNAs: Identification, biogenesis and function.

Circular RNAs are a novel class of non-coding RNA characterized by the presence of a covalent bond linking the 3' and 5' ends generated by backsplicing. Circular RNAs are widely expressed in a tissue and developmental-stage specific pattern and a subset displays conservation across species. Functional circRNAs have been shown to act as cytoplasmic microRNA sponges and RNA-binding protein sequestering agents as well as nuclear transcriptional regulators, illustrating the relevance of circular RNAs as participants in the regulatory networks governing gene expression. Here, we review the features that characterize circular RNAs, discuss putative circular RNA biogenesis pathways as well as review the uncovered functions of circular RNAs. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Insights into circular RNA biology.

Circular RNAs (circRNAs) are a novel class of non-coding RNA characterized by a covalently closed-loop structure generated through a special type of alternative splicing termed backsplicing. CircRNAs are emerging as a heterogeneous class of molecules involved in modulating gene expression by regulation of transcription, protein and miRNA functions. CircRNA expression is cell type and tissue specific and can be largely independent of the expression level of the linear host gene, indicating that regulation of expression might be an important aspect with regard to control of circRNA function. In this review, a brief introduction to the characteristics that define a circRNA will be given followed by a discussion of putative biogenesis pathways and modulators of circRNA expression as well as of the stage at which circRNA formation takes place. A brief summary of circRNA functions will also be provided and lastly, an outlook with a focus on unanswered questions regarding circRNA biology will be included.

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.

miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA.

MicroRNAs (miRNAs) are ∼22 nt non-coding RNAs that typically bind to the 3' UTR of target mRNAs in the cytoplasm, resulting in mRNA destabilization and translational repression. Here, we report that miRNAs can also regulate gene expression by targeting non-coding antisense transcripts in human cells. Specifically, we show that miR-671 directs cleavage of a circular antisense transcript of the Cerebellar Degeneration-Related protein 1 (CDR1) locus in an Ago2-slicer-dependent manner. The resulting downregulation of circular antisense has a concomitant decrease in CDR1 mRNA levels, independently of heterochromatin formation. This study provides the first evidence for non-coding antisense transcripts as functional miRNA targets, and a novel regulatory mechanism involving a positive correlation between mRNA and antisense circular RNA levels.

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.

Pilot Study of ONCOS-102 and Pembrolizumab: Remodeling of the Tumor Microenvironment and Clinical Outcomes in Anti-PD-1-Resistant Advanced Melanoma.

PURPOSE: Intratumoral oncolytic virotherapy may overcome anti-PD(L)-1 resistance by triggering pro-inflammatory remodeling of the tumor microenvironment. This pilot study investigated ONCOS-102 (oncolytic adenovirus expressing GM-CSF) plus anti-programmed cell death protein 1 (PD)-1 therapy in anti-PD-1-resistant melanoma. PATIENTS AND METHODS: Patients with advanced melanoma progressing after prior PD-1 blockade received intratumoral ONCOS-102 either as priming with 3 doses (3 × 1011 viral particles) during Week 1 [Part 1 (sequential treatment)] or as 4-dose priming and 8 booster doses every 3 weeks [Part 2 (combination treatment)]. From Week 3, all patients received pembrolizumab every 3 weeks (≤8 doses). The primary endpoint was safety. Objective response rate (ORR), progression-free survival, and immunologic activation in repeat biopsies were also investigated. RESULTS: In 21 patients (Part 1, n = 9; Part 2, n = 12) ONCOS-102 plus pembrolizumab was well tolerated: most adverse events (AE) were mild/moderate in severity. Pyrexia (43%), chills (43%), and nausea (28%) were the most common ONCOS-102-related AEs. There were no dose-limiting toxicities. ORR was 35% [response evaluation in solid tumors (RECIST) 1.1, irRECIST]. Reduction in size of ≥1 non-injected lesions observed in 53% patients indicated a systemic effect. In injected tumors, persistent immune-related gene expression and T-cell infiltration were associated with clinical benefit. Viral persistence and efficacy in injected and non-injected lesions without additional toxicity supported Part 2 dosing regimen in future studies. CONCLUSIONS: ONCOS-102 plus pembrolizumab was well tolerated and led to objective responses in patients with anti-PD-1-resistant advanced melanoma. ONCOS-102 promoted T-cell infiltration, particularly cytotoxic CD8+ T cells, which persisted at Week 9, driving clinical benefit. Further investigation of ONCOS-102 plus PD-1 blockade is warranted. See related commentary by Levi and Boland, p. 3.

Circular stable intronic RNAs possess distinct biological features and are deregulated in bladder cancer.

Until recently, intronic lariats were regarded as short-lasting splicing byproducts with no apparent function; however, increasing evidence of stable derivatives suggests regulatory roles. Yet little is known about their characteristics, functions, distribution, and expression in healthy and tumor tissue. Here, we profiled and characterized circular stable intronic sequence RNAs (sisRNAs) using total RNA-Seq data from bladder cancer (BC; n = 457, UROMOL cohort), healthy tissue (n = 46), and fractionated cell lines (n = 5). We found that the recently-discovered full-length intronic circles and the stable lariats formed distinct subclasses, with a surprisingly high intronic circle fraction in BC (∼45%) compared to healthy tissues (0-20%). The stable lariats and their host introns were characterized by small transcript sizes, highly conserved BP regions, enriched BP motifs, and localization in multiple cell fractions. Additionally, circular sisRNAs showed tissue-specific expression patterns. We found nine circular sisRNAs as differentially expressed across early-stage BC patients with different prognoses, and sisHNRNPK expression correlated with progression-free survival. In conclusion, we identify distinguishing biological features of circular sisRNAs and point to specific candidates (incl. sisHNRNPK, sisWDR13 and sisMBNL1) that were highly expressed, had evolutionary conserved sequences, or had clinical correlations, which may facilitate future studies and further insights into their functional roles.

Impact of U2AF1 mutations on circular RNA expression in myelodysplastic neoplasms.

Mutations in U2AF1 are relatively common in myelodysplastic neoplasms (MDS) and are associated with an inferior prognosis, but the molecular mechanisms underlying this are not fully elucidated. Circular RNAs (circRNAs) have been implicated in cancer, but it is unknown how mutations in splicing factors may impact on circRNA biogenesis. Here, we used RNA-sequencing to investigate the effects of U2AF1 mutations on circRNA expression in K562 cells with a doxycycline-inducible U2AF1S34 mutation, in a mouse model with a doxycycline-inducible U2AF1S34 mutation, and in FACS-sorted CD34+ bone marrow cells from MDS patients with either U2AF1S34 or U2AF1Q157 mutations. In all contexts, we found an increase in global circRNA levels in the U2AF1-mutated setting, which was independent of expression changes in the cognate linear host genes. In patients, the U2AF1S34 and U2AF1Q157 mutations were both associated with an overall increased expression of circRNAs. circRNAs generated by a non-Alu-mediated mechanism generally showed the largest increase in expression levels. Several well-described cancer-associated circRNAs, including circZNF609 and circCSNK1G3, were upregulated in MDS patients with U2AF1 mutations compared to U2AF1-wildtype MDS controls. In conclusion, high circRNA expression is observed in association with U2AF1 mutations in three biological systems, presenting an interesting possibility for biomarker and therapeutic investigation.

A screen of chemical modifications identifies position-specific modification by UNA to most potently reduce siRNA off-target effects.

Small interfering RNAs (siRNAs) are now established as the preferred tool to inhibit gene function in mammalian cells yet trigger unintended gene silencing due to their inherent miRNA-like behavior. Such off-target effects are primarily mediated by the sequence-specific interaction between the siRNA seed regions (position 2-8 of either siRNA strand counting from the 5'-end) and complementary sequences in the 3'UTR of (off-) targets. It was previously shown that chemical modification of siRNAs can reduce off-targeting but only very few modifications have been tested leaving more to be identified. Here we developed a luciferase reporter-based assay suitable to monitor siRNA off-targeting in a high throughput manner using stable cell lines. We investigated the impact of chemically modifying single nucleotide positions within the siRNA seed on siRNA function and off-targeting using 10 different types of chemical modifications, three different target sequences and three siRNA concentrations. We found several differently modified siRNAs to exercise reduced off-targeting yet incorporation of the strongly destabilizing unlocked nucleic acid (UNA) modification into position 7 of the siRNA most potently reduced off-targeting for all tested sequences. Notably, such position-specific destabilization of siRNA-target interactions did not significantly reduce siRNA potency and is therefore well suited for future siRNA designs especially for applications in vivo where siRNA concentrations, expectedly, will be low.

Coordinated epigenetic repression of the miR-200 family and miR-205 in invasive bladder cancer.

MicroRNAs (miRNA) are small noncoding RNAs commonly deregulated in cancer. The miR-200 family (miR-200a, -200b, -200c, -141 and -429) and miR-205 are frequently silenced in advanced cancer and have been implicated in epithelial to mesenchymal transition (EMT) and tumor invasion by targeting the transcriptional repressors of E-cadherin, ZEB1 and ZEB2. ZEB1 is also known to repress miR-200c-141 transcription in a negative feedback loop, but otherwise little is known about the transcriptional regulation of the miR-200 family and miR-205. Recently, miR-200 silencing was also reported in cancer stem cells, implying that miR-200 deregulation is a key event in multiple levels of tumor biology. However, what prevents miR-200 expression remains largely unanswered. Here we report concerted transcriptional regulation of the miR-200 and miR-205 loci in bladder tumors and bladder cell lines. Using a combination of miRNA expression arrays, qPCR assays and mass spectrometry DNA methylation analyses, we show that the miR-200 and miR-205 loci are specifically silenced and gain promoter hypermethylation and repressive chromatin marks in muscle invasive bladder tumors and undifferentiated bladder cell lines. Moreover, we report that miR-200c expression is significantly correlated with early stage T1 bladder tumor progression, and propose miR-200 and miR-205 silencing and DNA hypermethylation as possible prognostic markers in bladder cancer. In addition, we observe that the mesoderm transcription factor TWIST1 and miR-200 expression are inversely correlated in bladder tumor samples and cell lines. TWIST1 associates directly with the miR-200 and miR-205 promoters, and may act as a repressor of miR-200 and miR-205 expression.

Enhancing miRNA annotation confidence in miRBase by continuous cross dataset analysis.

The immaculate annotation of all microRNAs (miRNAs) is a prerequisite to study their biological function on a genome-wide scale. However, the original criteria for proper miRNA annotation seem unsuited for the automated analysis of the immense number of small RNA reads available in next generation sequencing (NGS) datasets. Here we analyze the confidence of past miRNA annotation in miRBase by cross-analyzing publicly available NGS datasets using strengthened annotation requirements. Our analysis highlights that a large number of annotated human miRNAs in miRBase seems to require more experimental validation to be confidently annotated. Notably, our dataset analysis also identified almost 300 currently non-annotated miRNA*s and 28 novel miRNAs. These observations hereby greatly increase the confidence of past miRNA annotation in miRBase but also illustrate the usefulness of continuous re-evaluating NGS datasets in the identification of novel miRNAs.