Prior to the advent of new, deep sequencing methods, small RNA (sRNA) discovery was dependent on Sanger sequencing, which was time-consuming and limited knowledge to only the most abundant sRNA. The innovation of large-scale, next-generation sequencing has exponentially increased knowledge of the biology, diversity and abundance of sRNA populations. In this review, we discuss issues involved in the design of sRNA sequencing experiments, including choosing a sequencing platform, inherent biases that affect sRNA measurements and replication. We outline the steps involved in preprocessing sRNA sequencing data and review both the principles behind and the current options for normalization. Finally, we discuss differential expression analysis in the absence and presence of biological replicates. While our focus is on sRNA sequencing experiments, many of the principles discussed are applicable to the sequencing of other RNA populations.

RNA-directed DNA methylation (RdDM) in plants requires two RNA polymerase (Pol) II-related RNA polymerases, namely Pol IV and Pol V. A genetic screen designed to reveal factors that are important for RdDM in a developmental context in Arabidopsis identified DEFECTIVE IN MERISTEM SILENCING 4 (DMS4). Unlike other mutants defective in RdDM, dms4 mutants have a pleiotropic developmental phenotype. The DMS4 protein is similar to yeast IWR1 (interacts with RNA polymerase II), a conserved putative transcription factor that interacts with Pol II subunits. The DMS4 complementary DNA partly complements the K1 killer toxin hypersensitivity of a yeast iwr1 mutant, suggesting some functional conservation. In the transgenic system studied, mutations in DMS4 directly or indirectly affect Pol IV-dependent secondary short interfering RNAs, Pol V-mediated RdDM, Pol V-dependent synthesis of intergenic non-coding RNA and expression of many Pol II-driven genes. These data suggest that DMS4 might be a regulatory factor for several RNA polymerases, thus explaining its diverse roles in the plant.

miRNAs are the most abundant class of small non-coding RNAs, and they are involved in post-transcriptional regulations, playing a crucial role in the refinement of genetic programming during plant development. Here we present a comprehensive picture of miRNA regulation in Vitis vinifera L. plant during its complete life cycle. Furthering our knowledge about the post-transcriptional regulation of plant development is fundamental to understand the biology of such an important crop.

3' uridylation is increasingly recognized as a conserved RNA modification process associated with RNA turnover in eukaryotes. 2'-O-methylation on the 3' terminal ribose protects micro(mi)RNAs from 3' truncation and 3' uridylation in Arabidopsis. Previously, we identified HESO1 as the nucleotidyl transferase that uridylates most unmethylated miRNAs in vivo, but substantial 3' tailing of miRNAs still remains in heso1 loss-of-function mutants. In this study, we found that among nine other potential nucleotidyl transferases, UTP:RNA uridylyltransferase 1 (URT1) is the single most predominant nucleotidyl transferase that tails miRNAs. URT1 and HESO1 prefer substrates with different 3' end nucleotides in vitro and act cooperatively to tail different forms of the same miRNAs in vivo. Moreover, both HESO1 and URT1 exhibit nucleotidyl transferase activity on AGO1-bound miRNAs. Although these enzymes are able to add long tails to AGO1-bound miRNAs, the tailed miRNAs remain associated with AGO1. Moreover, tailing of AGO1-bound miRNA165/6 drastically reduces the slicing activity of AGO1-miR165/6, suggesting that tailing reduces miRNA activity. However, monouridylation of miR171a by URT1 endows the miRNA the ability to trigger the biogenesis of secondary siRNAs. Therefore, 3' tailing could affect the activities of miRNAs in addition to leading to miRNA degradation.

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

Antibiotic use in humans has been associated with outgrowth of fungi. Here we used a murine model to investigate the gut microbiome over 76 days of treatment with vancomycin, ampicillin, neomycin, and metronidazole and subsequent recovery. Mouse stool was studied as a surrogate for the microbiota of the lower gastrointestinal tract. The abundance of fungi and bacteria was measured using quantitative PCR, and the proportional composition of the communities quantified using 454/Roche pyrosequencing of rRNA gene tags. Prior to treatment, bacteria outnumbered fungi by >3 orders of magnitude. Upon antibiotic treatment, bacteria dropped in abundance >3 orders of magnitude, so that the predominant 16S sequences detected became transients derived from food. Upon cessation of treatment, bacterial communities mostly returned to their previous numbers and types after 8 weeks, though communities remained detectably different from untreated controls. Fungal communities varied substantially over time, even in the untreated controls. Separate cages within the same treatment group showed radical differences, but mice within a cage generally behaved similarly. Fungi increased ∼40-fold in abundance upon antibiotic treatment but declined back to their original abundance after cessation of treatment. At the last time point, Candida remained more abundant than prior to treatment. These data show that 1) gut fungal populations change radically during normal mouse husbandry, 2) fungi grow out in the gut upon suppression of bacterial communities with antibiotics, and 3) perturbations due to antibiotics persist long term in both the fungal and bacterial microbiota.

MPSS (massively parallel signature sequencing) is a sequencing-based technology that uses a unique method to quantify gene expression level, generating millions of short sequence tags per library. We have created a series of databases for four species (Arabidopsis, rice, grape and Magnaporthe grisea, the rice blast fungus). Our MPSS databases measure the expression level of most genes under defined conditions and provide information about potentially novel transcripts (antisense transcripts, alternative splice isoforms and regulatory intergenic transcripts). A modified version of MPSS has been used to perform deep profiling of small RNAs from Arabidopsis, and we have recently adapted our database to display these data. Interpretation of the small RNA MPSS data is facilitated by the inclusion of extensive repeat data in our genome viewer. All the data and the tools introduced in this article are available at http://mpss.udel.edu.

We have further characterized floral organ-localized gene expression in the inflorescence of Arabidopsis thaliana by comparison of massively parallel signature sequencing (MPSS) data. Six libraries of RNA sequence tags from immature inflorescence tissues were constructed and matched to their respective loci in the annotated Arabidopsis genome. These signature libraries survey the floral transcriptome of wild-type tissue as well as the floral homeotic mutants, apetala1, apetala3, agamous, a superman/apetala1 double mutant, and differentiated ovules dissected from the gynoecia of wild-type inflorescences. Comparing and contrasting these MPSS floral expression libraries enabled demarcation of transcripts enriched in the petals, stamens, stigma-style, gynoecia, and those with predicted enrichment within the sepal/sepal-petals, petal-stamens, or gynoecia-stamens.

Progress in epigenetics has revealed mechanisms that can heritably regulate gene function independent of genetic alterations. Nevertheless, little is known about the role of epigenetics in evolution. This is due in part to scant data on epigenetic variation among natural populations. In plants, small interfering RNA (siRNA) is involved in both the initiation and maintenance of gene silencing by directing DNA methylation and/or histone methylation. Here, we report that, in the model plant Arabidopsis thaliana, a cluster of approximately 24 nt siRNAs found at high levels in the ecotype Landsberg erecta (Ler) could direct DNA methylation and heterochromatinization at a hAT element adjacent to the promoter of FLOWERING LOCUS C (FLC), a major repressor of flowering, whereas the same hAT element in ecotype Columbia (Col) with almost identical DNA sequence, generates a set of low abundance siRNAs that do not direct these activities. We have called this hAT element MPF for Methylated region near Promoter of FLC, although de novo methylation triggered by an inverted repeat transgene at this region in Col does not alter its FLC expression. DNA methylation of the Ler allele MPF is dependent on genes in known silencing pathways, and such methylation is transmissible to Col by genetic crosses, although with varying degrees of penetrance. A genome-wide comparison of Ler and Col small RNAs identified at least 68 loci matched by a significant level of approximately 24 nt siRNAs present specifically in Ler but not Col, where nearly half of the loci are related to repeat or TE sequences. Methylation analysis revealed that 88% of the examined loci (37 out of 42) were specifically methylated in Ler but not Col, suggesting that small RNA can direct epigenetic differences between two closely related Arabidopsis ecotypes.

Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), play a variety of important regulatory roles in many eukaryotes. Their small size has made it challenging to study them directly in live cells. Here we describe an RNA-based fluorescent sensor for small RNA detection both in vitro and in vivo, adaptable for any small RNA. It utilizes an sxRNA switch for detection of miRNA-mRNA interactions combined with a fluorophore-binding sequence 'Spinach', a GFP-like RNA aptamer for which the RNA-fluorophore complex exhibits strong and consistent fluorescence under an excitation wavelength. Two example sensors, FASTmiR171 and FASTmiR122, can rapidly detect and quantify the levels of miR171 and miR122 in vitro. The sensors can determine relative levels of miRNAs in total RNA extracts with sensitivity similar to small RNA sequencing and northern blots. FASTmiR sensors were also used to estimate the copy number range of miRNAs in total RNA extracts. To localize and analyze the spatial distribution of small RNAs in live, single cells, tandem copies of FASTmiR122 were expressed in different cell lines. FASTmiR122 was able to quantitatively detect the differences in miR122 levels in Huh7 and HEK293T cells demonstrating its potential for tracking miRNA expression and localization in vivo.

MicroRNAs (miRNAs) control expression of endogenous target genes through transcript cleavage or translational inhibition. Legume plants can form a specialized organ, the nodule, through interaction with nitrogen fixing soil bacteria. To understand the regulatory roles of miRNAs in the nodulation process, we functionally validated gma-miR171o and gma-miR171q and their target genes in soybean. These two miRNA sequences are identical in sequence but their miRNA genes are divergent and show unique, tissue-specific expression patterns. The expression levels of the two miRNAs are negatively correlated with that of their target genes. Ectopic expression of these miRNAs in transgenic hairy roots resulted in a significant reduction in nodule formation. Both gma-miR171o and gma-miR171q target members of the GRAS transcription factor superfamily, namely GmSCL-6 and GmNSP2. Transient interaction of miRNAs and their target genes in tobacco cells further confirmed their cleavage activity. The results suggest that gma-miR171o and gma-miR171q regulate GmSCL-6 and GmNSP2, which in turn, influence expression of the Nodule inception (NIN), Early Nodulin 40 (ENOD40), and Ethylene Response Factor Required for Nodulation (ERN) genes during the Bradyrhizobium japonicum-soybean nodulation process. Collectively, our data suggest a role for two miRNAs, gma-miR171o and gma-miR171q, in regulating the spatial and temporal aspects of soybean nodulation.

Sex chromosomes evolved from autosomes many times across the eukaryote phylogeny. Several models have been proposed to explain this transition, some involving male and female sterility mutations linked in a region of suppressed recombination between X and Y (or Z/W, U/V) chromosomes. Comparative and experimental analysis of a reference genome assembly for a double haploid YY male garden asparagus (Asparagus officinalis L.) individual implicates separate but linked genes as responsible for sex determination. Dioecy has evolved recently within Asparagus and sex chromosomes are cytogenetically identical with the Y, harboring a megabase segment that is missing from the X. We show that deletion of this entire region results in a male-to-female conversion, whereas loss of a single suppressor of female development drives male-to-hermaphrodite conversion. A single copy anther-specific gene with a male sterile Arabidopsis knockout phenotype is also in the Y-specific region, supporting a two-gene model for sex chromosome evolution.

Plastid transformation technology has been widely used to express traits of potential commercial importance, though the technology has been limited to traits that function while sequestered in the organelle. Prior research indicates that plastid contents can escape from the organelle, suggesting a possible mechanism for engineering plastid transgenes to function in other cellular locations. To test this hypothesis, we created tobacco (Nicotiana tabacum cv. Petit Havana) plastid transformants that express a fragment of the nuclear-encoded Phytoene desaturase (PDS) gene capable of catalyzing post-transcriptional gene silencing if RNA escapes into the cytoplasm. We found multiple lines of direct evidence that plastid-encoded PDS transgenes affect nuclear PDS gene silencing: knockdown of the nuclear-encoded PDS mRNA and/or its apparent translational inhibition, biogenesis of 21-nucleotide (nt) phased small interfering RNAs (phasiRNAs), and pigment-deficient plants. Furthermore, plastid-expressed dsRNA with no cognate nuclear-encoded pairing partner also produced abundant 21-nt phasiRNAs in the cytoplasm, demonstrating that a nuclear-encoded template is not required for siRNA biogenesis. Our results indicate that RNA escape from plastids to the cytoplasm occurs generally, with functional consequences that include entry into the gene silencing pathway. Furthermore, we uncover a method to produce plastid-encoded traits with functions outside of the organelle and open additional fields of study in plastid development, compartmentalization, and small RNA biogenesis.

In Arabidopsis thaliana, ARGONAUTE1 (AGO1) plays a central role in microRNA (miRNA) and small interfering RNA (siRNA)-mediated silencing. AGO1 associates to the rough endoplasmic reticulum to conduct miRNA-mediated translational repression, mRNA cleavage, and biogenesis of phased siRNAs. Here, we show that a 37°C heat stress (HS) promotes AGO1 protein accumulation in cytosolic condensates where it colocalizes with components of siRNA bodies and of stress granules. AGO1 contains a prion-like domain in its poorly characterized N-terminal Poly-Q domain, which is sufficient to undergo phase separation independently of the presence of SGS3. HS only moderately affects the small RNA repertoire, the loading of AGO1 by miRNAs, and the signatures of target cleavage, suggesting that its localization in condensates protects AGO1 rather than promoting or impairing its activity in reprogramming gene expression during stress. Collectively, our work sheds new light on the impact of high temperature on a main effector of RNA silencing in plants.

RNA-directed DNA methylation (RdDM) is a small interfering RNA (siRNA)-mediated epigenetic modification that contributes to transposon silencing in plants. RdDM requires a complex transcriptional machinery that includes specialized RNA polymerases, named Pol IV and Pol V, as well as chromatin remodelling proteins, transcription factors, RNA binding proteins, and other plant-specific proteins whose functions are not yet clarified. In Arabidopsis thaliana, DICER-LIKE3 and members of the ARGONAUTE4 group of ARGONAUTE (AGO) proteins are involved, respectively, in generating and using 24-nt siRNAs that trigger methylation and transcriptional gene silencing of homologous promoter sequences. AGO4 is the main AGO protein implicated in the RdDM pathway. Here we report the identification of the related AGO6 in a forward genetic screen for mutants defective in RdDM and transcriptional gene silencing in shoot and root apical meristems in Arabidopsis thaliana. The identification of AGO6, and not AGO4, in our screen is consistent with the primary expression of AGO6 in shoot and root growing points.

MicroRNAs are ubiquitous in plant genomes but vary greatly in their abundance within and conservation among plant lineages. To gain insight into the evolutionary birth/death dynamics of microRNA families, we sequenced small RNA and 5'-end PARE libraries generated from two closely related species of Gossypium. Here, we demonstrate that 33 microRNA families, with similar copy numbers and average evolutionary rates, are conserved in the two congeneric cottons. Analysis of the presence/absence of these microRNA families in other land plants sheds light on their depth of phylogenetic origin and lineage-specific loss/gain. Conserved microRNA families in Gossypium exhibit a striking interspecific asymmetry in expression, potentially connected to relative proximity to neighboring transposable elements. A complex correlated expression pattern of microRNA target genes with their controlling microRNAs indicates that possible functional divergence of conserved microRNA families can also exist even within a single plant genus.

In plants, post-transcriptional gene silencing (PTGS) is mediated by DICER-LIKE 1 (DCL1)-dependent microRNAs (miRNAs), which also trigger 21-nucleotide secondary short interfering RNAs (siRNAs) via RNA-DEPENDENT RNA POLYMERASE 6 (RDR6), DCL4 and ARGONAUTE 1 (AGO1), whereas transcriptional gene silencing (TGS) of transposons is mediated by 24-nucleotide heterochromatic (het)siRNAs, RDR2, DCL3 and AGO4 (ref. 4). Transposons can also give rise to abundant 21-nucleotide 'epigenetically activated' small interfering RNAs (easiRNAs) in DECREASED DNA METHYLATION 1 (ddm1) and DNA METHYLTRANSFERASE 1 (met1) mutants, as well as in the vegetative nucleus of pollen grains and in dedifferentiated plant cell cultures. Here we show that easiRNAs in Arabidopsis thaliana resemble secondary siRNAs, in that thousands of transposon transcripts are specifically targeted by more than 50 miRNAs for cleavage and processing by RDR6. Loss of RDR6, DCL4 or DCL1 in a ddm1 background results in loss of 21-nucleotide easiRNAs and severe infertility, but 24-nucleotide hetsiRNAs are partially restored, supporting an antagonistic relationship between PTGS and TGS. Thus miRNA-directed easiRNA biogenesis is a latent mechanism that specifically targets transposon transcripts, but only when they are epigenetically reactivated during reprogramming of the germ line. This ancient recognition mechanism may have been retained both by transposons to evade long-term heterochromatic silencing and by their hosts for genome defence.

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.

Vegetative phase change is regulated by a decrease in the abundance of the miRNAs, miR156 and miR157, and the resulting increase in the expression of their targets, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors. To determine how miR156/miR157 specify the quantitative and qualitative changes in leaf morphology that occur during vegetative phase change, we measured their abundance in successive leaves and characterized the phenotype of mutations in different MIR156 and MIR157 genes. miR156/miR157 decline rapidly between leaf 1&2 and leaf 3 and decrease more slowly after this point. The amount of miR156/miR157 in leaves 1&2 greatly exceeds the threshold required to specify their identity. Subsequent leaves have relatively low levels of miR156/miR157 and are sensitive to small changes in their abundance. In these later-formed leaves, the amount of miR156/miR157 is close to the threshold required to specify juvenile vs. adult identity; a relatively small decrease in the abundance of miR156/157 in these leaves produces a disproportionately large increase in SPL proteins and a significant change in leaf morphology. miR157 is more abundant than miR156 but has a smaller effect on shoot morphology and SPL gene expression than miR156. This may be attributable to the inefficiency with which miR157 is loaded onto AGO1, as well as to the presence of an extra nucleotide at the 5' end of miR157 that is mis-paired in the miR157:SPL13 duplex. miR156 represses different targets by different mechanisms: it regulates SPL9 by a combination of transcript cleavage and translational repression and regulates SPL13 primarily by translational repression. Our results offer a molecular explanation for the changes in leaf morphology that occur during shoot development in Arabidopsis and provide new insights into the mechanism by which miR156 and miR157 regulate gene expression.

Miscanthus x giganteus (Mxg) is a perennial grass that produces superior biomass yields in temperate environments. The essentially uncharacterized triploid genome (3n = 57, x = 19) of Mxg is likely critical for the rapid growth of this vegetatively propagated interspecific hybrid.