miRNA-target chimeras reveal miRNA 3'-end pairing as a major determinant of Argonaute target specificity

The dynamics of t1 adenosine binding on human Argonaute 2: Understanding recognition with conformational selection

The control of expression in genetic regulation is a fundamental process for cell life. In RNA-mediated silencing, human Argonaute-2 protein (hAgo2) uses sequence information encoded in small RNAs (guide) to identify complementary sites in messenger RNAs (target) for repression. The specificity of this molecular recognition lies at the basis of the mechanisms that control the expression of thousands of genes, which necessarily requires a fine tuning of complex events. Among these, the binding of the first nucleotide of the target RNA (t1) is emerging as an important modulator of hAgo2-mediated machinery. Using atomistic molecular dynamics-derived analyses, we address the mechanism behind t1-dependent regulation and study the impact of different t1 nucleotides (t1A, t1C, t1G, t1U) on the conformational dynamics of both hAgo2 and guide-target RNAs. Only when an adenine is found at this position, t1 directly interacts with a specific hAgo2 binding pocket, favoring the stabilization of target binding. Our findings show that hAgo2 exploits a dynamic recognition mechanism of the t1-target thanks to a modulation of RNA conformations. Here, t1-adenine is the only nucleobase endowed with a dual binding mode: a T-shape and a co-planar conformation, respectively, orthogonal and parallel to the following base-pairs of guide-target duplex. This triggers a composite set of molecular interactions that stabilizes distinctive conformational ensembles. Our comparative analyses show characteristic traits of local and global dynamic interplay between hAgo2 and the RNA molecules and highlight how t1A binding acts as a molecular switch for target recognition and complex stabilization. Implications for future mechanistic studies are discussed.

Single-molecule FRET uncovers hidden conformations and dynamics of human Argonaute 2

Human Argonaute 2 (hAgo2) constitutes the functional core of the RNA interference pathway. Guide RNAs direct hAgo2 to target mRNAs, which ultimately leads to hAgo2-mediated mRNA degradation or translational inhibition. Here, we combine site-specifically labeled hAgo2 with time-resolved single-molecule FRET measurements to monitor conformational states and dynamics of hAgo2 and hAgo2-RNA complexes in solution that remained elusive so far. We observe dynamic anchoring and release of the guide’s 3’-end from the PAZ domain during the stepwise target loading process even with a fully complementary target. We find differences in structure and dynamic behavior between partially and fully paired canonical hAgo2-guide/target complexes and the miRNA processing complex formed by hAgo2 and pre-miRNA451. Furthermore, we detect a hitherto unknown conformation of hAgo2-guide/target complexes that poises them for target-directed miRNA degradation. Taken together, our results show how the conformational flexibility of hAgo2-RNA complexes determines function and the fate of the ribonucleoprotein particle.

Single-molecule FRET measurements provide detailed insights into the conformational states and dynamics of human Argonaute 2 that are required for its function at the core of the eukaryotic RNA silencing pathway.

A deep learning method for miRNA/isomiR target detection

Accurate identification of microRNA (miRNA) targets at base-pair resolution has been an open problem for over a decade. The recent discovery of miRNA isoforms (isomiRs) adds more complexity to this problem. Despite the existence of many methods, none considers isomiRs, and their performance is still suboptimal. We hypothesize that by taking the isomiR-mRNA interactions into account and applying a deep learning model to study miRNA-mRNA interaction features, we may improve the accuracy of miRNA target predictions. We developed a deep learning tool called DMISO to capture the intricate features of miRNA/isomiR-mRNA interactions. Based on tenfold cross-validation, DMISO showed high precision (95%) and recall (90%). Evaluated on three independent datasets, DMISO had superior performance to five tools, including three popular conventional tools and two recently developed deep learning-based tools. By applying two popular feature interpretation strategies, we demonstrated the importance of the miRNA regions other than their seeds and the potential contribution of the RNA-binding motifs within miRNAs/isomiRs and mRNAs to the miRNA/isomiR-mRNA interactions.

Modulation of miR-181 influences dopaminergic neuronal degeneration in a mouse model of Parkinson's disease

Parkinson's disease (PD) is caused by the loss of dopaminergic (DA) neurons in the substantia nigra (SN). Although PD pathogenesis is not fully understood, studies implicate perturbations in gene regulation, mitochondrial function, and neuronal activity. MicroRNAs (miRs) are small gene regulatory RNAs that inhibit diverse subsets of target mRNAs, and several studies have noted miR expression alterations in PD brains. For example, miR-181a is abundant in the brain and is increased in PD patient brain samples; however, the disease relevance of this remains unclear. Here, we show that miR-181 target mRNAs are broadly downregulated in aging and PD brains. To address whether the miR-181 family plays a role in PD pathogenesis, we generated adeno-associated viruses (AAVs) to overexpress and inhibit the miR-181 isoforms. After co-injection with AAV overexpressing alpha-synuclein (aSyn) into mouse SN (PD model), we found that moderate miR-181a/b overexpression exacerbated aSyn-induced DA neuronal loss, whereas miR-181 inhibition was neuroprotective relative to controls (GFP alone and/or scrambled RNA). Also, prolonged miR-181 overexpression in SN alone elicited measurable neurotoxicity that is coincident with an increased immune response. mRNA-seq analyses revealed that miR-181a/b inhibits genes involved in synaptic transmission, neurite outgrowth, and mitochondrial respiration, along with several genes having known protective roles and genetic links in PD.

Target binding triggers hierarchical phosphorylation of human Argonaute-2 to promote target release

Argonaute (Ago) proteins play a central role in post-transcriptional gene regulation through RNA interference (RNAi). Agos bind small RNAs (sRNAs) including small interfering RNAs (siRNAs) and microRNAs (miRNAs) to form the functional core of the RNA-induced silencing complex (RISC). The sRNA is used as a guide to target mRNAs containing either partially or fully complementary sequences, ultimately leading to downregulation of the corresponding proteins. It was previously shown that the kinase CK1α phosphorylates a cluster of residues in the eukaryotic insertion (EI) of Ago, leading to the alleviation of miRNA-mediated repression through an undetermined mechanism. We show that binding of miRNA-loaded human Ago2 to target RNA with complementarity to the seed and 3’ supplementary regions of the miRNA primes the EI for hierarchical phosphorylation by CK1α. The added negative charges electrostatically promote target release, freeing Ago to seek out additional targets once it is dephosphorylated. The high conservation of potential phosphosites in the EI suggests that such a regulatory strategy may be a shared mechanism for regulating miRNA-mediated repression.

Proteins are the chemical ‘workhorses’ of the cell: some help maintain a cell’s shape or structure, while others carry out the chemical reactions necessary for life. Organisms therefore need to keep tight control over the production of proteins in their cells, so that the right amount of each protein is made at the right time, in the right place.

Instructions for making new proteins are encoded in a type of molecule called messenger RNA. Each messenger RNA contains the instructions for one protein, which are then ‘read’ and carried out by special cellular machinery called ribosomes. The cell can control how much protein it produces by regulating both the levels of different messenger RNA and the amount of protein ribosomes are allowed to make from those instructions.

The main way to regulate the levels of messenger RNA is through their transcription from the genome. However, this needs fine tuning. Cells can do this in a highly specific way using molecules called microRNAs. A microRNA works by directing a protein called Argonaute to the messenger RNA that it targets. Once Argonaute arrives, it can call in additional ‘helper proteins’ to shut down, or reduce, protein production from that messenger RNA, or alternatively to break down the messenger RNA altogether.

Cells can use an enzyme called CK1α to attach bulky chemical groups onto a specific part of the Argonaute protein, in a reaction termed phosphorylation. The ability to carry out this reaction (and to reverse it) also seems to be important for microRNAs to do their job properly, but why has remained unknown. Bibel et al. wanted to determine what triggers CK1α to phosphorylate Argonaute, and how this affects interactions between microRNAs, Argonaute and their target messenger RNAs.

A series of ‘test tube’ experiments looked at the interaction between purified CK1α and Argonaute under different conditions. These demonstrated that CK1α could only carry out its phosphorylation reaction when Argonaute was already interacting with a microRNA and its corresponding messenger RNA. Further measurements revealed that phosphorylation of Argonaute made it detach from the messenger RNA more quickly. This suggests that phosphorylation might be a way to let Argonaute seek out new messenger RNAs after blocking protein production at its first ‘target’.

These results shed new light on a fundamental mechanism that cells use to control protein production. Bibel et al. propose that this mechanism may be shared across many different species and could one day help guide the development of new medical therapies based on microRNAs.

microRNA-27b regulates hepatic lipase enzyme LIPC and reduces triglyceride degradation during hepatitis C virus infection

miRNAs are short, noncoding RNAs that negatively and specifically regulate protein expression, the cumulative effects of which can result in broad changes to cell systems and architecture. The miRNA miR-27b is known to regulate lipid regulatory pathways in the human liver and is also induced by the hepatitis C virus (HCV). However, the functional targets of miR-27b are not well established. Herein, an activity-based protein profiling method using a serine hydrolase probe, coupled with stable isotope labeling and mass spectrometry identified direct and indirect targets of miR-27b. The hepatic lipase C (LIPC) stood out as both highly dependent on miR-27b and as a major modulator of lipid pathway misregulation. Modulation of miR-27b using both exogenous miRNA mimics and inhibitors demonstrated that transcription factors Jun, PPARα, and HNF4α, all of which also influence LIPC levels and activity, are regulated by miR-27b. LIPC was furthermore shown to affect the progress of the life cycle of HCV and to decrease levels of intracellular triglycerides, upon which HCV is known to depend. In summary, this work has demonstrated that miR-27b mediates HCV infection by downregulating LIPC, thereby reducing triglyceride degradation, which in turn increases cellular lipid levels.

miR-182-5p Regulates Nogo-A Expression and Promotes Neurite Outgrowth of Hippocampal Neurons In Vitro

Nogo-A protein is a key myelin-associated inhibitor of axonal growth, regeneration, and plasticity in the central nervous system (CNS). Regulation of the Nogo-A/NgR1 pathway facilitates functional recovery and neural repair after spinal cord trauma and ischemic stroke. MicroRNAs are described as effective tools for the regulation of important processes in the CNS, such as neuronal differentiation, neuritogenesis, and plasticity. Our results show that miR-182-5p mimic specifically downregulates the expression of the luciferase reporter gene fused to the mouse Nogo-A 3′UTR, and Nogo-A protein expression in Neuro-2a and C6 cells. Finally, we observed that when rat primary hippocampal neurons are co-cultured with C6 cells transfected with miR-182-5p mimic, there is a promotion of the outgrowth of neuronal neurites in length. From all these data, we suggest that miR-182-5p may be a potential therapeutic tool for the promotion of axonal regeneration in different diseases of the CNS.

CRISPRi for specific inhibition of miRNA clusters and miRNAs with high sequence homology

miRNAs form a class of noncoding RNAs, involved in post-transcriptional regulation of gene expression, broadly studied for their involvement in physiological and pathological context. Inhibition of mature miRNA transcripts, commonly used in miRNA loss-of-function experiments, may not be specific in case of miRNAs with high sequence homology, e.g. miRNAs from the same seed family. Phenotypic effects of miRNA repression might be biased by the repression of highly similar miRNAs. Another challenge is simultaneous inhibition of multiple miRNAs encoded within policistronic clusters, potentially co-regulating common biological processes. To elucidate roles of miRNA clusters and miRNAs with high sequence homology, it is of key importance to selectively repress only the miRNAs of interest. Targeting miRNAs on genomic level with CRISPR/dCas9-based methods is an attractive alternative to blocking mature miRNAs. Yet, so far no clear guidelines on the design of CRISPR inhibition (CRISPRi) experiments, specifically for miRNA repression, have been proposed. To address this need, here we propose a strategy for effective inhibition of miRNAs and miRNA clusters using CRISPRi. We provide clues on how to approach the challenges in using CRISPR/dCas in miRNA studies, which include prediction of miRNA transcription start sites (TSSs) and the design of single guide RNAs (sgRNAs). The strategy implements three TSS prediction online tools, dedicated specifically for miRNAs: miRStart, FANTOM 5 miRNA atlas, DIANA-miRGen, and CRISPOR tool for sgRNAs design; it includes testing and selection of optimal sgRNAs. We demonstrate that compared to siRNA/shRNA-based miRNA silencing, CRISPRi improves the repression specificity for miRNAs with highly similar sequence and contribute to higher uniformity of the effects of silencing the whole miRNA clusters. This strategy may be adapted for CRISPR-mediated activation (CRISPRa) of miRNA expression.

Non-coding RNAs in Kawasaki disease: Molecular mechanisms and clinical implications

Kawasaki disease (KD) is an acute self-limiting vasculitis with coronary complications, usually occurring in children. The incidence of KD in children is increasing year by year, mainly in East Asian countries, but relatively stably in Europe and America. Although studies on KD have been reported, the pathogenesis of KD is unknown. With the development of high-throughput sequencing technology, growing number of regulatory noncoding RNAs (ncRNAs) including microRNA (miRNA), long noncoding RNA (lncRNA), and circular RNA (circRNA) have been identified to involved in KD. However, the role of ncRNAs in KD has not been comprehensively elucidated. Therefore, it is significative to study the regulatory role of ncRNA in KD, which might help to uncover new and effective therapeutic strategies for KD. In this review, we summarize recent studies on ncRNA in KD from the perspectives of immune disorders, inflammatory disorders, and endothelial dysfunction, and highlight the potential of ncRNAs as therapeutic targets for KD.

Circr, a Computational Tool to Identify miRNA:circRNA Associations

Circular RNAs (circRNAs) are known to act as important regulators of the microRNA (miRNA) activity. Yet, computational resources to identify miRNA:circRNA interactions are mostly limited to already annotated circRNAs or affected by high rates of false positive predictions. To overcome these limitations, we developed Circr, a computational tool for the prediction of associations between circRNAs and miRNAs. Circr combines three publicly available algorithms for de novo prediction of miRNA binding sites on target sequences (miRanda, RNAhybrid, and TargetScan) and annotates each identified miRNA:target pairs with experimentally validated miRNA:RNA interactions and binding sites for Argonaute proteins derived from either ChIPseq or CLIPseq data. The combination of multiple tools for the identification of a single miRNA recognition site with experimental data allows to efficiently prioritize candidate miRNA:circRNA interactions for functional studies in different organisms. Circr can use its internal annotation database or custom annotation tables to enhance the identification of novel and not previously annotated miRNA:circRNA sites in virtually any species. Circr is written in Python 3.6 and is released under the GNU GPL3.0 License at https://github.com/bicciatolab/Circr.

Identification of the global miR-130a targetome reveals a role for TBL1XR1 in hematopoietic stem cell self-renewal and t(8;21) AML

Gene expression profiling and proteome analysis of normal and malignant hematopoietic stem cells (HSCs) point to shared core stemness properties. However, discordance between mRNA and protein signatures highlights an important role for post-transcriptional regulation by microRNAs (miRNAs) in governing this critical nexus. Here, we identify miR-130a as a regulator of HSC self-renewal and differentiation. Enforced expression of miR-130a impairs B lymphoid differentiation and expands long-term HSCs. Integration of protein mass spectrometry and chimeric AGO2 crosslinking and immunoprecipitation (CLIP) identifies TBL1XR1 as a primary miR-130a target, whose loss of function phenocopies miR-130a overexpression. Moreover, we report that miR-130a is highly expressed in t(8;21) acute myeloid leukemia (AML), where it is critical for maintaining the oncogenic molecular program mediated by the AML1-ETO complex. Our study establishes that identification of the comprehensive miRNA targetome within primary cells enables discovery of genes and molecular networks underpinning stemness properties of normal and leukemic cells.

Human MicroRNAs Attenuate the Expression of Immediate Early Proteins and HCMV Replication during Lytic and Latent Infection in Connection with Enhancement of Phosphorylated RelA/p65 (Serine 536) That Binds to MIEP

Attenuating the expression of immediate early (IE) proteins is essential for controlling the lytic replication of human cytomegalovirus (HCMV). The human microRNAs (hsa-miRs), miR-200b-3p and miR-200c-3p, have been identified to bind the 3′-untranslated region (3′-UTR) of the mRNA encoding IE proteins. However, whether hsa-miRs can reduce IE72 expression and HCMV viral load or exhibit a crosstalk with the host cellular signaling machinery, most importantly the NF-κB cascade, has not been evaluated. In this study, argonaute-crosslinking and immunoprecipitation-seq revealed that miR-200b-3p and miR-200c-3p bind the 3′-UTR of UL123, which is a gene that encodes IE72. The binding of these miRNAs to the 3′-UTR of UL123 was verified in transfected cells stably expressing GFP. We used miR-200b-3p/miR-200c-3p mimics to counteract the downregulation of these miRNA after acute HCMV infection. This resulted in reduced IE72/IE86 expression and HCMV VL during lytic infection. We determined that IE72/IE86 alone can inhibit the phosphorylation of RelA/p65 at the Ser⁵³⁶ residue and that p-Ser⁵³⁶ RelA/p65 binds to the major IE promoter/enhancer (MIEP). The upregulation of miR-200b-3p and miR-200c-3p resulted in the phosphorylation of RelA/p65 at Ser⁵³⁶ through the downregulation of IE, and the binding of the resultant p-Ser⁵³⁶ RelA/p65 to MIEP resulted in a decreased production of pro-inflammatory cytokines. Overall, miR-200b-3p and miR-200c-3p—together with p-Ser⁵³⁶ RelA/p65—can prevent lytic HCMV replication during acute and latent infection

Prediction and pan-cancer analysis of mammalian transcripts involved in target directed miRNA degradation

It is currently unknown how many RNA transcripts are able to induce degradation of microRNAs (miRNA) via the mechanism known as target-directed miRNA degradation (TDMD). We developed TDMDfinder, a computational pipeline that identifies 'high confidence' TDMD interactions in the Human and Mouse transcriptomes by combining sequence alignment and feature selection approaches. Our predictions suggested that TDMD is widespread, with potentially every miRNA controlled by endogenous targets. We experimentally tested 37 TDMDfinder predictions, of which 17 showed TDMD effects as measured by RT-qPCR and small RNA sequencing, linking the miR-17, miR-19, miR-30, miR-221, miR-26 and miR-23 families to novel endogenous TDMDs. In some cases, TDMD was found to affect different members of the same miRNA family selectively. Features like complementarity to the miRNA 3' region, bulge size and hybridization energy appeared to be the main factors determining sensitivity. Computational analyses performed using the multiomic TCGA platform substantiated the involvement of many TDMD transcripts in human cancer and highlighted 36 highly significant interactions, suggesting TDMD as a new potential oncogenic mechanism. In conclusion, TDMDfinder provides the first inventory of bona fide human and mouse TDMDs. Available as a free webtool, TDMDfinder allows users to search for any TDMD interaction of interest by customizing its selection criteria.

New Insights into Mitochondrial-Nuclear Interactions Revealed through Analysis of Small RNAs

Mitochondrial sequence variants affect phenotypic function, often through interaction with the nuclear genome. These "mitonuclear" interactions have been linked both to evolutionary processes and human health. The study of these interactions has focused on mechanisms regulating communication between mitochondrial and nuclear proteins; the role of mitochondrial (mt) RNAs has received little attention. Here, we show that small mt-RNAs bind to the nuclear protein Argonaute 2, and that nuclear miRNAs bind to mt-mRNAs. We identify one small mt-RNA that binds to Argonaute 2 in human tissues whose expression and sequence remain unchanged across vertebrates. Although analyses of CLEAR-CLIP sequencing data sets of human and mouse did not reveal consistent interactions between small mt-RNAs and nuclear mRNAs, we found that MT-ND4 and MT-ATP6 mRNAs are bound by different nuclear miRNAs in humans and mice. Our work homes in on previously unknown interactions between nuclear and small mt-RNAs, which may play key roles in intergenomic communication.

Structure-Activity Relationships of Anti-microRNA Oligonucleotides Containing Cationic Guanidine-Modified Nucleic Acids

Anti-microRNA oligonucleotides (AMOs) are valuable tools for the treatment of diseases caused by the dysregulation of microRNA expression. However, the correlation between chemical modifications in AMO sequences and the microRNA-inhibitory activity has not been fully elucidated. In this study, we synthesized a series of AMOs containing cationic guanidine-bridged nucleic acids (GuNA) and evaluated their activities using a dual luciferase assay. We also optimized the site of GuNA substitution and found an effective design for the inhibition of microRNA-21, which was partially different from that of conventional nucleic acid derivatives. This study showed that GuNA-substituted AMOs are effective in inhibiting the function of microRNA.

Reexamining assumptions about miRNA-guided gene silencing

MicroRNAs (miRNAs) are short endogenously expressed RNAs that have the potential to regulate the expression of any RNA. This potential has led to the publication of several thousand papers each year connecting miRNAs to many different genes and human diseases. By contrast, relatively few papers appear that investigate the molecular mechanism used by miRNAs. There is a disconnect between rigorous understanding of mechanism and the extraordinary diversity of reported roles for miRNAs. Consequences of this disconnect include confusion about the assumptions underlying the basic science of human miRNAs and slow development of therapeutics that target miRNAs. Here, we present an overview of investigations into miRNAs and their impact on gene expression. Progress in our understanding of miRNAs would be aided by a greater focus on the mechanism of miRNAs and a higher burden of evidence on researchers who seek to link expression of a particular miRNA to a biological phenotype.

A nuclear function for an oncogenic microRNA as a modulator of snRNA and splicing

BACKGROUND

miRNAs are regulatory transcripts established as repressors of mRNA stability and translation that have been functionally implicated in carcinogenesis. miR-10b is one of the key onco-miRs associated with multiple forms of cancer. Malignant gliomas exhibit particularly striking dependence on miR-10b. However, despite the therapeutic potential of miR-10b targeting, this miRNA's poorly investigated and largely unconventional properties hamper the clinical translation.

METHODS

We utilized Covalent Ligation of Endogenous Argonaute-bound RNAs and their high-throughput RNA sequencing to identify miR-10b interactome and a combination of biochemical and imaging approaches for target validation. They included Crosslinking and RNA immunoprecipitation with spliceosomal proteins, a combination of miRNA FISH with protein immunofluorescence in glioma cells and patient-derived tumors, native Northern blotting, and the transcriptome-wide analysis of alternative splicing.

RESULTS

We demonstrate that miR-10b binds to U6 snRNA, a core component of the spliceosomal machinery. We provide evidence of the direct binding between miR-10b and U6, in situ imaging of miR-10b and U6 co-localization in glioma cells and tumors, and biochemical co-isolation of miR-10b with the components of the spliceosome. We further demonstrate that miR-10b modulates U6 N-6-adenosine methylation and pseudouridylation, U6 binding to splicing factors SART3 and PRPF8, and regulates U6 stability, conformation, and levels. These effects on U6 result in global splicing alterations, exemplified by the altered ratio of the isoforms of a small GTPase CDC42, reduced overall CDC42 levels, and downstream CDC42 -mediated effects on cell viability.

CONCLUSIONS

We identified U6 snRNA, the key RNA component of the spliceosome, as the top miR-10b target in glioblastoma. We, therefore, present an unexpected intersection of the miRNA and splicing machineries and a new nuclear function for a major cancer-associated miRNA.

CLASH Analyst: A Web Server to Identify In Vivo RNA–RNA Interactions from CLASH Data

Non-coding RNAs, such as miRNAs and piRNAs, play critical roles in gene regulation through base-pairing interactions with their target molecules. The recent development of the crosslinking, ligation, and sequencing of hybrids (CLASH) method has allowed scientists to map transcriptome-wide RNA–RNA interactions by identifying chimeric reads consisting of fragments from regulatory RNAs and their targets. However, analyzing CLASH data requires scientists to use advanced bioinformatics, and currently available tools are limited for users with little bioinformatic experience. In addition, many published CLASH studies do not show the full scope of RNA–RNA interactions that were captured, highlighting the importance of reanalyzing published data. Here, we present CLASH Analyst, a web server that can analyze raw CLASH data within a fully customizable and easy-to-use interface. CLASH Analyst accepts raw CLASH data as input and identifies the RNA chimeras containing the regulatory and target RNAs according to the user’s interest. Detailed annotation of the captured RNA–RNA interactions is then presented for the user to visualize within the server or download for further analysis. We demonstrate that CLASH Analyst can identify miRNA- and piRNA-targeting sites reported from published CLASH data and should be applicable to analyze other RNA–RNA interactions. CLASH Analyst is freely available for academic use.

miRNA Combinatorics and its Role in Cell State Control-A Probabilistic Approach

A hallmark of cancer evolution is that the tumor may change its cell identity and improve its survival and fitness. Drastic change in microRNA (miRNA) composition and quantities accompany such dynamic processes. Cancer samples are composed of cells’ mixtures of varying stages of cancerous progress. Therefore, cell-specific molecular profiling represents cellular averaging. In this study, we consider the degree to which altering miRNAs composition shifts cell behavior. We used COMICS, an iterative framework that simulates the stochastic events of miRNA-mRNA pairing, using a probabilistic approach. COMICS simulates the likelihood that cells change their transcriptome following many iterations (100 k). Results of COMICS from the human cell line (HeLa) confirmed that most genes are resistant to miRNA regulation. However, COMICS results suggest that the composition of the abundant miRNAs dictates the nature of the cells (across three cell lines) regardless of its actual mRNA steady-state. In silico perturbations of cell lines (i.e., by overexpressing miRNAs) allowed to classify genes according to their sensitivity and resilience to any combination of miRNA perturbations. Our results expose an overlooked quantitative dimension for a set of genes and miRNA regulation in living cells. The immediate implication is that even relatively modest overexpression of specific miRNAs may shift cell identity and impact cancer evolution.

Solid-Support Directional (SSD) RNA-Seq as a Companion Method to CLIP-Seq

CLIP-Seq (Deep Sequencing after in vivo Crosslinking and Immunoprecipitation, HITS-CLIP) has emerged as a key method for the study of RNA-binding proteins (RBPs), as it can scrutinize the RNAs bound by an RBP in vivo, with minimum manipulation of biological samples. CLIP-Seq is best used to reveal changes of the RNA cargo of an RBP and differences on binding patterns of the bound RNAs in living cells in different genetic backgrounds or after experimental treatment, rather than simply identifying RNA species. It is therefore crucial that a reference of the steady state levels of the RNAs present in the samples used for the CLIP-Seq experiment is included in the bioinformatic analysis. A simple directional RNA-Seq method was developed that uses the same oligonucleotides and the same PCR amplification steps as our CLIP-Seq method, which therefore can be analyzed using the same bioinformatic pipeline as the CLIP-Seq data. This greatly simplifies and streamlines the analysis process, and at the same time reduces the chances of protocol-specific artifacts and biases interfering with data interpretation. Some considerations on ways to integrate CLIP-Seq and RNA-Seq analyses are also provided herein.

Global profiling and annotation of templated isomiRs dynamics across Caenorhabditis elegans development

microRNAs (miRNAs) are small non-coding RNAs that regulate gene expression through translational repression and mRNA destabilization. During canonical miRNA biogenesis, several miRNA isoforms, or isomiRs, are produced from a single precursor miRNA. Templated isomiRs are generated through Drosha or Dicer cleavage at alternate positions on either the primary or the precursor miRNAs, generating truncated or extended 5' and/or 3' miRNA ends. As changes to the mature miRNA sequence can alter miRNA gene target repertoire, we investigated the extent of templated isomiR prevalence, providing a profiling map for templated isomiRs across stages of C. elegans development. While most miRNA loci did not produce abundant templated isomiRs, a substantial number of miRNA loci produced isomiRs were just as, or more, abundant than their annotated canonical mature miRNAs. 3' end miRNA alterations were more frequent than the seed-altering 5' end extensions or truncations. However, we identified several miRNA loci that produced a considerable amount of isomiRs with 5' end alterations, predicted to target new, distinct sets of genes. Overall, the presented annotation of templated isomiR dynamics across C. elegans developmental stages provides a basis for further studies into miRNA biogenesis and the intriguing potential of functional miRNA diversification through isomiR production.

Sequencing of Argonaute-bound microRNA/mRNA hybrids reveals regulation of the unfolded protein response by microRNA-320a

MicroRNAs (miRNA) are short non-coding RNAs widely implicated in gene regulation. Most metazoan miRNAs utilize the RNase III enzymes Drosha and Dicer for biogenesis. One notable exception is the RNA polymerase II transcription start sites (TSS) miRNAs whose biogenesis does not require Drosha. The functional importance of the TSS-miRNA biogenesis is uncertain. To better understand the function of TSS-miRNAs, we applied a modified Crosslinking, Ligation, and Sequencing of Hybrids on Argonaute (AGO-qCLASH) to identify the targets for TSS-miRNAs in HCT116 colorectal cancer cells with or without DROSHA knockout. We observed that miR-320a hybrids dominate in TSS-miRNA hybrids identified by AGO-qCLASH. Targets for miR-320a are enriched for the eIF2 signaling pathway, a downstream component of the unfolded protein response. Consistently, in miR-320a mimic- and antagomir- transfected cells, differentially expressed gene products are associated with eIF2 signaling. Within the AGO-qCLASH data, we identified the endoplasmic reticulum (ER) chaperone calnexin as a direct miR-320a down-regulated target, thus connecting miR-320a to the unfolded protein response. During ER stress, but not amino acid deprivation, miR-320a up-regulates ATF4, a critical transcription factor for resolving ER stress. In summary, our study investigates the targetome of the TSS-miRNAs in colorectal cancer cells and establishes miR-320a as a regulator of unfolded protein response.

Widespread microRNA degradation elements in target mRNAs can assist the encoded proteins

Binding of microRNAs (miRNAs) to mRNAs normally results in post-transcriptional repression of gene expression. However, extensive base-pairing between miRNAs and target RNAs can trigger miRNA degradation, a phenomenon called target RNA-directed miRNA degradation (TDMD). Here, we systematically analyzed Argonaute-CLASH (cross-linking, ligation, and sequencing of miRNA-target RNA hybrids) data and identified numerous candidate TDMD triggers, focusing on their ability to induce nontemplated nucleotide addition at the miRNA 3' end. When exogenously expressed in various cell lines, eight triggers induce degradation of corresponding miRNAs. Both the TDMD base-pairing and surrounding sequences are essential for TDMD. CRISPR knockout of endogenous trigger or ZSWIM8, a ubiquitin ligase essential for TDMD, reduced miRNA degradation. Furthermore, degradation of miR-221 and miR-222 by a trigger in BCL2L11, which encodes a proapoptotic protein, enhances apoptosis. Therefore, we uncovered widespread TDMD triggers in target RNAs and demonstrated an example that could functionally cooperate with the encoded protein.

Transgenic overexpression of the miR-200b/200a/429 cluster inhibits mammary tumor initiation

The miR-200 family consists of five members expressed as two clusters: miR-200c/141 cluster and miR-200b/200a/429 cluster. In the mammary gland, miR-200s maintain epithelial identity by decreasing the expression of mesenchymal markers leading to high expression of epithelial markers. While the loss of miR-200s is associated with breast cancer growth and metastasis the impact of miR-200 expression on mammary tumor initiation has not been investigated. Using mammary specific expression of the miR-200b/200a/429 cluster in transgenic mice, we found that elevated expression miR-200s could almost completely prevent mammary tumor development. Only 1 of 16 MTB-IGFIRba429 transgenic mice (expressing both the IGF-IR and miR-200b/200a/429 transgenes) developed a mammary tumor while 100% of MTB-IGFIR transgenic mice (expressing only the IGF-IR transgene) developed mammary tumors. RNA sequencing, qRT-PCR, and immunohistochemistry of mammary tissue from 55-day old mice found Spp1, Saa1, and Saa2 to be elevated in mammary tumors and inhibited by miR-200b/200a/429 overexpression. This study suggests that miR-200s could be used as a preventative strategy to protect women from developing breast cancer. One concern with this approach is the potential negative impact miR-200 overexpression may have on mammary function. However, transgenic overexpression of miR-200s, on their own, did not significantly impact mammary ductal development indicating the miR-200 overexpression should not significantly impact mammary function. Thus, this study provides the initial foundation for using miR-200s for breast cancer prevention and additional studies should be performed to identify strategies for increasing mammary miR-200 expression and determine whether miR-200s can prevent mammary tumor initiation by other genetic alterations.

Modulation of pancreatic cancer cell sensitivity to FOLFIRINOX through microRNA-mediated regulation of DNA damage

FOLFIRINOX, a combination of chemotherapy drugs (Fluorouracil, Oxaliplatin, Irinotecan -FOI), provides the best clinical benefit in pancreatic ductal adenocarcinoma (PDAC) patients. In this study we explore the role of miRNAs (MIR) as modulators of chemosensitivity to identify potential biomarkers of response. We find that 41 and 84 microRNA inhibitors enhance the sensitivity of Capan1 and MiaPaCa2 PDAC cells respectively. These include a MIR1307-inhibitor that we validate in further PDAC cell lines. Chemotherapy-induced apoptosis and DNA damage accumulation are higher in MIR1307 knock-out (MIR1307KO) versus control PDAC cells, while re-expression of MIR1307 in MIR1307KO cells rescues these effects. We identify binding of MIR1307 to CLIC5 mRNA through covalent ligation of endogenous Argonaute-bound RNAs cross-linking immunoprecipitation assay. We validate these findings in an in vivo model with MIR1307 disruption. In a pilot cohort of PDAC patients undergoing FOLFIRONX chemotherapy, circulating MIR1307 correlates with clinical outcome.

A regulatory network of microRNAs confers lineage commitment during early developmental trajectories of B and T lymphocytes

The commitment of hematopoietic multipotent progenitors (MPPs) toward a particular lineage involves activation of cell type-specific genes and silencing of genes that promote alternate cell fates. Although the gene expression programs of early-B and early-T lymphocyte development are mutually exclusive, we show that these cell types exhibit significantly correlated microRNA (miRNA) profiles. However, their corresponding miRNA targetomes are distinct and predominated by transcripts associated with natural killer, dendritic cell, and myeloid lineages, suggesting that miRNAs function in a cell-autonomous manner. The combinatorial expression of miRNAs miR-186-5p, miR-128-3p, and miR-330-5p in MPPs significantly attenuates their myeloid differentiation potential due to repression of myeloid-associated transcripts. Depletion of these miRNAs caused a pronounced de-repression of myeloid lineage targets in differentiating early-B and early-T cells, resulting in a mixed-lineage gene expression pattern. De novo motif analysis combined with an assay of promoter activities indicates that B as well as T lineage determinants drive the expression of these miRNAs in lymphoid lineages. Collectively, we present a paradigm that miRNAs are conserved between developing B and T lymphocytes, yet they target distinct sets of promiscuously expressed lineage-inappropriate genes to suppress the alternate cell-fate options. Thus, our studies provide a comprehensive compendium of miRNAs with functional implications for B and T lymphocyte development.

Prediction methods for microRNA targets in bilaterian animals: Toward a better understanding by biologists

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression at the posttranscriptional level. Because of their wide network of interactions, miRNAs have become the focus of many studies over the past decade, particularly in animal species. To streamline the number of potential wet lab experiments, the use of miRNA target prediction tools is currently the first step undertaken. However, the predictions made may vary considerably depending on the tool used, which is mostly due to the complex and still not fully understood mechanism of action of miRNAs. The discrepancies complicate the choice of the tool for miRNA target prediction. To provide a comprehensive view of this issue, we highlight in this review the main characteristics of miRNA-target interactions in bilaterian animals, describe the prediction models currently used, and provide some insights for the evaluation of predictor performance.

Molecular Mechanisms of Nutrient-Mediated Regulation of MicroRNAs in Pancreatic β-cells

Pancreatic β-cells within the islets of Langerhans respond to rising blood glucose levels by secreting insulin that stimulates glucose uptake by peripheral tissues to maintain whole body energy homeostasis. To different extents, failure of β-cell function and/or β-cell loss contribute to the development of Type 1 and Type 2 diabetes. Chronically elevated glycaemia and high circulating free fatty acids, as often seen in obese diabetics, accelerate β-cell failure and the development of the disease. MiRNAs are essential for endocrine development and for mature pancreatic β-cell function and are dysregulated in diabetes. In this review, we summarize the different molecular mechanisms that control miRNA expression and function, including transcription, stability, posttranscriptional modifications, and interaction with RNA binding proteins and other non-coding RNAs. We also discuss which of these mechanisms are responsible for the nutrient-mediated regulation of the activity of β-cell miRNAs and identify some of the more important knowledge gaps in the field.

Sequence determinants as key regulators in gene expression of T cells

T cell homeostasis, T cell differentiation, and T cell effector function rely on the constant fine-tuning of gene expression. To alter the T cell state, substantial remodeling of the proteome is required. This remodeling depends on the intricate interplay of regulatory mechanisms, including post-transcriptional gene regulation. In this review, we discuss how the sequence of a transcript influences these post-transcriptional events. In particular, we review how sequence determinants such as sequence conservation, GC content, and chemical modifications define the levels of the mRNA and the protein in a T cell. We describe the effect of different forms of alternative splicing on mRNA expression and protein production, and their effect on subcellular localization. In addition, we discuss the role of sequences and structures as binding hubs for miRNAs and RNA-binding proteins in T cells. The review thus highlights how the intimate interplay of post-transcriptional mechanisms dictate cellular fate decisions in T cells.

Modified Cross-Linking, Ligation, and Sequencing of Hybrids (qCLASH) to Identify MicroRNA Targets

This protocol was designed to identify microRNA (miRNA) targetomes from smaller‐input samples by performing a simplified workflow of the Cross‐Linking and Sequencing of Hybrids (CLASH) technique developed in the Tollervey group. In this ribonomics‐based technique, Cross‐Linking and Immunoprecipitation (CLIP) of Argonaute (Ago) is combined with an RNA ligase reaction that yields covalently bound “hybrids” between miRNAs and their target RNAs. While this iteration of CLIP identifies “high‐confidence” or “unambiguous” miRNA targets, the added ligation step is highly inefficient and therefore requires large numbers of cultured cells. To make this powerful approach applicable to smaller cell numbers, we created qCLASH, incorporating a workflow that performs all enzymatic reactions on bead‐bound complexes and omits gel purification of immunoprecipitated Ago complexes associated with major loss of RNA. At a sequencing depth of 100 million reads per library, which is highly feasible with rapidly decreasing sequencing costs, qCLASH, when used with three biological replicates, results in thousands of high‐confidence miRNA targets. qCLASH was first developed to identify viral miRNA targetomes of endothelial cells infected with Kaposi's sarcoma−associated herpesvirus. Since then, qCLASH has been applied to Epstein‐Barr virus− and MHV68‐infected cells, and more recently to metastatic melanoma and breast cancer cells. Currently, protocols are under development to apply qCLASH to human solid tumor specimens. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC.

Basic Protocol: Quick Cross‐Linking and Sequencing of Hybrids (qCLASH)

Support Protocol: Optimization of Ago immunoprecipitation

Emerging Role of isomiRs in Cancer: State of the Art and Recent Advances

The advent of Next Generation Sequencing technologies brought with it the discovery of several microRNA (miRNA) variants of heterogeneous lengths and/or sequences. Initially ascribed to sequencing errors/artifacts, these isoforms, named isomiRs, are now considered non-canonical variants that originate from physiological processes affecting the canonical miRNA biogenesis. To date, accurate IsomiRs abundance, biological activity, and functions are not completely understood; however, the study of isomiR biology is an area of great interest due to their high frequency in the human miRNome, their putative functions in cooperating with the canonical miRNAs, and potential for exhibiting novel functional roles. The discovery of isomiRs highlighted the complexity of the small RNA transcriptional landscape in several diseases, including cancer. In this field, the study of isomiRs could provide further insights into the miRNA biology and its implication in oncogenesis, possibly providing putative new cancer diagnostic, prognostic, and predictive biomarkers as well. In this review, a comprehensive overview of the state of research on isomiRs in different cancer types, including the most common tumors such as breast cancer, colorectal cancer, melanoma, and prostate cancer, as well as in the less frequent tumors, as for example brain tumors and hematological malignancies, will be summarized and discussed.

Making use of transcription factor enrichment to identify functional microRNA-regulons

microRNAs (miRNAs) are important modulators of messenger RNA stability and translation, controlling wide gene networks. Albeit generally modest on individual targets, the regulatory effect of miRNAs translates into meaningful pathway modulation through concurrent targeting of regulons with functional convergence. Identification of miRNA-regulons is therefore essential to understand the function of miRNAs and to help realise their therapeutic potential, but it remains challenging due to the large number of false positive target sites predicted per miRNA. In the current work, we investigated whether genes regulated by a given miRNA were under the transcriptional control of a predominant transcription factor (TF). Strikingly we found that for ~50% of the miRNAs analysed, their targets were significantly enriched in at least one common TF. We leveraged such miRNA-TF co-regulatory networks to identify pathways under miRNA control, and demonstrated that filtering predicted miRNA-target interactions (MTIs) relying on such pathways significantly enriched the proportion of predicted true MTIs. To our knowledge, this is the first description of an in- silico pipeline facilitating the identification of miRNA-regulons, to help understand miRNA function.

Exosomes as mediators of intercellular crosstalk in metabolism

Exosomes are nanoparticles secreted by all cell types and are a large component of the broader class of nanoparticles termed extracellular vesicles (EVs). Once secreted, exosomes gain access to the interstitial space and ultimately the circulation, where they exert local paracrine or distal systemic effects. Because of this, exosomes are important components of an intercellular and intraorgan communication system capable of carrying biologic signals from one cell type or tissue to another. The exosomal cargo consists of proteins, lipids, miRNAs, and other RNA species, and many of the biologic effects of exosomes have been attributed to miRNAs. Exosomal miRNAs have also been used as disease biomarkers. The field of exosome biology and metabolism is rapidly expanding, with new discoveries and reports appearing on a regular basis, and it is possible that potential therapeutic approaches for the use of exosomes or miRNAs in metabolic diseases will be initiated in the near future.

The long noncoding RNA Synage regulates synapse stability and neuronal function in the cerebellum

The brain is known to express many long noncoding RNAs (lncRNAs); however, whether and how these lncRNAs function in modulating synaptic stability remains unclear. Here, we report a cerebellum highly expressed lncRNA, Synage, regulating synaptic stability via at least two mechanisms. One is through the function of Synage as a sponge for the microRNA miR-325-3p, to regulate expression of the known cerebellar synapse organizer Cbln1. The other function is to serve as a scaffold for organizing the assembly of the LRP1-HSP90AA1-PSD-95 complex in PF-PC synapses. Although somewhat divergent in its mature mRNA sequence, the locus encoding Synage is positioned adjacent to the Cbln1 loci in mouse, rhesus macaque, and human, and Synage is highly expressed in the cerebella of all three species. Synage deletion causes a full-spectrum cerebellar ablation phenotype that proceeds from cerebellar atrophy, through neuron loss, on to synapse density reduction, synaptic vesicle loss, and finally to a reduction in synaptic activity during cerebellar development; these deficits are accompanied by motor dysfunction in adult mice, which can be rescued by AAV-mediated Synage overexpression from birth. Thus, our study demonstrates roles for the lncRNA Synage in regulating synaptic stability and function during cerebellar development.

MicroRNA Interactome Multiomics Characterization for Cancer Research and Personalized Medicine: An Expert Review

MicroRNAs (miRNAs) that are mutually modulated by their interacting partners (interactome) are being increasingly noted for their significant role in pathogenesis and treatment of various human cancers. Recently, miRNA interactome dissected with multiomics approaches has been the subject of focus since individual tools or methods failed to provide the necessary comprehensive clues on the complete interactome. Even though single-omics technologies such as proteomics can uncover part of the interactome, the biological and clinical understanding still remain incomplete. In this study, we present an expert review of studies involving multiomics approaches to identification of miRNA interactome and its application in mechanistic characterization, classification, and therapeutic target identification in a variety of cancers, and with a focus on proteomics. We also discuss individual or multiple miRNA-based interactome identification in various pathological conditions of relevance to clinical medicine. Various new single-omics methods that can be integrated into multiomics cancer research and the computational approaches to analyze and predict miRNA interactome are also highlighted in this review. In all, we contextulize the power of multiomics approaches and the importance of the miRNA interactome to achieve the vision and practice of predictive, preventive, and personalized medicine in cancer research and clinical oncology.

Effective tools for RNA-derived therapeutics: siRNA interference or miRNA mimicry

The approval of the first small interfering RNA (siRNA) drug Patisiran by FDA in 2018 marks a new era of RNA interference (RNAi) therapeutics. MicroRNAs (miRNA), an important post-transcriptional gene regulator, are also the subject of both basic research and clinical trials. Both siRNA and miRNA mimics are ~21 nucleotides RNA duplexes inducing mRNA silencing. Given the well performance of siRNA, researchers ask whether miRNA mimics are unnecessary or developed siRNA technology can pave the way for the emergence of miRNA mimic drugs. Through comprehensive comparison of siRNA and miRNA, we focus on (1) the common features and lessons learnt from the success of siRNAs; (2) the unique characteristics of miRNA that potentially offer additional therapeutic advantages and opportunities; (3) key areas of ongoing research that will contribute to clinical application of miRNA mimics. In conclusion, miRNA mimics have unique properties and advantages which cannot be fully matched by siRNA in clinical applications. MiRNAs are endogenous molecules and the gene silencing effects of miRNA mimics can be regulated or buffered to ameliorate or eliminate off-target effects. An in-depth understanding of the differences between siRNA and miRNA mimics will facilitate the development of miRNA mimic drugs.

Regulation of AMPAR expression by microRNAs

AMPA receptors (AMPARs) are the major excitatory neurotransmitter receptor in the brain, and their expression at synapses is a critical determinant of synaptic transmission and therefore brain function. Synaptic plasticity involves increases or decreases in synaptic strength, caused by changes in the number or subunit-specific subtype of AMPARs expressed at synapses, and resulting in modifications of functional connectivity of neuronal circuits, a process which is thought to underpin learning and the formation or loss of memories. Furthermore, numerous neurological disorders involve dysregulation of excitatory synaptic transmission or aberrant recruitment of plasticity processes. MicroRNAs (miRNAs) repress the translation of target genes by partial complementary base pairing with mRNAs, and are the core component of a mechanism widely used in a range of cell processes for regulating protein translation. MiRNA-dependent translational repression can occur locally in neuronal dendrites, close to synapses, and can also result in relatively rapid changes in protein expression. MiRNAs are therefore well-placed to regulate synaptic plasticity via the local control of AMPAR subunit synthesis, and can also result in synaptic dysfunction in the event of dysregulation in disease. Here, I will review the miRNAs that have been identified as playing a role in physiological or pathological changes in AMPAR subunit expression at synapses, focussing on miRNAs that target mRNAs encoding AMPAR subunits, and on miRNAs that target AMPAR accessory proteins involved in AMPAR trafficking and hence the regulation of AMPAR synaptic localisation.

Critical roles of microRNA-196 in normal physiology and non-malignant diseases: Diagnostic and therapeutic implications

MicroRNAs (miRNAs) have emerged as a critical component of regulatory networks that modulate and fine-tune gene expression in a post-transcriptional manner. The microRNA-196 family is encoded by three loci in the human genome, namely hsa-mir-196a-1, hsa-mir-196a-2, and hsa-mir-196b. Increasing evidence supports the roles of different components of this miRNA family in regulating key cellular processes during differentiation and development, ranging from inflammation and differentiation of stem cells to limb development and remodeling and structure of adipose tissue. This review first discusses about the genomic context and regulation of this miRNA family and then take a bird's eye view on the updated list of its target genes and their biological processes to obtain insights about various functions played by members of the microRNA-196 family. We then describe evidence supporting the involvement of the human microRNA-196 family in regulating critical cellular processes both in physiological and non-malignant inflammatory conditions, highlighting recent seminal findings that carry implications for developing novel therapeutic or diagnostic strategies.

MicroRNA Sequences Modulated by Beta Cell Lipid Metabolism: Implications for Type 2 Diabetes Mellitus

Alterations in lipid metabolism within beta cells and islets contributes to dysfunction and apoptosis of beta cells, leading to loss of insulin secretion and the onset of type 2 diabetes. Over the last decade, there has been an explosion of interest in understanding the landscape of gene expression which influences beta cell function, including the importance of small non-coding microRNA sequences in this context. This review sought to identify the microRNA sequences regulated by metabolic challenges in beta cells and islets, their targets, highlight their function and assess their possible relevance as biomarkers of disease progression in diabetic individuals. Predictive analysis was used to explore networks of genes targeted by these microRNA sequences, which may offer new therapeutic strategies to protect beta cell function and delay the onset of type 2 diabetes.

The MicroRNA Family Gets Wider: The IsomiRs Classification and Role

MicroRNAs (miRNAs or miRs) are the most characterized class of non-coding RNAs and are engaged in many cellular processes, including cell differentiation, development, and homeostasis. MicroRNA dysregulation was observed in several diseases, cancer included. Epitranscriptomics is a branch of epigenomics that embraces all RNA modifications occurring after DNA transcription and RNA synthesis and involving coding and non-coding RNAs. The development of new high-throughput technologies, especially deep RNA sequencing, has facilitated the discovery of miRNA isoforms (named isomiRs) resulting from RNA modifications mediated by enzymes, such as deaminases and exonucleases, and differing from the canonical ones in length, sequence, or both. In this review, we summarize the distinct classes of isomiRs, their regulation and biogenesis, and the active role of these newly discovered molecules in cancer and other diseases.

Biological differentiation of traditional Chinese medicine from excessive to deficient syndromes in AIDS: Comparative microRNA microarray profiling and syndrome-specific biomarker identification

Traditional Chinese medicine (TCM) has been widely applied as a supplementary therapy of human immunodeficiency virus infection and acquired immunodeficiency syndrome (HIV/AIDS) in China. TCM has a positive effect on improving the quality of life, prolonging life, and ameliorating the symptoms of HIV/AIDS patients. Yang deficiency of spleen and kidney (YDSK) syndrome is a typical deficient TCM syndrome in AIDS patients, and accumulation of heat-toxicity (AHT) syndrome is a common excessive syndrome in the earlier stage of AIDS. Thus, accurate diagnosis of these two syndromes can improve the targeted treatment effect, and predict the prognosis of the disease. However, the scientific basis of TCM syndromes remains lacking, greatly hindering the accuracy of diagnosis and effectiveness of treatment. In this research, microRNA (miRNA) microarray and quantitative real-time polymerase chain reaction combined with bioinformatics were used for comparative analysis between YDSK and AHT patients. Significantly differential expressed miRNAs (SDE-miRNAs) of each TCM syndrome were identified, including hsa-miR-766-3p and hsa-miR-1260a and so on, as well hsa-miR-6124, hsa-let-7g-5p and so on, for YDSK and AHT, respectively. Biological differences were found between their SDE-miRNAs based on bioinformatics analyses, for example, ErbB signaling pathway mainly linked to AHT, while focal adhesion dominated in YDSK. Syndrome-specific SDE-miRNAs were further identified as potential biomarkers, including hsa-miR-30e-5p, hsa-miR-144-5p for YDSK and hsa-let-7g-5p, hsa-miR-126-3p for AHT, respectively. All of them have laid biological and clinical bases for TCM diagnosis and treatment of AIDS syndrome at the miRNA level, offering potential diagnostic indicators of immune reconstitution.

A noncanonical microRNA derived from the snaR-A noncoding RNA targets a metastasis inhibitor

MicroRNAs (miRNAs) are small noncoding RNAs that function as critical posttranscriptional regulators in various biological processes. While most miRNAs are generated from processing of long primary transcripts via sequential Drosha and Dicer cleavage, some miRNAs that bypass Drosha cleavage can be transcribed as part of another small noncoding RNA. Here, we develop the target-oriented miRNA discovery (TOMiD) bioinformatic analysis method to identify Drosha-independent miRNAs from Argonaute crosslinking and sequencing of hybrids (Ago-CLASH) data sets. Using this technique, we discovered a novel miRNA derived from a primate specific noncoding RNA, the small NF90 associated RNA A (snaR-A). The miRNA derived from snaR-A (miR-snaR) arises independently of Drosha processing but requires Exportin-5 and Dicer for biogenesis. We identify that miR-snaR is concurrently up-regulated with the full snaR-A transcript in cancer cells. Functionally, miR-snaR associates with Ago proteins and targets NME1, a key metastasis inhibitor, contributing to snaR-A's role in promoting cancer cell migration. Our findings suggest a functional link between a novel miRNA and its precursor noncoding RNA.

Comprehensive machine-learning-based analysis of microRNA-target interactions reveals variable transferability of interaction rules across species

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression post-transcriptionally via base-pairing with complementary sequences on messenger RNAs (mRNAs). Due to the technical challenges involved in the application of high-throughput experimental methods, datasets of direct bona fide miRNA targets exist only for a few model organisms. Machine learning (ML)-based target prediction models were successfully trained and tested on some of these datasets. There is a need to further apply the trained models to organisms in which experimental training data are unavailable. However, it is largely unknown how the features of miRNA-target interactions evolve and whether some features have remained fixed during evolution, raising questions regarding the general, cross-species applicability of currently available ML methods.

RESULTS

We examined the evolution of miRNA-target interaction rules and used data science and ML approaches to investigate whether these rules are transferable between species. We analyzed eight datasets of direct miRNA-target interactions in four species (human, mouse, worm, cattle). Using ML classifiers, we achieved high accuracy for intra-dataset classification and found that the most influential features of all datasets overlap significantly. To explore the relationships between datasets, we measured the divergence of their miRNA seed sequences and evaluated the performance of cross-dataset classification. We found that both measures coincide with the evolutionary distance between the compared species.

CONCLUSIONS

The transferability of miRNA-targeting rules between species depends on several factors, the most associated factors being the composition of seed families and evolutionary distance. Furthermore, our feature-importance results suggest that some miRNA-target features have evolved while others remained fixed during the evolution of the species. Our findings lay the foundation for the future development of target prediction tools that could be applied to "non-model" organisms for which minimal experimental data are available.

AVAILABILITY AND IMPLEMENTATION

The code is freely available at https://github.com/gbenor/TPVOD .

Restoring miR-132 expression rescues adult hippocampal neurogenesis and memory deficits in Alzheimer's disease

Neural stem cells residing in the hippocampal neurogenic niche sustain lifelong neurogenesis in the adult brain. Adult hippocampal neurogenesis (AHN) is functionally linked to mnemonic and cognitive plasticity in humans and rodents. In Alzheimer's disease (AD), the process of generating new neurons at the hippocampal neurogenic niche is impeded, yet the mechanisms involved are unknown. Here we identify miR-132, one of the most consistently downregulated microRNAs in AD, as a potent regulator of AHN, exerting cell-autonomous proneurogenic effects in adult neural stem cells and their progeny. Using distinct AD mouse models, cultured human primary and established neural stem cells, and human patient material, we demonstrate that AHN is directly affected by AD pathology. miR-132 replacement in adult mouse AD hippocampus restores AHN and relevant memory deficits. Our findings corroborate the significance of AHN in mouse models of AD and reveal the possible therapeutic potential of targeting miR-132 in neurodegeneration.

TimiRGeN: R/Bioconductor package for time series microRNA-mRNA integration and analysis

Motivation

The analysis of longitudinal datasets and construction of gene regulatory networks (GRNs) provide a valuable means to disentangle the complexity of microRNA (miRNA)–mRNA interactions. However, there are no computational tools that can integrate, conduct functional analysis and generate detailed networks from longitudinal miRNA–mRNA datasets.

Results

We present TimiRGeN, an R package that uses time point-based differential expression results to identify miRNA–mRNA interactions influencing signaling pathways of interest. miRNA–mRNA interactions can be visualized in R or exported to PathVisio or Cytoscape. The output can be used for hypothesis generation and directing in vitro or further in silico work such as GRN construction.

Availability and implementation

TimiRGeN is available for download on Bioconductor (https://bioconductor.org/packages/TimiRGeN) and requires R v4.0.2 or newer and BiocManager v3.12 or newer.

Supplementary information

Supplementary data are available at Bioinformatics online.

A cohort of Caenorhabditis species lacking the highly conserved let-7 microRNA

The let-7 gene encodes a highly conserved microRNA with critical functions integral to cell fate specification and developmental progression in diverse animals. In Caenorhabditis elegans, let-7 is a component of the heterochronic (developmental timing) gene regulatory network, and loss-of-function mutations of let-7 result in lethality during the larval to adult transition due to misregulation of the conserved let-7 target, lin-41. To date, no bilaterian animal lacking let-7 has been characterized. In this study, we identify a cohort of nematode species within the genus Caenorhabditis, closely related to C. elegans, that lack the let-7 microRNA, owing to absence of the let-7 gene. Using Caenorhabditis sulstoni as a representative let-7-lacking species to characterize normal larval development in the absence of let-7, we demonstrate that, except for the lack of let-7, the heterochronic gene network is otherwise functionally conserved. We also report that species lacking let-7 contain a group of divergent let-7 paralogs-also known as the let-7-family of microRNAs-that have apparently assumed the role of targeting the LIN-41 mRNA.

Argonaute-CLIP delineates versatile, functional RNAi networks in Aedes aegypti, a major vector of human viruses

Argonaute (AGO) proteins bind small RNAs to silence complementary RNA transcripts, and they are central to RNA interference (RNAi). RNAi is critical for regulation of gene expression and antiviral defense in Aedes aegypti mosquitoes, which transmit Zika, chikungunya, dengue, and yellow fever viruses. In mosquitoes, AGO1 mediates miRNA interactions, while AGO2 mediates siRNA interactions. We applied AGO-crosslinking immunoprecipitation (AGO-CLIP) for both AGO1 and AGO2, and we developed a universal software package for CLIP analysis (CLIPflexR), identifying 230 small RNAs and 5,447 small RNA targets that comprise a comprehensive RNAi network map in mosquitoes. RNAi network maps predicted expression levels of small RNA targets in specific tissues. Additionally, this resource identified unexpected, context-dependent AGO2 target preferences, including endogenous viral elements and 3'UTRs. Finally, contrary to current thinking, mosquito AGO2 repressed imperfect targets. These findings expand our understanding of small RNA networks and have broad implications for the study of antiviral RNAi.

Network based multifactorial modelling of miRNA-target interactions

Competing endogenous RNA (ceRNA) regulations and crosstalk between various types of non-coding RNA in humans is an important and under-explored subject. Several studies have pointed out that an alteration in miRNA:target interaction can result in unexpected changes due to indirect and complex interactions. In this article, we defined a new network-based model that incorporates miRNA:ceRNA interactions with expression values. Our approach calculates network-wide effects of perturbations in the expression level of one or more nodes in the presence or absence of miRNA interaction factors such as seed type, binding energy. We carried out the analysis of large-scale miRNA:target networks from breast cancer patients. Highly perturbing genes identified by our approach coincide with breast cancer-associated genes and miRNAs. Our network-based approach takes the sponge effect into account and helps to unveil the crosstalk between nodes in miRNA:target network. The model has potential to reveal unforeseen regulations that are only evident in the network context. Our tool is scalable and can be plugged in with emerging miRNA effectors such as circRNAs, lncRNAs, and available as R package ceRNAnetsim: https://www.bioconductor.org/packages/release/bioc/html/ceRNAnetsim.html.

Cross-Linking Ligation and Sequencing of Hybrids (qCLASH) Reveals an Unpredicted miRNA Targetome in Melanoma Cells

MicroRNAs are key post-transcriptional gene regulators often displaying aberrant expression patterns in cancer. As microRNAs are promising disease-associated biomarkers and modulators of responsiveness to anti-cancer therapies, a solid understanding of their targetome is crucial. Despite enormous research efforts, the success rates of available tools to reliably predict microRNAs (miRNA)-target interactions remains limited. To investigate the disease-associated miRNA targetome, we have applied modified cross-linking ligation and sequencing of hybrids (qCLASH) to BRAF-mutant melanoma cells. The resulting RNA-RNA hybrid molecules provide a comprehensive and unbiased snapshot of direct miRNA-target interactions. The regulatory effects on selected miRNA target genes in predicted vs. non-predicted binding regions was validated by miRNA mimic experiments. Most miRNA-target interactions deviate from the central dogma of miRNA targeting up to 60% interactions occur via non-canonical seed pairing with a strong contribution of the 3' miRNA sequence, and over 50% display a clear bias towards the coding sequence of mRNAs. miRNAs targeting the coding sequence can directly reduce gene expression (miR-34a/CD68), while the majority of non-canonical miRNA interactions appear to have roles beyond target gene suppression (miR-100/AXL). Additionally, non-mRNA targets of miRNAs (lncRNAs) whose interactions mainly occur via non-canonical binding were identified in melanoma. This first application of CLASH sequencing to cancer cells identified over 8 K distinct miRNA-target interactions in melanoma cells. Our data highlight the importance non-canonical interactions, revealing further layers of complexity of post-transcriptional gene regulation in melanoma, thus expanding the pool of miRNA-target interactions, which have so far been omitted in the cancer field.

tRFTars: predicting the targets of tRNA-derived fragments

Background

tRNA-derived fragments (tRFs) are 14–40-nucleotide-long, small non-coding RNAs derived from specific tRNA cleavage events with key regulatory functions in many biological processes. Many studies have shown that tRFs are associated with Argonaute (AGO) complexes and inhibit gene expression in the same manner as miRNAs. However, there are currently no tools for accurately predicting tRF target genes.

Methods

We used tRF-mRNA pairs identified by crosslinking, ligation, and sequencing of hybrids (CLASH) and covalent ligation of endogenous AGO-bound RNAs (CLEAR)-CLIP to assess features that may participate in tRF targeting, including the sequence context of each site and tRF-mRNA interactions. We applied genetic algorithm (GA) to select key features and support vector machine (SVM) to construct tRF prediction models.

Results

We first identified features that globally influenced tRF targeting. Among these features, the most significant were the minimum free folding energy (MFE), position 8 match, number of bases paired in the tRF-mRNA duplex, and length of the tRF, which were consistent with previous findings. Our constructed model yielded an area under the receiver operating characteristic (ROC) curve (AUC) = 0.980 (0.977–0.983) in the training process and an AUC = 0.847 (0.83–0.861) in the test process. The model was applied to all the sites with perfect Watson–Crick complementarity to the seed in the 3′ untranslated region (3′-UTR) of the human genome. Seven of nine target/nontarget genes of tRFs confirmed by reporter assay were predicted. We also validated the predictions via quantitative real-time PCR (qRT-PCR). Thirteen potential target genes from the top of the predictions were significantly down-regulated at the mRNA levels by overexpression of the tRFs (tRF-3001a, tRF-3003a or tRF-3009a).

Conclusions

Predictions can be obtained online, tRFTars, freely available at http://trftars.cmuzhenninglab.org:3838/tar/, which is the first tool to predict targets of tRFs in humans with a user-friendly interface.