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.

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.

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.

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.

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.

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.

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.

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.

miR824-Regulated AGAMOUS-LIKE16 Contributes to Flowering Time Repression in Arabidopsis.

The timing of flowering is pivotal for maximizing reproductive success under fluctuating environmental conditions. Flowering time is tightly controlled by complex genetic networks that integrate endogenous and exogenous cues, such as light, temperature, photoperiod, and hormones. Here, we show that AGAMOUS-LIKE16 (AGL16) and its negative regulator microRNA824 (miR824) control flowering time in Arabidopsis thaliana. Knockout of AGL16 effectively accelerates flowering in nonvernalized Col-FRI, in which the floral inhibitor FLOWERING LOCUS C (FLC) is strongly expressed, but shows no effect if plants are vernalized or grown in short days. Alteration of AGL16 expression levels by manipulating miR824 abundance influences the timing of flowering quantitatively, depending on the expression level and number of functional FLC alleles. The effect of AGL16 is fully dependent on the presence of FLOWERING LOCUS T (FT). Further experiments show that AGL16 can interact directly with SHORT VEGETATIVE PHASE and indirectly with FLC, two proteins that form a complex to repress expression of FT. Our data reveal that miR824 and AGL16 modulate the extent of flowering time repression in a long-day photoperiod.

Development-related PcG target in the apex 4 controls leaf margin architecture in Arabidopsis thaliana.

In a reverse genetics screen based on a group of genes enriched for development-related Polycomb group targets in the apex (DPAs), we isolated DPA4 as a novel regulator of leaf margin shape. T-DNA insertion lines in the DPA4 locus display enhanced leaf margin serrations and enlarged petals, whereas overexpression of DPA4 results in smooth margins. DPA4 encodes a putative RAV (Related to ABI3/VP1) transcriptional repressor and is expressed in the lateral organ boundary region and in the sinus of leaf serrations. DPA4 expression domains overlap with those of the known leaf shape regulator CUP-SHAPED COTYLEDON 2 (CUC2) and we provide evidence that DPA4 negatively regulates CUC2 expression independently of MIR164A, an established regulator of CUC2. Taken together, the data suggest DPA4 as a newly identified player in the signalling network that controls leaf serrations in Arabidopsis thaliana.

PEP1 regulates perennial flowering in Arabis alpina.

Annual plants complete their life cycle in one year and initiate flowering only once, whereas perennials live for many years and flower repeatedly. How perennials undergo repeated cycles of vegetative growth and flowering that are synchronized to the changing seasons has not been extensively studied. Flowering is best understood in annual Arabidopsis thaliana, but many closely related species, such as Arabis alpina, are perennials. We identified the A. alpina mutant perpetual flowering 1 (pep1), and showed that PEP1 contributes to three perennial traits. It limits the duration of flowering, facilitating a return to vegetative development, prevents some branches from undergoing the floral transition allowing polycarpic growth habit, and confers a flowering response to winter temperatures that restricts flowering to spring. Here we show that PEP1 is the orthologue of the A. thaliana gene FLOWERING LOCUS C (FLC). The FLC transcription factor inhibits flowering until A. thaliana is exposed to winter temperatures, which trigger chromatin modifications that stably repress FLC transcription. In contrast, PEP1 is only transiently repressed by low temperatures, causing repeated seasonal cycles of repression and activation of PEP1 transcription that allow it to carry out functions characteristic of the cyclical life history of perennials. The patterns of chromatin modifications at FLC and PEP1 differ correlating with their distinct expression patterns. Thus we describe a critical mechanism by which flowering regulation differs between related perennial and annual species, and propose that differences in chromatin regulation contribute to this variation.

Tissue-specific expression of FLOWERING LOCUS T in Arabidopsis is maintained independently of polycomb group protein repression.

The Polycomb Group (PcG) pathway represses transcription through a mechanism conserved among plants and animals. PcG-mediated repression can determine spatial territories of gene expression, but it remains unclear whether PcG-mediated repression is a regulatory requirement for all targets. Here, we show the role of PcG proteins in the spatial regulation of FLOWERING LOCUS T (FT), a main activator of flowering in Arabidopsis thaliana exclusively expressed in the vasculature. Strikingly, the loss of PcG repression causes down-regulation of FT. In addition, our results show how the effect of PcG-mediated regulation differs for target genes and that, for FT expression, it relies primarily on tissue differentiation.

Measuring CO2 assimilation of Arabidopsis thaliana whole plants and seedlings.

Photosynthesis is an essential process in plants that synthesizes sugars used for growth and development, highlighting the importance of establishing robust methods to monitor photosynthetic activity. Infrared gas analysis (IRGA) can be used to track photosynthetic rates by measuring plant CO2 assimilation and release. Although much progress has been made in the development of IRGA technologies, challenges remain when using this technique on small herbaceous plants such as Arabidopsis thaliana. The use of whole plant chambers can overcome the difficulties associated with applying bulky leaf clamps to small delicate leaves. However, respiration from the roots and from soil-based microorganisms may skew these gas exchange measurements. Here, we present a simple method to efficiently perform IRGA on A. thaliana plants using a whole plant chamber that removes the confounding effects of respiration from roots and soil-based microorganisms from the measurements. We show that this method can be used to detect subtle changes in photosynthetic rates measured at different times of day, under different growth conditions, and between wild-type and plants with deficiencies in the photosynthetic machinery. Furthermore, we show that this method can be used to detect changes in photosynthetic rates even at very young developmental stages such as 10 d-old seedlings. This method contributes to the array of techniques currently used to perform IRGA on A. thaliana and can allow for the monitoring of photosynthetic rates of whole plants from young ages.

Genome-wide analysis of cis-regulatory divergence between species in the Arabidopsis genus.

Cis-regulatory DNA has been suspected to play a preeminent role in adaptive evolution, but understanding the role of cis-regulatory mutations in gene expression divergence first requires an accurate analysis of the functional differences associated with these regions. We analyzed allele-specific expression (ASE) in leaf and floral tissues of F1 interspecific hybrids generated between the two closely related species Arabidopsis thaliana and A. lyrata with a whole-genome SNP (single nucleotide polymorphism) tiling array. We observed 2,205 genes showing ASE pattern in at least one tissue. Nearly 90% of genes displaying ASE preferentially expressed the allele of A. lyrata. Genome-wide comparison of sequence divergence revealed that genes displaying ASE had a higher ratio of nonsynonymous to synonymous substitutions in coding regions. We further observe that the epigenetic landscape of histone methylation in A. thaliana genome associate with ASE. The asymmetry in the direction of allele-specific expression suggests interspecific differences in the efficiency of gene silencing in F1 hybrids.

Widespread interspecific divergence in cis-regulation of transposable elements in the Arabidopsis genus.

Transposable elements (TEs) are so abundant and variable that they count among the most important mutational sources in genomes. Nonetheless, little is known about the genetics of their variation in activity or silencing across closely related species. Here, we demonstrate that regulation of TE genes can differ dramatically between the two closely related Arabidopsis species A. thaliana and A. lyrata. In leaf and floral tissues of F1 interspecific hybrids, about 47% of TEs show allele-specific expression, with the A. lyrata copy being generally expressed at higher level. We confirm that TEs are generally expressed in A. lyrata but not in A. thaliana. Allele-specific differences in TE expression are associated with divergence in epigenetic modifications like DNA and histone methylation between species as well as with sequence divergence. Our data demonstrate that A. thaliana silences TEs much better than A. lyrata. For long terminal repeat retrotransposons, these differences are more pronounced for younger insertions. Interspecific differences in TE silencing may have a great impact on genome size changes.

Multiple loci and genetic interactions involving flowering time genes regulate stem branching among natural variants of Arabidopsis.

Shoot branching is a major determinant of plant architecture. Genetic variants for reduced stem branching in the axils of cauline leaves of Arabidopsis were found in some natural accessions and also at low frequency in the progeny of multiparent crosses. Detailed genetic analysis using segregating populations derived from backcrosses with the parental lines and bulked segregant analysis was used to identify the allelic variation controlling reduced stem branching. Eight quantitative trait loci (QTLs) contributing to natural variation for reduced stem branching were identified (REDUCED STEM BRANCHING 1-8 (RSB1-8)). Genetic analysis showed that RSB6 and RSB7, corresponding to flowering time genes FLOWERING LOCUS C (FLC) and FRIGIDA (FRI), epistatically regulate stem branching. Furthermore, FLOWERING LOCUS T (FT), which corresponds to RSB8 as demonstrated by fine-mapping, transgenic complementation and expression analysis, caused pleiotropic effects not only on flowering time, but, in the specific background of active FRI and FLC alleles, also on the RSB trait. The consequence of allelic variation only expressed in late-flowering genotypes revealed novel and thus far unsuspected roles of several genes well characterized for their roles in flowering time control.

Fast isogenic mapping-by-sequencing of ethyl methanesulfonate-induced mutant bulks.

Mapping-by-sequencing (or SHOREmapping) has revitalized the powerful concept of forward genetic screens in plants. However, as in conventional genetic mapping approaches, mapping-by-sequencing requires phenotyping of mapping populations established from crosses between two diverged accessions. In addition to the segregation of the focal phenotype, this introduces natural phenotypic variation, which can interfere with the recognition of quantitative phenotypes. Here, we demonstrate how mapping-by-sequencing and candidate gene identification can be performed within the same genetic background using only mutagen-induced changes as segregating markers. Using a previously unknown suppressor of mutants of like heterochromatin protein1 (lhp1), which in its functional form is involved in chromatin-mediated gene repression, we identified three closely linked ethyl methanesulfonate-induced changes as putative candidates. In order to assess allele frequency differences between such closely linked mutations, we introduced deep candidate resequencing using the new Ion Torrent Personal Genome Machine sequencing platform to our mutant identification pipeline and thereby reduced the number of causal candidate mutations to only one. Genetic analysis of two independent additional alleles confirmed that this mutation was causal for the suppression of lhp1.

cis-Regulatory elements and chromatin state coordinately control temporal and spatial expression of FLOWERING LOCUS T in Arabidopsis.

Flowering time of summer annual Arabidopsis thaliana accessions is largely determined by the timing of FLOWERING LOCUS T (FT) expression in the leaf vasculature. To understand the complex interplay between activating and repressive inputs controlling flowering through FT, cis-regulatory sequences of FT were identified in this study. A proximal and an approximately 5-kb upstream promoter region containing highly conserved sequence blocks were found to be essential for FT activation by CONSTANS (CO). Chromatin-associated protein complexes add another layer to FT regulation. In plants constitutively overexpressing CO, changes in chromatin status, such as a decrease in binding of LIKE HETEROCHROMATIN PROTEIN1 (LHP1) and increased acetylation of H3K9 and K14, were observed throughout the FT locus, although these changes appear to be a consequence of FT upregulation and not a prerequisite for activation. Binding of LHP1 was required to repress enhancer elements located between the CO-controlled regions. By contrast, the distal and proximal promoter sequences required for FT activation coincide with locally LHP1 and H3K27me3 depleted chromatin, indicating that chromatin status facilitates the accessibility of transcription factors to FT. Therefore, distant regulatory regions are required for FT transcription, reflecting the complexity of its control and differences in chromatin status delimit functionally important cis-regulatory regions.