Resources for the Comprehensive Discovery of Functional RNA Elements.

Transcriptome-wide maps of RNA binding protein (RBP)-RNA interactions by immunoprecipitation (IP)-based methods such as RNA IP (RIP) and crosslinking and IP (CLIP) are key starting points for evaluating the molecular roles of the thousands of human RBPs. A significant bottleneck to the application of these methods in diverse cell lines, tissues, and developmental stages is the availability of validated IP-quality antibodies. Using IP followed by immunoblot assays, we have developed a validated repository of 438 commercially available antibodies that interrogate 365 unique RBPs. In parallel, 362 short-hairpin RNA (shRNA) constructs against 276 unique RBPs were also used to confirm specificity of these antibodies. These antibodies can characterize subcellular RBP localization. With the burgeoning interest in the roles of RBPs in cancer, neurobiology, and development, these resources are invaluable to the broad scientific community. Detailed information about these resources is publicly available at the ENCODE portal (https://www.encodeproject.org/).

piRBase: integrating piRNA annotation in all aspects.

Piwi-interacting RNAs are a type of small noncoding RNA that have various functions. piRBase is a manually curated resource focused on assisting piRNA functional analysis. piRBase release v3.0 is committed to providing more comprehensive piRNA related information. The latest release covers >181 million unique piRNA sequences, including 440 datasets from 44 species. More disease-related piRNAs and piRNA targets have been collected and displayed. The regulatory relationships between piRNAs and targets have been visualized. In addition to the reuse and expansion of the content in the previous version, the latest version has additional new content, including gold standard piRNA sets, piRNA clusters, piRNA variants, splicing-junction piRNAs, and piRNA expression data. In addition, the entire web interface has been redesigned to provide a better experience for users. piRBase release v3.0 is free to access, browse, search, and download at http://bigdata.ibp.ac.cn/piRBase.

piRNAclusterDB 2.0: update and expansion of the piRNA cluster database.

PIWI-interacting RNAs (piRNAs) and their partnering PIWI proteins defend the animal germline against transposable elements and play a crucial role in fertility. Numerous studies in the past have uncovered many additional functions of the piRNA pathway, including gene regulation, anti-viral defense, and somatic transposon repression. Further, comparative analyses across phylogenetic groups showed that the PIWI/piRNA system evolves rapidly and exhibits great evolutionary plasticity. However, the presence of so-called piRNA clusters as the major source of piRNAs is common to nearly all metazoan species. These genomic piRNA-producing loci are highly divergent across taxa and critically influence piRNA populations in different evolutionary lineages. We launched the initial version of the piRNA cluster database to facilitate research on regulation and evolution of piRNA-producing loci across tissues und species. In recent years the amount of small RNA sequencing data that was generated and the abundance of species that were studied has grown rapidly. To keep up with this recent progress, we have released a major update for the piRNA cluster database (https://www.smallrnagroup.uni-mainz.de/piRNAclusterDB), expanding it from 12 to a total of 51 species with hundreds of new datasets, and revised its overall structure to enable easy navigation through this large amount of data.

miRNATissueAtlas2: an update to the human miRNA tissue atlas.

Small non-coding RNAs (sncRNAs) are pervasive regulators of physiological and pathological processes. We previously developed the human miRNA Tissue Atlas, detailing the expression of miRNAs across organs in the human body. Here, we present an updated resource containing sequencing data of 188 tissue samples comprising 21 organ types retrieved from six humans. Sampling the organs from the same bodies minimizes intra-individual variability and facilitates the making of a precise high-resolution body map of the non-coding transcriptome. The data allow shedding light on the organ- and organ system-specificity of piwi-interacting RNAs (piRNAs), transfer RNAs (tRNAs), microRNAs (miRNAs) and other non-coding RNAs. As use case of our resource, we describe the identification of highly specific ncRNAs in different organs. The update also contains 58 samples from six tissues of the Tabula Muris collection, allowing to check if the tissue specificity is evolutionary conserved between Homo sapiens and Mus musculus. The updated resource of 87 252 non-coding RNAs from nine non-coding RNA classes for all organs and organ systems is available online without any restrictions (https://www.ccb.uni-saarland.de/tissueatlas2).

cncRNAdb: a manually curated resource of experimentally supported RNAs with both protein-coding and noncoding function.

RNA endowed with both protein-coding and noncoding functions is referred to as 'dual-function RNA', 'binary functional RNA (bifunctional RNA)' or 'cncRNA (coding and noncoding RNA)'. Recently, an increasing number of cncRNAs have been identified, including both translated ncRNAs (ncRNAs with coding functions) and untranslated mRNAs (mRNAs with noncoding functions). However, an appropriate database for storing and organizing cncRNAs is still lacking. Here, we developed cncRNAdb, a manually curated database of experimentally supported cncRNAs, which aims to provide a resource for efficient manipulation, browsing and analysis of cncRNAs. The current version of cncRNAdb documents about 2600 manually curated entries of cncRNA functions with experimental evidence, involving more than 2,000 RNAs (including over 1300 translated ncRNAs and over 600 untranslated mRNAs) across over 20 species. In summary, we believe that cncRNAdb will help elucidate the functions and mechanisms of cncRNAs and develop new prediction methods. The database is available at http://www.rna-society.org/cncrnadb/.

Computational tools for plant small RNA detection and categorization.

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

LncRBase V.2: an updated resource for multispecies lncRNAs and ClinicLSNP hosting genetic variants in lncRNAs for cancer patients.

The recent discovery of long non-coding RNA as a regulatory molecule in the cellular system has altered the concept of the functional aptitude of the genome. Since our publication of the first version of LncRBase in 2014, there has been an enormous increase in the number of annotated lncRNAs of multiple species other than Human and Mouse. LncRBase V.2 hosts information of 549,648 lncRNAs corresponding to six additional species besides Human and Mouse, viz. Rat, Fruitfly, Zebrafish, Chicken, Cow and C.elegans. It provides additional distinct features such as (i) Transcription Factor Binding Site (TFBS) in the lncRNA promoter region, (ii) sub-cellular localization pattern of lncRNAs (iii) lnc-pri-miRNAs (iv) Possible small open reading frames (sORFs) within lncRNA. (v) Manually curated information of interacting target molecules and disease association of lncRNA genes (vi) Distribution of lncRNAs across multiple tissues of all species. Moreover, we have hosted ClinicLSNP within LncRBase V.2. ClinicLSNP has a comprehensive catalogue of lncRNA variants present within breast, ovarian, and cervical cancer inferred from 561 RNA-Seq data corresponding to these cancers. Further, we have checked whether these lncRNA variants overlap with (i)Repeat elements,(ii)CGI, (iii)TFBS within lncRNA loci (iv)SNP localization in trait-associated Linkage Disequilibrium(LD) region, (v)predicted the potentially pathogenic variants and (vi)effect of SNP on lncRNA secondary structure. Overall, LncRBaseV.2 is a user-friendly database to survey, search and retrieve information about multi-species lncRNAs. Further, ClinicLSNP will serve as a useful resource for cancer specific lncRNA variants and their related information. The database is freely accessible and available at http://dibresources.jcbose.ac.in/zhumur/lncrbase2/.

RNA interference in the nucleus: roles for small RNAs in transcription, epigenetics and beyond.

A growing number of functions are emerging for RNA interference (RNAi) in the nucleus, in addition to well-characterized roles in post-transcriptional gene silencing in the cytoplasm. Epigenetic modifications directed by small RNAs have been shown to cause transcriptional repression in plants, fungi and animals. Additionally, increasing evidence indicates that RNAi regulates transcription through interaction with transcriptional machinery. Nuclear small RNAs include small interfering RNAs (siRNAs) and PIWI-interacting RNAs (piRNAs) and are implicated in nuclear processes such as transposon regulation, heterochromatin formation, developmental gene regulation and genome stability.

Phased secondary small interfering RNAs in Panaxnotoginseng.

BACKGROUND: Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species. METHODS: Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5'-RACE to validate some of the identified phasiRNA targets. RESULTS: After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5'-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets. CONCLUSIONS: We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Comprehensive Annotation of Physcomitrella patens Small RNA Loci Reveals That the Heterochromatic Short Interfering RNA Pathway Is Largely Conserved in Land Plants.

Many plant small RNAs are sequence-specific negative regulators of target mRNAs and/or chromatin. In angiosperms, the two most abundant endogenous small RNA populations are usually 21-nucleotide microRNAs (miRNAs) and 24-nucleotide heterochromatic short interfering RNAs (siRNAs). Heterochromatic siRNAs are derived from repetitive regions and reinforce DNA methylation at targeted loci. The existence and extent of heterochromatic siRNAs in other land plant lineages has been unclear. Using small RNA-sequencing (RNA-seq) of the moss Physcomitrella patens, we identified 1090 loci that produce mostly 23- to 24-nucleotide siRNAs. These loci are mostly in intergenic regions with dense DNA methylation. Accumulation of siRNAs from these loci depends upon P. patens homologs of DICER-LIKE3 (DCL3), RNA-DEPENDENT RNA POLYMERASE2, and the largest subunit of DNA-DEPENDENT RNA POLYMERASE IV, with the largest subunit of a Pol V homolog contributing to expression at a smaller subset of the loci. A MINIMAL DICER-LIKE (mDCL) gene, which lacks the N-terminal helicase domain typical of DCL proteins, is specifically required for 23-nucleotide siRNA accumulation. We conclude that heterochromatic siRNAs, and their biogenesis pathways, are largely identical between angiosperms and P. patens, with the notable exception of the P. patens-specific use of mDCL to produce 23-nucleotide siRNAs.

A piece of the pi(e): The diverse roles of animal piRNAs and their PIWI partners.

Small non-coding RNAs are indispensable to many biological processes. A class of endogenous small RNAs, termed PIWI-interacting RNAs (piRNAs) because of their association with PIWI proteins, has known roles in safeguarding the genome against inordinate transposon mobilization, embryonic development, and stem cell regulation, among others. This review discusses the biogenesis of animal piRNAs and their diverse functions together with their PIWI protein partners, both in the germline and in somatic cells, and highlights the evolutionarily conserved aspects of these molecular players in animal biology.

From immunity to susceptibility: virus resistance induced in tomato by a silenced transgene is lost as TGS overcomes PTGS.

Tomato line 30.4 was obtained engineering the nucleocapsid (N) gene of tomato spotted wilt virus into plant genome, and immunity to tomato spotted wilt virus infection of its self-pollinated homozygous progeny was observed. Despite the presence of a high amount of transgenic transcripts, transgenic proteins have not been detected, suggesting a mechanism of resistance mediated by RNA. In the present study, we identify post-transcriptional gene silencing as the main mechanism of resistance, which is able to spread systemically through grafting, and show that the line 30.4 resistant plants produce both 24 and 21-22 nt N-gene specific siRNA classes. The transgenic locus in chromosome 4 shows complex multiple insertions of four T-DNA copies in various orientations, all with 3' end deletions in the terminator and part of the N gene. However, for three of them, polyadenylated transcripts are produced, due to flanking tomato genome sequences acting as alternative terminators. Interestingly, starting at the fifth generation after the transformation event, some individual plants show a tomato spotted wilt virus-susceptible phenotype. The change is associated with the disappearance of transgene-specific transcripts and siRNAs, and with hyper-methylation of the transgene, which proceeds gradually through the generations. Once it reaches a critical threshold, the shift from post-transcriptional gene silencing to transcriptional silencing of the transgene eliminates the previously well established virus resistance.

piRNAQuest: searching the piRNAome for silencers.

BACKGROUND: PIWI-interacting RNA (piRNA) is a novel and emerging class of small non-coding RNA (sncRNA). Ranging in length from 26-32 nucleotides, this sncRNA is a potent player in guiding the vital regulatory processes within a cellular system. Inspite of having such a wide role within cellular systems, piRNAs are not well organized and classified, so that a researcher can pool out the biologically relevant information concerning this class. DESCRIPTION: Here we present piRNAQuest- a unified and comprehensive database of 41749 human, 890078 mouse and 66758 rat piRNAs obtained from NCBI and different small RNA sequence experiments. This database provides piRNA annotation based on their localization in gene, intron, intergenic, CDS, 5/UTR, 3/UTR and repetitive regions which has not been done so far. We have also annotated piRNA clusters and have elucidated characteristic motifs within them. We have looked for the presence of piRNAs and piRNA clusters in pseudogenes, which are known to regulate the expression of protein coding transcripts by generating small RNAs. All these will help researchers progress towards solving the unanswered queries on piRNA biogenesis and their mode of action. Further, expression profile for piRNA in different tissues and from different developmental stages has been provided. In addition, we have provided several tools like 'homology search', 'dynamic cluster search' and 'pattern search'. Overall, piRNAQuest will serve as a useful resource for exploring human, mouse and rat piRNAome. The database is freely accessible and available at http://bicresources.jcbose.ac.in/zhumur/pirnaquest/. CONCLUSION: piRNAs play a remarkable role in stem cell self-renewal and various vital processes of developmental biology. Although researchers are mining different features on piRNAs, the exact regulatory mechanism is still fuzzy. Thus, understanding the true potential of these small regulatory molecules with respect to their origin, localization and mode of biogenesis is crucial. piRNAQuest will provide us with a better insight on piRNA origin and function which will help to explore the true potential of these sncRNAs.

mut-16 and other mutator class genes modulate 22G and 26G siRNA pathways in Caenorhabditis elegans.

Argonaute-associated siRNAs and Piwi-associated piRNAs have overlapping roles in silencing mobile genetic elements in animals. In Caenorhabditis elegans, mutator (mut) class genes mediate siRNA-guided repression of transposons as well as exogenous RNAi, but their roles in endogenous RNA silencing pathways are not well-understood. To characterize the endogenous small RNAs dependent on mut class genes, small RNA populations from a null allele of mut-16 as well as a regulatory mut-16(mg461) allele that disables only somatic RNAi were subjected to deep sequencing. Additionally, each of the mut class genes was tested for a requirement in 26G siRNA pathways. The results indicate that mut-16 is an essential factor in multiple endogenous germline and somatic siRNA pathways involving several distinct Argonautes and RNA-dependent RNA polymerases. The results also reveal essential roles for mut-2 and mut-7 in the ERGO-1 class 26G siRNA pathway and less critical roles for mut-8, mut-14, and mut-15. We show that transposons are hypersusceptible to mut-16-dependent silencing and identify a requirement for the siRNA machinery in piRNA biogenesis from Tc1 transposons. We also show that the soma-specific mut-16(mg461) mutant allele is present in multiple C. elegans laboratory strains.

Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway.

In flies, small silencing RNAs are sorted between Argonaute1 (Ago1), the central protein component of the microRNA (miRNA) pathway, and Argonaute2 (Ago2), which mediates RNA interference. Extensive double-stranded character-as is found in small interfering RNAs (siRNAs)-directs duplexes into Ago2, whereas central mismatches, like those found in miRNA/miRNA* duplexes, direct duplexes into Ago1. Central to this sorting decision is the affinity of the small RNA duplex for the Dcr-2/R2D2 heterodimer, which loads small RNAs into Ago2. Here, we show that while most Drosophila miRNAs are bound to Ago1, miRNA* strands accumulate bound to Ago2. Like siRNA loading, efficient loading of miRNA* strands in Ago2 favors duplexes with a paired central region and requires both Dcr-2 and R2D2. Those miRNA and miRNA* sequences bound to Ago2, like siRNAs diced in vivo from long double-stranded RNA, typically begin with cytidine, whereas Ago1-bound miRNA and miRNA* disproportionately begin with uridine. Consequently, some pre-miRNA generate two or more isoforms from the same side of the stem that differentially partition between Ago1 and Ago2. Our findings provide the first genome-wide test for the idea that Drosophila small RNAs are sorted between Ago1 and Ago2 according to their duplex structure and the identity of their first nucleotide.

Distinct RNA-dependent RNA polymerases are required for RNAi triggered by double-stranded RNA versus truncated transgenes in Paramecium tetraurelia.

In many eukaryotes, RNA-dependent RNA polymerases (RdRPs) play key roles in the RNAi pathway. They have been implicated in the recognition and processing of aberrant transcripts triggering the process, and in amplification of the silencing response. We have tested the functions of RdRP genes from the ciliate Paramecium tetraurelia in experimentally induced and endogenous mechanisms of gene silencing. In this organism, RNAi can be triggered either by high-copy, truncated transgenes or by directly feeding cells with double-stranded RNA (dsRNA). Surprisingly, dsRNA-induced silencing depends on the putatively functional RDR1 and RDR2 genes, which are required for the accumulation of both primary siRNAs and a distinct class of small RNAs suggestive of secondary siRNAs. In contrast, a third gene with a highly divergent catalytic domain, RDR3, is required for siRNA accumulation when RNAi is triggered by truncated transgenes. Our data further implicate RDR3 in the accumulation of previously described endogenous siRNAs and in the regulation of the surface antigen gene family. While only one of these genes is normally expressed in any clonal cell line, the knockdown of RDR3 leads to co-expression of multiple antigens. These results provide evidence for a functional specialization of Paramecium RdRP genes in distinct RNAi pathways operating during vegetative growth.

26G endo-siRNAs regulate spermatogenic and zygotic gene expression in Caenorhabditis elegans.

Endogenous small interfering RNAs (endo-siRNAs) regulate diverse gene expression programs in eukaryotes by either binding and cleaving mRNA targets or mediating heterochromatin formation; however, the mechanisms of endo-siRNA biogenesis, sorting, and target regulation remain poorly understood. Here we report the identification and function of a specific class of germline-generated endo-siRNAs in Caenorhabditis elegans that are 26 nt in length and contain a guanine at the first nucleotide position (i.e., 26G RNAs). 26G RNAs regulate gene expression during spermatogenesis and zygotic development, and their biogenesis requires the ERI-1 exonuclease and the RRF-3 RNA-dependent RNA polymerase (RdRP). Remarkably, we identified two nonoverlapping subclasses of 26G RNAs that sort into specific RNA-induced silencing complexes (RISCs) and differentially regulate distinct mRNA targets. Class I 26G RNAs target genes are expressed during spermatogenesis, whereas class II 26G RNAs are maternally inherited and silence gene expression during zygotic development. These findings implicate a class of endo-siRNAs in the global regulation of transcriptional programs required for fertility and development.

Silencing-associated and meiosis-specific small RNA pathways in Paramecium tetraurelia.

Distinct small RNA pathways are involved in the two types of homology-dependent effects described in Paramecium tetraurelia, as shown by a functional analysis of Dicer and Dicer-like genes and by the sequencing of small RNAs. The siRNAs that mediate post-transcriptional gene silencing when cells are fed with double-stranded RNA (dsRNA) were found to comprise two subclasses. DCR1-dependent cleavage of the inducing dsRNA generates approximately 23-nt primary siRNAs from both strands, while a different subclass of approximately 24-nt RNAs, characterized by a short untemplated poly-A tail, is strictly antisense to the targeted mRNA, suggestive of secondary siRNAs that depend on an RNA-dependent RNA polymerase. An entirely distinct pathway is responsible for homology-dependent regulation of developmental genome rearrangements after sexual reproduction. During early meiosis, the DCL2 and DCL3 genes are required for the production of a highly complex population of approximately 25-nt scnRNAs from all types of germline sequences, including both strands of exons, introns, intergenic regions, transposons and Internal Eliminated Sequences. A prominent 5'-UNG signature, and a minor fraction showing the complementary signature at positions 21-23, indicate that scnRNAs are cleaved from dsRNA precursors as duplexes with 2-nt 3' overhangs at both ends, followed by preferential stabilization of the 5'-UNG strand.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development.

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Predicting siRNA potency with random forests and support vector machines.

BACKGROUND: Short interfering RNAs (siRNAs) can be used to knockdown gene expression in functional genomics. For a target gene of interest, many siRNA molecules may be designed, whereas their efficiency of expression inhibition often varies. RESULTS: To facilitate gene functional studies, we have developed a new machine learning method to predict siRNA potency based on random forests and support vector machines. Since there were many potential sequence features, random forests were used to select the most relevant features affecting gene expression inhibition. Support vector machine classifiers were then constructed using the selected sequence features for predicting siRNA potency. Interestingly, gene expression inhibition is significantly affected by nucleotide dimer and trimer compositions of siRNA sequence. CONCLUSIONS: The findings in this study should help design potent siRNAs for functional genomics, and might also provide further insights into the molecular mechanism of RNA interference.