Search

RNA interference pathways can involve amplification of secondary siRNAs by RNA-dependent RNA polymerases. In plants, RDR6-dependent secondary siRNAs arise from transcripts targeted by some microRNAs (miRNAs). Here, Arabidopsis thaliana secondary siRNAs from mRNA as well as trans-acting siRNAs are shown to be triggered through initial targeting by a 22-nucleotide (nt) miRNA that associates with AGO1. In contrast to canonical 21-nt miRNAs, 22-nt miRNAs primarily arise from foldback precursors containing asymmetric bulges. Using artificial miRNA constructs, conversion of asymmetric foldbacks to symmetric foldbacks resulted in the production of 21-nt forms of miR173, miR472 and miR828. Both 21- and 22-nt forms associated with AGO1 and guided accurate slicer activity, but only 22-nt forms were competent to trigger RDR6-dependent siRNA production from target RNA. These data suggest that AGO1 functions differentially with 21- and 22-nt miRNAs to engage the RDR6-associated amplification apparatus.

MicroRNAs (miRNAs) are short regulatory RNAs processed from partially self-complementary foldbacks within longer MIRNA primary transcripts. Several MIRNA families are conserved deeply through land plants, but many are present only in closely related species or are species specific. The finding of numerous evolutionarily young MIRNA, many with low expression and few if any targets, supports a rapid birth-death model for MIRNA evolution. A systematic analysis of MIRNA genes and families in the close relatives, Arabidopsis thaliana and Arabidopsis lyrata, was conducted using both whole-genome comparisons and high-throughput sequencing of small RNAs. Orthologs of 143 A. thaliana MIRNA genes were identified in A. lyrata, with nine having significant sequence or processing changes that likely alter function. In addition, at least 13% of MIRNA genes in each species are unique, despite their relatively recent speciation (approximately 10 million years ago). Alignment of MIRNA foldbacks to the Arabidopsis genomes revealed evidence for recent origins of 32 families by inverted or direct duplication of mostly protein-coding gene sequences, but less than half of these yield miRNA that are predicted to target transcripts from the originating gene family. miRNA nucleotide divergence between A. lyrata and A. thaliana orthologs was higher for young MIRNA genes, consistent with reduced purifying selection compared with deeply conserved MIRNA genes. Additionally, target sites of younger miRNA were lost more frequently than for deeply conserved families. In summary, our systematic analyses emphasize the dynamic nature of the MIRNA complement of plant genomes.

Plants respond to virus infections by activation of RNA-based silencing, which limits infection at both the single-cell and system levels. Viruses encode RNA silencing suppressor proteins that interfere with this response. Wild-type Arabidopsis thaliana is immune to silencing suppressor (HC-Pro)-deficient Turnip mosaic virus, but immunity was lost in the absence of DICER-LIKE proteins DCL4 and DCL2. Systematic analysis of susceptibility and small RNA formation in Arabidopsis mutants lacking combinations of RNA-dependent RNA polymerase (RDR) and DCL proteins revealed that the vast majority of virus-derived small interfering RNAs (siRNAs) were dependent on DCL4 and RDR1, although full antiviral defense also required DCL2 and RDR6. Among the DCLs, DCL4 was sufficient for antiviral silencing in inoculated leaves, but DCL2 and DCL4 were both involved in silencing in systemic tissues (inflorescences). Basal levels of antiviral RNA silencing and siRNA biogenesis were detected in mutants lacking RDR1, RDR2, and RDR6, indicating an alternate route to form double-stranded RNA that does not depend on the three previously characterized RDR proteins.

Transacting siRNA (tasiRNA) biogenesis in Arabidopsis is initiated by microRNA (miRNA) -guided cleavage of primary transcripts. In the case of TAS3 tasiRNA formation, ARGONAUTE7 (AGO7)-miR390 complexes interact with primary transcripts at two sites, resulting in recruitment of RNA-DEPENDENT RNA POLYMERASE6 for dsRNA biosynthesis. An extensive screen for Arabidopsis mutants with specific defects in TAS3 tasiRNA biogenesis or function was done. This yielded numerous ago7 mutants, one dcl4 mutant, and two mutants that accumulated low levels of miR390. A direct genome sequencing-based approach to both map and rapidly identify one of the latter mutant alleles was developed. This revealed a G-to-A point mutation (mir390a-1) that was calculated to stabilize a relatively nonpaired region near the base of the MIR390a foldback, resulting in misprocessing of the miR390/miR390* duplex and subsequent reduced TAS3 tasiRNA levels. Directed substitutions, as well as analysis of variation at paralogous miR390-generating loci (MIR390a and MIR390b), indicated that base pair properties and nucleotide identity within a region 4-6 bases below the miR390/miR390* duplex region contributed to the efficiency and accuracy of precursor processing.

The advent of high-throughput sequencing (HTS) methods has enabled direct approaches to quantitatively profile small RNA populations. However, these methods have been limited by several factors, including representational artifacts and lack of established statistical methods of analysis. Furthermore, massive HTS data sets present new problems related to data processing and mapping to a reference genome. Here, we show that cluster-based sequencing-by-synthesis technology is highly reproducible as a quantitative profiling tool for several classes of small RNA from Arabidopsis thaliana. We introduce the use of synthetic RNA oligoribonucleotide standards to facilitate objective normalization between HTS data sets, and adapt microarray-type methods for statistical analysis of multiple samples. These methods were tested successfully using mutants with small RNA biogenesis (miRNA-defective dcl1 mutant and siRNA-defective dcl2 dcl3 dcl4 triple mutant) or effector protein (ago1 mutant) deficiencies. Computational methods were also developed to rapidly and accurately parse, quantify, and map small RNA data.

Correct daily phasing of transcription confers an adaptive advantage to almost all organisms, including higher plants. In this study, we describe a hypothesis-driven network discovery pipeline that identifies biologically relevant patterns in genome-scale data. To demonstrate its utility, we analyzed a comprehensive matrix of time courses interrogating the nuclear transcriptome of Arabidopsis thaliana plants grown under different thermocycles, photocycles, and circadian conditions. We show that 89% of Arabidopsis transcripts cycle in at least one condition and that most genes have peak expression at a particular time of day, which shifts depending on the environment. Thermocycles alone can drive at least half of all transcripts critical for synchronizing internal processes such as cell cycle and protein synthesis. We identified at least three distinct transcription modules controlling phase-specific expression, including a new midnight specific module, PBX/TBX/SBX. We validated the network discovery pipeline, as well as the midnight specific module, by demonstrating that the PBX element was sufficient to drive diurnal and circadian condition-dependent expression. Moreover, we show that the three transcription modules are conserved across Arabidopsis, poplar, and rice. These results confirm the complex interplay between thermocycles, photocycles, and the circadian clock on the daily transcription program, and provide a comprehensive view of the conserved genomic targets for a transcriptional network key to successful adaptation., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

Development of the Arabidopsis Small RNA Project (ASRP) Database, which provides information and tools for the analysis of microRNA, endogenous siRNA and other small RNA-related features, has been driven by the introduction of high-throughput sequencing technology. To accommodate the demands of increased data, numerous improvements and updates have been made to ASRP, including new ways to access data, more efficient algorithms for handling data, and increased integration with community-wide resources. New search and visualization tools have also been developed to improve access to small RNA classes and their targets. ASRP is publicly available through a web interface at http://asrp.cgrb.oregonstate.edu/db/.

Posttranscriptional RNA silencing of many endogenous transcripts, viruses, and transgenes involves the RNA-DEPENDENT RNA POLYMERASE6/DICER-LIKE4 (RDR6/DCL4)-dependent short interfering RNA (siRNA) biogenesis pathway. Arabidopsis thaliana contains several families of trans-acting siRNAs (tasiRNAs) that form in 21-nucleotide phased arrays through the RDR6/DCL4-dependent pathway and that negatively regulate target transcripts. Using deep sequencing technology and computational approaches, the phasing patterns of known tasiRNAs and tasiRNA-like loci from across the Arabidopsis genome were analyzed in wild-type plants and silencing-defective mutants. Several gene transcripts were found to be routed through the RDR6/DCL4-dependent pathway after initial targeting by one or multiple miRNAs or tasiRNAs, the most conspicuous example of which was an expanding clade of genes encoding pentatricopeptide repeat (PPR) proteins. Interestingly, phylogenetic analysis using Populus trichocarpa revealed evidence for small RNA-mediated regulatory mechanisms within a similarly expanded group of PPR genes. We suggest that posttranscriptional silencing mechanisms operate on an evolutionary scale to buffer the effects of rapidly expanding gene families.

Eukaryotes contain a diversified set of small RNA-guided pathways that control genes, repeated sequences, and viruses at the transcriptional and posttranscriptional levels. Genome-wide profiles and analyses of small RNAs, particularly the large class of 24-nucleotide (nt) short interfering RNAs (siRNAs), were done for wild-type Arabidopsis thaliana and silencing pathway mutants with defects in three RNA-dependent RNA polymerase (RDR) and four Dicer-like (DCL) genes. The profiling involved direct analysis using a multiplexed, parallel-sequencing strategy. Small RNA-generating loci, especially those producing predominantly 24-nt siRNAs, were found to be highly correlated with repetitive elements across the genome. These were found to be largely RDR2- and DCL3-dependent, although alternative DCL activities were detected on a widespread level in the absence of DCL3. In contrast, no evidence for RDR2-alternative activities was detected. Analysis of RDR2- and DCL3-dependent small RNA accumulation patterns in and around protein-coding genes revealed that upstream gene regulatory sequences systematically lack siRNA-generating activities. Further, expression profiling suggested that relatively few genes, proximal to abundant 24-nt siRNAs, are regulated directly by RDR2- and DCL3-dependent silencing. We conclude that the widespread accumulation patterns for RDR2- and DCL3-dependent siRNAs throughout the Arabidopsis genome largely reflect mechanisms to silence highly repeated sequences.

In plants, microRNAs (miRNAs) comprise one of two classes of small RNAs that function primarily as negative regulators at the posttranscriptional level. Several MIRNA genes in the plant kingdom are ancient, with conservation extending between angiosperms and the mosses, whereas many others are more recently evolved. Here, we use deep sequencing and computational methods to identify, profile and analyze non-conserved MIRNA genes in Arabidopsis thaliana. 48 non-conserved MIRNA families, nearly all of which were represented by single genes, were identified. Sequence similarity analyses of miRNA precursor foldback arms revealed evidence for recent evolutionary origin of 16 MIRNA loci through inverted duplication events from protein-coding gene sequences. Interestingly, these recently evolved MIRNA genes have taken distinct paths. Whereas some non-conserved miRNAs interact with and regulate target transcripts from gene families that donated parental sequences, others have drifted to the point of non-interaction with parental gene family transcripts. Some young MIRNA loci clearly originated from one gene family but form miRNAs that target transcripts in another family. We suggest that MIRNA genes are undergoing relatively frequent birth and death, with only a subset being stabilized by integration into regulatory networks.