Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Extensive natural Agrobacterium-induced transformation in the genus Camellia.

MAIN CONCLUSION: The genus Camellia underwent extensive natural transformation by Agrobacterium. Over a period of 15 million years, at least 12 different inserts accumulated in 72 investigated Camellia species. Like a wide variety of other wild and cultivated plants, Camellia species carry cellular T-DNA sequences (cT-DNAs) in their nuclear genomes, resulting from natural Agrobacterium-mediated transformation. Short and long DNA sequencing reads of 435 accessions belonging to 72 Camellia species (representing 12 out of 14 sections) were investigated for the occurrence of cT-DNA insertions. In all, 12 different cT-DNAs were recovered, either completely or partially, called CaTA to CaTL. Divergence analysis of internal cT-DNA repeats revealed that the insertion events span a period from 0.075 to 15 Mio years ago, and yielded an average transformation frequency of one event per 1.25 Mio years. The two oldest inserts, CaTA and CaTD, have been modified by spontaneous deletions and inversions, and by insertion of various plant sequences. In those cases where enough accessions were available (C. japonica, C. oleifera, C. chekiangoleosa, C. sasanqua and C. pitardii), the younger cT-DNA inserts showed a patchy distribution among different accessions of each species, indicating that they are not genetically fixed. It could be shown that Camellia breeding has led to intersectional transfer of cT-DNAs. Altogether, the cT-DNAs cover 374 kb, and carry 47 open reading frames (ORFs). Two Camellia cT-DNA genes, CaTH-orf358 and CaTK-orf8, represent new types of T-DNA genes. With its large number of cT-DNA sequences, the genus Camellia constitutes an interesting model for the study of natural Agrobacterium transformants.

CLSY docking to Pol IV requires a conserved domain critical for small RNA biogenesis and transposon silencing.

Eukaryotes must balance the need for gene transcription by RNA polymerase II (Pol II) against the danger of mutations caused by transposable element (TE) proliferation. In plants, these gene expression and TE silencing activities are divided between different RNA polymerases. Specifically, RNA polymerase IV (Pol IV), which evolved from Pol II, transcribes TEs to generate small interfering RNAs (siRNAs) that guide DNA methylation and block TE transcription by Pol II. While the Pol IV complex is recruited to TEs via SNF2-like CLASSY (CLSY) proteins, how Pol IV partners with the CLSYs remains unknown. Here we identified a conserved CYC-YPMF motif that is specific to Pol IV and is positioned on the complex exterior. Furthermore, we found that this motif is essential for the co-purification of all four CLSYs with Pol IV, but that only one CLSY is present in any given Pol IV complex. These findings support a "one CLSY per Pol IV" model where the CYC-YPMF motif acts as a CLSY-docking site. Indeed, mutations in and around this motif phenocopy pol iv null mutants. Together, these findings provide structural and functional insights into a critical protein feature that distinguishes Pol IV from other RNA polymerases, allowing it to promote genome stability by targeting TEs for silencing.

Long Noncoding RNAs in Plants.

Plants have an extraordinary diversity of transcription machineries, including five nuclear DNA-dependent RNA polymerases. Four of these enzymes are dedicated to the production of long noncoding RNAs (lncRNAs), which are ribonucleic acids with functions independent of their protein-coding potential. lncRNAs display a broad range of lengths and structures, but they are distinct from the small RNA guides of RNA interference (RNAi) pathways. lncRNAs frequently serve as structural, catalytic, or regulatory molecules for gene expression. They can affect all elements of genes, including promoters, untranslated regions, exons, introns, and terminators, controlling gene expression at various levels, including modifying chromatin accessibility, transcription, splicing, and translation. Certain lncRNAs protect genome integrity, while others respond to environmental cues like temperature, drought, nutrients, and pathogens. In this review, we explain the challenge of defining lncRNAs, introduce the machineries responsible for their production, and organize this knowledge by viewing the functions of lncRNAs throughout the structure of a typical plant gene.

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.

Integrated Genome-Scale Analysis and Northern Blot Detection of Retrotransposon siRNAs Across Plant Species.

Cells have sophisticated RNA-directed mechanisms to regulate genes, destroy viruses, or silence transposable elements (TEs). In terrestrial plants, a specialized non-coding RNA machinery involving RNA polymerase IV (Pol IV) and small interfering RNAs (siRNAs) targets DNA methylation and silencing to TEs. Here, we present a bioinformatics protocol for annotating and quantifying siRNAs that derive from long terminal repeat (LTR) retrotransposons. The approach was validated using small RNA northern blot analyses, comparing the species Arabidopsis thaliana and Brachypodium distachyon. To assist hybridization probe design, we configured a genome browser to show small RNA-seq mappings in distinct colors and shades according to their nucleotide lengths and abundances, respectively. Samples from wild-type and pol IV mutant plants, cross-species negative controls, and a conserved microRNA control validated the detected siRNA signals, confirming their origin from specific TEs and their Pol IV-dependent biogenesis. Moreover, an optimized labeling method yielded probes that could detect low-abundance siRNAs from B. distachyon TEs. The integration of de novo TE annotation, small RNA-seq profiling, and northern blotting, as outlined here, will facilitate the comparative genomic analysis of RNA silencing in crop plants and non-model species.

The NRPD1 N-terminus contains a Pol IV-specific motif that is critical for genome surveillance in Arabidopsis.

RNA-guided surveillance systems constrain the activity of transposable elements (TEs) in host genomes. In plants, RNA polymerase IV (Pol IV) transcribes TEs into primary transcripts from which RDR2 synthesizes double-stranded RNA precursors for small interfering RNAs (siRNAs) that guide TE methylation and silencing. How the core subunits of Pol IV, homologs of RNA polymerase II subunits, diverged to support siRNA biogenesis in a TE-rich, repressive chromatin context is not well understood. Here we studied the N-terminus of Pol IV's largest subunit, NRPD1. Arabidopsis lines harboring missense mutations in this N-terminus produce wild-type (WT) levels of NRPD1, which co-purifies with other Pol IV subunits and RDR2. Our in vitro transcription and genomic analyses reveal that the NRPD1 N-terminus is critical for robust Pol IV-dependent transcription, siRNA production and DNA methylation. However, residual RNA-directed DNA methylation observed in one mutant genotype indicates that Pol IV can operate uncoupled from the high siRNA levels typically observed in WT plants. This mutation disrupts a motif uniquely conserved in Pol IV, crippling the enzyme's ability to inhibit retrotransposon mobilization. We propose that the NRPD1 N-terminus motif evolved to regulate Pol IV function in genome surveillance.

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.

Mutation of Arabidopsis SMC4 identifies condensin as a corepressor of pericentromeric transposons and conditionally expressed genes.

In eukaryotes, transcriptionally inactive loci are enriched within highly condensed heterochromatin. In plants, as in mammals, the DNA of heterochromatin is densely methylated and wrapped by histones displaying a characteristic subset of post-translational modifications. Growing evidence indicates that these chromatin modifications are not sufficient for silencing. Instead, they are prerequisites for further assembly of higher-order chromatin structures that are refractory to transcription but not fully understood. We show that silencing of transposons in the pericentromeric heterochromatin of Arabidopsis thaliana requires SMC4, a core subunit of condensins I and II, acting in conjunction with CG methylation by MET1 (DNA METHYLTRANSFERASE 1), CHG methylation by CMT3 (CHROMOMETHYLASE 3), the chromatin remodeler DDM1 (DECREASE IN DNA METHYLATION 1), and histone modifications, including histone H3 Lys 27 monomethylation (H3K27me1), imparted by ATXR5 and ATXR6. SMC4/condensin also acts within the mostly euchromatic chromosome arms to suppress conditionally expressed genes involved in flowering or DNA repair, including the DNA glycosylase ROS1, which facilitates DNA demethylation. Collectively, our genome-wide analyses implicate condensin in the suppression of hundreds of loci, acting in both DNA methylation-dependent and methylation-independent pathways.

Hybrid incompatibility caused by an epiallele.

Hybrid incompatibility resulting from deleterious gene combinations is thought to be an important step toward reproductive isolation and speciation. Here, we demonstrate involvement of a silent epiallele in hybrid incompatibility. In Arabidopsis thaliana accession Cvi-0, one of the two copies of a duplicated histidine biosynthesis gene, HISN6A, is mutated, making HISN6B essential. In contrast, in accession Col-0, HISN6A is essential because HISN6B is not expressed. Owing to these differences, Cvi-0 × Col-0 hybrid progeny that are homozygous for both Cvi-0 HISN6A and Col-0 HISN6B do not survive. We show that HISN6B of Col-0 is not a defective pseudogene, but a stably silenced epiallele. Mutating HISTONE DEACETYLASE 6 (HDA6), or the cytosine methyltransferase genes MET1 or CMT3, erases HISN6B's silent locus identity, reanimating the gene to circumvent hisn6a lethality and hybrid incompatibility. These results show that HISN6-dependent hybrid lethality is a revertible epigenetic phenomenon and provide additional evidence that epigenetic variation has the potential to limit gene flow between diverging populations of a species.

Northern Blotting Techniques for Small RNAs.

Cells have evolved intricate RNA-directed mechanisms that destroy viruses, silence transposons, and regulate gene expression. These nucleic acid surveillance and gene silencing mechanisms rely upon the selective base-pairing of ~19-25 nt small RNAs to complementary RNA targets. This chapter describes northern blot hybridization techniques for the detection of such small RNAs. Blots spiked with synthetic standards are used to illustrate the detection specificity and sensitivity of DNA oligonucleotide probes. Known endogenous small RNAs are then analyzed in samples prepared from several model plants, including Arabidopsis thaliana, Nicotiana benthamiana, Oryza sativa, Zea mays, and Physcomitrella patens, as well as from the animals Drosophila melanogaster and Mus musculus. Finally, the value of northern blotting for dissecting small RNA biogenesis is shown using an example of virus infection in A. thaliana.

Chromosome-specific NOR inactivation explains selective rRNA gene silencing and dosage control in Arabidopsis.

In eukaryotes, scores of excess ribosomal RNA (rRNA) genes are silenced by repressive chromatin modifications. Given the near sequence identity of rRNA genes within a species, it is unclear how specific rRNA genes are reproducibly chosen for silencing. Using Arabidopsis thaliana ecotype (strain) Col-0, a systematic search identified sequence polymorphisms that differ between active and developmentally silenced rRNA gene subtypes. Recombinant inbred mapping populations derived from three different ecotype crosses were then used to map the chromosomal locations of silenced and active RNA gene subtypes. Importantly, silenced and active rRNA gene subtypes are not intermingled. All silenced rRNA gene subtypes mapped to the nucleolus organizer region (NOR) on chromosome 2 (NOR2). All active rRNA gene subtypes mapped to NOR4. Using an engineered A. thaliana line in which a portion of Col-0 chromosome 4 was replaced by sequences of another ecotype, we show that a major rRNA gene subtype silenced at NOR2 is active when introgressed into the genome at NOR4. Collectively, these results reveal that selective rRNA gene silencing is not regulated gene by gene based on mechanisms dependent on subtle gene sequence variation. Instead, we propose that a subchromosomal silencing mechanism operates on a multimegabase scale to inactivate NOR2.

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.

A two-step process for epigenetic inheritance in Arabidopsis.

In Arabidopsis, multisubunit RNA polymerases IV and V orchestrate RNA-directed DNA methylation (RdDM) and transcriptional silencing, but what identifies the loci to be silenced is unclear. We show that heritable silent locus identity at a specific subset of RdDM targets requires HISTONE DEACETYLASE 6 (HDA6) acting upstream of Pol IV recruitment and siRNA biogenesis. At these loci, epigenetic memory conferring silent locus identity is erased in hda6 mutants such that restoration of HDA6 activity cannot restore siRNA biogenesis or silencing. Silent locus identity is similarly lost in mutants for the cytosine maintenance methyltransferase, MET1. By contrast, pol IV or pol V mutants disrupt silencing without erasing silent locus identity, allowing restoration of Pol IV or Pol V function to restore silencing. Collectively, these observations indicate that silent locus specification and silencing are separable steps that together account for epigenetic inheritance of the silenced state.

Subnuclear partitioning of rRNA genes between the nucleolus and nucleoplasm reflects alternative epiallelic states.

Eukaryotes can have thousands of 45S ribosomal RNA (rRNA) genes, many of which are silenced during development. Using fluorescence-activated sorting techniques, we show that active rRNA genes in Arabidopsis thaliana are present within sorted nucleoli, whereas silenced rRNA genes are excluded. DNA methyltransferase (met1), histone deacetylase (hda6), or chromatin assembly (caf1) mutants that disrupt silencing abrogate this nucleoplasmic-nucleolar partitioning. Bisulfite sequencing data indicate that active nucleolar rRNA genes are nearly completely demethylated at promoter CGs, whereas silenced genes are nearly fully methylated. Collectively, the data reveal that rRNA genes occupy distinct but changeable nuclear territories according to their epigenetic state.

Primary and secondary siRNAs in geminivirus-induced gene silencing.

In plants, RNA silencing-based antiviral defense is mediated by Dicer-like (DCL) proteins producing short interfering (si)RNAs. In Arabidopsis infected with the bipartite circular DNA geminivirus Cabbage leaf curl virus (CaLCuV), four distinct DCLs produce 21, 22 and 24 nt viral siRNAs. Using deep sequencing and blot hybridization, we found that viral siRNAs of each size-class densely cover the entire viral genome sequences in both polarities, but highly abundant siRNAs correspond primarily to the leftward and rightward transcription units. Double-stranded RNA precursors of viral siRNAs can potentially be generated by host RDR-dependent RNA polymerase (RDR). However, genetic evidence revealed that CaLCuV siRNA biogenesis does not require RDR1, RDR2, or RDR6. By contrast, CaLCuV derivatives engineered to target 30 nt sequences of a GFP transgene by primary viral siRNAs trigger RDR6-dependent production of secondary siRNAs. Viral siRNAs targeting upstream of the GFP stop codon induce secondary siRNAs almost exclusively from sequences downstream of the target site. Conversely, viral siRNAs targeting the GFP 3'-untranslated region (UTR) induce secondary siRNAs mostly upstream of the target site. RDR6-dependent siRNA production is not necessary for robust GFP silencing, except when viral siRNAs targeted GFP 5'-UTR. Furthermore, viral siRNAs targeting the transgene enhancer region cause GFP silencing without secondary siRNA production. We conclude that the majority of viral siRNAs accumulating during geminiviral infection are RDR1/2/6-independent primary siRNAs. Double-stranded RNA precursors of these siRNAs are likely generated by bidirectional readthrough transcription of circular viral DNA by RNA polymerase II. Unlike transgenic mRNA, geminiviral mRNAs appear to be poor templates for RDR-dependent production of secondary siRNAs.

Histone methyltransferases regulating rRNA gene dose and dosage control in Arabidopsis.

Eukaryotes have hundreds of nearly identical 45S ribosomal RNA (rRNA) genes, each encoding the 18S, 5.8S, and 25S catalytic rRNAs. Because cellular demands for ribosomes and protein synthesis vary during development, the number of active rRNA genes is subject to dosage control. In genetic hybrids, one manifestation of dosage control is nucleolar dominance, an epigenetic phenomenon in which the rRNA genes of one progenitor are repressed. For instance, in Arabidopsis suecica, the allotetraploid hybrid of Arabidopsis thaliana and Arabidopsis arenosa, the A. thaliana-derived rRNA genes are selectively silenced. An analogous phenomenon occurs in nonhybrid A. thaliana, in which specific classes of rRNA gene variants are inactivated. An RNA-mediated knockdown screen identified SUVR4 {SUPPRESSOR OF VARIEGATION 3-9 [SU(VAR)3-9]-RELATED 4} as a histone H3 Lys 9 (H3K9) methyltransferase required for nucleolar dominance in A. suecica. H3K9 methyltransferases are also required for variant-specific silencing in A. thaliana, but SUVH5 [SU(VAR)3-9 HOMOLOG 5] and SUVH6, rather than SUVR4, are the key activities in this genomic context. Mutations disrupting the H3K27 methyltransferases ATXR5 or ATXR6 affect which rRNA gene variants are expressed or silenced, and in atxr5 atxr6 double mutants, dominance relationships among variants are reversed relative to wild type. Interestingly, these changes in gene expression are accompanied by changes in the relative abundance of the rRNA gene variants at the DNA level, including overreplication of the normally silenced class and decreased abundance of the normally dominant class. Collectively, our results indicate that histone methylation can affect both the doses of different variants and their differential silencing through the choice mechanisms that achieve dosage control.

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.

Massive production of small RNAs from a non-coding region of Cauliflower mosaic virus in plant defense and viral counter-defense.

To successfully infect plants, viruses must counteract small RNA-based host defense responses. During infection of Arabidopsis, Cauliflower mosaic pararetrovirus (CaMV) is transcribed into pregenomic 35S and subgenomic 19S RNAs. The 35S RNA is both reverse transcribed and also used as an mRNA with highly structured 600 nt leader. We found that this leader region is transcribed into long sense- and antisense-RNAs and spawns a massive quantity of 21, 22 and 24 nt viral small RNAs (vsRNAs), comparable to the entire complement of host-encoded small-interfering RNAs and microRNAs. Leader-derived vsRNAs were detected bound to the Argonaute 1 (AGO1) effector protein, unlike vsRNAs from other viral regions. Only negligible amounts of leader-derived vsRNAs were bound to AGO4. Genetic evidence showed that all four Dicer-like (DCL) proteins mediate vsRNA biogenesis, whereas the RNA polymerases Pol IV, Pol V, RDR1, RDR2 and RDR6 are not required for this process. Surprisingly, CaMV titers were not increased in dcl1/2/3/4 quadruple mutants that accumulate only residual amounts of vsRNAs. Ectopic expression of CaMV leader vsRNAs from an attenuated geminivirus led to increased accumulation of this chimeric virus. Thus, massive production of leader-derived vsRNAs does not restrict viral replication but may serve as a decoy diverting the silencing machinery from viral promoter and coding regions.

Nucleolin is required for DNA methylation state and the expression of rRNA gene variants in Arabidopsis thaliana.

In eukaryotes, 45S rRNA genes are arranged in tandem arrays in copy numbers ranging from several hundred to several thousand in plants. Although it is clear that not all copies are transcribed under normal growth conditions, the molecular basis controlling the expression of specific sets of rRNA genes remains unclear. Here, we report four major rRNA gene variants in Arabidopsis thaliana. Interestingly, while transcription of one of these rRNA variants is induced, the others are either repressed or remain unaltered in A. thaliana plants with a disrupted nucleolin-like protein gene (Atnuc-L1). Remarkably, the most highly represented rRNA gene variant, which is inactive in WT plants, is reactivated in Atnuc-L1 mutants. We show that accumulated pre-rRNAs originate from RNA Pol I transcription and are processed accurately. Moreover, we show that disruption of the AtNUC-L1 gene induces loss of symmetrical DNA methylation without affecting histone epigenetic marks at rRNA genes. Collectively, these data reveal a novel mechanism for rRNA gene transcriptional regulation in which the nucleolin protein plays a major role in controlling active and repressed rRNA gene variants in Arabidopsis.