Quality control and evaluation of plant epigenomics data.

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

Photoperiod-responsive changes in chromatin accessibility in phloem companion and epidermis cells of Arabidopsis leaves.

Photoperiod plays a key role in controlling the phase transition from vegetative to reproductive growth in flowering plants. Leaves are the major organs perceiving day-length signals, but how specific leaf cell types respond to photoperiod remains unknown. We integrated photoperiod-responsive chromatin accessibility and transcriptome data in leaf epidermis and vascular companion cells of Arabidopsis thaliana by combining isolation of nuclei tagged in specific cell/tissue types with assay for transposase-accessible chromatin using sequencing and RNA-sequencing. Despite a large overlap, vasculature and epidermis cells responded differently. Long-day predominantly induced accessible chromatin regions (ACRs); in the vasculature, more ACRs were induced and these were located at more distal gene regions, compared with the epidermis. Vascular ACRs induced by long days were highly enriched in binding sites for flowering-related transcription factors. Among the highly ranked genes (based on chromatin and expression signatures in the vasculature), we identified TREHALOSE-PHOSPHATASE/SYNTHASE 9 (TPS9) as a flowering activator, as shown by the late flowering phenotypes of T-DNA insertion mutants and transgenic lines with phloem-specific knockdown of TPS9. Our cell-type-specific analysis sheds light on how the long-day photoperiod stimulus impacts chromatin accessibility in a tissue-specific manner to regulate plant development.

Information integration and decision making in flowering time control.

In order to successfully reproduce, plants must sense changes in their environment and flower at the correct time. Many plants utilize day length and vernalization, a mechanism for verifying that winter has occurred, to determine when to flower. Our study used available temperature and day length data from different climates to provide a general understanding how this information processing of environmental signals could have evolved in plants. For climates where temperature fluctuation correlations decayed exponentially, a simple stochastic model characterizing vernalization was able to reconstruct the switch-like behavior of the core flowering regulatory genes. For these and other climates, artificial neural networks were used to predict flowering gene expression patterns. For temperate plants, long-term cold temperature and short-term day length measurements were sufficient to produce robust flowering time decisions from the neural networks. Additionally, evolutionary simulations on neural networks confirmed that the combined signal of temperature and day length achieved the highest fitness relative to neural networks with access to only one of those inputs. We suggest that winter temperature memory is a well-adapted strategy for plants' detection of seasonal changes, and absolute day length is useful for the subsequent triggering of flowering.

The domesticated transposase ALP2 mediates formation of a novel Polycomb protein complex by direct interaction with MSI1, a core subunit of Polycomb Repressive Complex 2 (PRC2).

A large fraction of plant genomes is composed of transposable elements (TE), which provide a potential source of novel genes through "domestication"-the process whereby the proteins encoded by TE diverge in sequence, lose their ability to catalyse transposition and instead acquire novel functions for their hosts. In Arabidopsis, ANTAGONIST OF LIKE HETEROCHROMATIN PROTEIN 1 (ALP1) arose by domestication of the nuclease component of Harbinger class TE and acquired a new function as a component of POLYCOMB REPRESSIVE COMPLEX 2 (PRC2), a histone H3K27me3 methyltransferase involved in regulation of host genes and in some cases TE. It was not clear how ALP1 associated with PRC2, nor what the functional consequence was. Here, we identify ALP2 genetically as a suppressor of Polycomb-group (PcG) mutant phenotypes and show that it arose from the second, DNA binding component of Harbinger transposases. Molecular analysis of PcG compromised backgrounds reveals that ALP genes oppose silencing and H3K27me3 deposition at key PcG target genes. Proteomic analysis reveals that ALP1 and ALP2 are components of a variant PRC2 complex that contains the four core components but lacks plant-specific accessory components such as the H3K27me3 reader LIKE HETEROCHROMATION PROTEIN 1 (LHP1). We show that the N-terminus of ALP2 interacts directly with ALP1, whereas the C-terminus of ALP2 interacts with MULTICOPY SUPPRESSOR OF IRA1 (MSI1), a core component of PRC2. Proteomic analysis reveals that in alp2 mutant backgrounds ALP1 protein no longer associates with PRC2, consistent with a role for ALP2 in recruitment of ALP1. We suggest that the propensity of Harbinger TE to insert in gene-rich regions of the genome, together with the modular two component nature of their transposases, has predisposed them for domestication and incorporation into chromatin modifying complexes.

Targeted DNA methylation represses two enhancers of FLOWERING LOCUS T in Arabidopsis thaliana.

FLOWERING LOCUS T (FT) plays a major role in regulating the floral transition in response to an inductive long day photoperiod in Arabidopsis thaliana. Expression of FT in leaves is dependent on the distal transcriptional enhancer Block C, located 5-kilobases (kb) upstream of the transcriptional start site (TSS). We expressed an inverted repeat of Block C to induce local DNA methylation and heterochromatin formation, which lead to FT downregulation in an inductive photoperiod. Using targeted DNA methylation as a tool to uncover further regulatory regions at the FT locus, we identified Block E, located 1 kb downstream of the gene, as a novel enhancer of FT. As Block C, Block E is conserved across Brassicaceae and located in accessible chromatin. Block C and E act as additive transcriptional enhancers that, in combination with the proximal FT promoter, control expression of FT in response to photoperiod in the leaf phloem.

AP2 transcription factor CBX1 with a specific function in symbiotic exchange of nutrients in mycorrhizal Lotus japonicus.

The arbuscular mycorrhizal (AM) symbiosis, a widespread mutualistic association between land plants and fungi, depends on reciprocal exchange of phosphorus driven by proton-coupled phosphate uptake into host plants and carbon supplied to AM fungi by host-dependent sugar and lipid biosynthesis. The molecular mechanisms and cis-regulatory modules underlying the control of phosphate uptake and de novo fatty acid synthesis in AM symbiosis are poorly understood. Here, we show that the AP2 family transcription factor CTTC MOTIF-BINDING TRANSCRIPTION FACTOR1 (CBX1), a WRINKLED1 (WRI1) homolog, directly binds the evolutionary conserved CTTC motif that is enriched in mycorrhiza-regulated genes and activates Lotus japonicus phosphate transporter 4 (LjPT4) in vivo and in vitro. Moreover, the mycorrhiza-inducible gene encoding H+-ATPase (LjHA1), implicated in energizing nutrient uptake at the symbiotic interface across the periarbuscular membrane, is coregulated with LjPT4 by CBX1. Accordingly, CBX1-defective mutants show reduced mycorrhizal colonization. Furthermore, genome-wide-binding profiles, DNA-binding studies, and heterologous expression reveal additional binding of CBX1 to AW box, the consensus DNA-binding motif for WRI1, that is enriched in promoters of glycolysis and fatty acid biosynthesis genes. We show that CBX1 activates expression of lipid metabolic genes including glycerol-3-phosphate acyltransferase RAM2 implicated in acylglycerol biosynthesis. Our finding defines the role of CBX1 as a regulator of host genes involved in phosphate uptake and lipid synthesis through binding to the CTTC/AW molecular module, and supports a model underlying bidirectional exchange of phosphorus and carbon, a fundamental trait in the mutualistic AM symbiosis.

Telobox motifs recruit CLF/SWN-PRC2 for H3K27me3 deposition via TRB factors in Arabidopsis.

Polycomb repressive complexes (PRCs) control organismic development in higher eukaryotes through epigenetic gene repression1-4. PRC proteins do not contain DNA-binding domains, thus prompting questions regarding how PRCs find their target loci 5 . Here we present genome-wide evidence of PRC2 recruitment by telomere-repeat-binding factors (TRBs) through telobox-related motifs in Arabidopsis. A triple trb1-2, trb2-1, and trb3-2 (trb1/2/3) mutant with a developmental phenotype and a transcriptome strikingly similar to those of strong PRC2 mutants showed redistribution of trimethyl histone H3 Lys27 (H3K27me3) marks and lower H3K27me3 levels, which were correlated with derepression of TRB1-target genes. TRB1-3 physically interacted with the PRC2 proteins CLF and SWN. A SEP3 reporter gene with a telobox mutation showed ectopic expression, which was correlated with H3K27me3 depletion, whereas tethering TRB1 to the mutated cis element partially restored repression. We propose that telobox-related motifs recruit PRC2 through the interaction between TRBs and CLF/SWN, a mechanism essential for H3K27me3 deposition at a subset of target genes.

Genome-wide mapping of transcriptional enhancer candidates using DNA and chromatin features in maize.

BACKGROUND: While most cells in multicellular organisms carry the same genetic information, in each cell type only a subset of genes is being transcribed. Such differentiation in gene expression depends, for a large part, on the activation and repression of regulatory sequences, including transcriptional enhancers. Transcriptional enhancers can be located tens of kilobases from their target genes, but display characteristic chromatin and DNA features, allowing their identification by genome-wide profiling. Here we show that integration of chromatin characteristics can be applied to predict distal enhancer candidates in Zea mays, thereby providing a basis for a better understanding of gene regulation in this important crop plant. RESULT: To predict transcriptional enhancers in the crop plant maize (Zea mays L. ssp. mays), we integrated available genome-wide DNA methylation data with newly generated maps for chromatin accessibility and histone 3 lysine 9 acetylation (H3K9ac) enrichment in young seedling and husk tissue. Approximately 1500 intergenic regions, displaying low DNA methylation, high chromatin accessibility and H3K9ac enrichment, were classified as enhancer candidates. Based on their chromatin profiles, candidate sequences can be classified into four subcategories. Tissue-specificity of enhancer candidates is defined based on the tissues in which they are identified and putative target genes are assigned based on tissue-specific expression patterns of flanking genes. CONCLUSIONS: Our method identifies three previously identified distal enhancers in maize, validating the new set of enhancer candidates and enlarging the toolbox for the functional characterization of gene regulation in the highly repetitive maize genome.

Ctf4-related protein recruits LHP1-PRC2 to maintain H3K27me3 levels in dividing cells in Arabidopsis thaliana.

Polycomb Repressive Complex (PRC) 2 catalyzes the H3K27me3 modification that warrants inheritance of a repressive chromatin structure during cell division, thereby assuring stable target gene repression in differentiated cells. It is still under investigation how H3K27me3 is passed on from maternal to filial strands during DNA replication; however, cell division can reinforce H3K27me3 coverage at target regions. To identify novel factors involved in the Polycomb pathway in plants, we performed a forward genetic screen for enhancers of the like heterochromatin protein 1 (lhp1) mutant, which shows relatively mild phenotypic alterations compared with other plant PRC mutants. We mapped enhancer of lhp1 (eol) 1 to a gene related to yeast Chromosome transmission fidelity 4 (Ctf4) based on phylogenetic analysis, structural similarities, physical interaction with the CMG helicase component SLD5, and an expression pattern confined to actively dividing cells. A combination of eol1 with the curly leaf (clf) allele, carrying a mutation in the catalytic core of PRC2, strongly enhanced the clf phenotype; furthermore, H3K27me3 coverage at target genes was strongly reduced in eol1 clf double mutants compared with clf single mutants. EOL1 physically interacted with CLF, its partially redundant paralog SWINGER (SWN), and LHP1. We propose that EOL1 interacts with LHP1-PRC2 complexes during replication and thereby participates in maintaining the H3K27me3 mark at target genes.

H2A monoubiquitination in Arabidopsis thaliana is generally independent of LHP1 and PRC2 activity.

BACKGROUND: Polycomb group complexes PRC1 and PRC2 repress gene expression at the chromatin level in eukaryotes. The classic recruitment model of Polycomb group complexes in which PRC2-mediated H3K27 trimethylation recruits PRC1 for H2A monoubiquitination was recently challenged by data showing that PRC1 activity can also recruit PRC2. However, the prevalence of these two mechanisms is unknown, especially in plants as H2AK121ub marks were examined at only a handful of Polycomb group targets. RESULTS: By using genome-wide analyses, we show that H2AK121ub marks are surprisingly widespread in Arabidopsis thaliana, often co-localizing with H3K27me3 but also occupying a set of transcriptionally active genes devoid of H3K27me3. Furthermore, by profiling H2AK121ub and H3K27me3 marks in atbmi1a/b/c, clf/swn, and lhp1 mutants we found that PRC2 activity is not required for H2AK121ub marking at most genes. In contrast, loss of AtBMI1 function impacts the incorporation of H3K27me3 marks at most Polycomb group targets. CONCLUSIONS: Our findings show the relationship between H2AK121ub and H3K27me3 marks across the A. thaliana genome and unveil that ubiquitination by PRC1 is largely independent of PRC2 activity in plants, while the inverse is true for H3K27 trimethylation.

The Arabidopsis Polycomb Repressive Complex 1 (PRC1) Components AtBMI1A, B, and C Impact Gene Networks throughout All Stages of Plant Development.

Polycomb Group regulation in Arabidopsis (Arabidopsis thaliana) is required to maintain cell differentiation and allow developmental phase transitions. This is achieved by the activity of three PcG repressive complex 2s (PRC2s) and the participation of a yet poorly defined PRC1. Previous results showed that apparent PRC1 components perform discrete roles during plant development, suggesting the existence of PRC1 variants; however, it is not clear in how many processes these components participate. We show that AtBMI1 proteins are required to promote all developmental phase transitions and to control cell proliferation during organ growth and development, expanding their proposed range of action. While AtBMI1 function during germination is closely linked to B3 domain transcription factors VAL1/2 possibly in combination with GT-box binding factors, other AtBMI1 regulatory networks require participation of different factor combinations. Conversely, EMF1 and LHP1 bind many H3K27me3 positive genes up-regulated in atbmi1a/b/c mutants; however, loss of their function affects expression of a different subset, suggesting that even if EMF1, LHP1, and AtBMI1 exist in a common PRC1 variant, their role in repression depends on the functional context.

Meta-analysis of Genome-Wide Chromatin Data.

Genome-wide analyses of chromatin factor-binding sites or histone modification localization generate lists of up to several thousand potential target genes. For many model organisms, large annotation databases are available to help with the characterization and classification of genomic datasets. The term meta-analysis has been coined for this type of multi-database comparison. In this chapter, we describe a workflow to perform a transcriptional and functional analysis of genome-wide target genes. Sources of transcription data and clustering tools to subdivide genes according to their expression pattern are described. For a functional analysis, we focus on the Gene Ontology (GO) vocabulary and methods to uncover over- or underrepresented functions among target genes. Genomic targets of the histone modification H3K27me3 are presented as a case study to demonstrate that meta-analysis can uncover functions that were hidden in genome-wide datasets.

The age of multiplexity: recruitment and interactions of Polycomb complexes in plants.

Polycomb group (PcG) proteins form distinct complexes that modify chromatin by histone H3 methylation and H2A mono-ubiquitination leading to chromatin compaction and epigenetic repression of target genes. A network of PcG protein complexes, associated partners and antagonistically acting chromatin modifiers is essential to regulate developmental transitions and cell fate in all multicellular eukaryotes. In this review, we discuss insights on the subfunctionalization of PcG complexes and their modes of recruitment to target sites based on data from the model organism Arabidopsis thaliana.

Complementary Activities of TELOMERE REPEAT BINDING Proteins and Polycomb Group Complexes in Transcriptional Regulation of Target Genes.

In multicellular organisms, Polycomb Repressive Complex 1 (PRC1) and PRC2 repress target genes through histone modification and chromatin compaction. Arabidopsis thaliana mutants strongly compromised in the pathway cannot develop differentiated organs. LIKE HETEROCHROMATIN PROTEIN1 (LHP1) is so far the only known plant PRC1 component that directly binds to H3K27me3, the histone modification set by PRC2, and also associates genome-wide with trimethylation of lysine 27 of histone H3 (H3K27me3). Surprisingly, lhp1 mutants show relatively mild phenotypic alterations. To explain this paradox, we screened for genetic enhancers of lhp1 mutants to identify novel components repressing target genes together with, or in parallel to, LHP1. Two enhancing mutations were mapped to TELOMERE REPEAT BINDING PROTEIN1 (TRB1) and its paralog TRB3. We show that TRB1 binds to thousands of genomic sites containing telobox or related cis-elements with a significant increase of sites and strength of binding in the lhp1 background. Furthermore, in combination with lhp1, but not alone, trb1 mutants show increased transcription of LHP1 targets, such as floral meristem identity genes, which are more likely to be bound by TRB1 in the lhp1 background. By contrast, expression of a subset of LHP1-independent TRB1 target genes, many involved in primary metabolism, is decreased in the absence of TRB1 alone. Thus, TRB1 is a bivalent transcriptional modulator that maintains downregulation of Polycomb Group (PcG) target genes in lhp1 mutants, while it sustains high expression of targets that are regulated independently of PcG.

Kicking against the PRCs - A Domesticated Transposase Antagonises Silencing Mediated by Polycomb Group Proteins and Is an Accessory Component of Polycomb Repressive Complex 2.

The Polycomb group (PcG) and trithorax group (trxG) genes play crucial roles in development by regulating expression of homeotic and other genes controlling cell fate. Both groups catalyse modifications of chromatin, particularly histone methylation, leading to epigenetic changes that affect gene activity. The trxG antagonizes the function of PcG genes by activating PcG target genes, and consequently trxG mutants suppress PcG mutant phenotypes. We previously identified the ANTAGONIST OF LIKE HETEROCHROMATIN PROTEIN1 (ALP1) gene as a genetic suppressor of mutants in the Arabidopsis PcG gene LIKE HETEROCHROMATIN PROTEIN1 (LHP1). Here, we show that ALP1 interacts genetically with several other PcG and trxG components and that it antagonizes PcG silencing. Transcriptional profiling reveals that when PcG activity is compromised numerous target genes are hyper-activated in seedlings and that in most cases this requires ALP1. Furthermore, when PcG activity is present ALP1 is needed for full activation of several floral homeotic genes that are repressed by the PcG. Strikingly, ALP1 does not encode a known chromatin protein but rather a protein related to PIF/Harbinger class transposases. Phylogenetic analysis indicates that ALP1 is broadly conserved in land plants and likely lost transposase activity and acquired a novel function during angiosperm evolution. Consistent with this, immunoprecipitation and mass spectrometry (IP-MS) show that ALP1 associates, in vivo, with core components of POLYCOMB REPRESSIVE COMPLEX 2 (PRC2), a widely conserved PcG protein complex which functions as a H3K27me3 histone methyltransferase. Furthermore, in reciprocal pulldowns using the histone methyltransferase CURLY LEAF (CLF), we identify not only ALP1 and the core PRC2 components but also plant-specific accessory components including EMBRYONIC FLOWER 1 (EMF1), a transcriptional repressor previously associated with PRC1-like complexes. Taken together our data suggest that ALP1 inhibits PcG silencing by blocking the interaction of the core PRC2 with accessory components that promote its HMTase activity or its role in inhibiting transcription. ALP1 is the first example of a domesticated transposase acquiring a novel function as a PcG component. The antagonistic interaction of a modified transposase with the PcG machinery is novel and may have arisen as a means for the cognate transposon to evade host surveillance or for the host to exploit features of the transposition machinery beneficial for epigenetic regulation of gene activity.

Molecular memories in the regulation of seasonal flowering: from competence to cessation.

Plants commit to flowering based on endogenous and exogenous information that they can remember across mitotic cell divisions. Here, we review how signal perception and epigenetic memory converge at key integrator genes, and we show how variation in their regulatory circuits supports the diversity of plant lifestyles.

The arabidopsis DNA polymerase δ has a role in the deposition of transcriptionally active epigenetic marks, development and flowering.

DNA replication is a key process in living organisms. DNA polymerase α (Polα) initiates strand synthesis, which is performed by Polε and Polδ in leading and lagging strands, respectively. Whereas loss of DNA polymerase activity is incompatible with life, viable mutants of Polα and Polε were isolated, allowing the identification of their functions beyond DNA replication. In contrast, no viable mutants in the Polδ polymerase-domain were reported in multicellular organisms. Here we identify such a mutant which is also thermosensitive. Mutant plants were unable to complete development at 28°C, looked normal at 18°C, but displayed increased expression of DNA replication-stress marker genes, homologous recombination and lysine 4 histone 3 trimethylation at the SEPALLATA3 (SEP3) locus at 24°C, which correlated with ectopic expression of SEP3. Surprisingly, high expression of SEP3 in vascular tissue promoted FLOWERING LOCUS T (FT) expression, forming a positive feedback loop with SEP3 and leading to early flowering and curly leaves phenotypes. These results strongly suggest that the DNA polymerase δ is required for the proper establishment of transcriptionally active epigenetic marks and that its failure might affect development by affecting the epigenetic control of master genes.

Induced and natural variation of promoter length modulates the photoperiodic response of FLOWERING LOCUS T.

FLOWERING LOCUS T (FT) regulates the floral transition in many plant species by integrating environmental seasonal signals and internal cues. Here we show that two interdependent regulatory regions are necessary and sufficient to convey photoperiod responsiveness to FT. While a minimal distance between the regulatory regions is required to fully suppress FT expression under short days, increased distance reduces promoter response to long days. Natural variation at FT creating promoter length differences is widespread, correlates with longitudinal and latitudinal clines and affects a promoter region physically interacting with both photoperiod control regions. Three major FT promoter variants correlate with differences in FT allele usage in F1 hybrids. We propose that FT variation in cis could be adaptive by conferring differences in FT transcriptional control ultimately translating to increased fitness.