Polycomb-dependent differential chromatin compartmentalization determines gene coregulation in Arabidopsis.

In animals, distant H3K27me3-marked Polycomb targets can establish physical interactions forming repressive chromatin hubs. In plants, growing evidence suggests that H3K27me3 acts directly or indirectly to regulate chromatin interactions, although how this histone modification modulates 3D chromatin architecture remains elusive. To decipher the impact of the dynamic deposition of H3K27me3 on the Arabidopsis thaliana nuclear interactome, we combined genetics, transcriptomics, and several 3D epigenomic approaches. By analyzing mutants defective for histone H3K27 methylation or demethylation, we uncovered the crucial role of this chromatin mark in short- and previously unnoticed long-range chromatin loop formation. We found that a reduction in H3K27me3 levels led to a decrease in the interactions within Polycomb-associated repressive domains. Regions with lower H3K27me3 levels in the H3K27 methyltransferase clf mutant established new interactions with regions marked with H3K9ac, a histone modification associated with active transcription, indicating that a reduction in H3K27me3 levels induces a global reconfiguration of chromatin architecture. Altogether, our results reveal that the 3D genome organization is tightly linked to reversible histone modifications that govern chromatin interactions. Consequently, nuclear organization dynamics shapes the transcriptional reprogramming during plant development and places H3K27me3 as a key feature in the coregulation of distant genes.

Hypometabolism to survive the long polar night and subsequent successful return to light in the diatom Fragilariopsis cylindrus.

Diatoms, the main eukaryotic phytoplankton of the polar marine regions, are essential for the maintenance of food chains specific to Arctic and Antarctic ecosystems, and are experiencing major disturbances under current climate change. As such, it is fundamental to understand the physiological mechanisms and associated molecular basis of their endurance during the long polar night. Here, using the polar diatom Fragilariopsis cylindrus, we report an integrative analysis combining transcriptomic, microscopic and biochemical approaches to shed light on the strategies used to survive the polar night. We reveal that in prolonged darkness, diatom cells enter a state of quiescence with reduced metabolic and transcriptional activity, during which no cell division occurs. We propose that minimal energy is provided by respiration and degradation of protein, carbohydrate and lipid stores and that homeostasis is maintained by autophagy in prolonged darkness. We also report internal structural changes that manifest the morphological acclimation of cells to darkness, including the appearance of a large vacuole. Our results further show that immediately following a return to light, diatom cells are able to use photoprotective mechanisms and rapidly resume photosynthesis, demonstrating the remarkable robustness of polar diatoms to prolonged darkness at low temperature.

Immunity onset alters plant chromatin and utilizes EDA16 to regulate oxidative homeostasis.

Perception of microbes by plants leads to dynamic reprogramming of the transcriptome, which is essential for plant health. The appropriate amplitude of this transcriptional response can be regulated at multiple levels, including chromatin. However, the mechanisms underlying the interplay between chromatin remodeling and transcription dynamics upon activation of plant immunity remain poorly understood. Here, we present evidence that activation of plant immunity by bacteria leads to nucleosome repositioning, which correlates with altered transcription. Nucleosome remodeling follows distinct patterns of nucleosome repositioning at different loci. Using a reverse genetic screen, we identify multiple chromatin remodeling ATPases with previously undescribed roles in immunity, including EMBRYO SAC DEVELOPMENT ARREST 16, EDA16. Functional characterization of the immune-inducible chromatin remodeling ATPase EDA16 revealed a mechanism to negatively regulate immunity activation and limit changes in redox homeostasis. Our transcriptomic data combined with MNase-seq data for EDA16 functional knock-out and over-expressor mutants show that EDA16 selectively regulates a defined subset of genes involved in redox signaling through nucleosome repositioning. Thus, collectively, chromatin remodeling ATPases fine-tune immune responses and provide a previously uncharacterized mechanism of immune regulation.

Arabidopsis Topless-related 1 mitigates physiological damage and growth penalties of induced immunity.

Transcriptional corepressors of the Topless (TPL) family regulate plant hormone and immunity signaling. The lack of a genome-wide profile of their chromatin associations limits understanding of the TPL family roles in transcriptional regulation. Chromatin immunoprecipitation with sequencing (ChIP-Seq) was performed on Arabidopsis thaliana lines expressing GFP-tagged Topless-related 1 (TPR1-GFP) with and without constitutive immunity via Enhanced Disease Susceptibility 1 (EDS1). RNA-Seq profiling of the TPR1-GFP lines and pathogen-infected tpl/tpr mutants, combined with measuring immunity, growth, and physiological parameters was employed to investigate TPL/TPR roles in immunity and defense homeostasis. TPR1 was enriched at promoter regions of c. 1400 genes and c. 10% of the detected binding required EDS1 immunity signaling. In a tpr1 tpl tpr4 (t3) mutant, resistance to bacteria was slightly compromised, and defense-related transcriptional reprogramming was weakly reduced or enhanced, respectively, at early (< 1 h) and late 24 h stages of bacterial infection. The t3 plants challenged with bacteria or pathogen-associated molecular pattern nlp24 displayed photosystem II dysfunctions. Also, t3 plants were hypersensitive to phytocytokine pep1 at the level of root growth inhibition. Transgenic expression of TPR1 rescued these t3 physiological defects. We propose that TPR1 and TPL family proteins function in Arabidopsis to reduce detrimental effects associated with activated transcriptional immunity.

GCN5 modulates salicylic acid homeostasis by regulating H3K14ac levels at the 5' and 3' ends of its target genes.

The modification of histones by acetyl groups has a key role in the regulation of chromatin structure and transcription. The Arabidopsis thaliana histone acetyltransferase GCN5 regulates histone modifications as part of the Spt-Ada-Gcn5 Acetyltransferase (SAGA) transcriptional coactivator complex. GCN5 was previously shown to acetylate lysine 14 of histone 3 (H3K14ac) in the promoter regions of its target genes even though GCN5 binding did not systematically correlate with gene activation. Here, we explored the mechanism through which GCN5 controls transcription. First, we fine-mapped its GCN5 binding sites genome-wide and then used several global methodologies (ATAC-seq, ChIP-seq and RNA-seq) to assess the effect of GCN5 loss-of-function on the expression and epigenetic regulation of its target genes. These analyses provided evidence that GCN5 has a dual role in the regulation of H3K14ac levels in their 5' and 3' ends of its target genes. While the gcn5 mutation led to a genome-wide decrease of H3K14ac in the 5' end of the GCN5 down-regulated targets, it also led to an increase of H3K14ac in the 3' ends of GCN5 up-regulated targets. Furthermore, genome-wide changes in H3K14ac levels in the gcn5 mutant correlated with changes in H3K9ac at both 5' and 3' ends, providing evidence for a molecular link between the depositions of these two histone modifications. To understand the biological relevance of these regulations, we showed that GCN5 participates in the responses to biotic stress by repressing salicylic acid (SA) accumulation and SA-mediated immunity, highlighting the role of this protein in the regulation of the crosstalk between diverse developmental and stress-responsive physiological programs. Hence, our results demonstrate that GCN5, through the modulation of H3K14ac levels on its targets, controls the balance between biotic and abiotic stress responses and is a master regulator of plant-environmental interactions.

Arabidopsis DNA Replication Initiates in Intergenic, AT-Rich Open Chromatin.

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2'-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips.

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Genome-Wide Analysis of the Arabidopsis Replication Timing Program.

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

Repliscan: a tool for classifying replication timing regions.

BACKGROUND: Replication timing experiments that use label incorporation and high throughput sequencing produce peaked data similar to ChIP-Seq experiments. However, the differences in experimental design, coverage density, and possible results make traditional ChIP-Seq analysis methods inappropriate for use with replication timing. RESULTS: To accurately detect and classify regions of replication across the genome, we present Repliscan. Repliscan robustly normalizes, automatically removes outlying and uninformative data points, and classifies Repli-seq signals into discrete combinations of replication signatures. The quality control steps and self-fitting methods make Repliscan generally applicable and more robust than previous methods that classify regions based on thresholds. CONCLUSIONS: Repliscan is simple and effective to use on organisms with different genome sizes. Even with analysis window sizes as small as 1 kilobase, reliable profiles can be generated with as little as 2.4x coverage.

Oxidative DNA damage bypass in Arabidopsis thaliana requires DNA polymerase λ and proliferating cell nuclear antigen 2.

The oxidized base 7,8-oxoguanine (8-oxo-G) is the most common DNA lesion generated by reactive oxygen species. This lesion is highly mutagenic due to the frequent misincorporation of A opposite 8-oxo-G during DNA replication. In mammalian cells, the DNA polymerase (pol) family X enzyme DNA pol λ catalyzes the correct incorporation of C opposite 8-oxo-G, together with the auxiliary factor proliferating cell nuclear antigen (PCNA). Here, we show that Arabidopsis thaliana DNA pol λ, the only member of the X family in plants, is as efficient in performing error-free translesion synthesis past 8-oxo-G as its mammalian homolog. Arabidopsis, in contrast with animal cells, possesses two genes for PCNA. Using in vitro and in vivo approaches, we observed that PCNA2, but not PCNA1, physically interacts with DNA pol λ, enhancing its fidelity and efficiency in translesion synthesis. The levels of DNA pol λ in transgenic plantlets characterized by overexpression or silencing of Arabidopsis POLL correlate with the ability of cell extracts to perform error-free translesion synthesis. The important role of DNA pol λ is corroborated by the observation that the promoter of POLL is activated by UV and that both overexpressing and silenced plants show altered growth phenotypes.

Histone H1 protects telomeric repeats from H3K27me3 invasion in Arabidopsis.

While the pivotal role of linker histone H1 in shaping nucleosome organization is well established, its functional interplays with chromatin factors along the epigenome are just starting to emerge. Here we show that, in Arabidopsis, as in mammals, H1 occupies Polycomb Repressive Complex 2 (PRC2) target genes where it favors chromatin condensation and H3K27me3 deposition. We further show that, contrasting with its conserved function in PRC2 activation at genes, H1 selectively prevents H3K27me3 accumulation at telomeres and large pericentromeric interstitial telomeric repeat (ITR) domains by restricting DNA accessibility to Telomere Repeat Binding (TRB) proteins, a group of H1-related Myb factors mediating PRC2 cis recruitment. This study provides a mechanistic framework by which H1 avoids the formation of gigantic H3K27me3-rich domains at telomeric sequences and contributes to safeguard nucleus architecture.

Functional analysis of the degradation of cellulosic substrates by a Chaetomium globosum endophytic isolate.

Most photosynthetically fixed carbon is contained in cell wall polymers present in plant biomasses, the largest organic carbon source in the biosphere. The degradation of these polymers for biotechnological purposes requires the combined action of several enzymes. To identify new activities, we examined which enzymes are activated by an endophytic strain of Chaetomium globosum to degrade cellulose-containing substrates. After biochemical analyses of the secretome of the fungus grown on cellulose or woody substrates, we took advantage of the available genomic data to identify potentially involved genes. After in silico identification of putative genes encoding either proteins able to bind to cellulose or glycohydrolases (GHs) of family 7, we investigated their transcript levels by reverse transcription-quantitative PCR (RT-qPCR). Our data suggest that eight genes compose the core of the cellulose-degrading system of C. globosum. Notably, the related enzymes belong structurally to the well-described GH families 5, 6, 7, 16, and 45, which are known to be the core of the cellulose degradation systems of several ascomycetes. The high expression levels of cellobiose dehydrogenase and two GH 61 enzymes suggest the involvement of this oxidoreductive synergic system in C. globosum. Transcript analysis along with relevant coding sequence (CDS) isolation and expression of recombinant proteins proved to be a key strategy for the determination of the features of two endoglucanases used by C. globosum for the first attack of crystalline cellulose. Finally, the possible involvement of transcriptional regulators described for other ascomycetes is discussed.

Molecular characterization of a Cyrtochilum loxense Somatic Embryogenesis Receptor-like Kinase (SERK) gene expressed during somatic embryogenesis.

Somatic embryogenesis is crucial for the propagation of endangered Ecuadorian orchid species, among them Cyrtochilum loxense, in view of the fact that their number in nature or in collections is quite reduced. One of the genes expressed during somatic and zygotic embryogenesis is Somatic Embryogenesis Receptor-like Kinase (SERK). Despite the development of somatic embryogenesis protocols for orchids, no SERK genes have been isolated from this family. This is the first report on the isolation of a full-length orchid SERK sequence, namely that of Cyrtochilum loxense (ClSERK). The identity of ClSERK was inferred by the presence of all domains typical of SERK proteins: a signal peptide, a leucine zipper domain, five LRRs, a serine proline-rich domain, a transmembrane domain, a kinase domain, and the C-terminal region. We have observed that the ClSERK gene is highly expressed in embryogenic calluses generated from protocorms at the time of appearance of embryonic morphological features. At later stages when embryos become well visible on calluses, ClSERK gene expression decreases. Compared to early stages of embryo formation on calluses, the expression detected in leaf tissue is far lower, thus suggesting a role of this gene during development.

New insights into homoeologous copy number variations in the hexaploid wheat genome.

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Wheat chromatin architecture is organized in genome territories and transcription factories.

BACKGROUND: Polyploidy is ubiquitous in eukaryotic plant and fungal lineages, and it leads to the co-existence of several copies of similar or related genomes in one nucleus. In plants, polyploidy is considered a major factor in successful domestication. However, polyploidy challenges chromosome folding architecture in the nucleus to establish functional structures. RESULTS: We examine the hexaploid wheat nuclear architecture by integrating RNA-seq, ChIP-seq, ATAC-seq, Hi-C, and Hi-ChIP data. Our results highlight the presence of three levels of large-scale spatial organization: the arrangement into genome territories, the diametrical separation between facultative and constitutive heterochromatin, and the organization of RNA polymerase II around transcription factories. We demonstrate the micro-compartmentalization of transcriptionally active genes determined by physical interactions between genes with specific euchromatic histone modifications. Both intra- and interchromosomal RNA polymerase-associated contacts involve multiple genes displaying similar expression levels. CONCLUSIONS: Our results provide new insights into the physical chromosome organization of a polyploid genome, as well as on the relationship between epigenetic marks and chromosome conformation to determine a 3D spatial organization of gene expression, a key factor governing gene transcription in polyploids.

Arabidopsis S2Lb links AtCOMPASS-like and SDG2 activity in H3K4me3 independently from histone H2B monoubiquitination.

BACKGROUND: The functional determinants of H3K4me3, their potential dependency on histone H2B monoubiquitination, and their contribution to defining transcriptional regimes are poorly defined in plant systems. Unlike in Saccharomyces cerevisiae, where a single SET1 protein catalyzes H3K4me3 as part of COMPlex of proteins ASsociated with Set1 (COMPASS), in Arabidopsis thaliana, this activity involves multiple histone methyltransferases. Among these, the plant-specific SET DOMAIN GROUP 2 (SDG2) has a prominent role. RESULTS: We report that SDG2 co-regulates hundreds of genes with SWD2-like b (S2Lb), a plant ortholog of the Swd2 axillary subunit of yeast COMPASS. We show that S2Lb co-purifies with the AtCOMPASS core subunit WDR5, and both S2Lb and SDG2 directly influence H3K4me3 enrichment over highly transcribed genes. S2Lb knockout triggers pleiotropic developmental phenotypes at the vegetative and reproductive stages, including reduced fertility and seed dormancy. However, s2lb seedlings display little transcriptomic defects as compared to the large repertoire of genes targeted by S2Lb, SDG2, or H3K4me3, suggesting that H3K4me3 enrichment is important for optimal gene induction during cellular transitions rather than for determining on/off transcriptional status. Moreover, unlike in budding yeast, most of the S2Lb and H3K4me3 genomic distribution does not rely on a trans-histone crosstalk with histone H2B monoubiquitination. CONCLUSIONS: Collectively, this study unveils that the evolutionarily conserved COMPASS-like complex has been co-opted by the plant-specific SDG2 histone methyltransferase and mediates H3K4me3 deposition through an H2B monoubiquitination-independent pathway in Arabidopsis.

Isolation of Plant Nuclei at Defined Cell Cycle Stages Using EdU Labeling and Flow Cytometry.

5-Ethynyl-2'-deoxyuridine (EdU) is a nucleoside analog of thymidine that can be rapidly incorporated into replicating DNA in vivo and, subsequently, detected by using "click" chemistry to couple its terminal alkyne group to fluorescent azides such as Alexa Fluor 488. Recently, EdU incorporation followed by coupling with a fluorophore has been used to visualize newly synthesized DNA in a wide range of plant species. One particularly useful application is in flow cytometry, where two-parameter sorting can be employed to analyze different phases of the cell cycle, as defined both by total DNA content and the amount of EdU pulse-labeled DNA. This approach allows analysis of the cell cycle without the need for synchronous cell populations, which can be difficult to obtain in many plant systems. The approach presented here, which was developed for fixed, EdU-labeled nuclei, can be used to prepare analytical profiles as well as to make highly purified preparations of G1, S, or G2/M phase nuclei for molecular or biochemical analysis. We present protocols for EdU pulse labeling, tissue fixation and harvesting, nuclei preparation, and flow sorting. Although developed for Arabidopsis suspension cells and maize root tips, these protocols should be modifiable to many other plant systems.