Global identification of mRNA-interacting circular RNAs by CLiPPR-Seq.

Although the functional role of circular RNA (circRNA) interaction with microRNAs and proteins has been studied extensively, circRNA interactions with the protein-coding mRNAs in intact cells remain largely unknown. Here, by employing AMT-mediated proximity ligation of RNA-RNA duplexes followed by circRNA enrichment and deep sequencing, we report a novel Cross-Linking Poly(A) Pulldown RNase R Sequencing (CLiPPR-seq) technology which identified hundreds of mRNA-interacting circRNAs in three different cell types, including ßTC6, C2C12 and HeLa cells. Furthermore, CLiPP-seq without RNase R treatment was also performed to identify the mRNA expression in these cells. BLAST analysis of circRNAs in CLiPPR-seq sample with the mRNAs in CLiPP-seq samples determined their potential complementary sequences for circRNA-mRNA interaction. Pulldown of circRNAs and poly(A) RNAs confirmed the direct interaction of circRNAs with target mRNAs. Silencing of mRNA-interacting circRNAs led to the altered expression of target mRNAs in ßTC6 cells, suggesting the role of direct interaction of circRNAs with mRNAs in gene expression regulation. CLiPPR-seq thus represents a novel method for illuminating the myriad of uncharacterized circRNA-mRNA hybrids that may regulate gene expression.

Subcellular localization of circular RNAs: Where and why.

Localization of RNAs at specific subcellular locations regulating various local cellular events has gained much attention recently. Like most other classes of RNAs, the function of newly discovered circular RNAs (circRNAs) is predominantly determined by their association with different cellular factors in the cell. CircRNAs function as transcriptional and posttranscriptional regulators of gene expression by interacting with transcription factors, splicing regulators, RNA-binding proteins, and microRNAs or by translating into functional polypeptides. Hence, studying their subcellular localization to assess their function is essential. The discovery of more than a million circRNA and increasing evidence of their involvement in development and diseases require a thorough analysis of their subcellular localization linking to their biological functions. Here, we summarize current knowledge of circRNA localization in cells and extracellular vesicles, factors regulating their subcellular localization, and the implications of circRNA localization on their cellular functions. Given the discovery of many circRNAs in all life forms and their implications in pathophysiology, we discuss the challenges in studying circRNA localization and the opportunities for unlocking the mystery of circRNA functions.

Synergistic Interaction between Polysaccharide-Based Extracellular Matrix and Mineralized Osteoblast-Derived EVs Promotes Bone Regeneration via miRNA-mRNA Regulatory Axis.

Extracellular vesicles (EVs) derived from bone progenitor cells are advantageous as cell-free and non-immunogenic cargo delivery vehicles. In this study, EVs are isolated from MC3T3-E1 cells before (GM-EVs) and after mineralization for 7 and 14 days (DM-EVs). It was observed that DM-EVs accelerate the process of differentiation in recipient cells more prominently. The small RNA sequencing of EVs revealed that miR-204-5p, miR-221-3p, and miR-148a-3p are among the highly upregulated miRNAs that have an inhibitory effect on the function of mRNAs, Sox11, Timp3, and Ccna2 in host cells, which is probably responsible for enhancing the activity of osteoblastic genes. To enhance the bioavailability of EVs, they are encapsulated in a chitosan-collagen composite hydrogel that serves as a bioresorbable extracellular matrix (ECM). The EVs-integrated scaffold (DM-EVs + Scaffold) enhances bone regeneration in critical-sized calvarial bone defects in rats within 8 weeks of implantation by providing the ECM cues. The shelf life of DM-EVs + Scaffold indicates that the bioactivity of EVs and their cargo in the polymer matrix remains intact for up to 30 days. Integrating mineralized cell-derived EVs into an ECM represents a bioresorbable matrix with a cell-free method for promoting new bone formation through the miRNA-mRNA regulatory axis.

HuR and GRSF1 modulate the nuclear export and mitochondrial localization of the lncRNA RMRP.

Some mitochondrial long noncoding RNAs (lncRNAs) are encoded by nuclear DNA, but the mechanisms that mediate their transport to mitochondria are poorly characterized. Using affinity RNA pull-down followed by mass spectrometry analysis, we found two RNA-binding proteins (RBPs), HuR (human antigen R) and GRSF1 (G-rich RNA sequence-binding factor 1), that associated with the nuclear DNA-encoded lncRNA RMRP and mobilized it to mitochondria. In cultured human cells, HuR bound RMRP in the nucleus and mediated its CRM1 (chromosome region maintenance 1)-dependent export to the cytosol. After RMRP was imported into mitochondria, GRSF1 bound RMRP and increased its abundance in the matrix. Loss of GRSF1 lowered the mitochondrial levels of RMRP, in turn suppressing oxygen consumption rates and modestly reducing mitochondrial DNA replication priming. Our findings indicate that RBPs HuR and GRSF1 govern the cytoplasmic and mitochondrial localization of the lncRNA RMRP, which is encoded by nuclear DNA but has key functions in mitochondria.

miR-431 promotes differentiation and regeneration of old skeletal muscle by targeting Smad4.

The myogenic capacity of myoblasts decreases in skeletal muscle with age. In addition to environmental factors, intrinsic factors are important for maintaining the regenerative potential of muscle progenitor cells, but their identities are largely unknown. Here, comparative analysis of microRNA (miRNA) expression profiles in young and old myoblasts uncovered miR-431 as a novel miRNA showing markedly reduced abundance in aged myoblasts. Importantly, elevating miR-431 improved the myogenic capacity of old myoblasts, while inhibiting endogenous miR-431 lowered myogenesis. Bioinformatic and biochemical analyses revealed that miR-431 directly interacted with the 3' untranslated region (UTR) of Smad4 mRNA, which encodes one of the downstream effectors of TGF-ß signaling. In keeping with the low levels of miR-431 in old myoblasts, SMAD4 levels increased in this myoblast population. Interestingly, in an in vivo model of muscle regeneration following cardiotoxin injury, ectopic miR-431 injection greatly improved muscle regeneration and reduced SMAD4 levels. Consistent with the finding that the mouse miR-431 seed sequence in the Smad4 3' UTR is conserved in the human SMAD4 3' UTR, inhibition of miR-431 also repressed the myogenic capacity of human skeletal myoblasts. Taken together, our results suggest that the age-associated miR-431 plays a key role in maintaining the myogenic ability of skeletal muscle with age.

circSamd4 represses myogenic transcriptional activity of PUR proteins.

By interacting with proteins and nucleic acids, the vast family of mammalian circRNAs is proposed to influence many biological processes. Here, RNA sequencing analysis of circRNAs differentially expressed during myogenesis revealed that circSamd4 expression increased robustly in mouse C2C12 myoblasts differentiating into myotubes. Moreover, silencing circSamd4, which is conserved between human and mouse, delayed myogenesis and lowered the expression of myogenic markers in cultured myoblasts from both species. Affinity pulldown followed by mass spectrometry revealed that circSamd4 associated with PURA and PURB, two repressors of myogenesis that inhibit transcription of the myosin heavy chain (MHC) protein family. Supporting the hypothesis that circSamd4 might complex with PUR proteins and thereby prevent their interaction with DNA, silencing circSamd4 enhanced the association of PUR proteins with the Mhc promoter, while overexpressing circSamd4 interfered with the binding of PUR proteins to the Mhc promoter. These effects were abrogated when using a mutant circSamd4 lacking the PUR binding site. Our results indicate that the association of PUR proteins with circSamd4 enhances myogenesis by contributing to the derepression of MHC transcription.

RPAD (RNase R treatment, polyadenylation, and poly(A)+ RNA depletion) method to isolate highly pure circular RNA.

Recent developments in high-throughput RNA sequencing methods coupled with innovative bioinformatic tools have uncovered thousands of circular (circ)RNAs. CircRNAs have emerged as a vast and novel class of regulatory RNAs with potential to modulate gene expression by acting as sponges for microRNAs (miRNAs) and RNA-binding proteins (RBPs). The biochemical enrichment of circRNAs by exoribonuclease treatment or by depletion of polyadenylated RNAs coupled with deep-sequencing is widely used for the systematic identification of circRNAs. Although these methods enrich circRNAs substantially, they do not eliminate efficiently non-polyadenylated and highly-structured RNAs. Here, we describe a method we termed RPAD, based on initial RNase R treatment followed by Polyadenylation and poly(A)+ RNA Depletion. These joint interventions drastically depleted linear RNAs leading to isolation of highly pure circRNAs from total RNA pools. By facilitating the isolation of highly pure circRNAs, RPAD enables the elucidation of circRNA biogenesis, sequence, and function.

Identification and Characterization of Circular Intronic RNAs Derived from Insulin Gene.

Circular RNAs (circRNAs) are a large family of noncoding RNAs that have emerged as novel regulators of gene expression. However, little is known about the function of circRNAs in pancreatic ß-cells. Here, transcriptomic analysis of mice pancreatic islet RNA-sequencing data identified 77 differentially expressed circRNAs between mice fed with a normal diet and a high-fat diet. Surprisingly, multiple circRNAs were derived from the intron 2 of the preproinsulin 2 (Ins2) gene and are termed as circular intronic (ci)-Ins2. The expression of ci-Ins2 transcripts in mouse pancreatic islets, and ßTC6 cells were confirmed by reverse transcription PCR, DNA sequencing, and RNase R treatment experiments. The level of ci-Ins2 was altered in ßTC6 cells upon exposure to elevated levels of palmitate and glucose. Computational analysis predicted the interaction of several RNA-binding proteins with ci-Ins2 and their flanking region, suggesting their role in the ci-Ins2 function or biogenesis. Additionally, bioinformatics analysis predicted the association of several microRNAs with ci-Ins2. Gene ontology and pathway analysis of genes targeted by miRNAs associated with ci-Ins2 suggested the regulation of several key biological processes. Together, our findings indicate that differential expression of circRNAs, especially ci-Ins2 transcripts, may regulate ß-cell function and may play a critical role in the development of diabetes.

High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs.

High-throughput RNA sequencing methods coupled with specialized bioinformatic analyses have recently uncovered tens of thousands of unique circular (circ)RNAs, but their complete sequences, genes of origin and functions are largely unknown. Given that circRNAs lack free ends and are thus relatively stable, their association with microRNAs (miRNAs) and RNA-binding proteins (RBPs) can influence gene expression programs. While exoribonuclease treatment is widely used to degrade linear RNAs and enrich circRNAs in RNA samples, it does not efficiently eliminate all linear RNAs. Here, we describe a novel method for the isolation of highly pure circRNA populations involving RNase R treatment followed by Polyadenylation and poly(A)+ RNA Depletion (RPAD), which removes linear RNA to near completion. High-throughput sequencing of RNA prepared using RPAD from human cervical carcinoma HeLa cells and mouse C2C12 myoblasts led to two surprising discoveries: (i) many exonic circRNA (EcircRNA) isoforms share an identical backsplice sequence but have different body sizes and sequences, and (ii) thousands of novel intronic circular RNAs (IcircRNAs) are expressed in cells. In sum, isolating high-purity circRNAs using the RPAD method can enable quantitative and qualitative analyses of circRNA types and sequence composition, paving the way for the elucidation of circRNA functions.

Identification of senescence-associated circular RNAs (SAC-RNAs) reveals senescence suppressor CircPVT1.

Using RNA sequencing (RNA-Seq), we compared the expression patterns of circular RNAs in proliferating (early-passage) and senescent (late-passage) human diploid WI-38 fibroblasts. Among the differentially expressed senescence-associated circRNAs (which we termed 'SAC-RNAs'), we identified CircPVT1, generated by circularization of an exon of the PVT1 gene, as a circular RNA showing markedly reduced levels in senescent fibroblasts. Reducing CircPVT1 levels in proliferating fibroblasts triggered senescence, as determined by a rise in senescence-associated ß-galactosidase activity, higher abundance of CDKN1A/P21 and TP53, and reduced cell proliferation. Although several microRNAs were predicted to bind CircPVT1, only let-7 was found enriched after pulldown of endogenous CircPVT1, suggesting that CircPVT1 might selectively modulate let-7 activity and hence expression of let-7-regulated mRNAs. Reporter analysis revealed that CircPVT1 decreased the cellular pool of available let-7, and antagonizing endogenous let-7 triggered cell proliferation. Importantly, silencing CircPVT1 promoted cell senescence and reversed the proliferative phenotype observed after let-7 function was impaired. Consequently, the levels of several proliferative proteins that prevent senescence, such as IGF2BP1, KRAS and HMGA2, encoded by let-7 target mRNAs, were reduced by silencing CircPVT1. Our findings indicate that the SAC-RNA CircPVT1, elevated in dividing cells and reduced in senescent cells, sequesters let-7 to enable a proliferative phenotype.

Rolling Circle cDNA Synthesis Uncovers Circular RNA Splice Variants.

High-throughput RNA sequencing and novel bioinformatic pipelines have identified thousands of circular (circ)RNAs containing backsplice junction sequences. However, circRNAs generated from multiple exons may contain different combinations of exons and/or introns arising from alternative splicing, while the backsplice junction sequence is the same. To be able to identify circRNA splice variants, we developed a method termed circRNA-Rolling Circle Amplification (circRNA-RCA). This method detects full-length circRNA sequences by performing reverse transcription (RT) in the absence of RNase H activity, followed by polymerase chain reaction (PCR) amplification of full-length circRNAs using a forward primer spanning the backsplice junction sequence and a reverse primer exactly upstream of the forward primer. By sequencing the PCR products, circRNA splice variants bearing the same backsplice junctions, which were otherwise only predicted computationally, could be experimentally validated. The splice variants were further predicted to associate with different subsets of target RNA-binding proteins and microRNAs, supporting the notion that different circRNA splice variants can have different biological impacts. In sum, the circRNA-RCA method allows the accurate identification of full-length circRNA sequences, offering unique insight into their individual function.

Novel RNA-binding activity of MYF5 enhances Ccnd1/Cyclin D1 mRNA translation during myogenesis.

Skeletal muscle contains long multinucleated and contractile structures known as muscle fibers, which arise from the fusion of myoblasts into multinucleated myotubes during myogenesis. The myogenic regulatory factor (MRF) MYF5 is the earliest to be expressed during myogenesis and functions as a transcription factor in muscle progenitor cells (satellite cells) and myocytes. In mouse C2C12 myocytes, MYF5 is implicated in the initial steps of myoblast differentiation into myotubes. Here, using ribonucleoprotein immunoprecipitation (RIP) analysis, we discovered a novel function for MYF5 as an RNA-binding protein which associated with a subset of myoblast mRNAs. One prominent MYF5 target was Ccnd1 mRNA, which encodes the key cell cycle regulator CCND1 (Cyclin D1). Biotin-RNA pulldown, UV-crosslinking and gel shift experiments indicated that MYF5 was capable of binding the 3' untranslated region (UTR) and the coding region (CR) of Ccnd1 mRNA. Silencing MYF5 expression in proliferating myoblasts revealed that MYF5 promoted CCND1 translation and modestly increased transcription of Ccnd1 mRNA. Accordingly, overexpressing MYF5 in C2C12 cells upregulated CCND1 expression while silencing MYF5 reduced myoblast proliferation as well as differentiation of myoblasts into myotubes. Moreover, MYF5 silencing reduced myogenesis, while ectopically restoring CCND1 abundance partially rescued the decrease in myogenesis seen after MYF5 silencing. We propose that MYF5 enhances early myogenesis in part by coordinately elevating Ccnd1 transcription and Ccnd1 mRNA translation.

Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1.

HuR influences gene expression programs and hence cellular phenotypes by binding to hundreds of coding and noncoding linear RNAs. However, whether HuR binds to circular RNAs (circRNAs) and impacts on their function is unknown. Here, we have identified en masse circRNAs binding HuR in human cervical carcinoma HeLa cells. One of the most prominent HuR target circRNAs was hsa_circ_0031288, renamed CircPABPN1 as it arises from the PABPN1 pre-mRNA. Further analysis revealed that HuR did not influence CircPABPN1 abundance; interestingly, however, high levels of CircPABPN1 suppressed HuR binding to PABPN1 mRNA. Evaluation of PABPN1 mRNA polysomes indicated that PABPN1 translation was modulated positively by HuR and hence negatively by CircPABPN1. We propose that the extensive binding of CircPABPN1 to HuR prevents HuR binding to PABPN1 mRNA and lowers PABPN1 translation, providing the first example of competition between a circRNA and its cognate mRNA for an RBP that affects translation.

CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs.

Circular RNAs (circRNAs) are widely expressed in animal cells, but their biogenesis and functions are poorly understood. CircRNAs have been shown to act as sponges for miRNAs and may also potentially sponge RNA-binding proteins (RBPs) and are thus predicted to function as robust posttranscriptional regulators of gene expression. The joint analysis of large-scale transcriptome data coupled with computational analyses represents a powerful approach to elucidate possible biological roles of ribonucleoprotein (RNP) complexes. Here, we present a new web tool, CircInteractome (circRNA interactome), for mapping RBP- and miRNA-binding sites on human circRNAs. CircInteractome searches public circRNA, miRNA, and RBP databases to provide bioinformatic analyses of binding sites on circRNAs and additionally analyzes miRNA and RBP sites on junction and junction-flanking sequences. CircInteractome also allows the user the ability to (1) identify potential circRNAs which can act as RBP sponges, (2) design junction-spanning primers for specific detection of circRNAs of interest, (3) design siRNAs for circRNA silencing, and (4) identify potential internal ribosomal entry sites (IRES). In sum, the web tool CircInteractome, freely accessible at http://circinteractome.nia.nih.gov, facilitates the analysis of circRNAs and circRNP biology.

7SL RNA represses p53 translation by competing with HuR.

Noncoding RNAs (ncRNAs) and RNA-binding proteins are potent post-transcriptional regulators of gene expression. The ncRNA 7SL is upregulated in cancer cells, but its impact upon the phenotype of cancer cells is unknown. Here, we present evidence that 7SL forms a partial hybrid with the 3'-untranslated region (UTR) of TP53 mRNA, which encodes the tumor suppressor p53. The interaction of 7SL with TP53 mRNA reduced p53 translation, as determined by analyzing p53 expression levels, nascent p53 translation and TP53 mRNA association with polysomes. Silencing 7SL led to increased binding of HuR to TP53 mRNA, an interaction that led to the promotion of p53 translation and increased p53 abundance. We propose that the competition between 7SL and HuR for binding to TP53 3'UTR contributes to determining the magnitude of p53 translation, in turn affecting p53 levels and the growth-suppressive function of p53. Our findings suggest that targeting 7SL may be effective in the treatment of cancers with reduced p53 levels.

Sanger Sequencing to Determine the Full-Length Sequence of Circular RNAs.

The pre-existing theory of pre-mRNA splicing into linear mature RNA was questioned with the introduction of circular RNAs (circRNAs). Hundreds of studies using high throughput RNA-sequencing (RNA-seq) techniques and novel computational programs reported the abundant and ubiquitous expression of circRNAs originating by pre-mRNA backsplicing. CircRNAs are mostly involved in gene expression by regulating functions of interacting microRNAs (miRNAs) and RNA-binding proteins (RBPs) or translating into functional polypeptides. Although all circRNA annotation tools identify circRNAs based on the backsplice junction (BSJ) sequences, only a few identify the internal sequences of circRNAs. However, the full-length sequence of circRNAs from RNA-seq data could be error-prone due to its similarity with the counterpart linear RNA. Since circRNA function depends on the mature sequence, validation of the mature sequence is the prerequisite for their further characterization. In this chapter, we discuss the validation of circRNA BSJ sequence by RT-PCR using divergent primer followed by Sanger sequencing. Furthermore, we describe the circRNA-rolling circle amplification (circRNA-RCA; circRNA enrichment by RNase R treatment, full-length cDNA synthesis, rolling circle PCR amplification using full-length primers, and Sanger sequencing of the PCR product) to validate the mature splice sequence of circRNAs. This chapter highlights the basic guidelines for designing divergent and full-length primers for PCR amplification and Sanger sequencing to validate circRNA sequences.

Posttranscriptional regulation of insulin family ligands and receptors.

Insulin system including ligands (insulin and IGFs) and their shared receptors (IR and IGFR) are critical regulators of insulin signaling and glucose homeostasis. Altered insulin system is associated with major pathological conditions like diabetes and cancer. The mRNAs encoding for these ligands and their receptors are posttranscriptionally controlled by three major groups of regulators; (i) alternative splicing regulatory factors; (ii) turnover and translation regulator RNA-binding proteins (TTR-RBPs); and (iii) non-coding RNAs including miRNAs and long non-coding RNAs (lncRNAs). In this review, we discuss the influence of these regulators on alternative splicing, mRNA stability and translation. Due to the pathological impacts of insulin system, we also discussed the possibilities of discovering new potential regulators which will improve understanding of insulin system and associated diseases.

Regulation of microRNA by circular RNA.

Circular (circ)RNAs have emerged as novel regulators of gene expression through various mechanisms. However, most publications focus on functional circRNAs regulating target gene expression by interacting with micro (mi)RNAs and acting as competing endogenous RNAs (ceRNAs). Although the theory of miRNA sponging by ceRNAs suggests the inhibition of miRNA activity, many studies are biased toward the selection of miRNAs showing a reverse expression pattern compared with circRNA expression. Although several computational tools and molecular assays have been used to predict and validate the interaction of miRNAs with circRNAs, the actual validation of functional in vivo interactions needs careful consideration of molecular experiments with specific controls. As extensive research is being performed on circRNA, many questions arise on the functional significance of circRNA-miRNA interactions. We hope the critical discussion on the criteria for selecting circRNA-miRNA pairs for functional analysis and providing standard methods for validating circRNA-miRNA interactions will advance our understanding of circRNAs as novel gene regulators. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs Translation > Regulation RNA Methods > RNA Analyses in Cells.

Identification of potential proteins translated from circular RNA splice variants.

Circular RNAs (circRNAs) are covalently closed RNA molecules generated from precursor RNAs by the head-to-tail backsplicing of exons. Hundreds of studies demonstrated that circRNAs are ubiquitously expressed and regulate cellular events by modulating microRNA (miRNA) and RNA-binding protein (RBP) activities. A few circRNAs are also known to translate into functional polypeptides regulating cellular physiology. All these functions primarily depend on the full-length sequence of the circRNAs. CircRNA backsplice junction sequence is the key to identifying circRNAs and their full-length mature sequence. However, some multi-exonic circRNAs exist in different isoforms sharing identical backsplice junction sequences and are termed circRNA splice variants. Here, we analyzed the previously published HeLa cell RNA-seq datasets to identify circRNA splice variants using the de novo module of the CIRCexplorer2 circRNA annotation pipeline. A subset of circRNAs with splice variants was validated by the circRNA-rolling circle amplification (circRNA-RCA) method. Interestingly, several validated circRNAs were predicted to translate into proteins by the riboCIRC database. Furthermore, polyribosome fractionation followed by quantitative PCR confirmed the association of a subset of circRNAs with polyribosome supporting their protein-coding potential. Finally, bioinformatics analysis of proteins derived from splice variants of circCORO1C and circASPH suggested altered protein sequences and structures that could affect their physiological functions. Together, our study identified novel circRNA splice variants and their potential translation into protein isoforms which may regulate various physiological processes.

Glucose-stimulated translation regulation of insulin by the 5' UTR-binding proteins.

Insulin is the key regulator of glucose homeostasis in mammals, and glucose-stimulated insulin biosynthesis is essential for maintaining glucose levels in a narrow range in mammals. Glucose specifically promotes the translation of insulin in pancreatic ß-islet, and the untranslated regions of insulin mRNA play a role in such regulation. Specific factors in the ß-islets bind to the insulin 5' UTR and regulate its translation. In the present study we identify protein-disulfide isomerase (PDI) as a key regulator of glucose-stimulated insulin biosynthesis. We show that both in vitro and in vivo PDI can specifically associate with the 5' UTR of insulin mRNA. Immunodepletion of PDI from the islet extract results in loss of glucose-stimulated translation indicating a critical role for PDI in insulin biosynthesis. Similarly, transient overexpression of PDI resulted in specific translation activation by glucose. We show that the RNA binding activity of PDI is mediated through PABP. PDI catalyzes the reduction of the PABP disulfide bond resulting in specific binding of PABP to the insulin 5' UTR. We also show that glucose stimulation of the islets results in activation of a specific kinase that can phosphorylate PDI. These findings identify PDI and PABP as important players in glucose homeostasis.