Bias-minimized quantification of microRNA reveals widespread alternative processing and 3' end modification

An ultrasensitive one-pot Cas13a-based microfluidic assay for rapid multiplexed detection of microRNAs

Aberrant microRNA expression is associated with tumor progression in various organs. Detecting microRNAs as clinical cancer biomarkers can facilitate early cancer diagnosis and monitoring. However, the rapid and accurate quantification of microRNAs from biological samples remains a significant challenge. Here we developed a one-pot isothermal assay utilizing a molecular circuit with CRISPR/Cas13a (CRISPR-circuit) to rapidly convert, amplify and report different microRNAs within 15 min at the attomolar (aM) level. Then the full process was performed on an active centrifugal microfluidic chip and its corresponding portable equipment for parallel detection of multiple microRNAs, including miR-21, miR-141, miR-196a, and miR-1246. We also demonstrated its application for identifying cell lines and clinical samples of cancer patients with varying microRNA levels, which showed a strong correlation with the RT-qPCR. The assay can be easily adapted for the detection of any microRNA by simply modifying the converter primer, thereby holding significant potential for accurate disease detection and clinical diagnosis.

Revealing microRNA regulation in single cells

MicroRNAs (miRNAs) are key regulators of gene expression and control cellular functions in physiological and pathophysiological states. miRNAs play important roles in disease, stress, and development, and are now being investigated for therapeutic approaches. Alternative processing of miRNAs during biogenesis results in the generation of miRNA isoforms (isomiRs) which further diversify miRNA gene regulation. Single-cell RNA-sequencing (scsRNA-seq) technologies, together with computational strategies, enable exploration of miRNAs, isomiRs, and interacting RNAs at the cellular level. By integration with other miRNA-associated single-cell modalities, miRNA roles can be resolved at different stages of processing and regulation. In this review we discuss (i) single-cell experimental assays that measure miRNA and isomiR abundances, and (ii) computational methods for their analysis to investigate the mechanisms of miRNA biogenesis and post-transcriptional regulation.

Temporal Expression Analysis to Unravel Gene Regulatory Dynamics by microRNAs

MicroRNAs (miRNAs) are a class of small non-coding RNAs (sncRNAs) of length 21-25 nucleotides. These sncRNAs hybridize to repress their target genes and inhibit protein translation, thereby controlling regulatory functions in the cell. Integration of time-series matched small and RNA-seq data enables investigation of dynamic gene regulation through miRNAs during development or in response to a stimulus, such as stress. Here we summarize analysis strategies, such as probabilistic and regression-based models, that take advantage of the temporal dimension to investigate the complexity of miRNA regulation.

The biogenesis and regulation of animal microRNAs

MicroRNAs (miRNAs) are small, yet profoundly influential, non-coding RNAs that base-pair with mRNAs to induce RNA silencing. Although the basic principles of miRNA biogenesis and function have been established, recent breakthroughs have yielded important new insights into the molecular mechanisms of miRNA biogenesis. In this Review, we discuss the metazoan miRNA biogenesis pathway step-by-step, focusing on the key biogenesis machinery, including the Drosha-DGCR8 complex (Microprocessor), exportin-5, Dicer and Argonaute. We also highlight newly identified cis-acting elements and their impact on miRNA maturation, informed by advanced high-throughput and structural studies, and discuss recently discovered mechanisms of clustered miRNA processing, target recognition and target-directed miRNA decay (TDMD). Lastly, we explore multiple regulatory layers of miRNA biogenesis, mediated by RNA-protein interactions, miRNA tailing (uridylation or adenylation) and RNA modifications.

An endogenous cluster of target-directed microRNA degradation sites induces decay of distinct microRNA families

While much is known about miRNA biogenesis and canonical miRNA targeting, relatively less is understood about miRNA decay. The major miRNA decay pathway in metazoans is mediated through target-directed miRNA degradation (TDMD), in which certain RNAs can "trigger" miRNA decay. All known triggers for TDMD base pair with the miRNA seed, and extensively base pair on the miRNA 3' end, a pattern that supposedly induces a TDMD-competent conformational change of Argonaute (Ago), allowing for miRNA turnover. Here, we utilized Ago1-CLASH to find that the Drosophila transcript Kah contains at least two triggers, a "trigger cluster", against miR-9b and the miR-279 family. One of these triggers contains minimal/non-canonical 3' end base pairing but is still sufficient to induce TDMD of the entire miR-279 family. We found that these clustered triggers likely lack cooperativity, the minimal 3' pairing is required for miR-279 family turnover, and probed the in-cell RNA structure of the Kah trigger cluster. Overall, this study expands the list of endogenous triggers and the unexpectedly complex regulatory network governing miRNA degradation.

Small and long non-coding RNAs: Past, present, and future

Since the introduction of the central dogma of molecular biology in 1958, various RNA species have been discovered. Messenger RNAs transmit genetic instructions from DNA to make proteins, a process facilitated by housekeeping non-coding RNAs (ncRNAs) such as small nuclear RNAs (snRNAs), ribosomal RNAs (rRNAs), and transfer RNAs (tRNAs). Over the past four decades, a wide array of regulatory ncRNAs have emerged as crucial players in gene regulation. In celebration of Cell's 50th anniversary, this Review explores our current understanding of the most extensively studied regulatory ncRNAs-small RNAs and long non-coding RNAs (lncRNAs)-which have profoundly shaped the field of RNA biology and beyond. While small RNA pathways have been well documented with clearly defined mechanisms, lncRNAs exhibit a greater diversity of mechanisms, many of which remain unknown. This Review covers pivotal events in their discovery, biogenesis pathways, evolutionary traits, action mechanisms, functions, and crosstalks among ncRNAs. We also highlight their roles in pathophysiological contexts and propose future research directions to decipher the unknowns of lncRNAs by leveraging lessons from small RNAs.

Biochemical and structural insights into a 5' to 3' RNA ligase reveal a potential role in tRNA ligation

ATP-grasp superfamily enzymes contain a hand-like ATP-binding fold and catalyze a variety of reactions using a similar catalytic mechanism. More than 30 protein families are categorized in this superfamily, and they are involved in a plethora of cellular processes and human diseases. Here, we identify C12orf29 (RLIG1) as an atypical ATP-grasp enzyme that ligates RNA. Human RLIG1 and its homologs autoadenylate on an active site Lys residue as part of a reaction intermediate that specifically ligates RNA halves containing a 5'-phosphate and a 3'-hydroxyl. RLIG1 binds tRNA in cells and can ligate tRNA within the anticodon loop in vitro. Transcriptomic analyses of Rlig1 knockout mice revealed significant alterations in global tRNA levels in the brains of female mice, but not in those of male mice. Furthermore, crystal structures of a RLIG1 homolog from Yasminevirus bound to nucleotides revealed a minimal and atypical RNA ligase fold with a conserved active site architecture that participates in catalysis. Collectively, our results identify RLIG1 as an RNA ligase and suggest its involvement in tRNA biology.

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

MicroRNA biogenesis is broadly disrupted by inhibition of the splicing factor SF3B1

In animals, microRNA (miRNA) biogenesis begins with cotranscriptional cleavage of the primary (pri-)miRNA by the Microprocessor complex. Cotranscriptional splicing has been shown to influence Microprocessor cleavage when miRNAs are hosted in introns of protein-coding pri-miRNAs, but the impact of splicing on production of miRNAs hosted in long non-coding (lnc)RNAs is largely unknown. Here, we investigated the role of splicing in the biogenesis of miR-122, an lncRNA-hosted, highly expressed, medically important, liver-specific miRNA. We found that splicing inhibition by the SF3B1 inhibitor pladienolide B (PlaB) led to strong and rapid reduction in transcription of endogenous, but not plasmid-encoded, pri-miR-122, resulting in reduced production of mature miR-122. To allow detection of rapid changes in miRNA biogenesis despite the high stability of mature miRNAs, we used SLAMseq to globally quantify the effects of short-term splicing inhibition on miRNA synthesis. We observed an overall decrease in biogenesis of mature miRNAs following PlaB treatment. Surprisingly, miRNAs hosted in exons and introns were similarly affected. Together, this study provides new insights into the emerging role of splicing in transcription, demonstrating novel biological importance in promotion of miR-122 biogenesis from an lncRNA, and shows that SF3B1 is important for global miRNA biogenesis.

Catalytic residues of microRNA Argonautes play a modest role in microRNA star strand destabilization in C. elegans

Many microRNA (miRNA)-guided Argonaute proteins can cleave RNA ('slicing'), even though miRNA-mediated target repression is generally cleavage-independent. Here we use Caenorhabditis elegans to examine the role of catalytic residues of miRNA Argonautes in organismal development. In contrast to previous work, mutations in presumed catalytic residues did not interfere with development when introduced by CRISPR. We find that unwinding and decay of miRNA star strands is weakly defective in the catalytic residue mutants, with the largest effect observed in embryos. Argonaute-Like Gene 2 (ALG-2) is more dependent on catalytic residues for unwinding than ALG-1. The miRNAs that displayed the greatest (albeit minor) dependence on catalytic residues for unwinding tend to form stable duplexes with their star strand, and in some cases, lowering duplex stability alleviates dependence on catalytic residues. While a few miRNA guide strands are reduced in the mutant background, the basis of this is unclear since changes were not dependent on EBAX-1, an effector of Target-Directed miRNA Degradation (TDMD). Overall, this work defines a role for the catalytic residues of miRNA Argonautes in star strand decay; future work should examine whether this role contributes to the selection pressure to conserve catalytic activity of miRNA Argonautes across the metazoan phylogeny.

Biochemical and structural insights into a 5' to 3' RNA ligase reveal a potential role in tRNA ligation

ATP-grasp superfamily enzymes contain a hand-like ATP-binding fold and catalyze a variety of reactions using a similar catalytic mechanism. More than 30 protein families are categorized in this superfamily, and they are involved in a plethora of cellular processes and human diseases. Here we identify C12orf29 as an atypical ATP-grasp enzyme that ligates RNA. Human C12orf29 and its homologs auto-adenylate on an active site Lys residue as part of a reaction intermediate that specifically ligates RNA halves containing a 5'-phosphate and a 3'-hydroxyl. C12orf29 binds tRNA in cells and can ligate tRNA within the anticodon loop in vitro. Genetic depletion of c12orf29 in female mice alters global tRNA levels in brain. Furthermore, crystal structures of a C12orf29 homolog from Yasminevirus bound to nucleotides reveal a minimal and atypical RNA ligase fold with a unique active site architecture that participates in catalysis. Collectively, our results identify C12orf29 as an RNA ligase and suggest its involvement in tRNA biology.

Structural atlas of human primary microRNAs generated by SHAPE-MaP

MicroRNA (miRNA) maturation is critically dependent on structural features of primary transcripts (pri-miRNAs). However, the scarcity of determined pri-miRNA structures has limited our understanding of miRNA maturation. Here, we employed selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), a high-throughput RNA structure probing method, to unravel the secondary structures of 476 high-confidence human pri-miRNAs. Our SHAPE-based structures diverge substantially from those inferred solely from computation, particularly in the apical loop and basal segments, underlining the need for experimental data in RNA structure prediction. By comparing the structures with high-throughput processing data, we determined the optimal structural features of pri-miRNAs. The sequence determinants are influenced substantially by their structural contexts. Moreover, we identified an element termed the bulged GWG motif (bGWG) with a 3' bulge in the lower stem, which promotes processing. Our structure-function mapping better annotates the determinants of pri-miRNA processing and offers practical implications for designing small hairpin RNAs and predicting the impacts of miRNA mutations.

Target-directed microRNA degradation: Mechanisms, significance, and functional implications

MicroRNAs (miRNAs) are small non-coding RNAs that play a fundamental role in enabling miRNA-mediated target repression, a post-transcriptional gene regulatory mechanism preserved across metazoans. Loss of certain animal miRNA genes can lead to developmental abnormalities, disease, and various degrees of embryonic lethality. These short RNAs normally guide Argonaute (AGO) proteins to target RNAs, which are in turn translationally repressed and destabilized, silencing the target to fine-tune gene expression and maintain cellular homeostasis. Delineating miRNA-mediated target decay has been thoroughly examined in thousands of studies, yet despite these exhaustive studies, comparatively less is known about how and why miRNAs are directed for decay. Several key observations over the years have noted instances of rapid miRNA turnover, suggesting endogenous means for animals to induce miRNA degradation. Recently, it was revealed that certain targets, so-called target-directed miRNA degradation (TDMD) triggers, can "trigger" miRNA decay through inducing proteolysis of AGO and thereby the bound miRNA. This process is mediated in animals via the ZSWIM8 ubiquitin ligase complex, which is recruited to AGO during engagement with triggers. Since its discovery, several studies have identified that ZSWIM8 and TDMD are indispensable for proper animal development. Given the rapid expansion of this field of study, here, we summarize the key findings that have led to and followed the discovery of ZSWIM8-dependent TDMD. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA in Disease and Development > RNA in Development.

The miRNA-target interactions: An underestimated intricacy

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Evaluation of microRNA variant maturation prior to genome edition

Assessment of the functionality of individual microRNA/target sites is a crucial issue. Genome editing techniques should theoretically permit a fine functional exploration of such interactions, allowing the mutation of microRNAs or individual binding sites in a complete in vivo setting, therefore abrogating or restoring individual interactions on demand. A major limitation to this experimental strategy is the influence of microRNA sequence on its accumulation level, which introduces a confounding effect when assessing phenotypic rescue by compensatorily mutated microRNA and target site. Here we describe a simple assay to identify microRNA variants most likely to accumulate at wild-type levels even though their sequence has been mutated. In this assay, quantification of a reporter construct in cultured cells predicts the efficiency of an early biogenesis step, the Drosha-dependent cleavage of microRNA precursors, which appears to be a major determinant of microRNA accumulation in our variant collection. This system allowed the generation of a mutant Drosophila strain expressing a bantam microRNA variant at wild-type levels.

The catalytic activity of microRNA Argonautes plays a modest role in microRNA star strand destabilization in C. elegans

Many Argonaute proteins can cleave RNA ("slicing") as part of the microRNA-induced silencing complex (miRISC), even though miRNA-mediated target repression is generally independent of target cleavage. Here we use genome editing in C. elegans to examine the role of miRNA-guided slicing in organismal development. In contrast to previous work, slicing-inactivating mutations did not interfere with normal development when introduced by CRISPR. We find that unwinding and decay of miRNA star strands is weakly defective in the absence of slicing, with the largest effect observed in embryos. Argonaute-Like Gene 2 (ALG-2) is more dependent on slicing for unwinding than ALG-1. The miRNAs that displayed the greatest (albeit minor) dependence on slicing for unwinding tend to form stable duplexes with their star strand, and in some cases, lowering duplex stability alleviates dependence on slicing. Gene expression changes were consistent with negligible to moderate loss of function for miRNA guides whose star strand was upregulated, suggesting a reduced proportion of mature miRISC in slicing mutants. While a few miRNA guide strands are reduced in the mutant background, the basis of this is unclear since changes were not dependent on EBAX-1, a factor in the Target-Directed miRNA Degradation (TDMD) pathway. Overall, this work defines a role for miRNA Argonaute slicing in star strand decay; future work should examine whether this role could have contributed to the selection pressure to conserve catalytic activity of miRNA Argonautes across the metazoan phylogeny.

IsoSeek for unbiased and UMI-informed sequencing of miRNAs from low input samples at single-nucleotide resolution

Besides canonical microRNAs (miRNAs), sequence-based variations called isomiRs have biological relevance and diagnostic potential; however, accurate calling of these post-transcriptional modifications is challenging, especially for low input samples. Here, we present IsoSeek, a sequencing protocol that reduces ligation and PCR amplification bias and improves the accuracy of miRNA detection in low input samples. We describe steps for using randomized adapters combined with unique molecular identifiers (UMI), library quantification, and sequencing, followed by detailed procedures for data processing and analysis. For complete details on the use and execution of this protocol, please refer to C. Gómez-Martín et al. (2023)1 and Van Eijndhoven et al. (2021).2.

Age and 17β-Estradiol (E2) Facilitate Nuclear Export and Argonaute Loading of microRNAs in the Female Brain

Aging in women is accompanied by a dramatic change in circulating sex steroid hormones. Specifically, the primary circulating estrogen, 17β-estradiol (E2), is nearly undetectable in post-menopausal women. This decline is associated with a variety of cognitive and mood disorders, yet hormone replacement therapy is only effective within a narrow window of time surrounding the menopausal transition. Our previous work identified microRNAs as a potential molecular substrate underlying the change in E2 efficacy associated with menopause in advanced age. Specifically, we showed that E2 regulated a small subset of mature miRNAs in the aging female brain. In this study, we hypothesized that E2 regulates the stability of mature miRNAs by altering their subcellular localization and their association with argonaute proteins. We also tested the hypothesis that the RNA binding protein, hnRNP A1, was an important regulator of mature miR-9-5p expression in neuronal cells. Our results demonstrated that E2 treatment affected miRNA subcellular localization and its association with argonaute proteins differently, depending on the length of time following E2 deprivation (i.e., ovariectomy). We also provide strong evidence that hnRNP A1 regulates the transcription of pri-miR-9 and likely plays a posttranscriptional role in mature miR-9-5p turnover. Taken together, these data have important implications for considering the optimal timing for hormone replacement therapy, which might be less dependent on age and more related to how long treatment is delayed following menopause.

Expression, not sequence, distinguishes miR-238 from its miR-239ab sister miRNAs in promoting longevity in Caenorhabditis elegans

MicroRNAs (miRNAs) regulate gene expression by base-pairing to target sequences in messenger RNAs (mRNAs) and recruiting factors that induce translational repression and mRNA decay. In animals, nucleotides 2-8 at the 5' end of the miRNA, called the seed region, are often necessary and sometimes sufficient for functional target interactions. MiRNAs that contain identical seed sequences are grouped into families where individual members have the potential to share targets and act redundantly. A rare exception seemed to be the miR-238/239ab family in Caenorhabditis elegans, as previous work indicated that loss of miR-238 reduced lifespan while deletion of the miR-239ab locus resulted in enhanced longevity and thermal stress resistance. Here, we re-examined these potentially opposing roles using new strains that individually disrupt each miRNA sister. We confirmed that loss of miR-238 is associated with a shortened lifespan but could detect no longevity or stress phenotypes in animals lacking miR-239a or miR-239b, individually or in combination. Additionally, dozens of genes were mis-regulated in miR-238 mutants but almost no gene expression changes were detected in either miR-239a or miR-239b mutants compared to wild type animals. We present evidence that the lack of redundancy between miR-238 and miR-239ab is independent of their sequence differences; miR-239a or miR-239b could substitute for the longevity role of miR-238 when expressed from the miR-238 locus. Altogether, these studies disqualify miR-239ab as negative regulators of aging and demonstrate that expression, not sequence, dictates the specific role of miR-238 in promoting longevity.

Consequences of depleting TNRC6, AGO, and DROSHA proteins on expression of microRNAs

The potential for microRNAs (miRNAs) to regulate gene expression remains incompletely understood. DROSHA initiates the biogenesis of miRNAs while variants of Argonaute (AGO) and trinucleotide repeat containing six (TNRC6) family proteins form complexes with miRNAs to facilitate RNA recognition and gene regulation. Here we investigate the fate of miRNAs in the absence of these critical RNAi protein factors. Knockout of DROSHA expression reduces levels of some miRNAs annotated in miRBase but not others. The identity of miRNAs with reduced expression matches the identity of miRNAs previously identified by experimental approaches. The MirGeneDB resource offers the closest alignment with experimental results. In contrast, the loss of TNRC6 proteins had much smaller effects on miRNA levels. Knocking out AGO proteins, which directly contact the mature miRNA, decreased expression of the miRNAs most strongly associated with AGO2 as determined from enhanced crosslinking immunoprecipitation (AGO2-eCLIP). Evaluation of miRNA binding to endogenously expressed AGO proteins revealed that miRNA:AGO association was similar for AGO1, AGO2, AGO3, and AGO4. Our data emphasize the need to evaluate annotated miRNAs based on approximate cellular abundance, DROSHA-dependence, and physical association with AGO when forming hypotheses related to their function.

A rapid protocol for ribosome profiling of low input samples

Ribosome profiling provides quantitative, comprehensive, and high-resolution snapshots of cellular translation by the high-throughput sequencing of short mRNA fragments that are protected by ribosomes from nucleolytic digestion. While the overall principle is simple, the workflow of ribosome profiling experiments is complex and challenging, and typically requires large amounts of sample, limiting its broad applicability. Here, we present a new protocol for ultra-rapid ribosome profiling from low-input samples. It features a robust strategy for sequencing library preparation within one day that employs solid phase purification of reaction intermediates, allowing to reduce the input to as little as 0.1 pmol of ∼30 nt RNA fragments. Hence, it is particularly suited for the analyses of small samples or targeted ribosome profiling. Its high sensitivity and its ease of implementation will foster the generation of higher quality data from small samples, which opens new opportunities in applying ribosome profiling.

Target-directed microRNA degradation regulates developmental microRNA expression and embryonic growth in mammals

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression that play critical roles in development and disease. Target-directed miRNA degradation (TDMD), a pathway in which miRNAs that bind to specialized targets with extensive complementarity are rapidly decayed, has emerged as a potent mechanism of controlling miRNA levels. Nevertheless, the biological role and scope of miRNA regulation by TDMD in mammals remains poorly understood. To address these questions, we generated mice with constitutive or conditional deletion of Zswim8, which encodes an essential TDMD factor. Loss of Zswim8 resulted in developmental defects in the heart and lungs, growth restriction, and perinatal lethality. Small RNA sequencing of embryonic tissues revealed widespread miRNA regulation by TDMD and greatly expanded the known catalog of miRNAs regulated by this pathway. These experiments also uncovered novel features of TDMD-regulated miRNAs, including their enrichment in cotranscribed clusters and examples in which TDMD underlies "arm switching," a phenomenon wherein the dominant strand of a miRNA precursor changes in different tissues or conditions. Importantly, deletion of two miRNAs, miR-322 and miR-503, rescued growth of Zswim8-null embryos, directly implicating the TDMD pathway as a regulator of mammalian body size. These data illuminate the broad landscape and developmental role of TDMD in mammals.

Target-directed microRNA degradation regulates developmental microRNA expression and embryonic growth in mammals

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression that play critical roles in development and disease. Target-directed miRNA degradation (TDMD), a pathway in which miRNAs that bind to specialized targets with extensive complementarity are rapidly decayed, has emerged as a potent mechanism of controlling miRNA levels. Nevertheless, the biological role and scope of miRNA regulation by TDMD in mammals remains poorly understood. To address these questions, we generated mice with constitutive or conditional deletion of Zswim8 , which encodes an essential TDMD factor. Loss of Zswim8 resulted in developmental defects in heart and lung, growth restriction, and perinatal lethality. Small RNA sequencing of embryonic tissues revealed widespread miRNA regulation by TDMD and greatly expanded the known catalog of miRNAs regulated by this pathway. These experiments also uncovered novel features of TDMD-regulated miRNAs, including their enrichment in co-transcribed clusters and examples in which TDMD underlies 'arm switching', a phenomenon wherein the dominant strand of a miRNA precursor changes in different tissues or conditions. Importantly, deletion of two miRNAs, miR-322 and miR-503, rescued growth of Zswim8 null embryos, directly implicating the TDMD pathway as a regulator of mammalian body size. These data illuminate the broad landscape and developmental role of TDMD in mammals.

The Exosome-like Vesicles of Giardia Assemblages A, B, and E Are Involved in the Delivering of Distinct Small RNA from Parasite to Parasite

The genetically related assemblages of the intestinal protozoa parasite Giardia lamblia are morphologically indistinguishable and are often derived from specific hosts. The Giardia assemblages are separated by large genetic distances, which might account for their relevant biological and pathogenic differences. In this work, we analyzed the RNAs cargo released into exosomal-like vesicles (ElVs) by the assemblages A and B, which differentially infect humans, and the assemblage E, which infects hoofed animals. The RNA sequencing analysis revealed that the ElVs of each assemblage contained distinct small RNA (sRNA) biotypes, suggesting a preference for specific packaging in each assemblage. These sRNAs were classified into three categories, ribosomal-small RNAs (rsRNAs), messenger-small RNAs (msRNAs), and transfer-small RNAs (tsRNAs), which may play a regulatory role in parasite communication and contribute to host-specificity and pathogenesis. Uptake experiments showed, for the first time, that ElVs were successfully internalized by the parasite trophozoites. Furthermore, we observed that the sRNAs contained inside these ElVs were first located below the plasma membrane but then distributed along the cytoplasm. Overall, the study provides new insights into the molecular mechanisms underlying the host-specificity and pathogenesis of G. lamblia and highlights the potential role of sRNAs in parasite communication and regulation.

NORMSEQ: a tool for evaluation, selection and visualization of RNA-Seq normalization methods

RNA-sequencing has become one of the most used high-throughput approaches to gain knowledge about the expression of all different RNA subpopulations. However, technical artifacts, either introduced during library preparation and/or data analysis, can influence the detected RNA expression levels. A critical step, especially in large and low input datasets or studies, is data normalization, which aims at eliminating the variability in data that is not related to biology. Many normalization methods have been developed, each of them relying on different assumptions, making the selection of the appropriate normalization strategy key to preserve biological information. To address this, we developed NormSeq, a free web-server tool to systematically assess the performance of normalization methods in a given dataset. A key feature of NormSeq is the implementation of information gain to guide the selection of the best normalization method, which is crucial to eliminate or at least reduce non-biological variability. Altogether, NormSeq provides an easy-to-use platform to explore different aspects of gene expression data with a special focus on data normalization to help researchers, even without bioinformatics expertise, to obtain reliable biological inference from their data. NormSeq is freely available at: https://arn.ugr.es/normSeq.

Reassessment of miRNA variant (isomiRs) composition by small RNA sequencing

IsomiRs, sequence variants of mature microRNAs, are usually detected and quantified using high-throughput sequencing. Many examples of their biological relevance have been reported, but sequencing artifacts identified as artificial variants might bias biological inference and therefore need to be ideally avoided. We conducted a comprehensive evaluation of 10 different small RNA sequencing protocols, exploring both a theoretically isomiR-free pool of synthetic miRNAs and HEK293T cells. We calculated that, with the exception of two protocols, less than 5% of miRNA reads can be attributed to library preparation artifacts. Randomized-end adapter protocols showed superior accuracy, with 40% of true biological isomiRs. Nevertheless, we demonstrate concordance across protocols for selected miRNAs in non-templated uridyl additions. Notably, NTA-U calling and isomiR target prediction can be inaccurate when using protocols with poor single-nucleotide resolution. Our results highlight the relevance of protocol choice for biological isomiRs detection and annotation, which has key potential implications for biomedical applications.

Role of the proline-rich disordered domain of DROSHA in intronic microRNA processing

DROSHA serves as a gatekeeper of the microRNA (miRNA) pathway by processing primary transcripts (pri-miRNAs). While the functions of structured domains of DROSHA have been well documented, the contribution of N-terminal proline-rich disordered domain (PRD) remains elusive. Here we show that the PRD promotes the processing of miRNA hairpins located within introns. We identified a DROSHA isoform (p140) lacking the PRD, which is produced by proteolytic cleavage. Small RNA sequencing revealed that p140 is significantly impaired in the maturation of intronic miRNAs. Consistently, our minigene constructs demonstrated that PRD enhances the processing of intronic hairpins, but not those in exons. Splice site mutations did not affect the PRD's enhancing effect on intronic constructs, suggesting that the PRD acts independently of splicing reaction by interacting with sequences residing within introns. The N-terminal regions from zebrafish and Xenopus DROSHA can replace the human counterpart, indicating functional conservation despite poor sequence alignment. Moreover, we found that rapidly evolving intronic miRNAs are generally more dependent on PRD than conserved ones, suggesting a role of PRD in miRNA evolution. Our study reveals a new layer of miRNA regulation mediated by a low-complexity disordered domain that senses the genomic contexts of miRNA loci.

Sequence determinant of small RNA production by DICER

RNA silencing relies on specific and efficient processing of double-stranded RNA by Dicer, which yields microRNAs (miRNAs) and small interfering RNAs (siRNAs)^1,2. However, our current knowledge of the specificity of Dicer is limited to the secondary structures of its substrates: a double-stranded RNA of approximately 22 base pairs with a 2-nucleotide 3' overhang and a terminal loop^3-11. Here we found evidence pointing to an additional sequence-dependent determinant beyond these structural properties. To systematically interrogate the features of precursor miRNAs (pre-miRNAs), we carried out massively parallel assays with pre-miRNA variants and human DICER (also known as DICER1). Our analyses revealed a deeply conserved cis-acting element, termed the 'GYM motif' (paired G, paired pyrimidine and mismatched C or A), near the cleavage site. The GYM motif promotes processing at a specific position and can override the previously identified 'ruler'-like counting mechanisms from the 5' and 3' ends of pre-miRNA^3-6. Consistently, integrating this motif into short hairpin RNA or Dicer-substrate siRNA potentiates RNA interference. Furthermore, we find that the C-terminal double-stranded RNA-binding domain (dsRBD) of DICER recognizes the GYM motif. Alterations in the dsRBD reduce processing and change cleavage sites in a motif-dependent fashion, affecting the miRNA repertoire in cells. In particular, the cancer-associated R1855L substitution in the dsRBD strongly impairs GYM motif recognition. This study uncovers an ancient principle of substrate recognition by metazoan Dicer and implicates its potential in the design of RNA therapeutics.

ALL-tRNAseq enables robust tRNA profiling in tissue samples

Transfer RNAs (tRNAs) are small adaptor RNAs essential for mRNA translation. Alterations in the cellular tRNA population can directly affect mRNA decoding rates and translational efficiency during cancer development and progression. To evaluate changes in the composition of the tRNA pool, multiple sequencing approaches have been developed to overcome reverse transcription blocks caused by the stable structures of these molecules and their numerous base modifications. However, it remains unclear whether current sequencing protocols faithfully capture tRNAs existing in cells or tissues. This is specifically challenging for clinical tissue samples that often present variable RNA qualities. For this reason, we developed ALL-tRNAseq, which combines the highly processive MarathonRT and RNA demethylation for the robust assessment of tRNA expression, together with a randomized adapter ligation strategy prior to reverse transcription to assess tRNA fragmentation levels in both cell lines and tissues. Incorporation of tRNA fragments not only informed on sample integrity but also significantly improved tRNA profiling of tissue samples. Our data showed that our profiling strategy effectively improves classification of oncogenic signatures in glioblastoma and diffuse large B-cell lymphoma tissues, particularly for samples presenting higher levels of RNA fragmentation, further highlighting the utility of ALL-tRNAseq for translational research.

Structure of the human DICER-pre-miRNA complex in a dicing state

Dicer has a key role in small RNA biogenesis, processing double-stranded RNAs (dsRNAs)^1,2. Human DICER (hDICER, also known as DICER1) is specialized for cleaving small hairpin structures such as precursor microRNAs (pre-miRNAs) and has limited activity towards long dsRNAs-unlike its homologues in lower eukaryotes and plants, which cleave long dsRNAs. Although the mechanism by which long dsRNAs are cleaved has been well documented, our understanding of pre-miRNA processing is incomplete because structures of hDICER in a catalytic state are lacking. Here we report the cryo-electron microscopy structure of hDICER bound to pre-miRNA in a dicing state and uncover the structural basis of pre-miRNA processing. hDICER undergoes large conformational changes to attain the active state. The helicase domain becomes flexible, which allows the binding of pre-miRNA to the catalytic valley. The double-stranded RNA-binding domain relocates and anchors pre-miRNA in a specific position through both sequence-independent and sequence-specific recognition of the newly identified 'GYM motif'³. The DICER-specific PAZ helix is also reoriented to accommodate the RNA. Furthermore, our structure identifies a configuration of the 5' end of pre-miRNA inserted into a basic pocket. In this pocket, a group of arginine residues recognize the 5' terminal base (disfavouring guanine) and terminal monophosphate; this explains the specificity of hDICER and how it determines the cleavage site. We identify cancer-associated mutations in the 5' pocket residues that impair miRNA biogenesis. Our study reveals how hDICER recognizes pre-miRNAs with stringent specificity and enables a mechanistic understanding of hDICER-related diseases.

isomiRdb: microRNA expression at isoform resolution

A significant fraction of mature miRNA transcripts carries sequence and/or length variations, termed isomiRs. IsomiRs are differentially abundant in cell types, tissues, body fluids or patients' samples. Not surprisingly, multiple studies describe a physiological and pathophysiological role. Despite their importance, systematically collected and annotated isomiR information available in databases remains limited. We thus developed isomiRdb, a comprehensive resource that compiles miRNA expression data at isomiR resolution from various sources. We processed 42 499 human miRNA-seq datasets (5.9 × 1011 sequencing reads) and consistently analyzed them using miRMaster and sRNAbench. Our database provides online access to the 90 483 most abundant isomiRs (>1 RPM in at least 1% of the samples) from 52 tissues and 188 cell types. Additionally, the full set of over 3 million detected isomiRs is available for download. Our resource can be queried at the sample, miRNA or isomiR level so users can quickly answer common questions about the presence/absence of a particular miRNA/isomiR in tissues of interest. Further, the database facilitates to identify whether a potentially interesting new isoform has been detected before and its frequency. In addition to expression tables, isomiRdb can generate multiple interactive visualisations including violin plots and heatmaps. isomiRdb is free to use and publicly available at: https://www.ccb.uni-saarland.de/isomirdb.

hnRNPC induces isoform shifts in miR-21-5p leading to cancer development

MicroRNA (miRNA) processing is a critical step in mature miRNA production. Its dysregulation leads to an increase in miRNA isoforms with heterogenous 5'-ends (isomiRs), which can recognize distinct target sites because of their shifted seed sequence. Although some miRNA genes display productive expression of their 5'-isomiRs in cancers, how their production is controlled and how 5'-isomiRs affect tumor progression have yet to be explored. In this study, based on integrative analyses of high-throughput sequencing data produced by our group and publicly available data, we demonstrate that primary miR-21 (pri-miR-21) is processed into the cancer-specific isomiR isomiR-21-5p | ±1, which suppresses growth hormone receptor (GHR) in liver cancer. Treatment with antagomirs against isomiR-21-5p | ±1 inhibited the in vitro tumorigenesis of liver cancer cells and allowed the recovery of GHR, whereas the introduction of isomiR-21-5p | ±1 mimics attenuated these effects. These effects were validated in a mouse model of spontaneous liver cancer. Heterogeneous nuclear ribonucleoprotein C and U2 small nuclear RNA auxiliary factor 2 were predicted to bind upstream of pre-miR-21 via a poly-(U) motif and influence Drosha processing to induce the production of isomiR-21-5p | ±1. Our findings suggest an oncogenic function for the non-canonical isomiR-21-5p | ±1 in liver cancer, and its production was shown to be regulated by hnRNPC.

High-throughput in vitro processing of human primary microRNA by the recombinant microprocessor

We describe a protocol to conduct a high-throughput in vitro processing assay, using 1,881 human primary microRNAs (pri-miRNAs) and recombinant Microprocessor complex, followed by deep sequencing library generation. This comprehensive approach allows the mapping of cleavage sites and the measurement of processing efficiency of a large number of substrates simultaneously. Our protocol is readily modifiable to investigate the effects of chemicals and regulatory proteins. Moreover, cis-acting elements can be examined by replacing the wild-type pri-miRNAs with mutant variants.

For complete details on the use and execution of this profile, please refer to Kim et al. (2021).

•

In vitro processing of human pri-miRNAs by the recombinant Microprocessor

•

Construction of cDNA libraries for high-throughput assay of pri-miRNA processing

•

This protocol allows mapping of cleavage sites and measuring processing efficiency

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.

chi-miR-487b-3p Inhibits Goat Myoblast Proliferation and Differentiation by Targeting IRS1 through the IRS1/PI3K/Akt Signaling Pathway

MicroRNAs (miRNAs) are endogenously expressed small noncoding RNAs and play critical roles in the regulation of post-transcriptional gene expression. Our previous study uncovered that chi-miR-487b-3p is widespread in different goat tissues, which is significantly higher in muscle, especially in lamb. Here, we demonstrate the role of chi-miR-487b-3p as a myogenic miRNA that regulates skeletal muscle development. chi-miR-487b-3p overexpression was demonstrated to significantly inhibit goat myoblast proliferation and differentiation, whereas chi-miR-487b-3p inhibition resulted in the opposite effects. Next, chi-miR-487b-3p was predicted to target the 3'UTR of insulin receptor substrate 1 (IRS1) gene by Target-Scan and miRDB. The results of dual-luciferase assay, RT-qPCR, and western blot all confirmed that IRS1 might be a direct target of chi-miR-487b-3p as its expression was negatively regulated by chi-miR-487b-3p. siRNA silencing of IRS1 further demonstrated significant inhibition on goat myoblast proliferation and differentiation, confirming the effect of IRS1 downregulation by chi-miR-487b-3p in myogenesis. In addition, chi-miR-487b-3p knockout goat myoblast clones were generated using CRISPR/Cas9 technology, and we further illustrated that chi-miR-487b-3p regulates goat myoblast growth through the PI3K/Akt signaling pathway by targeting IRS1. Collectively, our work demonstrated that chi-miR-487b-3p is a potent inhibitor of skeletal myogenesis and provided new insights into the mechanisms of miRNA on the regulation of goat growth.

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.

Modulating microRNAs in cancer: Next-generation therapies

MicroRNAs (miRNAs) are a class of endogenously expressed non-coding regulators of the genome with an ability to mediate a variety of biological and pathological processes. There is growing evidence demonstrating frequent dysregulation of microRNAs in cancer cells, which is associated with tumor initiation, development, migration, invasion, resisting cell death, and drug resistance. Studies have shown that modulation of these small RNAs is a novel and promising therapeutic tool in the treatment of a variety of diseases, especially cancer, due to their broad influence on multiple cellular processes. However, suboptimal delivery of the appropriate miRNA to the cancer sites, quick degradation by nucleases in the blood circulation, and off target effects have limited their research and clinical applications. Therefore, there is a pressing need to improve the therapeutic efficacy of miRNA modulators, while at the same time reducing their toxicities. Several delivery vehicles for miRNA modulators have been shown to be effective in vitro and in vivo. In this review, we will discuss the role and importance of miRNAs in cancer and provide perspectives on currently available carriers for miRNA modulation. We will also summarize the challenges and prospects for the clinical translation of miRNA-based therapeutic strategies.

The distinct RNA-interaction modes of a small ZnF domain underlay TUT4(7) diverse action in miRNA regulation

TUT4 and the closely related TUT7 are non-templated poly(U) polymerases required at different stages of development, and their mis-regulation or mutation has been linked to important cancer pathologies. While TUT4(7) interaction with its pre-miRNA targets has been characterized in detail, the molecular bases of the broader target recognition process are unclear. Here, we examine RNA binding by the ZnF domains of the protein. We show that TUT4(7) ZnF2 contains two distinct RNA binding surfaces that are used in the interaction with different RNA nucleobases in different targets, i.e that this small domain encodes diversity in TUT4(7) selectivity and molecular function. Interestingly and unlike other well-characterized CCHC ZnFs, ZnF2 is not physically coupled to the flanking ZnF3 and acts independently in miRNA recognition, while the remaining CCHC ZnF of TUT4(7), ZnF1, has lost its intrinsic RNA binding capability. Together, our data suggest that the ZnFs of TUT4(7) are independent units for RNA and, possibly, protein-protein interactions that underlay the protein's functional flexibility and are likely to play an important role in building its interaction network.

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.

The regulatory impact of RNA-binding proteins on microRNA targeting

Argonaute is the primary mediator of metazoan miRNA targeting (MT). Among the currently identified >1,500 human RNA-binding proteins (RBPs), there are only a handful of RBPs known to enhance MT and several others reported to suppress MT, leaving the global impact of RBPs on MT elusive. In this study, we have systematically analyzed transcriptome-wide binding sites for 150 human RBPs and evaluated the quantitative effect of individual RBPs on MT efficacy. In contrast to previous studies, we show that most RBPs significantly affect MT and that all of those MT-regulating RBPs function as MT enhancers rather than suppressors, by making the local secondary structure of the target site accessible to Argonaute. Our findings illuminate the unappreciated regulatory impact of human RBPs on MT, and as these RBPs may play key roles in the gene regulatory network governed by metazoan miRNAs, MT should be understood in the context of co-regulating RBPs.

miRNAs are loaded into Argonaute protein and repress complementary mRNA targets. Here the authors show the unappreciated role of RNA binding proteins for efficient miRNA targeting and expand the current understanding of miRNA targeting.

Quantitative mapping of the cellular small RNA landscape with AQRNA-seq

Current next-generation RNA-sequencing (RNA-seq) methods do not provide accurate quantification of small RNAs within a sample, due to sequence-dependent biases in capture, ligation and amplification during library preparation. We present a method, absolute quantification RNA-sequencing (AQRNA-seq), that minimizes biases and provides a direct, linear correlation between sequencing read count and copy number for all small RNAs in a sample. Library preparation and data processing were optimized and validated using a 963-member microRNA reference library, oligonucleotide standards of varying length, and RNA blots. Application of AQRNA-seq to a panel of human cancer cells revealed >800 detectable miRNAs that varied during cancer progression, while application to bacterial transfer RNA pools, with the challenges of secondary structure and abundant modifications, revealed 80-fold variation in tRNA isoacceptor levels, stress-induced site-specific tRNA fragmentation, quantitative modification maps, and evidence for stress-induced, tRNA-driven, codon-biased translation. AQRNA-seq thus provides a versatile means to quantitatively map the small RNA landscape in cells.

Bray P, Pliatsika V, Telonis AG, Rigoutsos I. IsoMiRmap: fast, deterministic and exhaustive mining of isomiRs from short RNA-seq datasets

MOTIVATION

MicroRNA (miRNA) precursor arms give rise to multiple isoforms simultaneously called 'isomiRs.' IsomiRs from the same arm typically differ by a few nucleotides at either their 5' or 3' termini or both. In humans, the identities and abundances of isomiRs depend on a person's sex and genetic ancestry as well as on tissue type, tissue state and disease type/subtype. Moreover, nearly half of the time the most abundant isomiR differs from the miRNA sequence found in public databases. Accurate mining of isomiRs from deep sequencing data is thus important.

RESULTS

We developed isoMiRmap, a fast, standalone, user-friendly mining tool that identifies and quantifies all isomiRs by directly processing short RNA-seq datasets. IsoMiRmap is a portable 'plug-and-play' tool, requires minimal setup, has modest computing and storage requirements, and can process an RNA-seq dataset with 50 million reads in just a few minutes on an average laptop. IsoMiRmap deterministically and exhaustively reports all isomiRs in a given deep sequencing dataset and quantifies them accurately (no double-counting). IsoMiRmap comprehensively reports all miRNA precursor locations from which an isomiR may be transcribed, tags as 'ambiguous' isomiRs whose sequences exist both inside and outside of the space of known miRNA sequences and reports the public identifiers of common single-nucleotide polymorphisms and documented somatic mutations that may be present in an isomiR. IsoMiRmap also identifies isomiRs with 3' non-templated post-transcriptional additions. Compared to similar tools, isoMiRmap is the fastest, reports more bona fide isomiRs, and provides the most comprehensive information related to an isomiR's transcriptional origin.

AVAILABILITY AND IMPLEMENTATION

The codes for isoMiRmap are freely available at https://cm.jefferson.edu/isoMiRmap/ and https://github.com/TJU-CMC-Org/isoMiRmap/. IsomiR profiles for the datasets of the 1000 Genomes Project, spanning five population groups, and The Cancer Genome Atlas (TCGA), spanning 33 cancer studies, are also available at https://cm.jefferson.edu/isoMiRmap/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A quantitative map of human primary microRNA processing sites

Maturation of canonical microRNA (miRNA) is initiated by DROSHA that cleaves the primary transcript (pri-miRNA). More than 1,800 miRNA loci are annotated in humans, but it remains largely unknown whether and at which sites pri-miRNAs are cleaved by DROSHA. Here, we performed in vitro processing on a full set of human pri-miRNAs (miRBase version 21) followed by sequencing. This comprehensive profiling enabled us to classify miRNAs on the basis of DROSHA dependence and map their cleavage sites with respective processing efficiency measures. Only 758 pri-miRNAs are confidently processed by DROSHA, while the majority may be non-canonical or false entries. Analyses of the DROSHA-dependent pri-miRNAs show key cis-elements for processing. We observe widespread alternative processing and unproductive cleavage events such as "nick" or "inverse" processing. SRSF3 is a broad-acting auxiliary factor modulating alternative processing and suppressing unproductive processing. The profiling data and methods developed in this study will allow systematic analyses of miRNA regulation.

Secretory microRNA Profiles of Third- and Fourth-Stage Dirofilaria immitis Larvae with Different Macrocyclic Lactone Susceptibility: In Search of Biomarkers for Early Detection of Infection

The canine heartworm, Dirofilaria immitis, is among the most important parasites of dogs in the United States and worldwide, and may cause severe and potentially fatal disease. Current diagnostic recommendations rely on serological detection of an adult female antigen, and visualization of microfilariae in the blood. Therefore, a reliable diagnosis can be only performed approximately six months post-infection. There is a growing need to characterize novel diagnostic markers that are capable of detecting the early stages of heartworm infection, in special markers associated with third-stage larvae (L3) and fourth-stage larvae (L4). The early detection of infection would guide medical interventions that could impede the development of patent infections and further parasite transmission. We cultured D. immitis L3 and L4 of two laboratorial strains with different susceptibility statuses to macrocyclic lactone drugs in vitro. Excretory/secretory microRNAs were sequenced and analyzed. We identified two miRNA novel candidates secreted abundantly by both L3 and L4 of both strains. These candidates were previously detected in the secretions of other D. immitis stages and one of them was found in the blood of D. immitis-infected dogs. These miRNAs have not been found in the secretions of other nematodes and could be D. immitis-specific diagnostic biomarkers, which could allow for the early detection of infection.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

The impact of microRNAs on myeloid-derived suppressor cells in cancer

Inflammation promotes cancer development. To a large extent, this can be attributed to the recruitment of myeloid-derived suppressor cells (MDSCs) to tumors. These cells are known for establishing an immunosuppressive tumor microenvironment by suppressing T cell activities. However, MDSCs also promote metastasis and angiogenesis. Critically, as small non-coding RNAs that regulate gene expression, microRNAs (miRNAs) control MDSC activities. In this review, we discuss how miRNA networks regulate key MDSC signaling pathways, how they shape MDSC development, differentiation and activation, and how this impacts tumor development. By targeting the expression of miRNAs in MDSCs, we can alter their main signaling pathways. In turn, this can compromise their ability to promote multiple hallmarks of cancer. Therefore, this may represent a new powerful strategy for cancer immunotherapy.

Paired-end small RNA sequencing reveals a possible overestimation in the isomiR sequence repertoire previously reported from conventional single read data analysis

BACKGROUND

Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real.

RESULTS

We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets.

CONCLUSIONS

We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account.

EZcount: An all-in-one software for microRNA expression quantification from NGS sequencing data

MicroRNAs (miRNAs) are short endogenous molecules of RNA that influence cell regulation by suppressing genes. Their ubiquity throughout all branches of the tree of life has suggested their central role in many cellular functions. Nowadays, several personalized medicine applications rely on miRNAs as biomarkers for diagnoses, prognoses, and prediction of drug response. The increasing ease of sequencing miRNAs contrasts with the difficulty of accurately quantifying their concentration. The use of general purpose aligners is only a partial solution as they have limited possibilities to accurately solve ambiguous mapping due to the short length of these sequences. We developed EZcount, an all-in-one software that, with a single command, performs the entire quantification process: from raw fastq files to read counts. Experiments show that EZcount is more sensitive and accurate than methods based on sequence alignment, independently of the library preparation protocol and sequencing machine. The parallel architecture of EZcount makes it fast enough to process a sample in minutes using a standard workstation. EZcount runs on all of the most common operating systems (Linux, Windows and MacOS) and is freely available for download at https://gitlab.com/BioAlgo/miR-pipe. A detailed description of the datasets, the raw experimental results, and all the scripts used for testing are available as supplementary material.

A ubiquitin ligase mediates target-directed microRNA decay independently of tailing and trimming

MicroRNAs (miRNAs) act in concert with Argonaute (AGO) proteins to repress target messenger RNAs. After AGO loading, miRNAs generally exhibit slow turnover. An important exception occurs when miRNAs encounter highly complementary targets, which can trigger a process called target-directed miRNA degradation (TDMD). During TDMD, miRNAs undergo tailing and trimming, suggesting that this is an important step in the decay mechanism. We identified a cullin-RING ubiquitin ligase (CRL), containing the substrate adaptor ZSWIM8, that mediates TDMD. The ZSWIM8 CRL interacts with AGO proteins, promotes TDMD in a tailing and trimming-independent manner, and regulates miRNA expression in multiple cell types. These findings suggest a model in which the ZSWIM8 ubiquitin ligase mediates TDMD by directing proteasomal decay of miRNA-containing complexes engaged with highly complementary targets.

An alternative miRISC targets a cancer-associated coding sequence mutation in FOXL2

Recent evidence suggests that animal microRNAs (miRNAs) can target coding sequences (CDSs); however, the pathophysiological importance of such targeting remains unknown. Here, we show that a somatic heterozygous missense mutation (c.402C>G; p.C134W) in FOXL2, a feature shared by virtually all adult-type granulosa cell tumors (AGCTs), introduces a target site for miR-1236, which causes haploinsufficiency of the tumor-suppressor FOXL2. This miR-1236-mediated selective degradation of the variant FOXL2 mRNA is preferentially conducted by a distinct miRNA-loaded RNA-induced silencing complex (miRISC) directed by the Argonaute3 (AGO3) and DHX9 proteins. In both patients and a mouse model of AGCT, abundance of the inversely regulated variant FOXL2 with miR-1236 levels is highly correlated with malignant features of AGCT. Our study provides a molecular basis for understanding the conserved FOXL2 CDS mutation-mediated etiology of AGCT, revealing the existence of a previously unidentified mechanism of miRNA-targeting disease-associated mutations in the CDS by forming a non-canonical miRISC.