Biogenesis and function of rice small RNAs from non-coding RNA precursors

Evolution, functions, and mysteries of plant ARGONAUTE proteins

ARGONAUTE (AGO) proteins bind small RNAs (sRNAs) to form RNA-induced silencing complexes for transcriptional and post-transcriptional gene silencing. Genomes of primitive plants encode only a few AGO proteins. The Arabidopsis thaliana genome encodes ten AGO proteins, designated AGO1 to AGO10. Most early studies focused on these ten proteins and their interacting sRNAs. AGOs in other flowering plant species have duplicated and diverged from this set, presumably corresponding to new, diverged or specific functions. Among these, the grass-specific AGO18 family has been discovered and implicated as playing important roles during plant reproduction and viral defense. This review covers our current knowledge about functions and features of AGO proteins in both eudicots and monocots and compares their similarities and differences. On the basis of these features, we propose a new nomenclature for some plant AGOs.

Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA precursors

Artificial microRNAs (amiRNAs) are used for selective gene silencing in plants. However, current methods to produce amiRNA constructs for silencing transcripts in monocot species are not suitable for simple, cost-effective and large-scale synthesis. Here, a series of expression vectors based on Oryza sativa MIR390 (OsMIR390) precursor was developed for high-throughput cloning and high expression of amiRNAs in monocots. Four different amiRNA sequences designed to target specifically endogenous genes and expressed from OsMIR390-based vectors were validated in transgenic Brachypodium distachyon plants. Surprisingly, amiRNAs accumulated to higher levels and were processed more accurately when expressed from chimeric OsMIR390-based precursors that include distal stem-loop sequences from Arabidopsis thaliana MIR390a (AtMIR390a). In all cases, transgenic plants displayed the predicted phenotypes induced by target gene repression, and accumulated high levels of amiRNAs and low levels of the corresponding target transcripts. Genome-wide transcriptome profiling combined with 5'-RLM-RACE analysis in transgenic plants confirmed that amiRNAs were highly specific.

Genetic basis of sRNA quantitative variation analyzed using an experimental population derived from an elite rice hybrid

We performed a genetic analysis of sRNA abundance in flag leaf from an immortalized F₂ (IMF2) population in rice. We identified 53,613,739 unique sRNAs and 165,797 sRNA expression traits (s-traits). A total of 66,649 s-traits mapped 40,049 local-sQTLs and 30,809 distant-sQTLs. By defining 80,362 sRNA clusters, 22,263 sRNA cluster QTLs (scQTLs) were recovered for 20,249 of all the 50,139 sRNA cluster expression traits (sc-traits). The expression levels for most of s-traits from the same genes or the same sRNA clusters were slightly positively correlated. While genetic co-regulation between sRNAs from the same mother genes and between sRNAs and their mother genes was observed for a portion of the sRNAs, most of the sRNAs and their mother genes showed little co-regulation. Some sRNA biogenesis genes were located in distant-sQTL hotspots and showed correspondence with specific length classes of sRNAs suggesting their important roles in the regulation and biogenesis of the sRNAs.

DOI:http://dx.doi.org/10.7554/eLife.03913.001

Genes within the DNA of a plant or animal contain instructions to make molecules called RNAs. Some RNA molecules can be decoded to make proteins, whereas others have different roles. A single gene often contains the instructions to make both protein-coding RNAs and non-coding RNAs.

Molecules called small RNAs (or sRNAs) do not code for proteins. Instead, sRNAs can control protein-coding RNA molecules or chemically alter the DNA itself; this allows them to perform many different roles in living organisms. In plants, for example, these molecules affect how the plant grows, the shapes and structures it forms, and how likely it is to survive challenges such as drought and diseases. Often different plants of the same species have different amounts of sRNAs, but the reasons for this remain unclear.

Now, Wang, Yao et al. have made use of a technique called ‘expression quantitative locus’ analysis to look at how sRNAs in rice plants are controlled by additional information encoded within DNA. The analysis identified over 53 million sRNA molecules from a population of rice plants. Many of these sRNAs varied in their abundance between different plants within the population. Wang, Yao et al. also found many thousands of individual instructions within the DNA of the rice that can either increase or reduce the abundance of their associated sRNA.

Some of the abundant sRNAs were influenced by instructions within their own genes; some were influenced by instructions from other genes; and some were influenced by both. Wang, Yao et al. also found that the control of protein-coding RNAs was not necessarily related to the control of sRNAs encoded by the same gene. Further work is now needed to identify which specific DNA sequences regulate the abundance of sRNA molecules in plants and other organisms.

DOI:http://dx.doi.org/10.7554/eLife.03913.002

Spatiotemporally dynamic, cell-type-dependent premeiotic and meiotic phasiRNAs in maize anthers

Maize anthers, the male reproductive floral organs, express two classes of phased small-interfering RNAs (phasiRNAs). PhasiRNA precursors are transcribed by RNA polymerase II and map to low-copy, intergenic regions similar to PIWI-interacting RNAs (piRNAs) in mammalian testis. From 10 sequential cohorts of staged maize anthers plus mature pollen we find that 21-nt phased siRNAs from 463 loci appear abruptly after germinal and initial somatic cell fate specification and then diminish, whereas 24-nt phasiRNAs from 176 loci coordinately accumulate during meiosis and persist as anther somatic cells mature and haploid gametophytes differentiate into pollen. Male-sterile ocl4 anthers defective in epidermal signaling lack 21-nt phasiRNAs. Male-sterile mutants with subepidermal defects--mac1 (excess meiocytes), ms23 (defective pretapetal cells), and msca1 (no normal soma or meiocytes)--lack 24-nt phasiRNAs. ameiotic1 mutants (normal soma, no meiosis) accumulate both 21-nt and 24-nt phasiRNAs, ruling out meiotic cells as a source or regulator of phasiRNA biogenesis. By in situ hybridization, miR2118 triggers of 21-nt phasiRNA biogenesis localize to epidermis; however, 21-PHAS precursors and 21-nt phasiRNAs are abundant subepidermally. The miR2275 trigger, 24-PHAS precursors, and 24-nt phasiRNAs all accumulate preferentially in tapetum and meiocytes. Therefore, each phasiRNA type exhibits independent spatiotemporal regulation with 21-nt premeiotic phasiRNAs dependent on epidermal and 24-nt meiotic phasiRNAs dependent on tapetal cell differentiation. Maize phasiRNAs and mammalian piRNAs illustrate putative convergent evolution of small RNAs in male reproduction.

Alterations in siRNA and miRNA expression profiles detected by deep sequencing of transgenic rice with siRNA-mediated viral resistance

RNA-mediated gene silencing has been demonstrated to serve as a defensive mechanism against viral pathogens by plants. It is known that specifically expressed endogenous siRNAs and miRNAs are involved in the self-defense process during viral infection. However, research has been rarely devoted to the endogenous siRNA and miRNA expression changes under viral infection if the resistance has already been genetically engineered in plants. Aiming to gain a deeper understanding of the RNA-mediated gene silencing defense process in plants, the expression profiles of siRNAs and miRNAs before and after viral infection in both wild type and transgenic anti-Rice stripe virus (RSV) rice plants were examined by small RNA high-throughput sequencing. Our research confirms that the newly generated siRNAs, which are derived from the engineered inverted repeat construct, is the major contributor of the viral resistance in rice. Further analysis suggests the accuracy of siRNA biogenesis might be affected when siRNAs machinery is excessively used in the transgenic plants. In addition, the expression levels of many known miRNAs are dramatically changed due to RSV infection on both wild type and transgenic rice plants, indicating potential function of those miRNAs involved in plant-virus interacting process.

Divergence of gene body DNA methylation and evolution of plant duplicate genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets, including new miRNA target-identification software

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.

GC content fluctuation around plant small RNA-generating sites

GC content of small RNA-generating sites and their flanking sequences in Arabidopsis thaliana and rice was systematically analyzed in silico. High GC content fluctuation (GCF) is observed at the borders of sRNA sites, while the GCF within sRNA sites is low. Furthermore, the GC content along sequences of some miniature inverted-repeat transposable element (MITE) families coincides with the abundance of MITE-derived small RNAs. The GCF within tasiRNA clusters is negatively correlated with its phasing score. We conclude that high GC content and low GCF could increase the expression of small RNA. Our results provide further insights on small RNA expression, which may be applied to improve the silencing efficiency of RNAi and virus-induced gene silencing.

Regulation of meiotic gene expression in plants

With the recent advances in genomics and sequencing technologies, databases of transcriptomes representing many cellular processes have been assembled. Meiotic transcriptomes in plants have been studied in Arabidopsis thaliana, rice (Oryza sativa), wheat (Triticum aestivum), petunia (Petunia hybrida), sunflower (Helianthus annuus), and maize (Zea mays). Studies in all organisms, but particularly in plants, indicate that a very large number of genes are expressed during meiosis, though relatively few of them seem to be required for the completion of meiosis. In this review, we focus on gene expression at the RNA level and analyze the meiotic transcriptome datasets and explore expression patterns of known meiotic genes to elucidate how gene expression could be regulated during meiosis. We also discuss mechanisms, such as chromatin organization and non-coding RNAs that might be involved in the regulation of meiotic transcription patterns.

Unresolved issues in pre-meiotic anther development

Compared to the diversity of other floral organs, the steps in anther ontogeny, final cell types, and overall organ shape are remarkably conserved among Angiosperms. Defects in pre-meiotic anthers that alter cellular composition or function typically result in male-sterility. Given the ease of identifying male-sterile mutants, dozens of genes with key roles in early anther development have been identified and cloned in model species, ordered by time of action and spatiotemporal expression, and used to propose explanatory models for critical steps in cell fate specification. Despite rapid progress, fundamental issues in anther development remain unresolved, and it is unclear if insights from one species can be applied to others. Here we construct a comparison of Arabidopsis, rice, and maize immature anthers to pinpoint distinctions in developmental pace. We analyze the mechanisms by which archesporial (pre-meiotic) cells are specified distinct from the soma, discuss what constitutes meiotic preparation, and review what is known about the secondary parietal layer and its terminal periclinal division that generates the tapetal and middle layers. Finally, roles for small RNAs are examined, focusing on the grass-specific phasiRNAs.