Sequence features of Drosha and Dicer cleavage sites affect the complexity of isomiRs

Substrate promiscuity of Dicer toward precursors of the let-7 family and their 3'-end modifications

The human let-7 miRNA family consists of thirteen members that play critical roles in many biological processes, including development timing and tumor suppression, and their levels are disrupted in several diseases. Dicer is the endoribonuclease responsible for processing the precursor miRNA (pre-miRNA) to yield the mature miRNA, and thereby plays a crucial role in controlling the cellular levels of let-7 miRNAs. It is well established that the sequence and structural features of pre-miRNA hairpins such as the 5'-phosphate, the apical loop, and the 2-nt 3'-overhang are important for the processing activity of Dicer. Exceptionally, nine precursors of the let-7 family (pre-let-7) contain a 1-nt 3'-overhang and get mono-uridylated in vivo, presumably to allow efficient processing by Dicer. Pre-let-7 are also oligo-uridylated in vivo to promote their degradation and likely prevent their efficient processing by Dicer. In this study, we systematically investigated the impact of sequence and structural features of all human let-7 pre-miRNAs, including their 3'-end modifications, on Dicer binding and processing. Through the combination of SHAPE structural probing, in vitro binding and kinetic studies using purified human Dicer, we show that despite structural discrepancies among pre-let-7 RNAs, Dicer exhibits remarkable promiscuity in binding and cleaving these substrates. Moreover, the 1- or 2-nt 3'-overhang, 3'-mono-uridylation, and 3'-oligo-uridylation of pre-let-7 substrates appear to have little effect on Dicer binding and cleavage rates. Thus, this study extends current knowledge regarding the broad substrate specificity of Dicer and provides novel insight regarding the effect of 3'-modifications on binding and cleavage by Dicer.

Association of microRNA Polymorphisms with Toxicities Induced by Methotrexate in Children with Acute Lymphoblastic Leukemia

Methotrexate (MTX), a structurally related substance to folic acid, is an important chemotherapeutic agent used for decades in the treatment of pediatric acute lymphoblastic leukemia (ALL) and other types of cancer as non-Hodgkin lymphomas and osteosarcomas. Despite the successful outcomes observed, the primary drawback is the variability in the pharmacokinetics and pharmacodynamics between patients. The main adverse events related to its use are nephrotoxicity, mucositis, and myelosuppression, especially when used in high doses. The potential adverse reactions and toxicities associated with MTX are a cause for concern and may lead to dose reduction or treatment interruption. Genetic variants in MTX transport genes have been linked to toxicity. Pharmacogenetic studies conducted in the past focused on single nucleotide polymorphisms (SNPs) in the coding and 5'-regulatory regions of genes. Recent studies have demonstrated a significant role of microRNAs (miRNAs) in the transport and metabolism of drugs and in the regulation of target genes. In the last few years, the number of annotated miRNAs has continually risen, in addition to the studies of miRNA polymorphisms and MTX toxicity. Therefore, the objective of the present study is to investigate the role of miRNA variants related to MTX adverse effects.

Interrogating the Role of miR-125b and Its 3'isomiRs in Protection against Hypoxia

MiR-125b has therapeutic potential in the amelioration of myocardial ischemic injury. MicroRNA isomiRs, with either 5' or 3' addition or deletion of nucleotide(s), have been reported from next-generation sequencing data (NGS). However, due to technical challenges, validation and functional studies of isomiRs are few. In this study, we discovered using NGS, four 3'isomiRs of miR-125b, i.e., addition of A (adenosine), along with deletions of A, AG (guanosine) and AGU (uridine) from rat and sheep heart. These findings were validated using RT-qPCR. Comprehensive functional studies were carried out in the H9C2 hypoxia model. After miR-125b, isomiRs of Plus A, Trim A, AG and AGU mimic transfection, the H9C2 cells were subjected to hypoxic challenge. As assessed using cell viability, apoptosis, CCK-8 and LDH release, miR-125b and isomiRs were all protective against hypoxia. However, Plus A and Trim A were more effective than miR-125b, whilst Trim AG and Trim AGU had far weaker effects than miR-125b. Interestingly, both the gene regulation profile and apoptotic gene validation indicated a major overlap among miR-125b, Plus A and Trim A, whilst Trims AG and AGU revealed a different profile compared to miR-125b. Conclusions: miR-125b and its 3' isomiRs are expressed stably in the heart. miR-125b and isomiRs with addition or deletion of A might function concurrently and concordantly under specific physiological and pathophysiological conditions. In-depth understanding of isomiRs' metabolism and function will contribute to better miRNA therapeutic drug design.

PLK1 Regulates MicroRNA Biogenesis through Drosha Phosphorylation

Polo-Like Kinase 1 (PLK1), a key mediator of cell-cycle progression, is associated with poor prognosis and is a therapeutic target in a number of malignancies. Putative phosphorylation sites for PLK1 have been identified on Drosha, the main catalytic component of the microprocessor responsible for miR biogenesis. Several kinases, including GSK3β, p70 S6 kinase, ABL, PAK5, p38 MAPK, CSNK1A1 and ANKRD52-PPP6C, have been shown to phosphorylate components of the miR biogenesis machinery, altering their activity and/or localisation, and therefore the biogenesis of distinct miR subsets. We hypothesised that PLK1 regulates miR biogenesis through Drosha phosphorylation. In vitro kinase assays confirmed PLK1 phosphorylation of Drosha at S300 and/or S302. PLK1 inhibition reduced serine-phosphorylated levels of Drosha and its RNA-dependent association with DGCR8. In contrast, a "phospho-mimic" Drosha mutant showed increased association with DGCR8. PLK1 phosphorylation of Drosha alters Drosha Microprocessor complex subcellular localisation, since PLK1 inhibition increased cytosolic protein levels of both DGCR8 and Drosha, whilst nuclear levels were decreased. Importantly, the above effects are independent of PLK1's cell cycle-regulatory role, since altered Drosha:DGCR8 localisation upon PLK1 inhibition occurred prior to significant accumulation of cells in M-phase, and PLK1-regulated miRs were not increased in M-phase-arrested cells. Small RNA sequencing and qPCR validation were used to assess downstream consequences of PLK1 activity on miR biogenesis, identifying a set of ten miRs (miR-1248, miR-1306-5p, miR-2277-5p, miR-29c-5p, miR-93-3p, miR-152-3p, miR-509-3-5p, miR-511-5p, miR-891a-5p and miR-892a) whose expression levels were statistically significantly downregulated by two pharmacological PLK1 kinase domain inhibitors, RO-5203280 and GSK461364. Opposingly, increased levels of these miRs were observed upon transfection of wild-type or constitutively active PLK1. Importantly, pre-miR levels were reduced upon PLK1 inhibition, and pri-miR levels decreased upon PLK1 activation, and hence, PLK1 Drosha phosphorylation regulates MiR biogenesis at the level of pri-miR-to-pre-miR processing. In combination with prior studies, this work identifies Drosha S300 and S302 as major integration points for signalling by several kinases, whose relative activities will determine the relative biogenesis efficiency of different miR subsets. Identified kinase-regulated miRs have potential for use as kinase inhibitor response-predictive biomarkers, in cancer and other diseases.

One locus, several functional RNAs-emerging roles of the mechanisms responsible for the sequence variability of microRNAs

With the development of modern molecular genetics, the original "one gene-one enzyme" hypothesis has been outdated. For protein coding genes, the discovery of alternative splicing and RNA editing provided the biochemical background for the RNA repertoire of a single locus, which also serves as an important pillar for the enormous protein variability of the genomes. Non-protein coding RNA genes were also revealed to produce several RNA species with distinct functions. The loci of microRNAs (miRNAs), encoding for small endogenous regulatory RNAs, were also found to produce a population of small RNAs, rather than a single defined product. This review aims to present the mechanisms contributing to the astonishing variability of miRNAs revealed by the new sequencing technologies. One important source is the careful balance of arm selection, producing sequentially different 5p- or 3p-miRNAs from the same pre-miRNA, thereby broadening the number of regulated target RNAs and the phenotypic response. In addition, the formation of 5', 3' and polymorphic isomiRs, with variable end and internal sequences also leads to a higher number of targeted sequences, and increases the regulatory output. These miRNA maturation processes, together with other known mechanisms such as RNA editing, further increase the potential outcome of this small RNA pathway. By discussing the subtle mechanisms behind the sequence diversity of miRNAs, this review intends to reveal this engaging aspect of the inherited "RNA world", how it contributes to the almost infinite molecular variability among living organisms, and how this variability can be exploited to treat human diseases.

The context signals of mitochondrial miRNAs (mitomiRs) of mammals

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression at the post-transcriptional level in the cytoplasm and play an important role in a wide range of biological processes. Recent studies have found that the miRNA sequences are presented not only in the cytoplasm, but also in the mitochondria. These miRNAs (the so-called mitomiRs) may be the sequences of nuclear or mitochondrial origin; some of them are involved in regulation of the mitochondrial gene functions, while the role of others is still unknown. The identification of nucleotide signals, which are unique to mitomiRs, may help to determine this role. We formed a dataset that combined the experimentally discovered mitomiRs in human, rat and mouse. To isolate signals that may be responsible for the mitomiRs' functions or for their translocation from or into mitochondria a context analysis was carried out for the sequences. For three species in the group mitomiRs/non-mitomiRs and the group of all miRNAs from the miRBase database statistically overrepresented 8-letter motifs were identified (p-value < 0.01 with Bonferroni correction for multiple comparisons), for these motifs the patterns of the localization in functionally important regions for different types of miRNAs were found. Also, for the group mitomiRs/non-mitomiRs we found the statistically significant features of the miRNA nucleotide context near the Dicer and Drosha cleavage sites (Pearson's χ2 test of independence for the first three positions of the miRNA, p-value < 0.05). The observed nucleotide frequencies may indicate a more homogeneous pri-miRNA cleavage by the Drosha complex during the formation of the 5' end of mitomiRs. The obtained results can help to determine the role of the nucleotide signals in the origin, processing, and functions of the mitomiRs.

Small RNAs and their protein partners in animal meiosis

Meiosis is characterized by highly regulated transitions in gene expression that require diverse mechanisms of gene regulation. For example, in male mammals, transcription undergoes a global shut-down in early prophase I of meiosis, followed by increasing transcriptional activity into pachynema. Later, as spermiogenesis proceeds, the histones bound to DNA are replaced with transition proteins, which are themselves replaced with protamines, resulting in a highly condensed nucleus with repressed transcriptional activity. In addition, two specialized gene silencing events take place during prophase I: meiotic silencing of unsynapsed chromatin (MSUC), and the sex chromatin specific mechanism, meiotic sex chromosome inactivation (MSCI). Notably, conserved roles for the RNA binding protein (RBP) machinery that functions with small non-coding RNAs have been described as participating in these meiosis-specific mechanisms, suggesting that RNA-mediated gene regulation is critical for fertility in many species. Here, we review roles of small RNAs and their associated RBPs in meiosis-related processes such as centromere function, silencing of unpaired chromatin and meiotic recombination. We will discuss the emerging evidence of non-canonical functions of these components in meiosis.

Predicting Drosha and Dicer Cleavage Sites with DeepMirCut

MicroRNAs are a class of small RNAs involved in post-transcriptional gene silencing with roles in disease and development. Many computational tools have been developed to identify novel microRNAs. However, there have been no attempts to predict cleavage sites for Drosha from primary sequence, or to identify cleavage sites using deep neural networks. Here, we present DeepMirCut, a recurrent neural network-based software that predicts both Dicer and Drosha cleavage sites. We built a microRNA primary sequence database including flanking genomic sequences for 34,713 microRNA annotations. We compare models trained on sequence data, sequence and secondary structure data, as well as input data with annotated structures. Our best model is able to predict cuts within closer average proximity than results reported for other methods. We show that a guanine nucleotide before and a uracil nucleotide after Dicer cleavage sites on the 3′ arm of the microRNA precursor had a positive effect on predictions while the opposite order (U before, G after) had a negative effect. Our analysis was also able to predict several positions where bulges had either positive or negative effects on the score. We expect that our approach and the data we have curated will enable several future studies.

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.

Hairpin sequence and structure is associated with features of isomiR biogenesis

MiRNA isoforms (isomiRs) are single stranded small RNAs originating from the same pri-miRNA hairpin as a result of cleavage by Drosha and Dicer enzymes. Variations at the 5'-end of a miRNA alter the seed region of the molecule, thus affecting the targetome of the miRNA. In this manuscript, we analysed the distribution of miRNA cleavage positions across 31 different cancers using miRNA sequencing data of TCGA project. As a result, we found that the processing positions are not tissue specific and that all miRNAs could be correctly classified as ones exhibiting homogeneous or heterogeneous cleavage at one of the four cleavage sites. In 42% of cases (42 out of 100 miRNAs), we observed imprecise 5'-end Dicer cleavage, while this fraction was only 14% for Drosha (14 out of 99). To the contrary, almost all cleavage sites of 3'-ends (either Drosha or Dicer) were heterogeneous. With the use of only four nucleotides surrounding a 5'-end Dicer cleavage position we built a model which allowed us to distinguish between homogeneous and heterogeneous cleavage with the reliable quality (ROC AUC = 0.68). Finally, we showed the possible applications of the study by the analysis of two 5'-end isoforms originating from the same exogeneous shRNA hairpin. It turned out that the less expressed shRNA variant was functionally active, which led to the increased off-targeting. Thus, the obtained results could be applied to the design of shRNAs whose processing will result in a single 5'-variant.

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.

In silico Analysis of Polymorphisms in microRNAs Deregulated in Alzheimer Disease

BACKGROUND

Alzheimer's disease (AD) is a degenerative condition characterized by progressive cognitive impairment and dementia. Findings have revolutionized current knowledge of miRNA in the neurological conditions. Two regulatory mechanisms determine the level of mature miRNA expression; one is miRNA precursor processing, and the other is gene expression regulation by transcription factors. This study is allocated to the in-silico investigation of miRNA's SNPs and their effect on other cell mechanisms.

METHODS

We used databases which annotate the functional effect of SNPs on mRNA-miRNA and miRNA-RBP interaction. Also, we investigated SNPs which are located on the promoter or UTR region.

RESULTS

miRNA SNP3.0 database indicated several SNPs in miR-339 and miR-34a in the upstream and downstream of pre-miRNA and mature miRNAs. While, for some miRNAs miR-124, and miR-125, no polymorphism was observed, and also miR-101 with ΔG -3.1 and mir-328 with ΔG 5.8 had the highest and lowest potencies to produce mature microRNA. SNP2TFBS web-server presented several SNPs which altered the Transcription Factor Binding Sites (TFBS) or generated novel TFBS in the promoter regions of related miRNA. At last, RBP-Var database provided a list of SNPs which alter miRNA-RBP interaction pattern and can also influence other miRNAs' expression.

DISCUSSION

The results indicated that SNPs microRNA affects both miRNA function and miRNA expression. Our study expands molecular insight into how SNPs in different parts of miRNA, including the regulatory (promoter), the precursor (pre-miRNA), functional regions (seed region of mature miRNA), and RBP-binding motifs, which theoretically may be correlated to the Alzheimer's disease.

A pan-cancer atlas of somatic mutations in miRNA biogenesis genes

It is a well-known and intensively studied phenomenon that the levels of many miRNAs are differentiated in cancer. miRNA biogenesis and functional expression are complex processes orchestrated by many proteins cumulatively called miRNA biogenesis proteins. To characterize cancer somatic mutations in the miRNA biogenesis genes and investigate their potential impact on the levels of miRNAs, we analyzed whole-exome sequencing datasets of over 10 000 cancer/normal sample pairs deposited within the TCGA repository. We identified and characterized over 3600 somatic mutations in 29 miRNA biogenesis genes and showed that some of the genes are overmutated in specific cancers and/or have recurrent hotspot mutations (e.g. SMAD4 in PAAD, COAD and READ; DICER1 in UCEC; PRKRA in OV and LIN28B in SKCM). We identified a list of miRNAs whose level is affected by particular types of mutations in either SMAD4, SMAD2 or DICER1 and showed that hotspot mutations in the RNase domains in DICER1 not only decrease the level of 5p-miRNAs but also increase the level of 3p-miRNAs, including many well-known cancer-related miRNAs. We also showed an association of the mutations with patient survival. Eventually, we created an atlas/compendium of miRNA biogenesis alterations providing a useful resource for different aspects of biomedical research.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

Genetic variability in noncoding RNAs: involvement of miRNAs and long noncoding RNAs in osteoporosis pathogenesis

The pathogenesis of osteoporosis is multifactorial and is the consequence of genetic, hormonal and lifestyle factors. Epigenetics, including noncoding RNA (ncRNA) deregulation, represents a link between susceptibility to develop the disease and environmental influences. The majority of studies investigated the expression of ncRNAs in osteoporosis patients; however, very little information is available on their genetic variability. In this review, we focus on two classes of ncRNAs: miRNAs and long noncoding RNAs (lncRNAs). We summarize recent findings on how polymorphisms in miRNAs and lncRNAs can perturb the lncRNA/miRNA/mRNA axis and may be involved in osteoporosis clinical outcome. We also provide a general overview on databases and bioinformatic tools useful for associating miRNAs and lncRNAs variability with complex genetic diseases.

Postprandial transfer of colostral extracellular vesicles and their protein and miRNA cargo in neonatal calves

Extracellular vesicles (EVs) such as exosomes are key regulators of intercellular communication that can be found in almost all bio fluids. Although studies in the last decade have made great headway in discerning the role of EVs in many physiological and pathophysiological processes, the bioavailability and impact of dietary EVs and their cargo still remain to be elucidated. Due to its widespread consumption and high content of EV-associated microRNAs and proteins, a major focus in this field has been set on EVs in bovine milk and colostrum. Despite promising in vitro studies in recent years that show high resiliency of milk EVs to degradation and uptake of milk EV cargo in a variety of intestinal and blood cell types, in vivo experiments continue to be inconclusive and sometimes outright contradictive. To resolve this discrepancy, we assessed the potential postprandial transfer of colostral EVs to the circulation of newborn calves by analysing colostrum-specific protein and miRNAs, including specific isoforms (isomiRs) in cells, EV isolations and unfractionated samples from blood and colostrum. Our findings reveal distinct populations of EVs in colostrum and blood from cows that can be clearly separated by density, particle concentration and protein content (BTN1A1, MFGE8). Postprandial blood samples of calves show a time-dependent increase in EVs that share morphological and protein characteristics of colostral EVs. Analysis of miRNA expression profiles by Next-Generation Sequencing gave a different picture however. Although significant postprandial expression changes could only be detected for calf EV samples, expression profiles show very limited overlap with highly expressed miRNAs in colostral EVs or colostrum in general. Taken together our results indicate a selective uptake of membrane-associated protein cargo but not luminal miRNAs from colostral EVs into the circulation of neonatal calves.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

Somatic Mutations in miRNA Genes in Lung Cancer-Potential Functional Consequences of Non-Coding Sequence Variants

A growing body of evidence indicates that miRNAs may either drive or suppress oncogenesis. However, little is known about somatic mutations in miRNA genes. To determine the frequency and potential consequences of miRNA gene mutations, we analyzed whole exome sequencing datasets of 569 lung adenocarcinoma (LUAD) and 597 lung squamous cell carcinoma (LUSC) samples generated in The Cancer Genome Atlas (TCGA) project. Altogether, we identified 1091 somatic sequence variants affecting 522 different miRNA genes and showed that half of all cancers had at least one such somatic variant/mutation. These sequence variants occurred in most crucial parts of miRNA precursors, including mature miRNA and seed sequences. Due to our findings, we hypothesize that seed mutations may affect miRNA:target interactions, drastically changing the pool of predicted targets. Mutations may also affect miRNA biogenesis by changing the structure of miRNA precursors, DROSHA and DICER cleavage sites, and regulatory sequence/structure motifs. We identified 10 significantly overmutated hotspot miRNA genes, including the miR-379 gene in LUAD enriched in mutations in the mature miRNA and regulatory sequences. The occurrence of mutations in the hotspot miRNA genes was also shown experimentally. We present a comprehensive analysis of somatic variants in miRNA genes and show that some of these genes are mutational hotspots, suggesting their potential role in cancer.

The in vivo dsRNA Cleavage Has Sequence Preference in Insects

Exogenous dsRNA enters the insect body and can induce the RNAi effect only when it is cleaved into siRNA. However, what kinds of base composition are easier to cut and what kinds of siRNA will be produced in vivo is largely unknown. In this study, we found that dsRNA processing into siRNA has sequence preference and regularity in insects. We injected 0.04 mg/g dsRNA into Asian corn borers or cotton bollworms according to their body weight, and then the siRNAs produced in vivo were analyzed by RNA-Seq. We discovered that a large number of siRNAs were produced with GGU nucleotide residues at the 5′- and 3′-ends and produced a siRNA peak on the sequence. Once the GGU site is mutated, the number of siRNAs will decrease significantly and the siRNA peak will also lost. However, in the red flour beetle, a member of Coleoptera, dsRNA was cut at more diverse sites, such as AAG, GUG, and GUU; more importantly, these enzyme restriction sites have a high conservation base of A/U. Our discovery regarding dsRNA in vivo cleavage preference and regularity will help us understand the RNAi mechanism and its application.

Evolutionary Divergence of Brain-specific Precursor miRNAs Drives Efficient Processing and Production of Mature miRNAs in Human

The hallmark of human evolution encompasses the dramatic increase in brain size and complexity. The intricate interplays of micro-RNAs (miRNAs) and their target genes are indispensable in brain development. Sequence divergence in distinct structural regions of Brain-specific precursor miRNAs (pre-miRNAs) and its consequence in the production of corresponding mature miRNAs in human are unknown. To address these questions, first we classified miRNAs into three categories based on tissue expression: Brain-specific (expressed exclusively in brain), Non-brain (expressed in Non-brain tissues) and Common (expressed in all tissues) and compared the sequence divergence of different structural regions (basal segment, lower and upper stem, internal and terminal loop) of categorized pre-miRNAs across human, non-human primates and rodents. Our analysis revealed that unpaired regions of Brain-specific pre-miRNAs in human bear traces of relatively high rate of evolutionary divergence compared to those in other species. Cross-tissue expression analysis unveiled the higher expression of the Brain-specific miRNAs in human compared to other species. Intriguingly, in human brain, expression levels of these miRNAs superseded the levels of the ubiquitously expressed "Common-miRNAs". Further analysis revealed that presence of certain motif and nucleotide preference in the Brain-specific pre-miRNAs may favor DROSHA and DICER to ameliorate miRNA processing. The higher processing efficiency of human Brain-specific miRNAs was reflected as an elevated production of corresponding mature miRNAs in the human brain. Finally, re-construction of gene-regulatory network uncovers different pathways driven by Brain-specific miRNAs that may contribute to the development of brain in human.

Small RNA-seq analysis of single porcine blastocysts revealed that maternal estradiol-17beta exposure does not affect miRNA isoform (isomiR) expression

BACKGROUND

The expression of microRNAs (miRNAs) is essential for the proper development of the mammalian embryo. A maternal exposure to endocrine disrupting chemicals during preimplantation bears the potential for transgenerational inheritance of disease through the epigenetic perturbation of the developing embryo. A comprehensive assembly of embryo-specific miRNAs and respective isoforms (isomiR) is lacking to date. We aimed at revealing the sex-specific miRNA expression profile of single porcine blastocysts developing in gilts orally exposed to exogenous estradiol-17 (E2). Therefore we analyzed the miRNA profile specifically focusing on isomiRs and potentially embryo-specific miRNAs.

RESULTS

Deep sequencing of small RNA (small RNA-seq) result in the detection of miRNA sequences mapping to known and predicted porcine miRNAs as well as novel miRNAs highly conserved in human and cattle. A set of highly abundant miRNAs and a large number of rarely expressed miRNAs were identified by using a small RNA analysis pipeline, which was integrated into a novel Galaxy workflow specifically benefits incompletely annotated species. In particular, orthologue species information increased the total number of annotated miRNAs, while mapping to other non-coding RNAs avoided falsely annotated miRNAs. Neither the low nor the high dose of E2 treatment (10 and 1000 µ E2/kg body weight daily, respectively) affected the miRNA profile in blastocysts despite the distinct differential mRNA expression and DNA methylation found in previous studies. The high number of generated sequence reads enabled a comprehensive analysis of the isomiR repertoire showing various templated and non-templated modifications. Furthermore, potentially blastocyst-specific miRNAs were identified.

CONCLUSIONS

In pre-implantation embryos, numerous distinct isomiRs were discovered indicating a high complexity of miRNA expression. Neither the sex of the embryo nor a maternal E2 exposure affected the miRNA expression profile of developing porcine blastocysts. The adaptation to the continuous duration of the E2 treatment might explain the lack of an effect.

3' RNA Uridylation in Epitranscriptomics, Gene Regulation, and Disease

Emerging evidence implicates a wide range of post-transcriptional RNA modifications that play crucial roles in fundamental biological processes including regulating gene expression. Collectively, they are known as epitranscriptomics. Recent studies implicate 3' RNA uridylation, the non-templated addition of uridine(s) to the terminal end of RNA, as a key player in epitranscriptomics. In this review, we describe the functional roles and significance of 3' terminal RNA uridylation that has diverse functions in regulating both mRNAs and non-coding RNAs. In mammals, three Terminal Uridylyl Transferases (TUTases) are primarily responsible for 3' RNA uridylation. These enzymes are also referred to as polyU polymerases. TUTase 1 (TUT1) is implicated in U6 snRNA maturation via uridylation. The TUTases TUT4 and/or TUT7 are the predominant mediators of all other cellular uridylation. Terminal uridylation promotes turnover for many polyadenylated mRNAs, replication-dependent histone mRNAs that lack polyA-tails, and aberrant structured noncoding RNAs. In addition, uridylation regulates biogenesis of a subset of microRNAs and generates isomiRs, sequent variant microRNAs that have altered function in specific cases. For example, the RNA binding protein and proto-oncogene LIN28A and TUT4 work together to polyuridylate pre-let-7, thereby blocking biogenesis and function of the tumor suppressor let-7 microRNA family. In contrast, monouridylation of Group II pre-miRNAs creates an optimal 3' overhang that promotes recognition and subsequent cleavage by the Dicer-TRBP complex that then yields the mature microRNA. Also, uridylation may play a role in non-canonical microRNA biogenesis. The overall significance of 3' RNA uridylation is discussed with an emphasis on mammalian development, gene regulation, and disease, including cancer and Perlman syndrome. We also introduce recent changes to the HUGO-approved gene names for multiple terminal nucleotidyl transferases that affects in part TUTase nomenclature (TUT1/TENT1, TENT2/PAPD4/GLD2, TUT4/ZCCHC11/TENT3A, TUT7/ZCCHC6/TENT3B, TENT4A/PAPD7, TENT4B/PAPD5, TENT5A/FAM46A, TENT5B/FAM46B, TENT5C/FAM46C, TENT5D/FAM46D, MTPAP/TENT6/PAPD1).

Classification of miRNA-related sequence variations

miRNA regulome is whole set of regulatory elements that regulate miRNA expression or are under control of miRNAs. Its understanding is vital for comprehension of miRNA functions. Classification of miRNA-related genetic variability is challenging because miRNA interact with different genomic elements and are studied at different omics levels. In the present study, miRNA-associated genetic variability is presented at three levels: miRNA genes and their upstream regulation, miRNA silencing machinery and miRNA targets. Several types of miRNA-associated genetic variations are known, including short and structural polymorphisms and epimutations. Differential expression can also affect miRNA regulome function. Classification of miRNA-associated genetic variability presents a baseline for complementing sequence variant nomenclature, planning of experiments, protocols for multi-omics data integration and development of biomarkers.

Comparing miRNA structure of mirtrons and non-mirtrons

BACKGROUND

MicroRNAs proceeds through the different canonical and non-canonical pathways; the most frequent of the non-canonical ones is the splicing-dependent biogenesis of mirtrons. We compare the mirtrons and non-mirtrons of human and mouse to explore how their maturation appears in the precursor structure around the miRNA.

RESULTS

We found the coherence of the overhang lengths what indicates the dependence between the cleavage sites. To explain this dependence we suggest the 2-lever model of the Dicer structure that couples the imprecisions in Drosha and Dicer. Considering the secondary structure of all animal pre-miRNAs we confirmed that single-stranded nucleotides tend to be located near the miRNA boundaries and in its center and are characterized by a higher mutation rate. The 5' end of the canonical 5' miRNA approaches the nearest single-stranded nucleotides what suggests the extension of the loop-counting rule from the Dicer to the Drosha cleavage site. A typical structure of the annotated mirtron pre-miRNAs differs from the canonical pre-miRNA structure and possesses the 1- and 2 nt hanging ends at the hairpin base. Together with the excessive variability of the mirtron Dicer cleavage site (that could be partially explained by guanine at its ends inherited from splicing) this is one more evidence for the 2-lever model. In contrast with the canonical miRNAs the mirtrons have higher snp densities and their pre-miRNAs are inversely associated with diseases. Therefore we supported the view that mirtrons are under positive selection while canonical miRNAs are under negative one and we suggested that mirtrons are an intrinsic source of silencing variability which produces the disease-promoting variants. Finally, we considered the interference of the pre-miRNA structure and the U2snRNA:pre-mRNA basepairing. We analyzed the location of the branchpoints and found that mirtron structure tends to expose the branchpoint site what suggests that the mirtrons can readily evolve from occasional hairpins in the immediate neighbourhood of the 3' splice site.

CONCLUSION

The miRNA biogenesis manifests itself in the footprints of the secondary structure. Close inspection of these structural properties can help to uncover new pathways of miRNA biogenesis and to refine the known miRNA data, in particular, new non-canonical miRNAs may be predicted or the known miRNAs can be re-classified.

Stress-induced changes in miRNA biogenesis and functioning

MicroRNAs (miRNAs) are small, noncoding RNAs that play key roles in the regulation of cellular homeostasis in eukaryotic organisms. There is emerging evidence that some of these processes are influenced by various forms of cellular stresses, including DNA damage, pathogen invasion or chronic stress associated with diseases. Many reports over the last decade demonstrate examples of stress-induced miRNA deregulation at the level of transcription, processing, subcellular localization and functioning. Moreover, core miRNA biogenesis proteins and their interactions with partners can be selectively regulated in response to stress signaling. However, little is known about the role of isomiRs and the interactions of miRNA with non-canonical targets in the context of the stress response. In this review, we summarize the current knowledge on miRNA functions under various stresses, including chronic stress and miRNA deregulation in the pathogenesis of age-associated neurodegenerative disorders.

Naturally existing isoforms of miR-222 have distinct functions

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

Expression of short hairpin RNAs using the compact architecture of retroviral microRNA genes

Short hairpin RNAs (shRNAs) are effective in generating stable repression of gene expression. RNA polymerase III (RNAP III) type III promoters (U6 or H1) are typically used to drive shRNA expression. While useful for some knockdown applications, the robust expression of U6/H1-driven shRNAs can induce toxicity and generate heterogeneous small RNAs with undesirable off-target effects. Additionally, typical U6/H1 promoters encompass the majority of the ∼270 base pairs (bp) of vector space required for shRNA expression. This can limit the efficacy and/or number of delivery vector options, particularly when delivery of multiple gene/shRNA combinations is required. Here, we develop a compact shRNA (cshRNA) expression system based on retroviral microRNA (miRNA) gene architecture that uses RNAP III type II promoters. We demonstrate that cshRNAs coded from as little as 100 bps of total coding space can precisely generate small interfering RNAs (siRNAs) that are active in the RNA-induced silencing complex (RISC). We provide an algorithm with a user-friendly interface to design cshRNAs for desired target genes. This cshRNA expression system reduces the coding space required for shRNA expression by >2-fold as compared to the typical U6/H1 promoters, which may facilitate therapeutic RNAi applications where delivery vector space is limiting.

Comprehensive miRNA expression profiling in human T-cell acute lymphoblastic leukemia by small RNA-sequencing

T-cell acute lymphoblastic leukemia (T-ALL) is a genetically heterogeneous disease that can be classified into different molecular genetic subtypes according to their mRNA gene expression profile. In this study, we applied RNA sequencing to investigate the full spectrum of miRNA expression in primary T-ALL patient samples, T-ALL leukemia cell lines and healthy donor thymocytes. Notably, this analysis revealed that genetic subtypes of human T-ALL also display unique miRNA expression signatures, which are largely conserved in human T-ALL cell lines with corresponding genetic background. Furthermore, small RNA-sequencing also unraveled the variety of isoforms that are expressed for each miRNA in T-ALL and showed that a significant number of miRNAs are actually represented by an alternative isomiR. Finally, comparison of CD34⁺ and CD4⁺CD8⁺ healthy donor thymocytes and T-ALL miRNA profiles allowed identifying several novel miRNAs with putative oncogenic or tumor suppressor functions in T-ALL. Altogether, this study provides a comprehensive overview of miRNA expression in normal and malignant T-cells and sets the stage for functional evaluation of novel miRNAs in T-ALL disease biology.

Knowledge about the presence or absence of miRNA isoforms (isomiRs) can successfully discriminate amongst 32 TCGA cancer types

Isoforms of human miRNAs (isomiRs) are constitutively expressed with tissue- and disease-subtype-dependencies. We studied 10 271 tumor datasets from The Cancer Genome Atlas (TCGA) to evaluate whether isomiRs can distinguish amongst 32 TCGA cancers. Unlike previous approaches, we built a classifier that relied solely on ‘binarized’ isomiR profiles: each isomiR is simply labeled as ‘present’ or ‘absent’. The resulting classifier successfully labeled tumor datasets with an average sensitivity of 90% and a false discovery rate (FDR) of 3%, surpassing the performance of expression-based classification. The classifier maintained its power even after a 15× reduction in the number of isomiRs that were used for training. Notably, the classifier could correctly predict the cancer type in non-TCGA datasets from diverse platforms. Our analysis revealed that the most discriminatory isomiRs happen to also be differentially expressed between normal tissue and cancer. Even so, we find that these highly discriminating isomiRs have not been attracting the most research attention in the literature. Given their ability to successfully classify datasets from 32 cancers, isomiRs and our resulting ‘Pan-cancer Atlas’ of isomiR expression could serve as a suitable framework to explore novel cancer biomarkers.

Shortcuts to a functional adipose tissue: The role of small non-coding RNAs

Metabolic diseases such as type 2 diabetes are a major public health issue worldwide. These diseases are often linked to a dysfunctional adipose tissue. Fat is a large, heterogenic, pleiotropic and rather complex tissue. It is found in virtually all cavities of the human body, shows unique plasticity among tissues, and harbors many cell types in addition to its main functional unit - the adipocyte. Adipose tissue function varies depending on the localization of the fat depot, the cell composition of the tissue and the energy status of the organism. While the white adipose tissue (WAT) serves as the main site for triglyceride storage and acts as an important endocrine organ, the brown adipose tissue (BAT) is responsible for thermogenesis. Beige adipocytes can also appear in WAT depots to sustain heat production upon certain conditions, and it is becoming clear that adipose tissue depots can switch phenotypes depending on cell autonomous and non-autonomous stimuli. To maintain such degree of plasticity and respond adequately to changes in the energy balance, three basic processes need to be properly functioning in the adipose tissue: i) adipogenesis and adipocyte turnover, ii) metabolism, and iii) signaling. Here we review the fundamental role of small non-coding RNAs (sncRNAs) in these processes, with focus on microRNAs, and demonstrate their importance in adipose tissue function and whole body metabolic control in mammals.

The variant of pri-mir-26a-1 polymorphism is associated with decreased risk of betel quid-related oral premalignant lesions and oral squamous cell carcinoma

OBJECTIVES

This case-control study evaluated the association of the single nucleotide polymorphism rs7372209 (T>C) in pri-mir-26a-1 with the risk and progression of betel quid (BQ)-related oral premalignant lesions (OPLs) and oral squamous cell carcinoma (OSCC).

STUDY DESIGN

In total, 597 BQ chewers were recruited: 196 healthy controls, 241 patients with OPLs, and 160 patients with OSCC. Genotypes were determined using the TaqMan real-time assay.

RESULTS

The C/T + T/T genotypes and T allele in pri-mir-26a-1 were correlated with a decreased risk of BQ-related OPLs (P = .038 and .005, respectively), oral leukoplakia (P = .01 and .001, respectively), and advanced-stage OSCC (P = .021 and .004, respectively). The effects of the C/T + T/T genotypes and T allele on the decreased risk of OPLs were potent in the older age group (both P_interaction < .001), heavy smokers (P_interaction ≤ .003 and .006, respectively) and alcohol drinkers (P_interaction ≤ .004 and .001, respectively). Furthermore, among patients with OSCC, the C/T + T/T genotypes and T allele were associated with a decreased risk of advanced pathologic stage (P = .032) and lymph node involvement (P = .017).

CONCLUSIONS

BQ chewers carrying the T allele or C/T + T/T genotypes in pri-mir-26a-1 may have a decreased risk of oral leukoplakia, OPLs, and advanced-stage OSCC.

Structural Characteristics of Simple RNA Repeats Associated with Disease and their Deleterious Protein Interactions

Short Tandem Repeats (STRs) are frequent entities in many transcripts, however, in some cases, pathological events occur when a critical repeat length is reached. This phenomenon is observed in various neurological disorders, such as myotonic dystrophy type 1 (DM1), fragile X-associated tremor/ataxia syndrome, C9orf72-related amyotrophic lateral sclerosis and frontotemporal dementia (C9ALS/FTD), and polyglutamine diseases, such as Huntington's disease (HD) and spinocerebellar ataxias (SCA). The pathological effects of these repeats are triggered by mutant RNA transcripts and/or encoded mutant proteins, which depend on the localization of the expanded repeats in non-coding or coding regions. A growing body of recent evidence revealed that the RNA structures formed by these mutant RNA repeat tracts exhibit toxic effects on cells. Therefore, in this review article, we present existing knowledge on the structural aspects of different RNA repeat tracts as revealed mainly using well-established biochemical and biophysical methods. Furthermore, in several cases, it was shown that these expanded RNA structures are potent traps for a variety of RNA-binding proteins and that the sequestration of these proteins from their normal intracellular environment causes alternative splicing aberration, inhibition of nuclear transport and export, or alteration of a microRNA biogenesis pathway. Therefore, in this review article, we also present the most studied examples of abnormal interactions that occur between mutant RNAs and their associated proteins.

QuickMIRSeq: a pipeline for quick and accurate quantification of both known miRNAs and isomiRs by jointly processing multiple samples from microRNA sequencing

BACKGROUND

Genome-wide miRNA expression data can be used to study miRNA dysregulation comprehensively. Although many open-source tools for microRNA (miRNA)-seq data analyses are available, challenges remain in accurate miRNA quantification from large-scale miRNA-seq dataset. We implemented a pipeline called QuickMIRSeq for accurate quantification of known miRNAs and miRNA isoforms (isomiRs) from multiple samples simultaneously.

RESULTS

QuickMIRSeq considers the unique nature of miRNAs and combines many important features into its implementation. First, it takes advantage of high redundancy of miRNA reads and introduces joint mapping of multiple samples to reduce computational time. Second, it incorporates the strand information in the alignment step for more accurate quantification. Third, reads potentially arising from background noise are filtered out to improve the reliability of miRNA detection. Fourth, sequences aligned to miRNAs with mismatches are remapped to a reference genome to further reduce false positives. Finally, QuickMIRSeq generates a rich set of QC metrics and publication-ready plots.

CONCLUSIONS

The rich visualization features implemented allow end users to interactively explore the results and gain more insights into miRNA-seq data analyses. The high degree of automation and interactivity in QuickMIRSeq leads to a substantial reduction in the time and effort required for miRNA-seq data analysis.

isomiR-SEA: an RNA-Seq analysis tool for miRNAs/isomiRs expression level profiling and miRNA-mRNA interaction sites evaluation

Background

Massive parallel sequencing of transcriptomes, revealed the presence of many miRNAs and miRNAs variants named isomiRs with a potential role in several cellular processes through their interaction with a target mRNA. Many methods and tools have been recently devised to detect and quantify miRNAs from sequencing data. However, all of them are implemented on top of general purpose alignment methods, thus providing poorly accurate results and no information concerning isomiRs and conserved miRNA-mRNA interaction sites.

Results

To overcome these limitations we present a novel algorithm named isomiR-SEA, that is able to provide users with very accurate miRNAs expression levels and both isomiRs and miRNA-mRNA interaction sites precise classifications. Tags are mapped on the known miRNAs sequences thanks to a specialized alignment algorithm developed on top of biological evidence concerning miRNAs structure. Specifically, isomiR-SEA checks for miRNA seed presence in the input tags and evaluates, during all the alignment phases, the positions of the encountered mismatches, thus allowing to distinguish among the different isomiRs and conserved miRNA-mRNA interaction sites.

Conclusions

isomiR-SEA performances have been assessed on two public RNA-Seq datasets proving that the implemented algorithm is able to account for more reliable and accurate miRNAs expression levels with respect to those provided by two compared state of the art tools. Moreover, differently from the few methods currently available to perform isomiRs detection, the proposed algorithm implements the evaluation of isomiRs and conserved miRNA-mRNA interaction sites already in the first alignment phases, thus avoiding any additional filtering stages potentially responsible for the loss of useful information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0958-0) contains supplementary material, which is available to authorized users.

Nucleotide sequence of miRNA precursor contributes to cleavage site selection by Dicer

The ribonuclease Dicer excises mature miRNAs from a diverse group of precursors (pre-miRNAs), most of which contain various secondary structure motifs in their hairpin stem. In this study, we analyzed Dicer cleavage in hairpin substrates deprived of such motifs. We searched for the factors other than the secondary structure, which may influence the length diversity and heterogeneity of miRNAs. We found that the nucleotide sequence at the Dicer cleavage site influences both of these miRNA characteristics. With regard to cleavage mechanism, we demonstrate that the Dicer RNase IIIA domain that cleaves within the 3′ arm of the pre-miRNA is more sensitive to the nucleotide sequence of its substrate than is the RNase IIIB domain. The RNase IIIA domain avoids releasing miRNAs with G nucleotide and prefers to generate miRNAs with a U nucleotide at the 5′ end. We also propose that the sequence restrictions at the Dicer cleavage site might be the factor that contributes to the generation of miRNA duplexes with 3′ overhangs of atypical lengths. This finding implies that the two RNase III domains forming the single processing center of Dicer may exhibit some degree of flexibility, which allows for the formation of these non-standard 3′ overhangs.

Viral Small-RNA Analysis of Bombyx mori Larval Midgut during Persistent and Pathogenic Cytoplasmic Polyhedrosis Virus Infection

The lepidopteran innate immune response against RNA viruses remains poorly understood, while in other insects several studies have highlighted an essential role for the exo-RNAi pathway in combating viral infection. Here, by using deep-sequencing technology for viral small-RNA (vsRNA) assessment, we provide evidence that exo-RNAi is operative in the silkworm Bombyx mori against both persistent and pathogenic infection of B. mori cytoplasmic polyhedrosis virus (BmCPV) which is characterized by a segmented double-stranded RNA (dsRNA) genome. Further, we show that Dicer-2 predominantly targets viral dsRNA and produces 20-nucleotide (nt) vsRNAs, whereas an additional pathway is responsive to viral mRNA derived from segment 10. Importantly, vsRNA distributions, which define specific hot and cold spot profiles for each viral segment, to a considerable degree overlap between Dicer-2-related (19 to 21 nt) and Dicer-2-unrelated vsRNAs, suggesting a common origin for these profiles. We found a degenerate motif significantly enriched at the cut sites of vsRNAs of various lengths which link an unknown RNase to the origins of vsRNAs biogenesis and distribution. Accordingly, the indicated RNase activity may be an important early factor for the host's antiviral defense in Lepidoptera.

IMPORTANCE

This work contributes to the elucidation of the lepidopteran antiviral response against infection of segmented double-stranded RNA (dsRNA) virus (CPV; Reoviridae) and highlights the importance of viral small-RNA (vsRNA) analysis for getting insights into host-pathogen interactions. Three vsRNA pathways are implicated in antiviral defense. For dsRNA, two pathways are proposed, either based on Dicer-2 cleavage to generate 20-nucleotide vsRNAs or based on the activity of an uncharacterized endo-RNase that cleaves the viral RNA substrate at a degenerate motif. The analysis also indicates the existence of a degradation pathway that targets the positive strand of segment 10.