Nuclear export of circular RNA.

Circular RNAs (circRNAs), which are increasingly being implicated in a variety of functions in normal and cancerous cells1-5, are formed by back-splicing of precursor mRNAs in the nucleus6-10. circRNAs are predominantly localized in the cytoplasm, indicating that they must be exported from the nucleus. Here we identify a pathway that is specific for the nuclear export of circular RNA. This pathway requires Ran-GTP, exportin-2 and IGF2BP1. Enhancing the nuclear Ran-GTP gradient by depletion or chemical inhibition of the major protein exporter CRM1 selectively increases the nuclear export of circRNAs, while reducing the nuclear Ran-GTP gradient selectively blocks circRNA export. Depletion or knockout of exportin-2 specifically inhibits nuclear export of circRNA. Analysis of nuclear circRNA-binding proteins reveals that interaction between IGF2BP1 and circRNA is enhanced by Ran-GTP. The formation of circRNA export complexes in the nucleus is promoted by Ran-GTP through its interactions with exportin-2, circRNA and IGF2BP1. Our findings demonstrate that adaptors such as IGF2BP1 that bind directly to circular RNAs recruit Ran-GTP and exportin-2 to export circRNAs in a mechanism that is analogous to protein export, rather than mRNA export.

ESRP1 controls biogenesis and function of a large abundant multiexon circRNA.

While the majority of circRNAs are formed from infrequent back-splicing of exons from protein coding genes, some can be produced at quite high level and in a regulated manner. We describe the regulation, biogenesis and function of circDOCK1(2-27), a large, abundant circular RNA that is highly regulated during epithelial-mesenchymal transition (EMT) and whose formation depends on the epithelial splicing regulator ESRP1. CircDOCK1(2-27) synthesis in epithelial cells represses cell motility both by diverting transcripts from DOCK1 mRNA production to circRNA formation and by direct inhibition of migration by the circRNA. HITS-CLIP analysis and CRISPR-mediated deletions indicate ESRP1 controls circDOCK1(2-27) biosynthesis by binding a GGU-containing repeat region in intron 1 and detaining its splicing until Pol II completes its 157 kb journey to exon 27. Proximity-dependent biotinylation (BioID) assay suggests ESRP1 may modify the RNP landscape of intron 1 in a way that disfavours communication of exon 1 with exon 2, rather than physically bridging exon 2 to exon 27. The X-ray crystal structure of RNA-bound ESRP1 qRRM2 domain reveals it binds to GGU motifs, with the guanines embedded in clamp-like aromatic pockets in the protein.

The landscape of alternative polyadenylation during EMT and its regulation by the RNA-binding protein Quaking.

Epithelial-mesenchymal transition (EMT) plays important roles in tumour progression and is orchestrated by dynamic changes in gene expression. While it is well established that post-transcriptional regulation plays a significant role in EMT, the extent of alternative polyadenylation (APA) during EMT has not yet been explored. Using 3' end anchored RNA sequencing, we mapped the alternative polyadenylation (APA) landscape following Transforming Growth Factor (TGF)-ß-mediated induction of EMT in human mammary epithelial cells and found APA generally causes 3'UTR lengthening during this cell state transition. Investigation of potential mediators of APA indicated the RNA-binding protein Quaking (QKI), a splicing factor induced during EMT, regulates a subset of events including the length of its own transcript. Analysis of QKI crosslinked immunoprecipitation (CLIP)-sequencing data identified the binding of QKI within 3' untranslated regions (UTRs) was enriched near cleavage and polyadenylation sites. Following QKI knockdown, APA of many transcripts is altered to produce predominantly shorter 3'UTRs associated with reduced gene expression. These findings reveal the changes in APA that occur during EMT and identify a potential role for QKI in this process.

Basonuclin-2 regulates extracellular matrix production and degradation.

Epithelial-mesenchymal transition is essential for tissue patterning and organization. It involves both regulation of cell motility and alterations in the composition and organization of the ECM-a complex environment of proteoglycans and fibrous proteins essential for tissue homeostasis, signaling in response to chemical and biomechanical stimuli, and is often dysregulated under conditions such as cancer, fibrosis, and chronic wounds. Here, we demonstrate that basonuclin-2 (BNC2), a mesenchymal-expressed gene, that is, strongly associated with cancer and developmental defects across genome-wide association studies, is a novel regulator of ECM composition and degradation. We find that at endogenous levels, BNC2 controls the expression of specific collagens, matrix metalloproteases, and other matrisomal components in breast cancer cells, and in fibroblasts that are primarily responsible for the production and processing of the ECM within the tumour microenvironment. In so doing, BNC2 modulates the motile and invasive properties of cancers, which likely explains the association of high BNC2 expression with increasing cancer grade and poor patient prognosis.

On the rules of engagement for microRNAs targeting protein coding regions.

MiRNAs post-transcriptionally repress gene expression by binding to mRNA 3'UTRs, but the extent to which they act through protein coding regions (CDS regions) is less well established. MiRNA interaction studies show a substantial proportion of binding occurs in CDS regions, however sequencing studies show much weaker effects on mRNA levels than from 3'UTR interactions, presumably due to competition from the translating ribosome. Consequently, most target prediction algorithms consider only 3'UTR interactions. However, the consequences of CDS interactions may have been underestimated, with the reporting of a novel mode of miRNA-CDS interaction requiring base pairing of the miRNA 3' end, but not the canonical seed site, leading to repression of translation with little effect on mRNA turnover. Using extensive reporter, western blotting and bioinformatic analyses, we confirm that miRNAs can indeed suppress genes through CDS-interaction in special circumstances. However, in contrast to that previously reported, we find repression requires extensive base-pairing, including of the canonical seed, but does not strictly require base pairing of the 3' miRNA terminus and is mediated through reducing mRNA levels. We conclude that suppression of endogenous genes can occur through miRNAs binding to CDS, but the requirement for extensive base-pairing likely limits the regulatory impacts to modest effects on a small subset of targets.

DraculR: A Web-Based Application for In Silico Haemolysis Detection in High-Throughput microRNA Sequencing Data.

The search for novel microRNA (miRNA) biomarkers in plasma is hampered by haemolysis, the lysis and subsequent release of red blood cell contents, including miRNAs, into surrounding fluid. The biomarker potential of miRNAs comes in part from their multicompartment origin and the long-lived nature of miRNA transcripts in plasma, giving researchers a functional window for tissues that are otherwise difficult or disadvantageous to sample. The inclusion of red-blood-cell-derived miRNA transcripts in downstream analysis introduces a source of error that is difficult to identify posthoc and may lead to spurious results. Where access to a physical specimen is not possible, our tool will provide an in silico approach to haemolysis prediction. We present DraculR, an interactive Shiny/R application that enables a user to upload miRNA expression data from a short-read sequencing of human plasma as a raw read counts table and interactively calculate a metric that indicates the degree of haemolysis contamination. The code, DraculR web tool and its tutorial are freely available as detailed herein.

Haemolysis Detection in MicroRNA-Seq from Clinical Plasma Samples.

The abundance of cell-free microRNA (miRNA) has been measured in blood plasma and proposed as a source of novel, minimally invasive biomarkers for several diseases. Despite improvements in quantification methods, there is no consensus regarding how haemolysis affects plasma miRNA content. We propose a method for haemolysis detection in miRNA high-throughput sequencing (HTS) data from libraries prepared using human plasma. To establish a miRNA haemolysis signature we tested differential miRNA abundance between plasma samples with known haemolysis status. Using these miRNAs with statistically significant higher abundance in our haemolysed group, we further refined the set to reveal high-confidence haemolysis association. Given our specific context, i.e., women of reproductive age, we also tested for significant differences between pregnant and non-pregnant groups. We report a novel 20-miRNA signature used to identify the presence of haemolysis in silico in HTS miRNA-sequencing data. Further, we validated the signature set using firstly an all-male cohort (prostate cancer) and secondly a mixed male and female cohort (radiographic knee osteoarthritis). Conclusion: Given the potential for haemolysis contamination, we recommend that assays for haemolysis detection become standard pre-analytical practice and provide here a simple method for haemolysis detection.

Differentially expressed miRNAs in bone after methotrexate treatment.

Previous studies have shown that administration of antimetabolite methotrexate (MTX) caused a reduced trabecular bone volume and increased marrow adiposity (bone/fat switch), for which the underlying molecular mechanisms and recovery potential are unclear. Altered expression of microRNAs (miRNAs) has been shown to be associated with dysregulation of osteogenic and/or adipogenic differentiation by disrupting target gene expression. First, the current study confirmed the bone/fat switch following MTX treatment in precursor cell culture models in vitro. Then, using a rat intensive 5-once daily MTX treatment model, this study aimed to identify miRNAs associated with bone damage and recovery (in a time course over Days 3, 6, 9, and 14 after the first MTX treatment). RNA isolated from bone samples of treated and control rats were subjected to miRNA array and reverse transcription-polymerase chain reaction validation, which identified five upregulated miRNA candidates, namely, miR-155-5p, miR-154-5p, miR-344g, miR-6215, and miR-6315. Target genes of these miRNAs were predicted using TargetScan and miRDB. Then, the protein-protein network was established via STRING database, after which the miRNA-key messenger RNA (mRNA) network was constructed by Cytoscape. Functional annotation and pathway enrichment analyses for miR-6315 were performed by DAVID database. We found that TGF-ß signaling was the most significantly enriched pathway and subsequent dual-luciferase assays suggested that Smad2 was the direct target of miR-6315. Our current study showed that miR-6315 might be a vital regulator involved in bone and marrow fat formation. Also, this study constructed a comprehensive miRNA-mRNA regulatory network, which may contribute to the pathogenesis/prognosis of MTX-associated bone loss and bone marrow adiposity.

Post-transcriptional Gene Regulation by MicroRNA-194 Promotes Neuroendocrine Transdifferentiation in Prostate Cancer.

Potent therapeutic inhibition of the androgen receptor (AR) in prostate adenocarcinoma can lead to the emergence of neuroendocrine prostate cancer (NEPC), a phenomenon associated with enhanced cell plasticity. Here, we show that microRNA-194 (miR-194) is a regulator of epithelial-neuroendocrine transdifferentiation. In clinical prostate cancer samples, miR-194 expression and activity were elevated in NEPC and inversely correlated with AR signaling. miR-194 facilitated the emergence of neuroendocrine features in prostate cancer cells, a process mediated by its ability to directly target a suite of genes involved in cell plasticity. One such target was FOXA1, which encodes a transcription factor with a vital role in maintaining the prostate epithelial lineage. Importantly, a miR-194 inhibitor blocked epithelial-neuroendocrine transdifferentiation and inhibited the growth of cell lines and patient-derived organoids possessing neuroendocrine features. Overall, our study reveals a post-transcriptional mechanism regulating the plasticity of prostate cancer cells and provides a rationale for targeting miR-194 in NEPC.

Insufficiently complex unique-molecular identifiers (UMIs) distort small RNA sequencing.

The attachment of unique molecular identifiers (UMIs) to RNA molecules prior to PCR amplification and sequencing, makes it possible to amplify libraries to a level that is sufficient to identify rare molecules, whilst simultaneously eliminating PCR bias through the identification of duplicated reads. Accurate de-duplication is dependent upon a sufficiently complex pool of UMIs to allow unique labelling. In applications dealing with complex libraries, such as total RNA-seq, only a limited variety of UMIs are required as the variation in molecules to be sequenced is enormous. However, when sequencing a less complex library, such as small RNAs for which there is a more limited range of possible sequences, we find increased variation in UMIs are required, even beyond that provided in a commercial kit specifically designed for the preparation of small RNA libraries for sequencing. We show that a pool of UMIs randomly varying across eight nucleotides is not of sufficient depth to uniquely tag the microRNAs to be sequenced. This results in over de-duplication of reads and the marked under-estimation of expression of the more abundant microRNAs. Whilst still arguing for the utility of UMIs, this work demonstrates the importance of their considered design to avoid errors in the estimation of gene expression in libraries derived from select regions of the transcriptome or small genomes.

Extensive transcriptional responses are co-ordinated by microRNAs as revealed by Exon-Intron Split Analysis (EISA).

Epithelial-mesenchymal transition (EMT) has been a subject of intense scrutiny as it facilitates metastasis and alters drug sensitivity. Although EMT-regulatory roles for numerous miRNAs and transcription factors are known, their functions can be difficult to disentangle, in part due to the difficulty in identifying direct miRNA targets from complex datasets and in deciding how to incorporate 'indirect' miRNA effects that may, or may not, represent biologically relevant information. To better understand how miRNAs exert effects throughout the transcriptome during EMT, we employed Exon-Intron Split Analysis (EISA), a bioinformatic technique that separates transcriptional and post-transcriptional effects through the separate analysis of RNA-Seq reads mapping to exons and introns. We find that in response to the manipulation of miRNAs, a major effect on gene expression is transcriptional. We also find extensive co-ordination of transcriptional and post-transcriptional regulatory mechanisms during both EMT and mesenchymal to epithelial transition (MET) in response to TGF-ß or miR-200c respectively. The prominent transcriptional influence of miRNAs was also observed in other datasets where miRNA levels were perturbed. This work cautions against a narrow approach that is limited to the analysis of direct targets, and demonstrates the utility of EISA to examine complex regulatory networks involving both transcriptional and post-transcriptional mechanisms.

miRNA length variation during macrophage stimulation confounds the interpretation of results: implications for miRNA quantification by RT-qPCR.

Most microRNAs (miRNAs) are expressed as a mix of length isoforms (referred to as isomiRs). IsomiR stoichiometry can be differentially impacted upon cell stimulation, as recently evidenced by our group in the context of immune responses induced by type-I interferon (IFN). Here, we revisit published RNA-seq data sets of human and mouse macrophages stimulated with bacterial products at the isomiR level. We demonstrate that for several miRNAs, macrophage stimulation induces changes in isomiR stoichiometry. Critically, we find that changes in miRNA expression can be misinterpreted when miRNAs are quantified by RT-qPCR, as primers directed against canonical miRNA sequences may not equally target the different isomiRs that are regulated endogenously. Beyond the case of phagocyte stimulation, our analyses reinforce the concept that analysis of miRNA expression at the isoform level should become standard practice.

Combinatorial Targeting by MicroRNAs Co-ordinates Post-transcriptional Control of EMT.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression, functioning in part by facilitating the degradation of target mRNAs. They have an established role in controlling epithelial-mesenchymal transition (EMT), a reversible phenotypic program underlying normal and pathological processes. Many studies demonstrate the role of individual miRNAs using overexpression at levels greatly exceeding physiological abundance. This can influence transcripts with relatively poor targeting and may in part explain why over 130 different miRNAs are directly implicated as EMT regulators. Analyzing a human mammary cell model of EMT we found evidence that a set of miRNAs, including the miR-200 and miR-182/183 family members, co-operate in post-transcriptional regulation, both reinforcing and buffering transcriptional output. Investigating this, we demonstrate that combinatorial treatment altered cellular phenotype with miRNA concentrations much closer to endogenous levels and with less off-target effects. This suggests that co-operative targeting by miRNAs is important for their physiological function and future work classifying miRNAs should consider such combinatorial effects.

miR-200/375 control epithelial plasticity-associated alternative splicing by repressing the RNA-binding protein Quaking.

Members of the miR-200 family are critical gatekeepers of the epithelial state, restraining expression of pro-mesenchymal genes that drive epithelial-mesenchymal transition (EMT) and contribute to metastatic cancer progression. Here, we show that miR-200c and another epithelial-enriched miRNA, miR-375, exert widespread control of alternative splicing in cancer cells by suppressing the RNA-binding protein Quaking (QKI). During EMT, QKI-5 directly binds to and regulates hundreds of alternative splicing targets and exerts pleiotropic effects, such as increasing cell migration and invasion and restraining tumour growth, without appreciably affecting mRNA levels. QKI-5 is both necessary and sufficient to direct EMT-associated alternative splicing changes, and this splicing signature is broadly conserved across many epithelial-derived cancer types. Importantly, several actin cytoskeleton-associated genes are directly targeted by both QKI and miR-200c, revealing coordinated control of alternative splicing and mRNA abundance during EMT These findings demonstrate the existence of a miR-200/miR-375/QKI axis that impacts cancer-associated epithelial cell plasticity through widespread control of alternative splicing.

miR-222 isoforms are differentially regulated by type-I interferon.

Endogenous microRNAs (miRNAs) often exist as multiple isoforms (known as "isomiRs") with predominant variation around their 3'-end. Increasing evidence suggests that different isomiRs of the same family can have diverse functional roles, as recently demonstrated with the example of miR-222-3p 3'-end variants. While isomiR levels from a same miRNA family can vary between tissues and cell types, change of templated isomiR stoichiometry to stimulation has not been reported to date. Relying on small RNA-sequencing analyses, we demonstrate here that miR-222-3p 3'-end variants >23 nt are specifically decreased upon interferon (IFN) ß stimulation of human fibroblasts, while shorter isoforms are spared. This length-dependent dynamic regulation of long miR-222-3p 3'-isoforms and >40 other miRNA families was confirmed in human monocyte-derived dendritic cells following infection with Salmonella Typhimurium, underlining the breadth of 3'-length regulation by infection, beyond the example of miR-222-3p. We further show that stem-loop miRNA Taqman RT-qPCR exhibits selectivity between 3'-isoforms, according to their length, and that this can lead to misinterpretation of results when these isoforms are differentially regulated. Collectively, and to our knowledge, this work constitutes the first demonstration that the stoichiometry of highly abundant templated 3'-isoforms of a same miRNA family can be dynamically regulated by a stimulus. Given that such 3'-isomiRs can have different functions, our study underlines the need to consider isomiRs when investigating miRNA-based regulation.

Naturally existing isoforms of miR-222 have distinct functions.

Deep-sequencing reveals extensive variation in the sequence of endogenously expressed microRNAs (termed 'isomiRs') in human cell lines and tissues, especially in relation to the 3' end. From the immunoprecipitation of the microRNA-binding protein Argonaute and the sequencing of associated small RNAs, we observe extensive 3'-isomiR variation, including for miR-222 where the majority of endogenously expressed miR-222 is extended by 1-5 nt compared to the canonical sequence. We demonstrate this 3' heterogeneity has dramatic implications for the phenotype of miR-222 transfected cells, with longer isoforms promoting apoptosis in a size (but not 3' sequence)-dependent manner. The transfection of longer miR-222 isomiRs did not induce an interferon response, but did downregulate the expression of many components of the pro-survival PI3K-AKT pathway including PIK3R3, a regulatory subunit whose knockdown phenocopied the expression of longer 222 isoforms in terms of apoptosis and the inhibition of other PI3K-AKT genes. As this work demonstrates the capacity for 3' isomiRs to mediate differential functions, we contend more attention needs to be given to 3' variance given the prevalence of this class of isomiR.

p53 Represses the Oncogenic Sno-MiR-28 Derived from a SnoRNA.

p53 is a master tumour repressor that participates in vast regulatory networks, including feedback loops involving microRNAs (miRNAs) that regulate p53 and that themselves are direct p53 transcriptional targets. We show here that a group of polycistronic miRNA-like non-coding RNAs derived from small nucleolar RNAs (sno-miRNAs) are transcriptionally repressed by p53 through their host gene, SNHG1. The most abundant of these, sno-miR-28, directly targets the p53-stabilizing gene, TAF9B. Collectively, p53, SNHG1, sno-miR-28 and TAF9B form a regulatory loop which affects p53 stability and downstream p53-regulated pathways. In addition, SNHG1, SNORD28 and sno-miR-28 are all significantly upregulated in breast tumours and the overexpression of sno-miR-28 promotes breast epithelial cell proliferation. This research has broadened our knowledge of the crosstalk between small non-coding RNA pathways and roles of sno-miRNAs in p53 regulation.

The RNA binding protein quaking regulates formation of circRNAs.

Circular RNAs (circRNAs), formed by non-sequential back-splicing of pre-mRNA transcripts, are a widespread form of non-coding RNA in animal cells. However, it is unclear whether the majority of circRNAs represent splicing by-products without function or are produced in a regulated manner to carry out specific cellular functions. We show that hundreds of circRNAs are regulated during human epithelial-mesenchymal transition (EMT) and find that the production of over one-third of abundant circRNAs is dynamically regulated by the alternative splicing factor, Quaking (QKI), which itself is regulated during EMT. Furthermore, by modulating QKI levels, we show the effect on circRNA abundance is dependent on intronic QKI binding motifs. Critically, the addition of QKI motifs is sufficient to induce de novo circRNA formation from transcripts that are normally linearly spliced. These findings demonstrate circRNAs are both purposefully synthesized and regulated by cell-type specific mechanisms, suggesting they play specific biological roles in EMT.

Assessing the gene regulatory properties of Argonaute-bound small RNAs of diverse genomic origin.

High-throughput sequencing reveals an abundance of microRNA-sized fragments derived from larger non-coding RNAs. Roles for these small RNAs in gene silencing are suggested by their co-precipitation with Argonaute, the microRNA effector protein, though the extent to which they suppress gene expression endogenously remains unclear. To address this, we used luciferase reporters to determine the endogenous functionality of small RNAs from a diverse range of sources. We demonstrate small RNAs derived from snoRNAs have the capacity to act in a microRNA-like manner, though we note the vast majority of these are bound to Argonaute at levels below that required for detectable silencing activity. We show Argonaute exhibits a high degree of selectivity for the small RNAs with which it interacts and note that measuring Argonaute-associated levels is a better indicator of function than measuring total expression. Although binding to Argonaute at sufficient levels is necessary for demonstrating microRNA functionality in our reporter assay, this alone is not enough as some small RNAs derived from other non-coding RNAs (tRNAs, rRNAs, Y-RNAs) are associated with Argonaute at very high levels yet do not serve microRNA-like roles.