Heterochromatic siRNAs and DDM1 independently silence aberrant 5S rDNA transcripts in Arabidopsis

Atypical epigenetic and small RNA control of degenerated transposons and their fragments in clonally reproducing Spirodela polyrhiza

A handful of model plants have provided insight into silencing of transposable elements (TEs) through RNA-directed DNA methylation (RdDM). Guided by 24 nt long small-interfering RNAs (siRNAs), this epigenetic regulation installs DNA methylation and histone modifications like H3K9me2, which can be subsequently maintained independently of siRNAs. However, the genome of the clonally propagating duckweed Spirodela polyrhiza (Lemnaceae) has low levels of DNA methylation, very low expression of RdDM components, and near absence of 24 nt siRNAs. Moreover, some genes encoding RdDM factors, DNA methylation maintenance, and RNA silencing mechanisms are missing from the genome. Here, we investigated the distribution of TEs and their epigenetic marks in the Spirodela genome. Although abundant degenerated TEs have largely lost DNA methylation and H3K9me2 is low, they remain marked by the heterochromatin-associated H3K9me1 and H3K27me1 modifications. In contrast, we find high levels of DNA methylation and H3K9me2 in the relatively few intact TEs, which are source of 24 nt siRNAs, like RdDM-controlled TEs in other angiosperms. The data suggest that, potentially as adaptation to vegetative propagation, RdDM extent, silencing components, and targets are different from other angiosperms, preferentially focused on potentially intact TEs. It also provides evidence for heterochromatin maintenance independently of DNA methylation in flowering plants. These discoveries highlight the diversity of silencing mechanisms that exist in plants and the importance of using disparate model species to discover these mechanisms.

Functions of RNAi Pathways in Ribosomal RNA Regulation

Argonaute proteins, guided by small RNAs, play crucial roles in gene regulation and genome protection through RNA interference (RNAi)-related mechanisms. Ribosomal RNAs (rRNAs), encoded by repeated rDNA units, constitute the core of the ribosome being the most abundant cellular transcripts. rDNA clusters also serve as sources of small RNAs, which are loaded into Argonaute proteins and are able to regulate rDNA itself or affect other gene targets. In this review, we consider the impact of small RNA pathways, specifically siRNAs and piRNAs, on rRNA gene regulation. Data from diverse eukaryotic organisms suggest the potential involvement of small RNAs in various molecular processes related to the rDNA transcription and rRNA fate. Endogenous siRNAs are integral to the chromatin-based silencing of rDNA loci in plants and have been shown to repress rDNA transcription in animals. Small RNAs also play a role in maintaining the integrity of rDNA clusters and may function in the cellular response to rDNA damage. Studies on the impact of RNAi and small RNAs on rRNA provide vast opportunities for future exploration.

Advances in the Study of the Transcriptional Regulation Mechanism of Plant miRNAs

MicroRNAs (miRNA) are a class of endogenous, non-coding, small RNAs with about 22 nucleotides (nt), that are widespread in plants and are involved in various biological processes, such as development, flowering phase transition, hormone signal transduction, and stress response. The transcriptional regulation of miRNAs is an important process of miRNA gene regulation, and it is essential for miRNA biosynthesis and function. Like mRNAs, miRNAs are transcribed by RNA polymerase II, and these transcription processes are regulated by various transcription factors and other proteins. Consequently, the upstream genes regulating miRNA transcription, their specific expression, and the regulating mechanism were reviewed to provide more information for further research on the miRNA regulatory mechanism and help to further understand the regulatory networks of plant miRNAs.

Enzymatic reactions of AGO4 in RNA-directed DNA methylation: siRNA duplex loading, passenger strand elimination, target RNA slicing, and sliced target retention

RNA-directed DNA methylation in plants is guided by 24-nt siRNAs generated in parallel with 23-nt RNAs of unknown function. We show that 23-nt RNAs function as passenger strands during 24-nt siRNA incorporation into AGO4. The 23-nt RNAs are then sliced into 11- and 12-nt fragments, with 12-nt fragments remaining associated with AGO4. Slicing recapitulated with recombinant AGO4 and synthetic RNAs reveals that siRNAs of 21-24 nt, with any 5'-terminal nucleotide, can guide slicing, with sliced RNAs then retained by AGO4. In vivo, RdDM target locus RNAs that copurify with AGO4 also display a sequence signature of slicing. Comparing plants expressing slicing-competent versus slicing-defective AGO4 shows that slicing elevates cytosine methylation levels at virtually all RdDM loci. We propose that siRNA passenger strand elimination and AGO4 tethering to sliced target RNAs are distinct modes by which AGO4 slicing enhances RNA-directed DNA methylation.

DDM1-Mediated TE Silencing in Plants

Epigenetic modifications are indispensable for regulating gene bodies and TE silencing. DECREASE IN DNA METHYLATION 1 (DDM1) is a chromatin remodeller involved in histone modifications and DNA methylation. Apart from maintaining the epigenome, DDM1 also maintains key plant traits such as flowering time and heterosis. The role of DDM1 in epigenetic regulation is best characterised in plants, especially arabidopsis, rice, maize and tomato. The epigenetic changes induced by DDM1 establish the stable inheritance of many plant traits for at least eight generations, yet DDM1 does not methylate protein-coding genes. The DDM1 TE silencing mechanism is distinct and has evolved independently of other silencing pathways. Unlike the RNA-directed DNA Methylation (RdDM) pathway, DDM1 does not depend on siRNAs to enforce the heterochromatic state of TEs. Here, we review DDM1 TE silencing activity in the RdDM and non-RdDM contexts. The DDM1 TE silencing machinery is strongly associated with the histone linker H1 and histone H2A.W. While the linker histone H1 excludes the RdDM factors from methylating the heterochromatin, the histone H2A.W variant prevents TE mobility. The DDM1-H2A.W strategy alone silences nearly all the mobile TEs in the arabidopsis genome. Thus, the DDM1-directed TE silencing essentially preserves heterochromatic features and abolishes mobile threats to genome stability.

Thousands of high-quality sequencing samples fail to show meaningful correlation between 5S and 45S ribosomal DNA arrays in humans

The ribosomal RNA genes (rDNA) are tandemly arrayed in most eukaryotes and exhibit vast copy number variation. There is growing interest in integrating this variation into genotype-phenotype associations. Here, we explored a possible association of rDNA copy number variation with autism spectrum disorder and found no difference between probands and unaffected siblings. Because short-read sequencing estimates of rDNA copy number are error prone, we sought to validate our 45S estimates. Previous studies reported tightly correlated, concerted copy number variation between the 45S and 5S arrays, which should enable the validation of 45S copy number estimates with pulsed-field gel-verified 5S copy numbers. Here, we show that the previously reported strong concerted copy number variation may be an artifact of variable data quality in the earlier published 1000 Genomes Project sequences. We failed to detect a meaningful correlation between 45S and 5S copy numbers in thousands of samples from the high-coverage Simons Simplex Collection dataset as well as in the recent high-coverage 1000 Genomes Project sequences. Our findings illustrate the challenge of genotyping repetitive DNA regions accurately and call into question the accuracy of recently published studies of rDNA copy number variation in cancer that relied on diverse publicly available resources for sequence data.

Non-coding RNA polymerases that silence transposable elements and reprogram gene expression in plants

Multisubunit RNA polymerase (Pol) complexes are the core machinery for gene expression in eukaryotes. The enzymes Pol I, Pol II and Pol III transcribe distinct subsets of nuclear genes. This family of nuclear RNA polymerases expanded in terrestrial plants by the duplication of Pol II subunit genes. Two Pol II-related enzymes, Pol IV and Pol V, are highly specialized in the production of regulatory, non-coding RNAs. Pol IV and Pol V are the central players of RNA-directed DNA methylation (RdDM), an RNA interference pathway that represses transposable elements (TEs) and selected genes. Genetic and biochemical analyses of Pol IV/V subunits are now revealing how these enzymes evolved from ancestral Pol II to sustain non-coding RNA biogenesis in silent chromatin. Intriguingly, Pol IV-RdDM regulates genes that influence flowering time, reproductive development, stress responses and plant–pathogen interactions. Pol IV target genes vary among closely related taxa, indicating that these regulatory circuits are often species-specific. Data from crops like maize, rice, tomato and Brassicarapa suggest that dynamic repositioning of TEs, accompanied by Pol IV targeting to TE-proximal genes, leads to the reprogramming of plant gene expression over short evolutionary timescales.

Systematic Review of Plant Ribosome Heterogeneity and Specialization

Plants dedicate a high amount of energy and resources to the production of ribosomes. Historically, these multi-protein ribosome complexes have been considered static protein synthesis machines that are not subject to extensive regulation but only read mRNA and produce polypeptides accordingly. New and increasing evidence across various model organisms demonstrated the heterogeneous nature of ribosomes. This heterogeneity can constitute specialized ribosomes that regulate mRNA translation and control protein synthesis. A prominent example of ribosome heterogeneity is seen in the model plant, Arabidopsis thaliana, which, due to genome duplications, has multiple paralogs of each ribosomal protein (RP) gene. We support the notion of plant evolution directing high RP paralog divergence toward functional heterogeneity, underpinned in part by a vast resource of ribosome mutants that suggest specialization extends beyond the pleiotropic effects of single structural RPs or RP paralogs. Thus, Arabidopsis is a highly suitable model to study this phenomenon. Arabidopsis enables reverse genetics approaches that could provide evidence of ribosome specialization. In this review, we critically assess evidence of plant ribosome specialization and highlight steps along ribosome biogenesis in which heterogeneity may arise, filling the knowledge gaps in plant science by providing advanced insights from the human or yeast fields. We propose a data analysis pipeline that infers the heterogeneity of ribosome complexes and deviations from canonical structural compositions linked to stress events. This analysis pipeline can be extrapolated and enhanced by combination with other high-throughput methodologies, such as proteomics. Technologies, such as kinetic mass spectrometry and ribosome profiling, will be necessary to resolve the temporal and spatial aspects of translational regulation while the functional features of ribosomal subpopulations will become clear with the combination of reverse genetics and systems biology approaches.

3-Hydroxybutyrate Is Active Compound in Flax That Upregulates Genes Involved in DNA Methylation

In mammalian cells, 3-hydroxybutyrate (3-HB) is not only an intermediate metabolite during the oxidation of fatty acids, but also an important signaling molecule. On the other hand, the information about the metabolism or function of this compound in plants is scarce. In our study, we show for the first time that this compound naturally occurs in flax. The expression of bacterial β-ketothiolase in flax affects expression of endogenous genes of the 3-HB biosynthesis pathway and the compound content. The increase in 3-HB content in transgenic plants or after control plants treatment with 3-HB resulted in upregulation of genes involved in chromatin remodeling. The observation that 3-HB is an endogenous activator of methyltransferase 3 (CMT3), decreased DNA methylation I (DDM1), DEMETER DNA glycosylase (DME), and an inhibitor of sirtuin 1 (SRT1) provides an example of integration of different genes in chromatin remodeling. The changes in chromatin remodeling gene expression concomitant with those involved in phenolics and the lignin biosynthesis pathway suggest potential integration of secondary metabolic status with epigenetic changes.

Genetic and epigenetic variation in 5S ribosomal RNA genes reveals genome dynamics in Arabidopsis thaliana

Organized in tandem repeat arrays in most eukaryotes and transcribed by RNA polymerase III, expression of 5S rRNA genes is under epigenetic control. To unveil mechanisms of transcriptional regulation, we obtained here in depth sequence information on 5S rRNA genes from the Arabidopsis thaliana genome and identified differential enrichment in epigenetic marks between the three 5S rDNA loci situated on chromosomes 3, 4 and 5. We reveal the chromosome 5 locus as the major source of an atypical, long 5S rRNA transcript characteristic of an open chromatin structure. 5S rRNA genes from this locus translocated in the Landsberg erecta ecotype as shown by linkage mapping and chromosome-specific FISH analysis. These variations in 5S rDNA locus organization cause changes in the spatial arrangement of chromosomes in the nucleus. Furthermore, 5S rRNA gene arrangements are highly dynamic with alterations in chromosomal positions through translocations in certain mutants of the RNA-directed DNA methylation pathway and important copy number variations among ecotypes. Finally, variations in 5S rRNA gene sequence, chromatin organization and transcripts indicate differential usage of 5S rDNA loci in distinct ecotypes. We suggest that both the usage of existing and new 5S rDNA loci resulting from translocations may impact neighboring chromatin organization.

Analysis of siRNA Precursors Generated by RNA Polymerase IV and RNA-Dependent RNA Polymerase 2 in Arabidopsis

Noncoding RNAs perform diverse regulatory functions in living cells. In plants, two RNA polymerase II-related enzymes, RNA polymerases IV and V (Pol IV and V), specialize in the synthesis of noncoding RNAs that silence a subset of transposable elements and genes via RNA-directed DNA methylation (RdDM). In this process, Pol IV partners with RNA-dependent RNA polymerase 2 (RDR2) to produce double-stranded RNAs that are then cut by an RNase III enzyme, Dicer-like 3 (DCL3), into 24 nt small interfering RNAs (siRNAs). The siRNAs are loaded into an Argonaute family protein, primarily AGO4, and guide the complex to complementary DNA target sequences where RdDM and repressive chromatin modifications ensue. The dependence of 24 nt siRNA biogenesis on Pol IV and RDR2 has been known for more than a decade, but the elusive pre-siRNA transcripts synthesized by Pol IV and RDR2 have only recently been identified. This chapter describes the approaches that enabled our identification of Pol IV/RDR2-dependent RNAs (P4R2 RNAs) in Arabidopsis thaliana. These included the use of a triple Dicer mutant (dcl2 dcl3 dcl4) to cause P4R2 RNAs to accumulate, genome-wide identification and mapping of P4R2 RNAs using a modified Illumina small RNA-Seq protocol, and multiple bioinformatic pipelines for data analysis and displaying results.

The PEAT protein complexes are required for histone deacetylation and heterochromatin silencing

In eukaryotes, heterochromatin regions are typically subjected to transcriptional silencing. DNA methylation has an important role in such silencing and has been studied extensively. However, little is known about how methylated heterochromatin regions are subjected to silencing. We conducted a genetic screen and identified an epcr (enhancer of polycomb-related) mutant that releases heterochromatin silencing in Arabidopsis thaliana We demonstrated that EPCR1 functions redundantly with its paralog EPCR2 and interacts with PWWP domain-containing proteins (PWWPs), AT-rich interaction domain-containing proteins (ARIDs), and telomere repeat binding proteins (TRBs), thus forming multiple functionally redundant protein complexes named PEAT (PWWPs-EPCRs-ARIDs-TRBs). The PEAT complexes mediate histone deacetylation and heterochromatin condensation and thereby facilitate heterochromatin silencing. In heterochromatin regions, the production of small interfering RNAs (siRNAs) and DNA methylation is repressed by the PEAT complexes. The study reveals how histone deacetylation, heterochromatin condensation, siRNA production, and DNA methylation interplay with each other and thereby maintain heterochromatin silencing.

Inhibition of RNA polymerase II allows controlled mobilisation of retrotransposons for plant breeding

Background

Retrotransposons play a central role in plant evolution and could be a powerful endogenous source of genetic and epigenetic variability for crop breeding. To ensure genome integrity several silencing mechanisms have evolved to repress retrotransposon mobility. Even though retrotransposons fully depend on transcriptional activity of the host RNA polymerase II (Pol II) for their mobility, it was so far unclear whether Pol II is directly involved in repressing their activity.

Results

Here we show that plants defective in Pol II activity lose DNA methylation at repeat sequences and produce more extrachromosomal retrotransposon DNA upon stress in Arabidopsis and rice. We demonstrate that combined inhibition of both DNA methylation and Pol II activity leads to a strong stress-dependent mobilization of the heat responsive ONSEN retrotransposon in Arabidopsis seedlings. The progenies of these treated plants contain up to 75 new ONSEN insertions in their genome which are stably inherited over three generations of selfing. Repeated application of heat stress in progeny plants containing increased numbers of ONSEN copies does not result in increased activation of this transposon compared to control lines. Progenies with additional ONSEN copies show a broad panel of environment-dependent phenotypic diversity.

Conclusions

We demonstrate that Pol II acts at the root of transposon silencing. This is important because it suggests that Pol II can regulate the speed of plant evolution by fine-tuning the amplitude of transposon mobility. Our findings show that it is now possible to study induced transposon bursts in plants and unlock their use to induce epigenetic and genetic diversity for crop breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1265-4) contains supplementary material, which is available to authorized users.

Mobile small RNAs and their role in regulating cytosine methylation of DNA

Small (s)RNAs of 21 to 24 nucleotides are associated with RNA silencing and methylation of DNA cytosine residues. All sizes can move from cell-to-cell and long distance in plants, directing RNA silencing in destination cells. Twenty-four nucleotide sRNAs are the predominant long-distance mobile species. Thousands move from shoot to root, where they target methylation of transposable elements both directly and indirectly. We derive several classes of interaction between small RNAs and methylation and use these to explore the mechanisms of methylation and gene expression that associate with mobile sRNA signaling.

Identification of Pol IV and RDR2-dependent precursors of 24 nt siRNAs guiding de novo DNA methylation in Arabidopsis

In Arabidopsis thaliana, abundant 24 nucleotide small interfering RNAs (24 nt siRNA) guide the cytosine methylation and silencing of transposons and a subset of genes. 24 nt siRNA biogenesis requires nuclear RNA polymerase IV (Pol IV), RNA-dependent RNA polymerase 2 (RDR2) and DICER-like 3 (DCL3). However, siRNA precursors are mostly undefined. We identified Pol IV and RDR2-dependent RNAs (P4R2 RNAs) that accumulate in dcl3 mutants and are diced into 24 nt RNAs by DCL3 in vitro. P4R2 RNAs are mostly 26-45 nt and initiate with a purine adjacent to a pyrimidine, characteristics shared by Pol IV transcripts generated in vitro. RDR2 terminal transferase activity, also demonstrated in vitro, may account for occasional non-templated nucleotides at P4R2 RNA 3' termini. The 24 nt siRNAs primarily correspond to the 5' or 3' ends of P4R2 RNAs, suggesting a model whereby siRNAs are generated from either end of P4R2 duplexes by single dicing events.

Non-coding RNA transcription and RNA-directed DNA methylation in Arabidopsis

RNA-directed DNA methylation (RdDM) is responsible for transcriptional silencing of endogenous transposable elements and introduced transgenes. This process requires non-coding RNAs produced by DNA-dependent RNA polymerases IV and V (Pol IV and Pol V). Pol IV-produced non-coding RNAs are precursors of 24-nt small interfering RNAs, whereas Pol V-produced ncRNAs directly act as scaffold RNAs. In this review, we summarize recent advances in the understanding of RdDM. In particular, we focus on the mechanisms underlying the recruitment of Pol IV and Pol V to chromatin and the targeting of RdDM.

Cytosine hypomethylation at CHG and CHH sites in the pleiotropic mutants of Mendelian inheritance in Catharanthus roseus

The 5S and 18S rDNA sequences of Catharanthus roseus cv 'Nirmal' (wild type) and its leafless inflorescence (lli), evergreen dwarf (egd) and irregular leaf lamina (ill) single mutants and lli egd, lli ill and egd ill double mutants were characterized. The lli, egd and ill mutants of Mendelian inheritance bore the names after their most conspicuous morphological feature(s). They had been chemically induced and isolated for their salt tolerance. The double mutants were isolated as morphological segregants from crosses between single mutants. The morphological features of the two parents accompanied salt tolerance in the double mutants. All the six mutants were hypomethylated at repeat sequences, upregulated and downregulated for many genes and carried pleiotropic alterations for several traits. Here the 5S and 18S rDNAs of C. roseus were found to be relatively low in cytosine content. Cytosines were preponderantly in CG context (53%) and almost all of them were methylated (97%). The cytosines in CHH and CHG (where H = A, T or C) contexts were largely demethylated (92%) in mutants. The demethylation was attributable to reduced expression of RDR2 and DRM2 led RNA dependant DNA methylation and CMT3 led maintenance methylation pathways. Mutants had gained some cytosines by substitution of C at T sites. These perhaps arose on account of errors in DNA replication, mediated by widespread cytosine demethylation at CHG and CHH sites. It was concluded that the regulation of cytosine ethylation mechanisms was disturbed in the mutants. ILL, EGD and LLI genes were identified as the positive regulators of other genes mediating the RdDM and CMT3 pathways, for establishment and maintenance of cytosine methylation in C. roseus.

Functional consequences of subunit diversity in RNA polymerases II and V

Multisubunit RNA polymerases IV and V (Pol IV and Pol V) evolved as specialized forms of Pol II that mediate RNA-directed DNA methylation (RdDM) and transcriptional silencing of transposons, viruses, and endogenous repeats in plants. Among the subunits common to Arabidopsis thaliana Pols II, IV, and V are 93% identical alternative ninth subunits, NRP(B/D/E)9a and NRP(B/D/E)9b. The 9a and 9b subunit variants are incompletely redundant with respect to Pol II; whereas double mutants are embryo lethal, single mutants are viable, yet phenotypically distinct. Likewise, 9a or 9b can associate with Pols IV or V but RNA-directed DNA methylation is impaired only in 9b mutants. Based on genetic and molecular tests, we attribute the defect in RdDM to impaired Pol V function. Collectively, our results reveal a role for the ninth subunit in RNA silencing and demonstrate that subunit diversity generates functionally distinct subtypes of RNA polymerases II and V.

Poplars with a PtDDM1-RNAi transgene have reduced DNA methylation and show aberrant post-dormancy morphology

The Arabidopsis thaliana DDM1 (Decreased DNA Methylation) gene is necessary for the maintenance of DNA methylation and heterochromatin assembly. In Arabidopsis, ddm1 mutants exhibit strong but delayed morphological phenotypes. We used RNA interference (RNAi) to suppress transcripts of two orthologous DDM1 paralogs in Populus trichocarpa and examined effects on whole plant phenotypes during perennial growth and seasonal dormancy. The RNAi-PtDDM1 transgenic poplars showed a wide range of DDM1 transcript suppression; the most strongly suppressed line had 37.5 % of the expression of the non-transgenic control. Genomic cytosine methylation (mC %) was 11.1 % in the non-transgenic control, compared with 9.1 % for the transgenic event with lowest mC %, a reduction of 18.1 %. An evaluation of greenhouse growth directly after acclimation of in vitro grown plants showed no developmental or growth rate abnormalities associated with the decrease in PtDDM1 expression. However, after a dormancy cycle and growth outdoors, a mottled leaf phenotype appeared in some of the transgenic insertion events that had strongly reduced PtDDM1 expression and DNA methylation. The phenotypic consequences of reduced DDM1 activity and DNA methylation appears to increase with cumulative plant propagation and growth.

The Arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin

Nucleosome remodelers of the DDM1/Lsh family are required for DNA methylation of transposable elements, but the reason for this is unknown. How DDM1 interacts with other methylation pathways, such as small-RNA-directed DNA methylation (RdDM), which is thought to mediate plant asymmetric methylation through DRM enzymes, is also unclear. Here, we show that most asymmetric methylation is facilitated by DDM1 and mediated by the methyltransferase CMT2 separately from RdDM. We find that heterochromatic sequences preferentially require DDM1 for DNA methylation and that this preference depends on linker histone H1. RdDM is instead inhibited by heterochromatin and absolutely requires the nucleosome remodeler DRD1. Together, DDM1 and RdDM mediate nearly all transposon methylation and collaborate to repress transposition and regulate the methylation and expression of genes. Our results indicate that DDM1 provides DNA methyltransferases access to H1-containing heterochromatin to allow stable silencing of transposable elements in cooperation with the RdDM pathway.

Heterochromatin dynamics during developmental transitions in Arabidopsis - a focus on ribosomal DNA loci

The Arabidopsis chromosomes contain conspicuous heterochromatin domains comprising the repetitive 45S and 5S ribosomal DNA loci as well as centromeric and pericentromeric repeats that organize into chromocenters during interphase. During developmental phase transitions such as seed maturation, germination, seedling growth and flowering that require large-scale reprogramming of gene expression patterns, the organization of repetitive sequences into chromocenters dynamically changes. Here we illustrate recent studies that shed light on the heterochromatin dynamics in cotyledons, the first aerial tissues preformed in the embryo, and in true leaves. We will summarize available data for the 5S rDNA repeat loci, in particular their chromatin organization and expression dynamics during the first days of post-germination development, and discuss how the plant accommodates 5S rRNA transcription during large-scale chromatin reorganization events.

Asymmetric purine-pyrimidine distribution in cellular small RNA population of papaya

Background

The small RNAs (sRNA) are a regulatory class of RNA mainly represented by the 21 and 24-nucleotide size classes. The cellular sRNAs are processed by RNase III family enzyme dicer (Dicer like in plant) from a self-complementary hairpin loop or other type of RNA duplexes. The papaya genome has been sequenced, but its microRNAs and other regulatory RNAs are yet to be analyzed.

Results

We analyzed the genomic features of the papaya sRNA population from three sRNA deep sequencing libraries made from leaves, flowers, and leaves infected with Papaya Ringspot Virus (PRSV). We also used the deep sequencing data to annotate the micro RNA (miRNA) in papaya. We identified 60 miRNAs, 24 of which were conserved in other species, and 36 of which were novel miRNAs specific to papaya. In contrast to the Chargaff’s purine-pyrimidine equilibrium, cellular sRNA was significantly biased towards a purine rich population. Of the two purine bases, higher frequency of adenine was present in 23nt or longer sRNAs, while 22nt or shorter sRNAs were over represented by guanine bases. However, this bias was not observed in the annotated miRNAs in plants. The 21nt species were expressed from fewer loci but expressed at higher levels relative to the 24nt species. The highly expressed 21nt species were clustered in a few isolated locations of the genome. The PRSV infected leaves showed higher accumulation of 21 and 22nt sRNA compared to uninfected leaves. We observed higher accumulation of miRNA* of seven annotated miRNAs in virus-infected tissue, indicating the potential function of miRNA* under stressed conditions.

Conclusions

We have identified 60 miRNAs in papaya. Our study revealed the asymmetric purine-pyrimidine distribution in cellular sRNA population. The 21nt species of sRNAs have higher expression levels than 24nt sRNA. The miRNA* of some miRNAs shows higher accumulation in PRSV infected tissues, suggesting that these strands are not totally functionally redundant. The findings open a new avenue for further investigation of the sRNA silencing pathway in plants.

Expression of 5 S rRNA genes linked to 35 S rDNA in plants, their epigenetic modification and regulatory element divergence

BACKGROUND

In plants, the 5 S rRNA genes usually occur as separate tandems (S-type arrangement) or, less commonly, linked to 35 S rDNA units (L-type). The activity of linked genes remains unknown so far. We studied the homogeneity and expression of 5 S genes in several species from family Asteraceae known to contain linked 35 S-5 S units. Additionally, their methylation status was determined using bisulfite sequencing. Fluorescence in situ hybridization was applied to reveal the sub-nuclear positions of rDNA arrays.

RESULTS

We found that homogenization of L-type units went to completion in most (4/6) but not all species. Two species contained major L-type and minor S-type units (termed L(s)-type). The linked genes dominate 5 S rDNA expression while the separate tandems do not seem to be expressed. Members of tribe Anthemideae evolved functional variants of the polymerase III promoter in which a residing C-box element differs from the canonical angiosperm motif by as much as 30%. On this basis, a more relaxed consensus sequence of a plant C-box: (5'-RGSWTGGGTG-3') is proposed. The 5 S paralogs display heavy DNA methylation similarly as to their unlinked counterparts. FISH revealed the close association of 35 S-5 S arrays with nucleolar periphery indicating that transcription of 5 S genes may occur in this territory.

CONCLUSIONS

We show that the unusual linked arrangement of 5 S genes, occurring in several plant species, is fully compatible with their expression and functionality. This extraordinary 5 S gene dynamics is manifested at different levels, such as variation in intrachromosomal positions, unit structure, epigenetic modification and considerable divergence of regulatory motifs.

RNAi-independent de novo DNA methylation revealed in Arabidopsis mutants of chromatin remodeling gene DDM1

Methylation of histone H3 lysine 9 (H3K9me) and small RNAs are associated with constitutively silent chromatin in diverse eukaryotes including plants. In plants, silent transposons are also marked by cytosine methylation, especially at non-CpG sites. Transposon-specific non-CpG methylation in plants is controlled by small RNAs and H3K9me. Although it is often assumed that small RNA directs H3K9me, interaction between small RNA and H3K9me has not been directly demonstrated in plants. We have previously shown that a mutation in the chromatin remodeling gene DDM1 (DECREASE IN DNA METHYLATION 1) induces a global decrease but a local increase of cytosine methylation and accumulation of small RNA at a locus called BONSAI. Here we show that de novo BONSAI methylation does not depend on RNAi but does depend on H3K9me. In mutants of H3K9 methyltransferase gene KRYPTONITE or the H3K9me-dependent DNA methyltransferase gene CHROMOMETHYALSE3, the ddm1-induced de novo cytosine methylation was abolished for all three contexts (CpG, CpHpG and CpHpH). Furthermore, RNAi mutants showed strong developmental defects when combined with the ddm1 mutation. Our results revealed unexpected interactions of epigenetic modifications that may be conserved among diverse eukaryotes.

Deficiency in DNA methylation increases meiotic crossover rates in euchromatic but not in heterochromatic regions in Arabidopsis

Meiotic recombination is tightly regulated by cis- and trans-acting factors. Although DNA methylation and chromatin remodeling affect chromosome structure, their impact on meiotic recombination is not well understood. To study the effect of DNA methylation on the landscape of chromosomal recombination, we analyzed meiotic recombination in the decreased DNA methylation 1 (ddm1) mutant. DDM1 is a SWI2/SNF2-like chromatin-remodeling protein necessary for DNA methylation and heterochromatin maintenance in Arabidopsis thaliana. The rate of meiotic recombination between markers located in euchromatic regions was significantly higher in both heterozygous (DDM1/ddm1) and homozygous (ddm1/ddm1) backgrounds than in WT plants. The effect on recombination was similar for both male and female meiocytes. Contrary to expectations, ddm1 had no effect on the number of crossovers between markers in heterochromatic pericentric regions that underwent demethylation. These results are surprising, because the pericentromeric regions are hypermethylated and were expected to be the regions most affected by demethylation. Thus, DDM1 loss of function may trigger changes that enhance meiotic recombination in euchromatin regions but are not sufficient to induce the same events in heterochromatic segments. This work uncovers the repressive role of methylation on meiotic recombination in euchromatic regions and suggests that additional factors may have a role in controlling the suppression of recombination in heterochromatin.

Regulation of Pol I-transcribed 45S rDNA and Pol III-transcribed 5S rDNA in Arabidopsis

The 18S, 5.8S and 25S rRNAs, which result from the 45S precursor, together with 5S rRNAs, are central components of the ribosome. The integration of one molecule of each rRNA per ribosome necessitates an elaborate coordination between transcriptions of the two ribosomal DNA (rDNA) families. Even though 5S rDNA is transcribed by RNA polymerase III and 45S rDNA by RNA polymerase I, the two rDNA families present certain similarities in their transcriptional regulation. This review aims to compare 5S and 45S rRNA genes in the plant model Arabidopsis thaliana in terms of organization, transcription and regulation, and draws parallels between the two rDNA families.

Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing

BACKGROUND

Melon (Cucumis melo L.) is a commercially important fruit crop that is cultivated worldwide. The melon research community has recently benefited from the determination of a complete draft genome sequence and the development of associated genomic tools, which have allowed us to focus on small RNAs (sRNAs). These are short, non-coding RNAs 21-24 nucleotides in length with diverse physiological roles. In plants, they regulate gene expression and heterochromatin assembly, and control protection against virus infection. Much remains to be learned about the role of sRNAs in melon.

RESULTS

We constructed 10 sRNA libraries from two stages of developing ovaries, fruits and photosynthetic cotyledons infected with viruses, and carried out high-throughput pyrosequencing. We catalogued and analysed the melon sRNAs, resulting in the identification of 26 known miRNA families (many conserved with other species), the prediction of 84 melon-specific miRNA candidates, the identification of trans-acting siRNAs, and the identification of chloroplast, mitochondrion and transposon-derived sRNAs. In silico analysis revealed more than 400 potential targets for the conserved and novel miRNAs.

CONCLUSION

We have discovered and analysed a large number of conserved and melon-specific sRNAs, including miRNAs and their potential target genes. This provides insight into the composition and function of the melon small RNAome, and paves the way towards an understanding of sRNA-mediated processes that regulate melon fruit development and melon-virus interactions.

Small RNA-mediated epigenetic modifications in plants

Epigenetic modifications in plants can be directed and mediated by small RNAs (sRNAs). This regulation is composed of a highly interactive network of sRNA-directed DNA methylation, histone, and chromatin modifications, all of which control transcription. Identification and functional characterization of components of the siRNA-directed DNA methylation pathway have provided insights into epigenetic pathways that form heterochromatin and into chromatin-based pathways for gene silencing, paramutation, genetic imprinting, and epigenetic reprogramming. Next-generation sequencing technologies have facilitated new discoveries and have helped create a basic blueprint of the plant epigenome. As the multiple layers of epigenetic regulation in plants are dissected, a more comprehensive understanding of the biological importance of epigenetic marks and states has been developed.

Incomplete homogenization of 18 S ribosomal DNA coding regions in Arabidopsis thaliana

Background

As a result of concerted evolution, coding regions of ribosomal DNA sequences are highly conserved within species and variation is generally thought to be limited to a few nucleotides. However, rDNA sequence variation has not been systematically examined in plant genomes, including that of the model plant Arabidopsis thaliana whose genome was the first to be sequenced.

Findings

Both genomic and transcribed 18 S sequences were sampled and revealed that most deviation from the consensus sequence was limited to single nucleotide substitutions except for a variant with a 270 bp deletion from position 456 to 725 in Arabidopsis numbering. The deletion maps to the functionally important and highly conserved 530 loop or helix18 in the structure of E. coli 16 S. The expression of the deletion variant is tightly controlled during developmental growth stages. Transcripts were not detectable in young seedlings but could be amplified from RNA extracts of mature leaves, stems, flowers and roots of Arabidopsis thaliana ecotype Columbia. We also show polymorphism for the deletion variant among four Arabidopsis ecotypes examined.

Conclusion

Despite a strong purifying selection that might be expected against functionally impaired rDNAs, the newly identified variant is maintained in the Arabidopsis genome. The expression of the variant and the polymorphism displayed by Arabidopsis ecotypes suggest a transition state in concerted evolution.

Regulation and function of DNA methylation in plants and animals

DNA methylation is an important epigenetic mark involved in diverse biological processes. In plants, DNA methylation can be established through the RNA-directed DNA methylation pathway, an RNA interference pathway for transcriptional gene silencing (TGS), which requires 24-nt small interfering RNAs. In mammals, de novo DNA methylation occurs primarily at two developmental stages: during early embryogenesis and during gametogenesis. While it is not clear whether establishment of DNA methylation patterns in mammals involves RNA interference in general, de novo DNA methylation and suppression of transposons in germ cells require 24-32-nt piwi-interacting small RNAs. DNA methylation status is dynamically regulated by DNA methylation and demethylation reactions. In plants, active DNA demethylation relies on the repressor of silencing 1 family of bifunctional DNA glycosylases, which remove the 5-methylcytosine base and then cleave the DNA backbone at the abasic site, initiating a base excision repair (BER) pathway. In animals, multiple mechanisms of active DNA demethylation have been proposed, including a deaminase- and DNA glycosylase-initiated BER pathway. New information concerning the effects of various histone modifications on the establishment and maintenance of DNA methylation has broadened our understanding of the regulation of DNA methylation. The function of DNA methylation in plants and animals is also discussed in this review.

Differential sensitivity of Arabidopsis siRNA biogenesis mutants to genotoxic stress

Plant response to stress has been linked to different RNA-silencing processes and epigenetic mechanisms. Our recent results showed that Arabidopsis thaliana Dicer-like (DCL) mutants were impaired in transgenerational changes, recombination frequency and stress tolerance. We also found that transgenerational changes were dependent on changes in DNA methylation. Here, we hypothesized that plants deficient in the production of small RNAs would show an impaired abiotic stress response. To test this, we exposed A. thaliana dcl2, dcl3, dcl4, dcl2 dcl3 (d2d3), dcl2 dcl4 (d2d4), dcl2 dcl3 dcl4 (d2d3d4), nrpd1a, rdr2 and rdr6 mutants to methyl methane sulfonate (MMS). We found dcl4 and rdr6 to be more sensitive and dcl2, dcl3, d2d3 and rdr2 plants more resistant to MMS, as shown by fresh weight, root length and survival rate. The in vitro repair assay showed the lower ability of dcl2 and dcl3 to repair UV-damaged DNA. To summarize, we found that whereas mutants impaired in transactivating siRNA biogenesis were more sensitive to MMS, mutants impaired in natural antisense siRNA and heterochromatic siRNA biogeneses were more tolerant. Our data suggest that plant response to MMS is in part regulated through biogenesis of various siRNAs.

Overexpression of a natural chloroplast-encoded antisense RNA in tobacco destabilizes 5S rRNA and retards plant growth

BACKGROUND

The roles of non-coding RNAs in regulating gene expression have been extensively studied in both prokaryotes and eukaryotes, however few reports exist as to their roles in organellar gene regulation. Evidence for accumulation of natural antisense RNAs (asRNAs) in chloroplasts comes from the expressed sequence tag database and cDNA libraries, while functional data have been largely obtained from artificial asRNAs. In this study, we used Nicotiana tabacum to investigate the effect on sense strand transcripts of overexpressing a natural chloroplast asRNA, AS5, which is complementary to the region which encodes the 5S rRNA and tRNAArg.

RESULTS

AS5-overexpressing (AS5ox) plants obtained by chloroplast transformation exhibited slower growth and slightly pale green leaves. Analysis of AS5 transcripts revealed four distinct species in wild-type (WT) and AS5ox plants, and additional AS5ox-specific products. Of the corresponding sense strand transcripts, tRNAArg overaccumulated several-fold in transgenic plants whereas 5S rRNA was unaffected. However, run-on transcription showed that the 5S-trnR region was transcribed four-fold more in the AS5ox plants compared to WT, indicating that overexpression of AS5 was associated with decreased stability of 5S rRNA. In addition, polysome analysis of the transformants showed less 5S rRNA and rbcL mRNA associated with ribosomes.

CONCLUSIONS

Our results suggest that AS5 can modulate 5S rRNA levels, giving it the potential to affect Chloroplast translation and plant growth. More globally, overexpression of asRNAs via chloroplast transformation may be a useful strategy for defining their functions.

Genome organization and function: a view from yeast and Arabidopsis

Recent progress in understanding higher-order chromatin organization in the nucleus has been considerable. From single gene to chromosome territory, realistic biophysical models can now accurately predict some of the structural feature of cell nuclei. Despite growing evidence of a deterministic nuclear organization, the physiological consequence of spatial genome organization is still unclear. In the simple eukaryotic model, Saccharomyces cerevisiae, clear correlation between gene position and transcription has been established. In this review, we will focus on higher-order chromatin organization in yeast with respect to the nuclear envelope and nucleolus. In Arabidopsis thaliana, a model plant for which we have a complete genome sequence, chromosome territory (CT) arrangement and somatic homologous pairing in interphase nuclei seem to occur randomly. Since chromosomes containing nucleolar organizer regions associate more frequently to form a single nucleolar structure, as in yeast, the nucleolus seems to play a major role in organizing nuclear space. Recent findings have begun to elucidate how plant regulatory factors, such as chromatin remodeling or histone chaperones, affect the chromatin state of ribosomal DNA genes located in two distinct CT arrangements in the nucleus. The functional outcome of yeast nuclear organization allowed us to propose how nuclear organization might contribute to a novel type of epigenetic regulation: the spatial regulation of transcription.

Repeat elements and the Arabidopsis DNA methylation landscape

DNA methylation is an epigenetic mark that has key roles in the control of genome activity in plants and mammals. It is critical for the stable silencing of repeat elements and is also involved in the epigenetic regulation of some genes. Despite similarities in the controlling functions of DNA methylation, its dynamics and deposition patterns differ in several respects between plants and mammals. One of the most striking differences is that plants tend to propagate pre-existing DNA methylation states across generations, whereas mammals re-establish them genome wide at every generation. Here, we review our current understanding of DNA methylation in the flowering plant Arabidopsis. We discuss in particular the role of RNAi in the incremental methylation and silencing of repeat elements over successive generations. We argue that paramutation, an epigenetic phenomenon first described in maize, is an extreme manifestation of this RNAi-dependent pathway.

Deep sequencing of viroid-derived small RNAs from grapevine provides new insights on the role of RNA silencing in plant-viroid interaction

BACKGROUND

Viroids are circular, highly structured, non-protein-coding RNAs that, usurping cellular enzymes and escaping host defense mechanisms, are able to replicate and move through infected plants. Similarly to viruses, viroid infections are associated with the accumulation of viroid-derived 21-24 nt small RNAs (vd-sRNAs) with the typical features of the small interfering RNAs characteristic of RNA silencing, a sequence-specific mechanism involved in defense against invading nucleic acids and in regulation of gene expression in most eukaryotic organisms.

METHODOLOGY/PRINCIPAL FINDINGS

To gain further insights on the genesis and possible role of vd-sRNAs in plant-viroid interaction, sRNAs isolated from Vitis vinifera infected by Hop stunt viroid (HSVd) and Grapevine yellow speckle viroid 1 (GYSVd1) were sequenced by the high-throughput platform Solexa-Illumina, and the vd-sRNAs were analyzed. The large majority of HSVd- and GYSVd1-sRNAs derived from a few specific regions (hotspots) of the genomic (+) and (-) viroid RNAs, with a prevalence of those from the (-) strands of both viroids. When grouped according to their sizes, vd-sRNAs always assumed a distribution with prominent 21-, 22- and 24-nt peaks, which, interestingly, mapped at the same hotspots.

CONCLUSIONS/SIGNIFICANCE

These findings show that different Dicer-like enzymes (DCLs) target viroid RNAs, preferentially accessing to the same viroid domains. Interestingly, our results also suggest that viroid RNAs may interact with host enzymes involved in the RNA-directed DNA methylation pathway, indicating more complex scenarios than previously thought for both vd-sRNAs genesis and possible interference with host gene expression.

A Pol V-mediated silencing, independent of RNA-directed DNA methylation, applies to 5S rDNA

The plant-specific RNA polymerases Pol IV and Pol V are essential to RNA–directed DNA methylation (RdDM), which also requires activities from RDR2 (RNA–Dependent RNA Polymerase 2), DCL3 (Dicer-Like 3), AGO4 (Argonaute), and DRM2 (Domains Rearranged Methyltransferase 2). RdDM is dedicated to the methylation of target sequences which include transposable elements, regulatory regions of several protein-coding genes, and 5S rRNA–encoding DNA (rDNA) arrays. In this paper, we have studied the expression of the 5S-210 transcript, a marker of silencing release at 5S RNA genes, to show a differential impact of RNA polymerases IV and V on 5S rDNA arrays during early development of the plant. Using a combination of molecular and cytological assays, we show that Pol IV, RDR2, DRM2, and Pol V, actors of the RdDM, are required to maintain a transcriptional silencing of 5S RNA genes at chromosomes 4 and 5. Moreover, we have shown a derepression associated to chromatin decondensation specific to the 5S array from chromosome 4 and restricted to the Pol V–loss of function. In conclusion, our results highlight a new role for Pol V on 5S rDNA, which is RdDM–independent and comes specifically at chromosome 4, in addition to the RdDM pathway.

In plant genomes, the RNA–directed DNA methylation (RdDM) process induces de novo methylation of cytosines at repeated sequences. The RNA polymerases Pol IV and Pol V are two key components of the RdDM pathway. Pol IV acts with RDR2 (RNA–dependent RNA polymerase 2) and DCL3 (Dicer-Like protein 3) to generate short interfering RNAs (siRNAs). Pol V, in a partnership including AGO4 (Argonaute4) and DRM2 (Domains Rearranged Methyltransferase 2), drives DNA methylation at the targeted sequence. Changes in 5S (ribosomal DNA) rDNA methylation, 5S rDNA chromatin compaction, and 5S siRNA accumulation in Pol IV/V mutants have been reported. However, 5S rDNA arrays were considered together. In the present study, we observed an overexpression of the atypic 5S-210 transcript, restricted to the 5S rDNA array from chomosome 4. This derepression is specific to the Pol V–loss of function (and not to Pol IV) and comes in addition to the RdDM pathway. The Pol V–loss of function induces also the chromatin decondensation of the derepressed 5S locus at chomosome 4. Our results highlight a new role for Pol V which, suprisingly, appears to be Pol IV– and RdDM–independent.