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 Arabidopsis thaliana plants with each of the 3 human H2A.Z variants 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 significant divergence in the N-terminal tail sequences of Arabidopsis and human H2A.Zs. In contrast, the brain-specific splice variant H2A.Z.2.2 has a dominant-negative effect in wild-type plants, mimicking an H2A.Z deficiency phenotype. Furthermore, human H2A.Z.1 almost completely reestablishes normal H2A.Z chromatin occupancy in h2a.z plants and restores the expression of more than 84% of misexpressed genes. Finally, we used a series of N-terminal tail truncations of Arabidopsis HTA11 to reveal that the N-terminal tail of Arabidopsis H2A.Z is not necessary for normal plant development under optimal growth conditions but does play an important role in mounting proper abiotic stress responses.

Dynamic protein post-translational methylation is essential for cellular function, highlighted by the essential role of methylation in transcriptional regulation and its aberrant dysregulation in diseases, including cancer. This underscores the importance of cataloging the cellular methylproteome. However, comprehensive analysis of the methylproteome remains elusive due to limitations in current enrichment and analysis pipelines. Here, we employ an l-methionine analogue, ProSeMet, that is chemoenzymatically converted to the SAM analogue ProSeAM in cells and in vivo to tag proteins with a biorthogonal alkyne that can be directly detected via liquid chromatography and tandem mass spectrometry (LC-MS/MS), or functionalized for subsequent selective enrichment and LC-MS/MS identification. Without enrichment, we identify known and novel lysine mono-, di-, and tripargylation, histidine propargylation, and arginine propargylation with site-specific resolution on proteins including heat shock protein HSPA8, the translational elongation factor eEF1A1, and the metabolic enzyme phosphoglycerate mutase 1, or PGAM1, for which methylation has been implicated in human disease. With enrichment, we identify 486 proteins known to be methylated and 221 proteins with novel propargylation sites encompassing diverse cellular functions. Systemic ProSeMet delivery in mice propargylates proteins across organ systems with blood-brain barrier penetrance and identifies site-specific propargylation in vivo with LC-MS/MS. Leveraging these pipelines to define the cellular methylproteome may have broad applications for understanding the methylproteome in the context of disease.

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.

The histone variant H2A.Z is important for transcriptional regulation across eukaryotes, where it can alternately promote or repress transcription. In plants, actively transcribed genes show H2A.Z enrichment in nucleosomes immediately downstream of the transcription start site (TSS), while silent genes show H2A.Z enrichment across the gene body. Previous work showed that silent genes responsive to temperature and far-red light lose gene body H2A.Z upon activation, but whether H2A.Z loss is generally required for transcription is not clear. We profiled H2A.Z and components of its deposition complex, SWR1, before and after treating Arabidopsis thaliana with the hormone abscisic acid (ABA). Our results show that transcribed genes with TSS-enriched H2A.Z have high SWR1 binding at steady-state, indicating continuous replacement of H2A.Z, while silent genes with gene body H2A.Z show lower SWR1 binding. Surprisingly, upon ABA treatment, thousands of previously silent genes activate, coincident with recruitment of SWR1 and retention of gene body H2A.Z enrichment. We also found that the SWR1-interacting protein MBD9 is not required for SWR1 recruitment to activated genes. These results provide new insights into the relationship between H2A.Z and transcription and the mechanics of H2A.Z targeting to chromatin. ### Competing Interest Statement The authors have declared no competing interest.

Altered nutrient conditions can trigger massive transcriptional reprogramming in plants, leading to the activation and silencing of thousands of genes. To gain a deeper understanding of the phosphate starvation response and the relationships between transcriptional and epigenomic changes that occur during this reprogramming, we conducted a time-resolved analysis of transcriptome and chromatin alterations in root hair cells of Arabidopsis thaliana during phosphate (P) starvation and subsequent resupply. We found that 96 hours of P starvation causes induction or repression of thousands of transcripts, and most of these recover to pre-starvation levels within 4 hours of P resupply. Among the phosphate starvation-induced genes are many polycomb targets with high levels of H3K27me3 and histone variant H2A.Z. When induced, these genes often show increased H3K4me3 consistent with active transcription, but surprisingly minimal loss of H3K27me3 or H2A.Z. These results indicate that the removal of silencing marks is not a prerequisite for activation of these genes. Our data provide a cell type- and time-resolved resource for studying the dynamics of a systemic nutrient stress and recovery and suggest that our current understanding of the mechanisms for switching between silent and active transcriptional states is incomplete. ### Competing Interest Statement The authors have declared no competing interest.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Isolated nuclei provide access to early steps in gene regulation involving chromatin as well as transcript production and processing. Here, we describe transfer of the isolation of nuclei from tagged specific cell types (INTACT) to the monocot rice (Oryza sativa L.). The purification of biotinylated nuclei was redesigned by replacing the outer nuclear-envelope-targeting domain of the nuclear tagging fusion (NTF) protein with an outer nuclear-envelope-anchored domain. This modified NTF was combined with codon-optimized Escherichia coli BirA in a single T-DNA construct. We also developed inexpensive methods for INTACT, T-DNA insertion mapping, and profiling of the complete nuclear transcriptome, including a ribosomal RNA degradation procedure that minimizes pre-ribosomal RNA (pre-rRNA) transcripts. A high-resolution comparison of nuclear and steady-state poly(A)+ transcript populations of seedling root tips confirmed the capture of pre-messenger RNA (pre-mRNA) and exposed distinctions in diversity and abundance of the nuclear and total transcriptomes. This retooled INTACT can enable high-resolution monitoring of the nuclear transcriptome and chromatin in specific cell types of rice and other species.

Gene expression is dynamically regulated on many levels, including chromatin accessibility and transcription. In order to study these nuclear regulatory events, we describe our method to purify nuclei with Isolation of Nuclei in TAgged Cell Types (INTACT). As nuclear RNA is low in polyadenylated transcripts and conventional pulldown methods would not capture non-polyadenylated pre-mRNA, we also present our method to remove ribosomal RNA from the total nuclear RNA in preparation for nuclear RNA-Seq.

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.

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.