Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements.

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Replacement of Arabidopsis H2A.Z with human H2A.Z orthologs reveals extensive functional conservation and limited importance of the N-terminal tail sequence for Arabidopsis development.

The incorporation of histone variants, distinct paralogs of core histones, into chromatin affects all DNA-templated processes in the cell, including the regulation of transcription. In recent years, much research has been focused on H2A.Z, an evolutionarily conserved H2A variant found in all eukaryotes. In order to investigate the functional conservation of H2A.Z histones during eukaryotic evolution we transformed h2a.z deficient plants with three human H2A.Z proteins to assess their ability to rescue the mutant defects. We discovered that human H2A.Z.1 and H2A.Z.2.1 fully complement the phenotypic abnormalities of h2a.z plants despite the fact that Arabidopsis and human H2A.Z N-terminal tail sequences are quite divergent. In contrast, the brain-specific splice variant H2A.Z.2.2 has a dominant-negative effect in wild-type plants. Furthermore, H2A.Z.1 almost completely re-establishes normal H2A.Z chromatin occupancy in h2a.z plants and restores the transcript levels of more than 84 % of misexpressed genes. Finally, our hypothesis that the N-terminal tail of Arabidopsis H2A.Z is not crucial for its developmental functions was supported by the ability of N-terminal end truncations of Arabidopsis HTA11 to largely rescue the defects of h2a.z mutants.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements.

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Gene regulatory networks shape developmental plasticity of root cell types under water extremes in rice.

Understanding how roots modulate development under varied irrigation or rainfall is crucial for development of climate-resilient crops. We established a toolbox of tagged rice lines to profile translating mRNAs and chromatin accessibility within specific cell populations. We used these to study roots in a range of environments: plates in the lab, controlled greenhouse stress and recovery conditions, and outdoors in a paddy. Integration of chromatin and mRNA data resolves regulatory networks of the following: cycle genes in proliferating cells that attenuate DNA synthesis under submergence; genes involved in auxin signaling, the circadian clock, and small RNA regulation in ground tissue; and suberin biosynthesis, iron transporters, and nitrogen assimilation in endodermal/exodermal cells modulated with water availability. By applying a systems approach, we identify known and candidate driver transcription factors of water-deficit responses and xylem development plasticity. Collectively, this resource will facilitate genetic improvements in root systems for optimal climate resilience.

The VIL gene CRAWLING ELEPHANT controls maturation and differentiation in tomato via polycomb silencing.

VERNALIZATION INSENSITIVE 3-LIKE (VIL) proteins are PHD-finger proteins that recruit the repressor complex Polycomb Repressive Complex 2 (PRC2) to the promoters of target genes. Most known VIL targets are flowering repressor genes. Here, we show that the tomato VIL gene CRAWLING ELEPHANT (CREL) promotes differentiation throughout plant development by facilitating the trimethylation of Histone H3 on lysine 27 (H3K27me3). We identified the crel mutant in a screen for suppressors of the simple-leaf phenotype of entire (e), a mutant in the AUX/IAA gene ENTIRE/SlIAA9, involved in compound-leaf development in tomato. crel mutants have increased leaf complexity, and suppress the ectopic blade growth of e mutants. In addition, crel mutants are late flowering, and have delayed and aberrant stem, root and flower development. Consistent with a role for CREL in recruiting PRC2, crel mutants show drastically reduced H3K27me3 enrichment at approximately half of the 14,789 sites enriched in wild-type plants, along with upregulation of many underlying genes. Interestingly, this reduction in H3K27me3 across the genome in crel is also associated with gains in H3K27me3 at a smaller number of sites that normally have modest levels of the mark in wild-type plants, suggesting that PRC2 activity is no longer limiting in the absence of CREL. Our results uncover a wide role for CREL in plant and organ differentiation in tomato and suggest that CREL is required for targeting PRC2 activity to, and thus silencing, a specific subset of polycomb targets.

Histone Variants in the Specialization of Plant Chromatin.

The basic unit of chromatin, the nucleosome, is an octamer of four core histone proteins (H2A, H2B, H3, and H4) and serves as a fundamental regulatory unit in all DNA-templated processes. The majority of nucleosome assembly occurs during DNA replication when these core histones are produced en masse to accommodate the nascent genome. In addition, there are a number of nonallelic sequence variants of H2A and H3 in particular, known as histone variants, that can be incorporated into nucleosomes in a targeted and replication-independent manner. By virtue of their sequence divergence from the replication-coupled histones, these histone variants can impart unique properties onto the nucleosomes they occupy and thereby influence transcription and epigenetic states, DNA repair, chromosome segregation, and other nuclear processes in ways that profoundly affect plant biology. In this review, we discuss the evolutionary origins of these variants in plants, their known roles in chromatin, and their impacts on plant development and stress responses. We focus on the individual and combined roles of histone variants in transcriptional regulation within euchromatic and heterochromatic genome regions. Finally, we highlight gaps in our understanding of plant variants at the molecular, cellular, and organismal levels, and we propose new directions for study in the field of plant histone variants.

Innovation, conservation, and repurposing of gene function in root cell type development.

Plant species have evolved myriads of solutions, including complex cell type development and regulation, to adapt to dynamic environments. To understand this cellular diversity, we profiled tomato root cell type translatomes. Using xylem differentiation in tomato, examples of functional innovation, repurposing, and conservation of transcription factors are described, relative to the model plant Arabidopsis. Repurposing and innovation of genes are further observed within an exodermis regulatory network and illustrate its function. Comparative translatome analyses of rice, tomato, and Arabidopsis cell populations suggest increased expression conservation of root meristems compared with other homologous populations. In addition, the functions of constitutively expressed genes are more conserved than those of cell type/tissue-enriched genes. These observations suggest that higher order properties of cell type and pan-cell type regulation are evolutionarily conserved between plants and animals.

Considerations in the analysis of plant chromatin accessibility data.

Transcriptional control is exerted primarily through the binding of transcription factor proteins to regulatory elements in DNA. By virtue of eukaryotic DNA being complexed with histones, transcription factor binding to DNA alters or eliminates histone-DNA contacts, leading to increased accessibility of the DNA region to nuclease enzymes. This hypersensitivity to nuclease digestion has been used to define DNA binding events and regulatory elements across genomes, and to compare these attributes between cell types or conditions. These approaches make it possible to define the regulatory elements in a genome as well as to predict the regulatory networks of transcription factors and their target genes in a given cell state. As these chromatin accessibility assays are increasingly used, it is important to consider how to analyze the resulting data to avoid artifactual results or misinterpretation. In this review, we focus on some of the key technical and computational caveats associated with plant chromatin accessibility data, including strategies for sample preparation, sequencing, read mapping, and downstream analyses.

Evolutionary flexibility in flooding response circuitry in angiosperms.

Flooding due to extreme weather threatens crops and ecosystems. To understand variation in gene regulatory networks activated by submergence, we conducted a high-resolution analysis of chromatin accessibility and gene expression at three scales of transcript control in four angiosperms, ranging from a dryland-adapted wild species to a wetland crop. The data define a cohort of conserved submergence-activated genes with signatures of overlapping cis regulation by four transcription factor families. Syntenic genes are more highly expressed than nonsyntenic genes, yet both can have the cis motifs and chromatin accessibility associated with submergence up-regulation. Whereas the flexible circuitry spans the eudicot-monocot divide, the frequency of specific cis motifs, extent of chromatin accessibility, and degree of submergence activation are more prevalent in the wetland crop and may have adaptive importance.

Methyl-CpG-binding domain 9 (MBD9) is required for H2A.Z incorporation into chromatin at a subset of H2A.Z-enriched regions in the Arabidopsis genome.

The SWR1 chromatin remodeling complex, which deposits the histone variant H2A.Z into nucleosomes, has been well characterized in yeast and animals, but its composition in plants has remained uncertain. We used the conserved SWR1 subunit ACTIN RELATED PROTEIN 6 (ARP6) as bait in tandem affinity purification experiments to isolate associated proteins from Arabidopsis thaliana. We identified all 11 subunits found in yeast SWR1 and the homologous mammalian SRCAP complexes, demonstrating that this complex is conserved in plants. We also identified several additional proteins not previously associated with SWR1, including Methyl-CpG-BINDING DOMAIN 9 (MBD9) and three members of the Alfin1-like protein family, all of which have been shown to bind modified histone tails. Since mbd9 mutant plants were phenotypically similar to arp6 mutants, we explored a potential role for MBD9 in H2A.Z deposition. We found that MBD9 is required for proper H2A.Z incorporation at thousands of discrete sites, which represent a subset of the genomic regions normally enriched with H2A.Z. We also discovered that MBD9 preferentially interacts with acetylated histone H4 peptides, as well as those carrying mono- or dimethylated H3 lysine 4, or dimethylated H3 arginine 2 or 8. Considering that MBD9-dependent H2A.Z sites show a distinct histone modification profile, we propose that MBD9 recognizes particular nucleosome modifications via its PHD- and Bromo-domains and thereby guides SWR1 to these sites for H2A.Z deposition. Our data establish the SWR1 complex as being conserved across eukaryotes and suggest that MBD9 may be involved in targeting the complex to specific genomic sites through nucleosomal interactions. The finding that MBD9 does not appear to be a core subunit of the Arabidopsis SWR1 complex, along with the synergistic phenotype of arp6;mbd9 double mutants, suggests that MBD9 also has important roles beyond H2A.Z deposition.

The histone variant H2A.Z and chromatin remodeler BRAHMA act coordinately and antagonistically to regulate transcription and nucleosome dynamics in Arabidopsis.

Plants adapt to environmental changes by regulating transcription and chromatin organization. The histone H2A variant H2A.Z and the SWI2/SNF2 ATPase BRAHMA (BRM) have overlapping roles in positively and negatively regulating environmentally responsive genes in Arabidopsis, but the extent of this overlap was uncharacterized. Both factors have been associated with various changes in nucleosome positioning and stability in different contexts, but their specific roles in transcriptional regulation and chromatin organization need further characterization. We show that H2A.Z and BRM co-localize at thousands of sites, where they interact both cooperatively and antagonistically in transcriptional repression and activation of genes involved in development and responses to environmental stimuli. We identified eight classes of genes that show distinct relationships between H2A.Z and BRM with respect to their roles in transcription. These include activating and silencing transcription both redundantly and antagonistically. We found that H2A.Z contributes to a range of different nucleosome properties, while BRM stabilizes nucleosomes where it binds and destabilizes or repositions flanking nucleosomes. We also found that, at many genes regulated by both BRM and H2A.Z, both factors overlap with binding sites of the light-regulated transcription factor FAR1-Related Sequence 9 (FRS9) and that a subset of these FRS9 binding sites are dependent on H2A.Z and BRM for accessibility. Collectively, we comprehensively characterized the antagonistic and cooperative contributions of H2A.Z and BRM to transcriptional regulation, and illuminated several interrelated roles in chromatin organization. The variability observed in their individual functions implies that both BRM and H2A.Z have more context-dependent roles than previously assumed.

The Chromatin Remodelers PKL and PIE1 Act in an Epigenetic Pathway That Determines H3K27me3 Homeostasis in Arabidopsis.

Selective, tissue-specific gene expression is facilitated by the epigenetic modification H3K27me3 (trimethylation of lysine 27 on histone H3) in plants and animals. Much remains to be learned about how H3K27me3-enriched chromatin states are constructed and maintained. Here, we identify a genetic interaction in Arabidopsis thaliana between the chromodomain helicase DNA binding chromatin remodeler PICKLE (PKL), which promotes H3K27me3 enrichment, and the SWR1-family remodeler PHOTOPERIOD INDEPENDENT EARLY FLOWERING1 (PIE1), which incorporates the histone variant H2A.Z. Chromatin immunoprecipitation-sequencing and RNA-sequencing reveal that PKL, PIE1, and the H3K27 methyltransferase CURLY LEAF act in a common gene expression pathway and are required for H3K27me3 levels genome-wide. Additionally, H3K27me3-enriched genes are largely a subset of H2A.Z-enriched genes, further supporting the functional linkage between these marks. We also found that recombinant PKL acts as a prenucleosome maturation factor, indicating that it promotes retention of H3K27me3. These data support the existence of an epigenetic pathway in which PIE1 promotes H2A.Z, which in turn promotes H3K27me3 deposition. After deposition, PKL promotes retention of H3K27me3 after DNA replication and/or transcription. Our analyses thus reveal roles for H2A.Z and ATP-dependent remodelers in construction and maintenance of H3K27me3-enriched chromatin in plants.

Changes in chromatin accessibility between Arabidopsis stem cells and mesophyll cells illuminate cell type-specific transcription factor networks.

Cell differentiation is driven by changes in the activity of transcription factors (TFs) and subsequent alterations in transcription. To study this process, differences in TF binding between cell types can be deduced by probing chromatin accessibility. We used cell type-specific nuclear purification followed by the assay for transposase-accessible chromatin (ATAC-seq) to delineate differences in chromatin accessibility and TF regulatory networks between stem cells of the shoot apical meristem (SAM) and differentiated leaf mesophyll cells in Arabidopsis thaliana. Chromatin accessibility profiles of SAM stem cells and leaf mesophyll cells were very similar at a qualitative level, yet thousands of regions having quantitatively different chromatin accessibility were also identified. Analysis of the genomic regions preferentially accessible in each cell type identified hundreds of overrepresented TF-binding motifs, highlighting sets of TFs that are probably important for each cell type. Within these sets, we found evidence for extensive co-regulation of target genes by multiple TFs that are preferentially expressed in each cell type. Interestingly, the TFs within each of these cell type-enriched sets also showed evidence of extensively co-regulating each other. We further found that preferentially accessible chromatin regions in mesophyll cells tended to also be substantially accessible in the stem cells, whereas the converse was not true. This observation suggests that the generally higher accessibility of regulatory elements in stem cells might contribute to their developmental plasticity. This work demonstrates the utility of cell type-specific chromatin accessibility profiling for the rapid development of testable models of regulatory control differences between cell types.

Profiling of Accessible Chromatin Regions across Multiple Plant Species and Cell Types Reveals Common Gene Regulatory Principles and New Control Modules.

The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Identification of Open Chromatin Regions in Plant Genomes Using ATAC-Seq.

Identifying and characterizing highly accessible chromatin regions assists in determining the location of genomic regulatory elements and understanding transcriptional regulation. In this chapter, we describe an approach to map accessible chromatin features in plants using the Assay for Transposase-Accessible Chromatin, combined with high-throughput sequencing (ATAC-seq), which was originally developed for cultured animal cells. This technique utilizes a hyperactive Tn5 transposase to cause DNA cleavage and simultaneous insertion of sequencing adapters into open chromatin regions of the input nuclei. The application of ATAC-seq to plant tissue has been challenging due to the difficulty of isolating nuclei sufficiently free of interfering organellar DNA. Here we present two different approaches to purify plant nuclei for ATAC-seq: the INTACT method (Isolation of Nuclei TAgged in specific Cell Types) to isolate nuclei from individual cell types of the plant, and tissue lysis followed by sucrose sedimentation to isolate sufficiently pure total nuclei. We provide detailed instructions for transposase treatment of nuclei isolated using either approach, as well as subsequent preparation of ATAC-seq libraries. Sequencing-ready ATAC-seq libraries can be prepared from plant tissue in as little as one day. The procedures described here are optimized for Arabidopsis thaliana but can also be applied to other plant species.

A Global View of RNA-Protein Interactions Identifies Post-transcriptional Regulators of Root Hair Cell Fate.

The Arabidopsis thaliana root epidermis is comprised of two cell types, hair and nonhair cells, which differentiate from the same precursor. Although the transcriptional programs regulating these events are well studied, post-transcriptional factors functioning in this cell fate decision are mostly unknown. Here, we globally identify RNA-protein interactions and RNA secondary structure in hair and nonhair cell nuclei. This analysis reveals distinct structural and protein binding patterns across both transcriptomes, allowing identification of differential RNA binding protein (RBP) recognition sites. Using these sequences, we identify two RBPs that regulate hair cell development. Specifically, we find that SERRATE functions in a microRNA-dependent manner to inhibit hair cell fate, while also terminating growth of root hairs mostly independent of microRNA biogenesis. In addition, we show that GLYCINE-RICH PROTEIN 8 promotes hair cell fate while alleviating phosphate starvation stress. In total, this global analysis reveals post-transcriptional regulators of plant root epidermal cell fate.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Epigenome profiling of specific plant cell types using a streamlined INTACT protocol and ChIP-seq.

Plants consist of many functionally specialized cell types, each with its own unique epigenome, transcriptome, and proteome. Characterization of these cell type-specific properties is essential to understanding cell fate specification and the responses of individual cell types to the environment. In this chapter we describe an approach to map chromatin features in specific cell types of Arabidopsis thaliana using nuclei purification from individual cell types with the INTACT method (isolation of nuclei tagged in specific cell types) followed by chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq). The INTACT system employs two transgenes to generate affinity-labeled nuclei in the cell type of interest, and these tagged nuclei can then be selectively purified from tissue homogenates. The primary transgene encodes the nuclear tagging fusion protein (NTF), which consists of a nuclear envelope-targeting domain, the green fluorescent protein, and a biotin ligase recognition peptide, while the second transgene encodes the E. coli biotin ligase (BirA), which selectively biotinylates NTF. Expression of NTF and BirA in a specific cell type thus yields nuclei that are coated with biotin and can be purified by virtue of their affinity for streptavidin-coated magnetic beads. Compared with the original INTACT nuclei purification protocol, the procedure presented here is greatly simplified and shortened. After nuclei purification, we provide detailed instructions for chromatin isolation, shearing, and immunoprecipitation. Finally, we present a low input ChIP-seq library preparation protocol based on the nano-ChIP-seq method of Adli and Bernstein, and we describe multiplex Illumina sequencing of these libraries to produce high quality, cell type-specific epigenome profiles at a relatively low cost. The procedures given here are optimized for Arabidopsis but should be easily adaptable to other plant species.

Global analysis of the RNA-protein interaction and RNA secondary structure landscapes of the Arabidopsis nucleus.

Posttranscriptional regulation in eukaryotes requires cis- and trans-acting features and factors including RNA secondary structure and RNA-binding proteins (RBPs). However, a comprehensive view of the structural and RBP interaction landscape of nuclear RNAs has yet to be compiled for any organism. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approaches to globally profile these features in Arabidopsis seedling nuclei in vivo. We reveal anticorrelated patterns of secondary structure and RBP binding throughout nuclear mRNAs that demarcate sites of alternative splicing and polyadenylation. We also uncover a collection of protein-bound sequence motifs, and identify their structural contexts, co-occurrences in transcripts encoding functionally related proteins, and interactions with putative RBPs. Finally, using these motifs, we find that the chloroplast RBP CP29A also interacts with nuclear mRNAs. In total, we provide a simultaneous view of the RNA secondary structure and RBP interaction landscapes in a eukaryotic nucleus.