scanMiR: a biochemically-based toolkit for versatile and efficient microRNA target prediction

MicroRNA-138-5p suppresses excitatory synaptic strength at the cerebellar input layer

MicroRNAs are small, highly conserved non-coding RNAs that negatively regulate mRNA translation and stability. In the brain, miRNAs contribute to neuronal development, synaptogenesis, and synaptic plasticity. MicroRNA 138-5p (miR-138-5p) controls inhibitory synaptic transmission in the hippocampus and is highly expressed in cerebellar excitatory neurons. However, its specific role in cerebellar synaptic transmission remains unknown. Here, we investigated excitatory transmission in the cerebellum of mice expressing a sponge construct that sequesters endogenous miR-138-5p. Mossy fibre stimulation-evoked EPSCs in granule cells were ∼40% larger in miR-138-5p sponge mice compared to controls. Furthermore, we observed larger miniature EPSC amplitudes, suggesting an increased number of functional postsynaptic AMPA receptors. High-frequency train stimulation revealed enhanced short-term depression following miR-138-5p downregulation. Together with computational modelling, this suggests a negative regulation of presynaptic release probability. Overall, our results demonstrate that miR-138-5p suppresses synaptic strength through pre- and postsynaptic mechanisms, providing a potentially powerful mechanism for tuning excitatory synaptic input into the cerebellum. KEY POINTS: MicroRNAs are powerful regulators of mRNA translation and control key cell biological processes including synaptic transmission, but their role in regulating synaptic function in the cerebellum has remained elusive. In this study, we investigated how microRNA-138-5p (miR-138-5p) modulates excitatory transmission at adult murine cerebellar mossy fibre to granule cell synapses. Downregulation of miR-138-5p enhances excitatory synaptic strength at the cerebellar input layer and increases short-term depression. miR-138-5p exerts its regulatory function through both pre- and postsynaptic mechanisms by negatively regulating release probability at mossy fibre boutons, as well as functional AMPA receptor numbers in granule cells. These findings provide insights into the role of miR-138-5p in the cerebellum and expand our understanding of microRNA-dependent control of excitatory synaptic transmission and short-term plasticity.

microRNA-218-5p coordinates scaling of excitatory and inhibitory synapses during homeostatic synaptic plasticity

Homeostatic synaptic plasticity (HSP) is a neuronal mechanism that allows networks to compensate for prolonged changes in activity by adjusting synaptic strength. This process is crucial for maintaining stable brain function and has been implicated in memory consolidation during sleep. While scaling of both excitatory and inhibitory synapses plays an important role during homeostatic synaptic plasticity, molecules coordinating these processes are unknown. In this study, we investigate the role of miR-218-5p as a regulator of inhibitory and excitatory synapses in the context of picrotoxin (PTX)-induced homeostatic synaptic downscaling (HSD) in rat hippocampal neurons. Using enrichment analysis of microRNA-binding sites in genes changing upon PTX-induced HSD, we bioinformatically predict and experimentally validate increased miR-218-5p activity upon PTX treatment. By electrophysiological recordings and confocal microscopy, we demonstrate that inhibiting miR-218-5p activity exerts a dual effect during HSD: It occludes the downscaling of excitatory synapses and dendritic spines, while at the same time attenuating inhibitory synapse upscaling. Furthermore, we identify the Neuroligin2 interacting molecule Mdga1 as a direct miR-218-5p target which potentially mediates the effect of miR-218-5p on homeostatic upscaling of inhibitory synapses. By performing long-term electroencephalographic recordings, we further reveal that local inhibition of miR-218-5p in the somatosensory cortex reduces local slow-wave activity during non-rapid-eye-movement sleep. In summary, this study uncovers miR-218-5p as a key player in coordinating inhibitory and excitatory synapses during homeostatic plasticity and sleep. Our findings contribute to a deeper understanding of how neural circuits maintain stability in the face of activity-induced perturbations, with implications for pathophysiology.

A functional screen uncovers circular RNAs regulating excitatory synaptogenesis in hippocampal neurons

Circular RNAs (circRNAs) are an expanding class of largely unexplored RNAs which are prominently enriched in the mammalian brain. Here, we systematically interrogate their role in excitatory synaptogenesis of rat hippocampal neurons using RNA interference. Thereby, we identify seven circRNAs as negative regulators of excitatory synapse formation, many of which contain high-affinity microRNA binding sites. Knockdown of one of these candidates, circRERE, promotes the formation of electrophysiologically silent synapses. Mechanistically, circRERE knockdown results in a preferential upregulation of synaptic mRNAs containing binding sites for miR-128-3p. Overexpression of circRERE stabilizes miR-128-3p and rescues exaggerated synapse formation upon circRERE knockdown in a miR-128-3p binding site-specific manner. Overall, our results uncover circRERE-mediated stabilization of miR-128-3p as a means to restrict the formation of silent excitatory synaptic co-clusters and more generally implicate circRNA-dependent microRNA regulation in the control of synapse development and function.

Large-scale single-nuclei profiling identifies role for ATRNL1 in atrial fibrillation

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans, yet the molecular basis of AF remains incompletely understood. To determine the cell type-specific transcriptional changes underlying AF, we perform single-nucleus RNA-seq (snRNA-seq) on left atrial (LA) samples from patients with AF and controls. From more than 175,000 nuclei we find that only cardiomyocytes (CMs) and macrophages (MΦs) have a significant number of differentially expressed genes in patients with AF. Attractin Like 1 (ATRNL1) was overexpressed in CMs among patients with AF and localized to the intercalated disks. Further, in both knockdown and overexpression experiments we identify a potent role for ATRNL1 in cell stress response, and in the modulation of the cardiac action potential. Finally, we detect an unexpected expression pattern for a leading AF candidate gene, KCNN3. In sum, we uncover a role for ATRNL1 which may serve as potential therapeutic target for this common arrhythmia.

Semen sEV tRF-Based Models Increase Non-Invasive Prediction Accuracy of Clinically Significant Prostate Cancer among Patients with Moderately Altered PSA Levels

PSA screening has led to an over-diagnosis of prostate cancer (PCa) and unnecessary biopsies of benign conditions due to its low cancer specificity. Consequently, more accurate, preferentially non-invasive, tests are needed. We aim to evaluate the potential of semen sEV (small extracellular vesicles) tsRNAs (tRNA-derived small RNAs) as PCa indicators. Initially, following a literature review in the OncotRF database and high-throughput small RNA-sequencing studies in PCa tissue together with the sncRNA profile in semen sEVs, we selected four candidate 5'tRF tsRNAs for validation as PCa biomarkers. RT-qPCR analysis in semen sEVs from men with moderately elevated serum PSA levels successfully shows that the differential expression of the four tRFs between PCa and healthy control groups can be detected in a non-invasive manner. The combined model incorporating PSA and specific tRFs (5'-tRNA-Glu-TTC-9-1_L30 and 5'-tRNA-Val-CAC-3-1_L30) achieved high predictive accuracy in identifying samples with a Gleason score ≥ 7 and staging disease beyond IIA, supporting that the 5'tRF fingerprint in semen sEV can improve the PSA predictive value to discriminate between malignant and indolent prostate conditions. The in silico study allowed us to map target genes for the four 5'tRFs possibly involved in PCa. Our findings highlight the synergistic use of multiple biomarkers as an efficient approach to improve PCa screening and prognosis.

Deletion of miR-146a enhances therapeutic protein restoration in model of dystrophin exon skipping

Duchenne muscular dystrophy (DMD) is a progressive muscle disease caused by the absence of dystrophin protein. One current DMD therapeutic strategy, exon skipping, produces a truncated dystrophin isoform using phosphorodiamidate morpholino oligomers (PMOs). However, the potential of exon skipping therapeutics has not been fully realized as increases in dystrophin protein have been minimal in clinical trials. Here, we investigate how miR-146a-5p, which is highly elevated in dystrophic muscle, impacts dystrophin protein levels. We find inflammation strongly induces miR-146a in dystrophic, but not wild-type myotubes. Bioinformatics analysis reveals that the dystrophin 3' UTR harbors a miR-146a binding site, and subsequent luciferase assays demonstrate miR-146a binding inhibits dystrophin translation. In dystrophin-null mdx52 mice, co-injection of miR-146a reduces dystrophin restoration by an exon 51 skipping PMO. To directly investigate how miR-146a impacts therapeutic dystrophin rescue, we generated mdx52 with body-wide miR-146a deletion (146aX). Administration of an exon skipping PMO via intramuscular or intravenous injection markedly increases dystrophin protein levels in 146aX vs. mdx52 muscles while skipped dystrophin transcript levels are unchanged supporting a post-transcriptional mechanism of action. Together, these data show that miR-146a expression opposes therapeutic dystrophin restoration, suggesting miR-146a inhibition warrants further research as a potential DMD exon skipping co-therapy.

microRNA-1 Regulates Metabolic Flexibility in Skeletal Muscle via Pyruvate Metabolism

MicroRNA-1 (miR-1) is the most abundant miRNA in adult skeletal muscle. To determine the function of miR-1 in adult skeletal muscle, we generated an inducible, skeletal muscle-specific miR-1 knockout (KO) mouse. Integration of RNA-sequencing (RNA-seq) data from miR-1 KO muscle with Argonaute 2 enhanced crosslinking and immunoprecipitation sequencing (AGO2 eCLIP-seq) from human skeletal muscle identified miR-1 target genes involved with glycolysis and pyruvate metabolism. The loss of miR-1 in skeletal muscle induced cancer-like metabolic reprogramming, as shown by higher pyruvate kinase muscle isozyme M2 (PKM2) protein levels, which promoted glycolysis. Comprehensive bioenergetic and metabolic phenotyping combined with skeletal muscle proteomics and metabolomics further demonstrated that miR-1 KO induced metabolic inflexibility as a result of pyruvate oxidation resistance. While the genetic loss of miR-1 reduced endurance exercise performance in mice and in C. elegans, the physiological down-regulation of miR-1 expression in response to a hypertrophic stimulus in both humans and mice causes a similar metabolic reprogramming that supports muscle cell growth. Taken together, these data identify a novel post-translational mechanism of adult skeletal muscle metabolism regulation mediated by miR-1.

miRNA-mediated inhibition of an actomyosin network in hippocampal pyramidal neurons restricts sociability in adult male mice

Social deficits are frequently observed in patients suffering from neurodevelopmental disorders, but the molecular mechanisms regulating sociability are still poorly understood. We recently reported that the loss of the microRNA (miRNA) cluster miR-379-410 leads to hypersocial behavior and anxiety in mice. Here, we show that ablating miR-379-410 in excitatory neurons of the postnatal mouse hippocampus recapitulates hypersociability, but not anxiety. At the cellular level, miR-379-410 loss in excitatory neurons leads to larger dendritic spines, increased excitatory synaptic transmission, and upregulation of an actomyosin gene network. Re-expression of three cluster miRNAs, as well as pharmacological inhibition of the actomyosin activator ROCK, is sufficient to reinstate normal sociability in miR-379-410 knockout mice. Several actomyosin genes and miR-379-410 family members are reciprocally dysregulated in isogenic human induced pluripotent stem cell (iPSC)-derived neurons harboring a deletion present in patients with Williams-Beuren syndrome, characterized by hypersocial behavior. Together, our results show an miRNA-actomyosin pathway involved in social behavior regulation.

Influence of RNA circularity on Target RNA-Directed MicroRNA Degradation

A subset of circular RNAs (circRNAs) and linear RNAs have been proposed to 'sponge' or block microRNA activity. Additionally, certain RNAs induce microRNA destruction through the process of Target RNA-Directed MicroRNA Degradation (TDMD), but whether both linear and circular transcripts are equivalent in driving TDMD is unknown. Here, we studied whether circular/linear topology of endogenous and artificial RNA targets affects TDMD. Consistent with previous knowledge that Cdr1as (ciRS-7) circular RNA protects miR-7 from Cyrano-mediated TDMD, we demonstrate that depletion of Cdr1as reduces miR-7 abundance. In contrast, overexpression of an artificial linear version of Cdr1as drives miR-7 degradation. Using plasmids that express a circRNA with minimal co-expressed cognate linear RNA, we show differential effects on TDMD that cannot be attributed to the nucleotide sequence, as the TDMD properties of a sequence often differ when in a circular versus linear form. By analysing RNA sequencing data of a neuron differentiation system, we further detect potential effects of circRNAs on microRNA stability. Our results support the view that RNA circularity influences TDMD, either enhancing or inhibiting it on specific microRNAs.

Cirscan: a shiny application to identify differentially active sponge mechanisms and visualize circRNA-miRNA-mRNA networks

BACKGROUND

Non-coding RNAs represent a large part of the human transcriptome and have been shown to play an important role in disease such as cancer. However, their biological functions are still incompletely understood. Among non-coding RNAs, circular RNAs (circRNAs) have recently been identified for their microRNA (miRNA) sponge function which allows them to modulate the expression of miRNA target genes by taking on the role of competitive endogenous RNAs (ce-circRNAs). Today, most computational tools are not adapted to the search for ce-circRNAs or have not been developed for the search for ce-circRNAs from user's transcriptomic data.

RESULTS

In this study, we present Cirscan (CIRcular RNA Sponge CANdidates), an interactive Shiny application that automatically infers circRNA-miRNA-mRNA networks from human multi-level transcript expression data from two biological conditions (e.g. tumor versus normal conditions in the case of cancer study) in order to identify on a large scale, potential sponge mechanisms active in a specific condition. Cirscan ranks each circRNA-miRNA-mRNA subnetwork according to a sponge score that integrates multiple criteria based on interaction reliability and expression level. Finally, the top ranked sponge mechanisms can be visualized as networks and an enrichment analysis is performed to help its biological interpretation. We showed on two real case studies that Cirscan is capable of retrieving sponge mechanisms previously described, as well as identifying potential novel circRNA sponge candidates.

CONCLUSIONS

Cirscan can be considered as a companion tool for biologists, facilitating their ability to prioritize sponge mechanisms for experimental validations and identifying potential therapeutic targets. Cirscan is implemented in R, released under the license GPL-3 and accessible on GitLab ( https://gitlab.com/geobioinfo/cirscan_Rshiny ). The scripts used in this paper are also provided on Gitlab ( https://gitlab.com/geobioinfo/cirscan_paper ).

Transcription readthrough is prevalent in healthy human tissues and associated with inherent genomic features

Transcription termination is a crucial step in the production of conforming mRNAs and functional proteins. Under cellular stress conditions, the transcription machinery fails to identify the termination site and continues transcribing beyond gene boundaries, a phenomenon designated as transcription readthrough. However, the prevalence and impact of this phenomenon in healthy human tissues remain unexplored. Here, we assessed transcription readthrough in almost 3000 transcriptome profiles representing 23 human tissues and found that 34% of the expressed protein-coding genes produced readthrough transcripts. The production of readthrough transcripts was restricted in genomic regions with high transcriptional activity and was associated with inefficient splicing and increased chromatin accessibility in terminal regions. In addition, we showed that these transcripts contained several binding sites for the same miRNA, unravelling a potential role as miRNA sponges. Overall, this work provides evidence that transcription readthrough is pervasive and non-stochastic, not only in abnormal conditions but also in healthy tissues. This suggests a potential role for such transcripts in modulating normal cellular functions.

SeedMatchR: identify off-target effects mediated by siRNA seed regions in RNA-seq experiments

MOTIVATION

On-target gene knockdown, using siRNA, ideally results from binding fully complementary regions in mRNA transcripts to induce direct cleavage. Off-target siRNA gene knockdown can occur through several modes, one being a seed-mediated mechanism mimicking miRNA gene regulation. Seed-mediated off-target effects occur when the ∼8 nucleotides at the 5' end of the guide strand, called a seed region, bind the 3' untranslated regions of mRNA, causing reduced translation. Experiments using siRNA knockdown paired with RNA-seq can be used to detect siRNA sequences with off-target effects driven by the seed region. However, there are limited computational tools designed specifically for detecting siRNA off-target effects mediated by the seed region in differential gene expression experiments.

RESULTS

SeedMatchR is an R package developed to provide users a single, unified resource for detecting and visualizing seed-mediated off-target effects of siRNA using RNA-seq experiments. SeedMatchR is designed to extend current differential expression analysis tools, such as DESeq2, by annotating results with predicted seed matches. Using publicly available data, we demonstrate the ability of SeedMatchR to detect cumulative changes in differential gene expression attributed to siRNA seed region activity.

AVAILABILITY

SeedMatchR is available on CRAN. Documentation and example workflows are available through the SeedMatchR GitHub page at https://github.com/tacazares/SeedMatchR.

Screening and Characterization of Functional circRNAs in Neuronal Cultures

This chapter describes a methodology for the screening and characterization of functional circRNAs, particularly in the context of neural circuit development. Taking advantage of a primary rat neuron culture model of synaptogenesis, we propose a means of selecting from the plethora of circRNA species based on their expression levels, dendritic localization, conservation, and activity regulation. These candidates are then knocked down with RNAi approaches in a functional screen for their potential role in the formation and maturation of excitatory synapses.Upon identification of top candidates regulating synaptogenesis, we tie together different "Omics" approaches to explore the molecular mechanisms underlying the phenotypes observed upon circRNA knockdown. We utilized our EnrichMir algorithm to identify overrepresented miRNA binding sites in differentially expressed genes from polyA-RNA-seq following circRNA knockdown. Furthermore, our ScanMiR web tool allows for the miRNA binding prediction of reconstructed internal circular RNA sequences. Small-RNA sequencing is used to monitor changes in miRNA levels in the circRNA knockdown to complement results obtained from EnrichMiR. Finally, the experimental validation of promising miRNA-circRNA pairs sets the stage for in-depth biochemical exploration of the circRNA interactome and mechanism of action.

Efficient 3'-pairing renders microRNA targeting less sensitive to mRNA seed accessibility

MicroRNAs (miRNAs) are short RNAs that post-transcriptionally regulate gene expression by binding to specific sites in mRNAs. Site recognition is primarily mediated by the seed region (nucleotides g2-g8 in the miRNA), but pairing beyond the seed (3'-pairing) is important for some miRNA:target interactions. Here, we use SHAPE, luciferase reporter assays and transcriptomics analyses to study the combined effect of 3'-pairing and secondary structures in mRNAs on repression efficiency. Using the interaction between miR-34a and its SIRT1 binding site as a model, we provide structural and functional evidence that 3'-pairing can compensate for low seed-binding site accessibility, enabling repression of sites that would otherwise be ineffective. We show that miRNA 3'-pairing regions can productively base-pair with nucleotides far upstream of the seed-binding site and that both hairpins and unstructured bulges within the target site are tolerated. We use SHAPE to show that sequences that overcome inaccessible seed-binding sites by strong 3'-pairing adopt the predicted structures and corroborate the model using luciferase assays and high-throughput modelling of 8177 3'-UTR targets for six miRNAs. Finally, we demonstrate that PHB2, a target of miR-141, is an inaccessible target rescued by efficient 3'-pairing. We propose that these results could refine predictions of effective target sites.

Deletion of miR-146a enhances therapeutic protein restoration in model of dystrophin exon skipping

Duchenne muscular dystrophy (DMD) is a progressive muscle disease caused by the absence of dystrophin protein. One current DMD therapeutic strategy, exon skipping, produces a truncated dystrophin isoform using phosphorodiamidate morpholino oligomers (PMOs). However, the potential of exon skipping therapeutics has not been fully realized as increases in dystrophin protein have been minimal in clinical trials. Here, we investigate how miR-146a-5p, which is highly elevated in dystrophic muscle, impacts dystrophin protein levels. We find inflammation strongly induces miR-146a in dystrophic, but not wild-type myotubes. Bioinformatics analysis reveals that the dystrophin 3'UTR harbors a miR-146a binding site, and subsequent luciferase assays demonstrate miR-146a binding inhibits dystrophin translation. In dystrophin-null mdx52 mice, co-injection of miR-146a reduces dystrophin restoration by an exon 51 skipping PMO. To directly investigate how miR-146a impacts therapeutic dystrophin rescue, we generated mdx52 with body-wide miR-146a deletion (146aX). Administration of an exon skipping PMO via intramuscular or intravenous injection markedly increases dystrophin protein levels in 146aX versus mdx52 muscles; skipped dystrophin transcript levels are unchanged, suggesting a post-transcriptional mechanism-of-action. Together, these data show that miR-146a expression opposes therapeutic dystrophin restoration, suggesting miR-146a inhibition warrants further research as a potential DMD exon skipping co-therapy.

Full-length circular RNA profiling by nanopore sequencing with CIRI-long

Circular RNAs (circRNAs) have important roles in regulating developmental processes and disease progression. As most circRNA sequences are highly similar to their cognate linear transcripts, the current short-read sequencing-based methods rely on the back-spliced junction signal for distinguishing circular and linear reads, which does not allow circRNAs' full-length structure to be effectively reconstructed. Here we describe a long-read sequencing-based protocol, CIRI-long, for the detection of full-length circular RNAs. The CIRI-long protocol combines rolling circular reverse transcription and nanopore sequencing to capture full-length circRNA sequences. After poly(A) tailing, RNase R treatment, and size selection of polymerase chain reaction products, CIRI-long achieves an increased percentage (6%) of circular reads in the constructed library, which is 20-fold higher compared with previous Illumina-based strategies. This method can be applied in cell lines or tissue samples, enabling accurate detection of full-length circRNAs in the range of 100-3,000 bp. The entire protocol can be completed in 1 d, and can be scaled up for large-scale analysis using the nanopore barcoding kit and PromethION sequencing device. CIRI-long can serve as an effective and user-friendly protocol for characterizing full-length circRNAs, generating direct and convincing evidence for the existence of detected circRNAs. The analytical pipeline offers convenient functions for identification of full-length circRNA isoforms and integration of multiple datasets. The assembled full-length transcripts and their splicing patterns provide indispensable information to explore the biological function of circRNAs.

Regulation of the Key Epithelial Cancer Suppressor miR-124 Function by Competing Endogenous RNAs

A decrease in the miR-124 expression was observed in various epithelial cancers. Like a classical suppressor, miR-124 can inhibit the translation of multiple oncogenic proteins. Epigenetic mechanisms play a significant role in the regulation of miR-124 expression and involve hypermethylation of the MIR-124-1/-2/-3 genes and the effects of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) according to the model of competing endogenous RNAs (ceRNAs). More than 40 interactomes (lncRNA/miR-124/mRNA) based on competition between lncRNAs and mRNAs for miR-124 binding have been identified in various epithelial cancers. LncRNAs MALAT1, NEAT1, HOXA11-AS, and XIST are the most represented in these axes. Fourteen axes (e.g., SND1-IT1/miR-124/COL4A1) are involved in EMT and/or metastasis. Moreover, eight axes (e.g., OIP5-AS1/miR-124-5p/IDH2) are involved in key pathways, such as Wnt/b-catenin, E2F1, TGF-β, SMAD, ERK/MAPK, HIF-1α, Notch, PI3K/Akt signaling, and cancer cell stemness. Additionally, 15 axes impaired patient survival and three axes reduced chemo- or radiosensitivity. To date, 14 cases of miR-124 regulation by circRNAs have been identified. Half of them involve circHIPK3, which belongs to the exonic ecircRNAs and stimulates cell proliferation, EMT, autophagy, angiogenesis, and multidrug resistance. Thus, miR-124 and its interacting partners may be considered promising targets for cancer therapy.

enrichMiR predicts functionally relevant microRNAs based on target collections

MicroRNAs (miRNAs) are small non-coding RNAs that are among the main post-transcriptional regulators of gene expression. A number of data collections and prediction tools have gathered putative or confirmed targets of these regulators. It is often useful, for discovery and validation, to harness such collections to perform target enrichment analysis in given transcriptional signatures or gene-sets in order to predict involved miRNAs. While several methods have been proposed to this end, a flexible and user-friendly interface for such analyses using various approaches and collections is lacking. enrichMiR (https://ethz-ins.org/enrichMiR/) addresses this gap by enabling users to perform a series of enrichment tests, based on several target collections, to rank miRNAs according to their likely involvement in the control of a given transcriptional signature or gene-set. enrichMiR results can furthermore be visualised through interactive and publication-ready plots. To guide the choice of the appropriate analysis method, we benchmarked various tests across a panel of experiments involving the perturbation of known miRNAs. Finally, we showcase enrichMiR functionalities in a pair of use cases.