Enhancers associated with unstable RNAs are rare in plants.

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Single-nucleus RNA-seq reveals that MBD5, MBD6, and SILENZIO maintain silencing in the vegetative cell of developing pollen.

Silencing of transposable elements (TEs) drives the evolution of numerous redundant mechanisms of transcriptional regulation. Arabidopsis MBD5, MBD6, and SILENZIO act as TE repressors downstream of DNA methylation. Here, we show, via single-nucleus RNA-seq of developing male gametophytes, that these repressors are critical for TE silencing in the pollen vegetative cell, a companion cell important for fertilization that undergoes chromatin decompaction. Instead, other silencing mutants (met1, ddm1, mom1, morc) show loss of silencing in all pollen nucleus types and somatic cells. We show that TEs repressed by MBD5/6 gain chromatin accessibility in wild-type vegetative nuclei despite remaining silent, suggesting that loss of DNA compaction makes them sensitive to loss of MBD5/6. Consistently, crossing mbd5/6 to histone 1 mutants, which have decondensed chromatin in leaves, reveals derepression of MBD5/6-dependent TEs in leaves. MBD5/6 and SILENZIO thus act as a silencing system particularly important when chromatin compaction is compromised.

The role of ATXR6 expression in modulating genome stability and transposable element repression in Arabidopsis.

ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) AND ATXR6 are required for the deposition of H3K27me1 and for maintaining genomic stability in Arabidopsis Reduction of ATXR5/6 activity results in activation of DNA damage response genes, along with tissue-specific derepression of transposable elements (TEs), chromocenter decompaction, and genomic instability characterized by accumulation of excess DNA from heterochromatin. How loss of ATXR5/6 and H3K27me1 leads to these phenotypes remains unclear. Here we provide extensive characterization of the atxr5/6 hypomorphic mutant by comprehensively examining gene expression and epigenetic changes in the mutant. We found that the tissue-specific phenotypes of TE derepression and excessive DNA in this atxr5/6 mutant correlated with residual ATXR6 expression from the hypomorphic ATXR6 allele. However, up-regulation of DNA damage genes occurred regardless of ATXR6 levels and thus appears to be a separable process. We also isolated an atxr6-null allele which showed that ATXR5 and ATXR6 are required for female germline development. Finally, we characterize three previously reported suppressors of the hypomorphic atxr5/6 mutant and show that these rescue atxr5/6 via distinct mechanisms, two of which involve increasing H3K27me1 levels.

Arabidopsis MORC proteins function in the efficient establishment of RNA directed DNA methylation.

The Microrchidia (MORC) family of ATPases are required for transposable element (TE) silencing and heterochromatin condensation in plants and animals, and C. elegans MORC-1 has been shown to topologically entrap and condense DNA. In Arabidopsis thaliana, mutation of MORCs has been shown to reactivate silent methylated genes and transposons and to decondense heterochromatic chromocenters, despite only minor changes in the maintenance of DNA methylation. Here we provide the first evidence localizing Arabidopsis MORC proteins to specific regions of chromatin and find that MORC4 and MORC7 are closely co-localized with sites of RNA-directed DNA methylation (RdDM). We further show that MORC7, when tethered to DNA by an artificial zinc finger, can facilitate the establishment of RdDM. Finally, we show that MORCs are required for the efficient RdDM mediated establishment of DNA methylation and silencing of a newly integrated FWA transgene, even though morc mutations have no effect on the maintenance of preexisting methylation at the endogenous FWA gene. We propose that MORCs function as a molecular tether in RdDM complexes to reinforce RdDM activity for methylation establishment. These findings have implications for MORC protein function in a variety of other eukaryotic organisms.

Transcriptional and imprinting complexity in Arabidopsis seeds at single-nucleus resolution.

Seeds are a key life cycle stage for many plants. Seeds are also the basis of agriculture and the primary source of calories consumed by humans1. Here, we employ single-nucleus RNA-sequencing to generate a transcriptional atlas of developing Arabidopsis thaliana seeds, with a focus on endosperm. Endosperm, the primary site of gene imprinting in flowering plants, mediates the relationship between the maternal parent and the embryo2. We identify transcriptionally uncharacterized nuclei types in the chalazal endosperm, which interfaces with maternal tissue for nutrient unloading3,4. We demonstrate that the extent of parental bias of maternally expressed imprinted genes varies with cell-cycle phase, and that imprinting of paternally expressed imprinted genes is strongest in chalazal endosperm. Thus, imprinting is spatially and temporally heterogeneous. Increased paternal expression in the chalazal region suggests that parental conflict, which is proposed to drive imprinting evolution, is fiercest at the boundary between filial and maternal tissues.

CRISPR-based targeting of DNA methylation in Arabidopsis thaliana by a bacterial CG-specific DNA methyltransferase.

CRISPR-based targeted modification of epigenetic marks such as DNA cytosine methylation is an important strategy to regulate the expression of genes and their associated phenotypes. Although plants have DNA methylation in all sequence contexts (CG, CHG, CHH, where H = A, T, C), methylation in the symmetric CG context is particularly important for gene silencing and is very efficiently maintained through mitotic and meiotic cell divisions. Tools that can directly add CG methylation to specific loci are therefore highly desirable but are currently lacking in plants. Here we have developed two CRISPR-based CG-specific targeted DNA methylation systems for plants using a variant of the bacterial CG-specific DNA methyltransferase MQ1 with reduced activity but high specificity. We demonstrate that the methylation added by MQ1 is highly target specific and can be heritably maintained in the absence of the effector. These tools should be valuable both in crop engineering and in plant genetic research.

The histone variant H2A.W and linker histone H1 co-regulate heterochromatin accessibility and DNA methylation.

In flowering plants, heterochromatin is demarcated by the histone variant H2A.W, elevated levels of the linker histone H1, and specific epigenetic modifications, such as high levels of DNA methylation at both CG and non-CG sites. How H2A.W regulates heterochromatin organization and interacts with other heterochromatic features is unclear. Here, we create a h2a.w null mutant via CRISPR-Cas9, h2a.w-2, to analyze the in vivo function of H2A.W. We find that H2A.W antagonizes deposition of H1 at heterochromatin and that non-CG methylation and accessibility are moderately decreased in h2a.w-2 heterochromatin. Compared to H1 loss alone, combined loss of H1 and H2A.W greatly increases accessibility and facilitates non-CG DNA methylation in heterochromatin, suggesting co-regulation of heterochromatic features by H2A.W and H1. Our results suggest that H2A.W helps maintain optimal heterochromatin accessibility and DNA methylation by promoting chromatin compaction together with H1, while also inhibiting excessive H1 incorporation.

A viral guide RNA delivery system for CRISPR-based transcriptional activation and heritable targeted DNA demethylation in Arabidopsis thaliana.

Plant RNA viruses are used as delivery vectors for their high level of accumulation and efficient spread during virus multiplication and movement. Utilizing this concept, several viral-based guide RNA delivery platforms for CRISPR-Cas9 genome editing have been developed. The CRISPR-Cas9 system has also been adapted for epigenome editing. While systems have been developed for CRISPR-Cas9 based gene activation or site-specific DNA demethylation, viral delivery of guide RNAs remains to be developed for these purposes. To address this gap we have developed a tobacco rattle virus (TRV)-based single guide RNA delivery system for epigenome editing in Arabidopsis thaliana. Because tRNA-like sequences have been shown to facilitate the cell-to-cell movement of RNAs in plants, we used the tRNA-guide RNA expression system to express guide RNAs from the viral genome to promote heritable epigenome editing. We demonstrate that the tRNA-gRNA system with TRV can be used for both transcriptional activation and targeted DNA demethylation of the FLOWERING WAGENINGEN gene in Arabidopsis. We achieved up to ~8% heritability of the induced demethylation phenotype in the progeny of virus inoculated plants. We did not detect the virus in the next generation, indicating effective clearance of the virus from plant tissues. Thus, TRV delivery, combined with a specific tRNA-gRNA architecture, provides for fast and effective epigenome editing.

RNA-directed DNA Methylation.

RNA-directed DNA methylation (RdDM) is a biological process in which non-coding RNA molecules direct the addition of DNA methylation to specific DNA sequences. The RdDM pathway is unique to plants, although other mechanisms of RNA-directed chromatin modification have also been described in fungi and animals. To date, the RdDM pathway is best characterized within angiosperms (flowering plants), and particularly within the model plant Arabidopsis thaliana. However, conserved RdDM pathway components and associated small RNAs (sRNAs) have also been found in other groups of plants, such as gymnosperms and ferns. The RdDM pathway closely resembles other sRNA pathways, particularly the highly conserved RNAi pathway found in fungi, plants, and animals. Both the RdDM and RNAi pathways produce sRNAs and involve conserved Argonaute, Dicer and RNA-dependent RNA polymerase proteins. RdDM has been implicated in a number of regulatory processes in plants. The DNA methylation added by RdDM is generally associated with transcriptional repression of the genetic sequences targeted by the pathway. Since DNA methylation patterns in plants are heritable, these changes can often be stably transmitted to progeny. As a result, one prominent role of RdDM is the stable, transgenerational suppression of transposable element (TE) activity. RdDM has also been linked to pathogen defense, abiotic stress responses, and the regulation of several key developmental transitions. Although the RdDM pathway has a number of important functions, RdDM-defective mutants in Arabidopsis thaliana are viable and can reproduce, which has enabled detailed genetic studies of the pathway. However, RdDM mutants can have a range of defects in different plant species, including lethality, altered reproductive phenotypes, TE upregulation and genome instability, and increased pathogen sensitivity. Overall, RdDM is an important pathway in plants that regulates a number of processes by establishing and reinforcing specific DNA methylation patterns, which can lead to transgenerational epigenetic effects on gene expression and phenotype.

Identification and Comparison of Imprinted Genes Across Plant Species.

Genomic imprinting is a phenomenon that occurs in flowering plants and mammals, whereby a gene is expressed in a parent-of-origin-specific manner. Although imprinting has now been examined genome-wide in a number of species using RNA-seq, the analyses used to assess imprinting vary between studies, making consistent comparisons between species difficult. Here we present a simple, easy-to-use bioinformatic pipeline for imprinting analyses suitable for any tissue, including plant endosperm. All relevant scripts can be downloaded. As an illustrative example, we reanalyze published data from A. thaliana and Z. mays endosperm using the pipeline and then demonstrate how to use the results to assess the conservation of imprinting between these species. We also introduce the Plant Imprinting Database, a repository for published imprinting datasets in plants that can be used to view, compare, and download data.

Proximal methylation features associated with nonrandom changes in gene body methylation.

BACKGROUND: Gene body methylation at CG dinucleotides is a widely conserved feature of methylated genomes but remains poorly understood. The Arabidopsis thaliana strain Cvi has depleted gene body methylation relative to the reference strain Col. Here, we leverage this natural epigenetic difference to investigate gene body methylation stability. RESULTS: Recombinant inbred lines derived from Col and Cvi were used to examine the transmission of distinct gene body methylation states. The vast majority of genic CG methylation patterns are faithfully transmitted over nine generations according to parental genotype, with only 1-4% of CGs either losing or gaining methylation relative to the parent. Genic CGs that fail to maintain the parental methylation state are shared among independent lines, suggesting that these are not random occurrences. We use a logistic regression framework to identify features that best predict sites that fail to maintain parental methylation state. Intermediate levels of CG methylation around a dynamic CG site and high methylation variability across many A. thaliana strains at that site are the strongest predictors. These data suggest that the dynamic CGs we identify are not specific to the Col-Cvi recombinant inbred lines but have an epigenetic state that is inherently less stable within the A. thaliana species. Extending this, variably methylated genic CGs in maize and Brachypodium distachyon are also associated with intermediate local CG methylation. CONCLUSIONS: These results provide new insights into the features determining the inheritance of gene body methylation and demonstrate that two different methylation equilibria can be maintained within single individuals.

Conserved imprinting associated with unique epigenetic signatures in the Arabidopsis genus.

In plants, imprinted gene expression occurs in endosperm seed tissue and is sometimes associated with differential DNA methylation between maternal and paternal alleles(1). Imprinting is theorized to have been selected for because of conflict between parental genomes in offspring(2), but most studies of imprinting have been conducted in Arabidopsis thaliana, an inbred primarily self-fertilizing species that should have limited parental conflict. We examined embryo and endosperm allele-specific expression and DNA methylation genome-wide in the wild outcrossing species Arabidopsis lyrata. Here we show that the majority of A. lyrata imprinted genes also exhibit parentally biased expression in A. thaliana, suggesting that there is evolutionary conservation in gene imprinting. Surprisingly, we discovered substantial interspecies differences in methylation features associated with paternally expressed imprinted genes (PEGs). Unlike in A. thaliana, the maternal allele of many A. lyrata PEGs was hypermethylated in the CHG context. Increased maternal allele CHG methylation was associated with increased expression bias in favour of the paternal allele. We propose that CHG methylation maintains or reinforces repression of maternal alleles of PEGs. These data suggest that the genes subject to imprinting are largely conserved, but there is flexibility in the epigenetic mechanisms employed between closely related species to maintain monoallelic expression. This supports the idea that imprinting of specific genes is a functional phenomenon, and not simply a byproduct of seed epigenomic reprogramming.

Natural epigenetic polymorphisms lead to intraspecific variation in Arabidopsis gene imprinting.

Imprinted gene expression occurs during seed development in plants and is associated with differential DNA methylation of parental alleles, particularly at proximal transposable elements (TEs). Imprinting variability could contribute to observed parent-of-origin effects on seed development. We investigated intraspecific variation in imprinting, coupled with analysis of DNA methylation and small RNAs, among three Arabidopsis strains with diverse seed phenotypes. The majority of imprinted genes were parentally biased in the same manner among all strains. However, we identified several examples of allele-specific imprinting correlated with intraspecific epigenetic variation at a TE. We successfully predicted imprinting in additional strains based on methylation variability. We conclude that there is standing variation in imprinting even in recently diverged genotypes due to intraspecific epiallelic variation. Our data demonstrate that epiallelic variation and genomic imprinting intersect to produce novel gene expression patterns in seeds.

Host cell factor 1 inhibits SKN-1 to modulate oxidative stress responses in Caenorhabditis elegans.

Host cell factor-1 (HCF-1) is a conserved regulator of the longevity and stress response functions of DAF-16/FOXO. SKN-1 transcription factor is an evolutionarily conserved xenobiotic stress regulator and a pro-longevity factor. Here, we demonstrate that SKN-1 contributes to the enhanced oxidative stress resistance incurred by hcf-1 mutation in C. elegans. HCF-1 prevents the nuclear accumulation of SKN-1 and represses the transcriptional activation of SKN-1 specifically at target genes involved in cellular detoxification pathways. Our findings reveal a novel and context-specific regulatory relationship between two highly conserved longevity and stress response factors HCF-1 and SKN-1.

The evolutionarily conserved longevity determinants HCF-1 and SIR-2.1/SIRT1 collaborate to regulate DAF-16/FOXO.

The conserved DAF-16/FOXO transcription factors and SIR-2.1/SIRT1 deacetylases are critical for diverse biological processes, particularly longevity and stress response; and complex regulation of DAF-16/FOXO by SIR-2.1/SIRT1 is central to appropriate biological outcomes. Caenorhabditis elegans Host Cell Factor 1 (HCF-1) is a longevity determinant previously shown to act as a co-repressor of DAF-16. We report here that HCF-1 represents an integral player in the regulatory loop linking SIR-2.1/SIRT1 and DAF-16/FOXO in both worms and mammals. Genetic analyses showed that hcf-1 acts downstream of sir-2.1 to influence lifespan and oxidative stress response in C. elegans. Gene expression profiling revealed a striking 80% overlap between the DAF-16 target genes responsive to hcf-1 mutation and sir-2.1 overexpression. Subsequent GO-term analyses of HCF-1 and SIR-2.1-coregulated DAF-16 targets suggested that HCF-1 and SIR-2.1 together regulate specific aspects of DAF-16-mediated transcription particularly important for aging and stress responses. Analogous to its role in regulating DAF-16/SIR-2.1 target genes in C. elegans, the mammalian HCF-1 also repressed the expression of several FOXO/SIRT1 target genes. Protein-protein association studies demonstrated that SIR-2.1/SIRT1 and HCF-1 form protein complexes in worms and mammalian cells, highlighting the conservation of their regulatory relationship. Our findings uncover a conserved interaction between the key longevity determinants SIR-2.1/SIRT1 and HCF-1, and they provide new insights into the complex regulation of FOXO proteins.