Elisabeth J. Chapman

High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes.

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.

P1/HC-Pro, a Viral Suppressor of RNA Silencing, Interferes with Arabidopsis Development and miRNA Function

The molecular basis for virus-induced disease in plants has been a long-standing mystery. Infection of Arabidopsis by Turnip mosaic virus (TuMV) induces a number of developmental defects in vegetative and reproductive organs. We found that these defects, many of which resemble those in miRNA-deficient dicer-like1 (dcl1) mutants, were due to the TuMV-encoded RNA-silencing suppressor, P1/HC-Pro. Suppression of RNA silencing is a counterdefensive mechanism that enables systemic infection by TuMV. The suppressor interfered with the activity of miR171 (also known as miRNA39), which directs cleavage of several mRNAs coding for Scarecrow-like transcription factors, by inhibiting miR171-guided nucleolytic function. Out of ten other mRNAs that were validated as miRNA-guided cleavage targets, eight accumulated to elevated levels in the presence of P1/HC-Pro. The basis for TuMV- and other virus-induced disease in plants may be explained, at least partly, by interference with miRNA-controlled developmental pathways that share components with the antiviral RNA-silencing pathway.

Genome-Wide Profiling and Analysis of Arabidopsis siRNAs

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.

Small RNA binding is a common strategy to suppress RNA silencing by several viral suppressors.

RNA silencing is an evolutionarily conserved system that functions as an antiviral mechanism in higher plants and insects. To counteract RNA silencing, viruses express silencing suppressors that interfere with both siRNA‐ and microRNA‐guided silencing pathways. We used comparative in vitro and in vivo approaches to analyse the molecular mechanism of suppression by three well‐studied silencing suppressors. We found that silencing suppressors p19, p21 and HC‐Pro each inhibit the intermediate step of RNA silencing via binding to siRNAs, although the molecular features required for duplex siRNA binding differ among the three proteins. None of the suppressors affected the activity of preassembled RISC complexes. In contrast, each suppressor uniformly inhibited the siRNA‐initiated RISC assembly pathway by preventing RNA silencing initiator complex formation.

Genome-Wide Analysis of the RNA-DEPENDENT RNA POLYMERASE6/DICER-LIKE4 Pathway in Arabidopsis Reveals Dependency on miRNA- and tasiRNA-Directed Targeting

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.

Arabidopsis RNA-Dependent RNA Polymerases and Dicer-Like Proteins in Antiviral Defense and Small Interfering RNA Biogenesis during Turnip Mosaic Virus Infection

Viruses promote infection by suppressing antiviral RNA silencing responses. Using suppressor-deficient and suppressor-competent Turnip mosaic virus, a genetic analysis of antiviral RNA silencing showed that Arabidopsis Dicers and RNA-dependent RNA polymerases have modular roles in antiviral defense. 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.

Small RNAs serve as a genetic buffer against genomic shock in Arabidopsis interspecific hybrids and allopolyploids.

Small RNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), and trans-acting siRNAs (tasiRNAs), control gene expression and epigenetic regulation. Although the roles of miRNAs and siRNAs have been extensively studied, their expression diversity and evolution in closely related species and interspecific hybrids are poorly understood. Here, we show comprehensive analyses of miRNA expression and siRNA distributions in two closely related species Arabidopsis thaliana and Arabidopsis arenosa, a natural allotetraploid Arabidopsis suecica, and two resynthesized allotetraploid lines (F1 and F7) derived from A. thaliana and A. arenosa. We found that repeat- and transposon-associated siRNAs were highly divergent between A. thaliana and A. arenosa. A. thaliana siRNA populations underwent rapid changes in F1 but were stably maintained in F7 and A. suecica. The correlation between siRNAs and nonadditive gene expression in allopolyploids is insignificant. In contrast, miRNA and tasiRNA sequences were conserved between species, but their expression patterns were highly variable between the allotetraploids and their progenitors. Many miRNAs tested were nonadditively expressed (deviating from the mid-parent value, MPV) in the allotetraploids and triggered unequal degradation of A. thaliana or A. arenosa targets. The data suggest that small RNAs produced during interspecific hybridization or polyploidization serve as a buffer against the genomic shock in interspecific hybrids and allopolyploids: Stable inheritance of repeat-associated siRNAs maintains chromatin and genome stability, whereas expression variation of miRNAs leads to changes in gene expression, growth vigor, and adaptation.

MicroRNA Gene Evolution in Arabidopsis lyrata and Arabidopsis thaliana

A whole-genome analysis of MIRNA from Arabidopsis thaliana and close relative Arabidopsis lyrata suggests that evolutionarily young MIRNA are diverging in sequence and function more rapidly than are more deeply conserved MIRNA. These and other results shed light on the birth, divergence, and death of MIRNA genes in plants. 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 (∼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.

Transgenically expressed viral RNA silencing suppressors interfere with microRNA methylation in Arabidopsis.

HEN1‐dependent methylation of the 3′‐terminal nucleotide is a crucial step in plant microRNA (miRNA) biogenesis. Here we report that several viral RNA silencing suppressors (P1/HC‐Pro, p21 and p19) inhibit miRNA methylation. These suppressors have distinct effects on different miRNAs. We also show that miRNA∗ is methylated in vivo in a suppressor‐sensitive manner, suggesting that the viral proteins interfere with miRNA/miRNA∗ duplexes. p19 and p21 bind both methylated and unmethylated miRNA/miRNA∗ duplexes in vivo. These findings suggest miRNA/miRNA∗ as the in vivo substrates for the HEN1 miRNA methyltransferase and raise intriguing possibilities regarding the cellular location of miRNA methylation.

Computational and analytical framework for small RNA profiling by high-throughput sequencing

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.