Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

Background

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

Results

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

Conclusions

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.     

Efficient reduction of target RNAs by small interfering RNA and RNase H-dependent antisense agents. A comparative analysis

RNA interference can be considered as an antisense mechanism of action that utilizes a double-stranded RNase to promote hydrolysis of the target RNA. We have performed a comparative study of optimized antisense oligonucleotides designed to work by an RNA interference mechanism to oligonucleotides designed to work by an RNase H-dependent mechanism in human cells. The potency, maximal effectiveness, duration of action, and sequence specificity of optimized RNase H-dependent oligonucleotides and small interfering RNA (siRNA) oligonucleotide duplexes were evaluated and found to be comparable. Effects of base mismatches on activity were determined to be position-dependent for both siRNA oligonucleotides and RNase H-dependent oligonucleotides. In addition, we determined that the activity of both siRNA oligonucleotides and RNase H-dependent oligonucleotides is affected by the secondary structure of the target mRNA. To determine whether positions on target RNA identified as being susceptible for RNase H-mediated degradation would be coincident with siRNA target sites, we evaluated the effectiveness of siRNAs designed to bind the same position on the target mRNA as RNase H-dependent oligonucleotides. Examination of 80 siRNA oligonucleotide duplexes designed to bind to RNA from four distinct human genes revealed that, in general, activity correlated with the activity to RNase H-dependent oligonucleotides designed to the same site, although some exceptions were noted. The one major difference between the two strategies is that RNase H-dependent oligonucleotides were determined to be active when directed against targets in the pre-mRNA, whereas siRNAs were not. These results demonstrate that siRNA oligonucleotide- and RNase H-dependent antisense strategies are both valid strategies for evaluating function of genes in cell-based assays.     

Small RNA deep sequencing identifies microRNAs and other small noncoding RNAs from human herpesvirus 6B

Roseolovirus, or human herpesvirus 6 (HHV-6), is a ubiquitous human pathogen infecting over 95% of the population by the age of 2 years. As with other herpesviruses, reactivation of HHV-6 can present with severe complications in immunocompromised individuals. Recent studies have highlighted the importance of herpesvirus-derived microRNAs (miRNAs) in modulating both cellular and viral gene expression. An initial report which computed the likelihood of various viruses to encode miRNAs did not predict HHV-6 miRNAs. To experimentally screen for small HHV-6-encoded RNAs, we conducted large-scale sequencing of Sup-T-1 cells lytically infected with a laboratory strain of HHV-6B. This revealed an abundant, 60- to 65-nucleotide RNA of unknown function derived from the lytic origin of replication (OriLyt) that gave rise to smaller RNA species of 18 or 19 nucleotides. In addition, we identified four pre-miRNAs whose mature forms accumulated in Argonaute 2. In contrast to the case for other betaherpesviruses, HHV-6B miRNAs are expressed from direct repeat regions (DR(L) and DR(R)) located at either side of the genome. All miRNAs are conserved in the closely related HHV-6A variant, and one of them is a seed ortholog of the human miRNA miR-582-5p. Similar to alphaherpesvirus miRNAs, they are expressed in antisense orientation relative to immediate-early open reading frames (ORFs) and thus have the potential to regulate key viral genes.     

Basis for regulated RNA cleavage by functional analysis of RNase L and Ire1p

RNase L and Ire1p are members of a superfamily of regulated endoribonucleases that play essential roles in mediating diverse types of cellular stress responses. 2'-5' oligoadenylates, produced in response to interferon treatment and viral double-stranded RNA, are necessary to activate RNase L. In contrast, unfolded proteins in the endoplasmic reticulum activate Ire1p, a transmembrane serine/threonine kinase and endoribonuclease. To probe their similarities and differences, molecular properties of wild-type and mutant forms of human RNase L and yeast Ire1p were compared. Surprisingly, RNase L and Ire1p showed mutually exclusive RNA substrate specificity and partially overlapping but not identical requirements for phylogenetically conserved amino acid residues in their nuclease domains. A functional model for RNase L was generated based on the comparative analysis with Ire1p that assigns novel roles for ankyrin repeats and kinase-like domains.     

Structural motifs in ribosomal RNAs: implications for RNA design and genomics

The various motifs of RNA molecules are closely related to their structural and functional properties. To better understand the nature and distributions of such structural motifs (i.e., paired and unpaired bases in stems, junctions, hairpin loops, bulges, and internal loops) and uncover characteristic features, we analyze the large 16S and 23S ribosomal RNAs of Escherichia coli. We find that the paired and unpaired bases in structural motifs have characteristic distribution shapes and ranges; for example, the frequency distribution of paired bases in stems declines linearly with the number of bases, whereas that for unpaired bases in junctions has a pronounced peak. Significantly, our survey reveals that the ratio of total (over the entire molecule) unpaired to paired bases (0.75) and the fraction of bases in stems (0.6), junctions (0.16), hairpin loops (0.12), and bulges/internal loops (0.12) are shared by 16S and 23S ribosomal RNAs, suggesting that natural RNAs may maintain certain proportions of bases in various motifs to ensure structural integrity. These findings may help in the design of novel RNAs and in the search (via constraints) for RNA-coding motifs in genomes, problems of intense current focus.     

Synthesis of chimeric RNAs between U6 small nuclear RNA and (-)sTRSV and analysis of their cleavage activities against the substrate RNA.

U6 small nuclear RNA (U6 snRNA) is one of the spliceosomal RNAs essential for pre-mRNA splicing. Highly conserved region of U6 snRNA shows a structural similarity with the catalytic center of the negative strand of the satellite RNA of tobacco ring spot virus [(-)sTRSV], supporting the hypothesis that U6 snRNA has a catalytic role in pre-mRNA splicing. To test this hypothesis, we examined in vitro whether synthetic RNAs consisting of the sequence of the highly conserved region of U6 snRNA or various chimeric RNAs between the U6 region and the catalytic center of (-)sTRSV could cleave a substrate RNA that can partially base-pair with them and has a GU sequence between the pairing regions. Chimeric RNAs with 70 to 83% sequence identity with the conserved region of S. pombe U6 snRNA cleaved the substrate RNA at the 5' side of the GU sequence. In addition, we found that the highly conserved region of U6 snRNA is similar in structure to the catalytic core region of the group I self-splicing intron in cyanobacteria. These results support the hypothesis that U6 snRNA catalyzes the pre-mRNA splicing reaction and U6 snRNA may originate from the catalytic domain of an ancient self-splicing intron.     

Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing

Distinct classes of small RNAs, 20-32 nucleotides long, play important regulatory roles for diverse cellular processes. It is therefore important to identify and quantify small RNAs as a function of development, tissue and cell type, in normal and disease states. Here we describe methods to prepare cDNA libraries from pools of small RNAs isolated from organisms, tissues or cells. These methods enable the identification of new members or new classes of small RNAs, and they are also suitable to obtain miRNA expression profiles based on clone count frequencies. This protocol includes the use of new deep sequencing methods (454/Roche and Solexa) to facilitate the characterization of diverse sequence pools of small RNAs.     

Molecular basis for target RNA recognition and cleavage by human RISC

The RNA-Induced Silencing Complex (RISC) is a ribonucleoprotein particle composed of a single-stranded short interfering RNA (siRNA) and an endonucleolytically active Argonaute protein, capable of cleaving mRNAs complementary to the siRNA. The mechanism by which RISC cleaves a target RNA is well understood, however it remains enigmatic how RISC finds its target RNA. Here, we show, both in vitro and in vivo, that the accessibility of the target site correlates directly with the efficiency of cleavage, demonstrating that RISC is unable to unfold structured RNA. In the course of target recognition, RISC transiently contacts single-stranded RNA nonspecifically and promotes siRNA-target RNA annealing. Furthermore, the 5' part of the siRNA within RISC creates a thermodynamic threshold that determines the stable association of RISC and the target RNA. We therefore provide mechanistic insights by revealing features of RISC and target RNAs that are crucial to achieve efficiency and specificity in RNA interference.     

MicroRNAs in Amoebozoa: Deep sequencing of the small RNA population in the social amoeba Dictyostelium discoideum reveals developmentally regulated microRNAs

The RNA interference machinery has served as a guardian of eukaryotic genomes since the divergence from prokaryotes. Although the basic components have a shared origin, silencing pathways directed by small RNAs have evolved in diverse directions in different eukaryotic lineages. Micro (mi)RNAs regulate protein-coding genes and play vital roles in plants and animals, but less is known about their functions in other organisms. Here, we report, for the first time, deep sequencing of small RNAs from the social amoeba Dictyostelium discoideum. RNA from growing single-cell amoebae as well as from two multicellular developmental stages was sequenced. Computational analyses combined with experimental data reveal the expression of miRNAs, several of them exhibiting distinct expression patterns during development. To our knowledge, this is the first report of miRNAs in the Amoebozoa supergroup. We also show that overexpressed miRNA precursors generate miRNAs and, in most cases, miRNA* sequences, whose biogenesis is dependent on the Dicer-like protein DrnB, further supporting the presence of miRNAs in D. discoideum. In addition, we find miRNAs processed from hairpin structures originating from an intron as well as from a class of repetitive elements. We believe that these repetitive elements are sources for newly evolved miRNAs.     

Deep sequencing analysis of small non-coding RNAs reveals the diversity of microRNAs and piRNAs in the human epididymis

The epididymis plays a crucial role in regulating the development of sperm motility and fertilizing capacity. Small non-coding RNAs (sncRNAs), especially microRNAs (miRNAs), can participate in the regulation of various physiological pathways. However, their abundance and whether they are involved in the regulation of gene expression in the human epididymis are unknown. By adopting the Solexa deep sequencing approach, we systematically investigated the sncRNAs in the adult human epididymis. A total of 4903 unique sequences representing 527 known miRNA were discovered. Eighteen novel miRNA genes encoding 23 mature miRNAs were also identified and the expression of some of them was confirmed by qRT-PCR. The presence of Piwi-interacting RNAs (piRNAs) in the library also adds to the diversity of the sncRNA population in the human epididymis. This research will contribute to a preliminary database for their functional study in male reproductive system.