Search

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants., (© 2023 Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: Genetic variation in regulatory sequences that alter transcription factor (TF) binding is a major cause of phenotypic diversity. Brassinosteroid is a growth hormone that has major effects on plant phenotypes. Genetic variation in brassinosteroid-responsive cis-elements likely contributes to trait variation. Pinpointing such regulatory variations and quantitative genomic analysis of the variation in TF-target binding, however, remains challenging. How variation in transcriptional targets of signaling pathways such as the brassinosteroid pathway contributes to phenotypic variation is an important question to be investigated with innovative approaches., Results: Here, we use a hybrid allele-specific chromatin binding sequencing (HASCh-seq) approach and identify variations in target binding of the brassinosteroid-responsive TF ZmBZR1 in maize. HASCh-seq in the B73xMo17 F1s identifies thousands of target genes of ZmBZR1. Allele-specific ZmBZR1 binding (ASB) has been observed for 18.3% of target genes and is enriched in promoter and enhancer regions. About a quarter of the ASB sites correlate with sequence variation in BZR1-binding motifs and another quarter correlate with haplotype-specific DNA methylation, suggesting that both genetic and epigenetic variations contribute to the high level of variation in ZmBZR1 occupancy. Comparison with GWAS data shows linkage of hundreds of ASB loci to important yield and disease-related traits., Conclusion: Our study provides a robust method for analyzing genome-wide variations of TF occupancy and identifies genetic and epigenetic variations of the brassinosteroid response transcription network in maize., (© 2023. The Author(s).)

The Heat Shock Factor (HSF) transcription factor family is a central and required component of plant heat stress responses and acquired thermotolerance. The HSF family has dramatically expanded in plant lineages, often including a repertoire of 20 or more genes. Here we assess and compare the composition, heat responsiveness, and chromatin profiles of the HSF families in maize and Setaria viridis ( Setaria ), two model C4 panicoid grasses. Both species encode a similar number of HSFs, and examples of both conserved and variable expression responses to a heat stress event were observed between the two species. Chromatin accessibility and genome-wide DNA-binding profiles were generated to assess the chromatin of HSF family members with distinct responses to heat stress. We observed significant variability for both chromatin accessibility and promoter occupancy within similarly regulated sets of HSFs between Setaria and maize, as well as between syntenic pairs of maize HSFs retained following its most recent genome duplication event. Additionally, we observed the widespread presence of TF binding at HSF promoters in control conditions, even at HSFs that are only expressed in response to heat stress. TF-binding peaks were typically near putative HSF-binding sites in HSFs upregulated in response to heat stress, but not in stable or not expressed HSFs. These observations collectively support a complex scenario of expansion and subfunctionalization within this transcription factor family and suggest that within-family HSF transcriptional regulation is a conserved, defining feature of the family., Competing Interests: The authors declare that they have no conflicts of interest., (© 2023 The Authors. Plant Direct published by American Society of Plant Biologists and the Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: Many plant species exhibit genetic variation for coping with environmental stress. However, there are still limited approaches to effectively uncover the genomic region that regulates distinct responsive patterns of the gene across multiple varieties within the same species under abiotic stress., Results: By analyzing the transcriptomes of more than 100 maize inbreds, we reveal many cis- and trans-acting eQTLs that influence the expression response to heat stress. The cis-acting eQTLs in response to heat stress are identified in genes with differential responses to heat stress between genotypes as well as genes that are only expressed under heat stress. The cis-acting variants for heat stress-responsive expression likely result from distinct promoter activities, and the differential heat responses of the alleles are confirmed for selected genes using transient expression assays. Global footprinting of transcription factor binding is performed in control and heat stress conditions to document regions with heat-enriched transcription factor binding occupancies., Conclusions: Footprints enriched near proximal regions of characterized heat-responsive genes in a large association panel can be utilized for prioritizing functional genomic regions that regulate genotype-specific responses under heat stress., (© 2022. The Author(s).)

Main Conclusion: The Arabidopsis DORNRÖSCHEN-LIKE enhancer 2 comprises a high-occupancy target region in the IM periphery that integrates signals for the spiral phyllotactic pattern and cruciferous arrangement of sepals. Transcription of the DORNRÖSCHEN-LIKE (DRNL) gene marks lateral organ founder cells (LOFCs) in the peripheral zone of the inflorescence meristem (IM) and enhancer 2 (En2) in the DRNL promoter upstream region essentially contributes to this phyllotactic transcription pattern. Further analysis focused on the phylogenetically highly conserved 100-bp En2 core element, which was sufficient to promote the phyllotactic pattern, but was recalcitrant to further shortening. Here, we show that En2 core functions independent of orientation and create a series of mutations to study consequences on the transcription pattern. Their analysis shows that, first, in addition to in the inflorescence apex, En2 core acts in the embryo; second, cis-regulatory target sequences are distributed throughout the 100-bp element, although substantial differences exist in their function between embryo and IM. Third, putative core auxin response elements (AuxREs) spatially activate or restrict DRNL expression, and fourth, according to chromatin configuration data, En2 core enhancer activity in LOFCs correlates with an open chromatin structure at the DRNL transcription start. In combination, mutational and chromatin analyses imply that En2 core comprises a high-occupancy target (HOT) region for transcription factors, which implements phyllotactic information for the spiral LOFC pattern in the IM periphery and coordinates the cruciferous array of floral sepals. Our data disfavor a contribution of activating auxin response factors (ARFs) but do not exclude auxin as a morphogenetic signal.

In June 2019, more than a hundred plant researchers met in Cologne, Germany, for the 6 th European Workshop on Plant Chromatin (EWPC). This conference brought together a highly dynamic community of researchers with the common aim to understand how chromatin organization controls gene expression, development, and plant responses to the environment. New evidence showing how epigenetic states are set, perpetuated, and inherited were presented, and novel data related to the three-dimensional organization of chromatin within the nucleus were discussed. At the level of the nucleosome, its composition by different histone variants and their specialized histone deposition complexes were addressed as well as the mechanisms involved in histone post-translational modifications and their role in gene expression. The keynote lecture on plant DNA methylation by Julie Law (SALK Institute) and the tribute session to Lars Hennig, honoring the memory of one of the founders of the EWPC who contributed to promote the plant chromatin and epigenetic field in Europe, added a very special note to this gathering. In this perspective article we summarize some of the most outstanding data and advances on plant chromatin research presented at this workshop., (Copyright © 2020 Moreno-Romero, Probst, Trindade, Kalyanikrishna, Engelhorn and Farrona.)

Background: GMZ2 is a recombinant malaria vaccine inducing immune responses against Plasmodium falciparum (Pf) merozoite surface protein-3 and glutamate-rich protein. We used standardized controlled human malaria infection (CHMI) to assess the efficacy of this asexual blood-stage vaccine., Methods: We vaccinated 50 healthy, adult volunteers with lifelong exposure to Pf 3 times, at 4-week intervals, with 30 or 100 µg GMZ2 formulated in CAF01, a liposome-based adjuvant; 100 µg GMZ2, formulated in Alhydrogel; or a control vaccine (Verorab). Approximately 13 weeks after the last vaccination, 35/50 volunteers underwent CHMI by direct venous inoculation of 3200 Pf sporozoites (Sanaria® PfSPZ Challenge)., Results: Adverse events were similarly distributed between GMZ2 and control vaccinees. Baseline-corrected anti-GMZ2 antibody concentrations 4 weeks after the last vaccination were higher in all 3 GMZ2-vaccinated arms, compared to the control group. All GMZ2 formulations induced similar antibody levels. CHMI resulted in 29/34 (85%) volunteers with Pf parasitemia and 15/34 (44%) with malaria (parasitemia and symptoms). The proportion of participants with malaria (2/5 control, 6/10 GMZ2-Alhydrogel, 2/8 30 µg GMZ2-CAF01, and 5/11 100 µg GMZ2-CAF01) and the time it took them to develop malaria were similar in all groups. Baseline, vaccine-specific antibody concentrations were associated with protection against malaria., Conclusions: GMZ2 is well tolerated and immunogenic in lifelong-Pf-exposed adults from Gabon, with similar antibody responses regardless of formulation. CHMI showed no protective effect of prior vaccination with GMZ2, although baseline, vaccine-specific antibody concentrations were associated with protection. CHMI with the PfSPZ Challenge is a potent new tool to validate asexual, blood-stage malaria vaccines in Africa., Clinical Trials Registration: Pan-African Clinical Trials: PACTR201503001038304., (© The Author(s) 2018. Published by Oxford University Press for the Infectious Diseases Society of America. All rights reserved. For permissions, e-mail: journals.permissions@oup.com.)

Fluorescence-activated cell sorting (FACS) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were combined to analyse the chromatin state of lateral organ founder cells (LOFCs) in the peripheral zone of the Arabidopsis apetala1-1 cauliflower-1 double mutant inflorescence meristem. On a genome-wide level, we observed a striking correlation between transposase hypersensitive sites (THSs) detected by ATAC-seq and DNase I hypersensitive sites (DHSs). The mostly expanded DHSs were often substructured into several individual THSs, which correlated with phylogenetically conserved DNA sequences or enhancer elements. Comparing chromatin accessibility with available RNA-seq data, THS change configuration was reflected by gene activation or repression and chromatin regions acquired or lost transposase accessibility in direct correlation with gene expression levels in LOFCs. This was most pronounced immediately upstream of the transcription start, where genome-wide THSs were abundant in a complementary pattern to established H3K4me3 activation or H3K27me3 repression marks. At this resolution, the combined application of FACS/ATAC-seq is widely applicable to detect chromatin changes during cell-type specification and facilitates the detection of regulatory elements in plant promoters., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Experimental Biology.)

Enhancers are critical for developmental stage-specific gene expression, but their dynamic regulation in plants remains poorly understood. Here we compare genome-wide localization of H3K27ac, chromatin accessibility and transcriptomic changes during flower development in Arabidopsis. H3K27ac prevalently marks promoter-proximal regions, suggesting that H3K27ac is not a hallmark for enhancers in Arabidopsis. We provide computational and experimental evidence to confirm that distal DNase І hypersensitive sites are predictive of enhancers. The predicted enhancers are highly stage-specific across flower development, significantly associated with SNPs for flowering-related phenotypes, and conserved across crucifer species. Through the integration of genome-wide transcription factor (TF) binding datasets, we find that floral master regulators and stage-specific TFs are largely enriched at developmentally dynamic enhancers. Finally, we show that enhancer clusters and intronic enhancers significantly associate with stage-specific gene regulation by floral master TFs. Our study provides insights into the functional flexibility of enhancers during plant development, as well as hints to annotate plant enhancers.

Targeted changes in chromatin state at thousands of genes are central to eukaryotic development. RELATIVE OF EARLY FLOWERING 6 (REF6) is a Jumonji-type histone demethylase that counteracts Polycomb repressive complex 2 (PRC2)-mediated gene silencing in plants and was reported to select its binding sites in a direct, sequence-specific manner 1-3 . Here we show that REF6 and its two close paralogues determine spatial 'boundaries' of the repressive histone H3K27me3 mark in the genome and control the tissue-specific release from PRC2-mediated gene repression. Targeted mutagenesis revealed that these histone demethylases display pleiotropic, redundant functions in plant development, several of which depend on trans factor-mediated recruitment. Thus, Jumonji-type histone demethylases restrict repressive chromatin domains and contribute to tissue-specific gene activation via complementary targeting mechanisms.

Chromatin modifiers and remodelers are involved in generating dynamic changes at the chromatin, which allow differential and specific readouts of the genome. While genetic evidence indicates that several chromatin factors play a key role in controlling basic developmental programs for inflorescence and flower morphogenesis, it remained unknown until recently how they exert their specificity toward gene expression, both temporally and spatially. An emerging topic is the recruitment or eviction of chromatin factors through the activity of sequence-specific DNA-binding domains, present in the chromatin factors themselves or in partnering transcription factors. Here we summarize recent progress that has been made in this regard in the model plant Arabidopsis thaliana. We further outline the different possible modes through which chromatin complexes specifically target genes involved in flower development.

Covalent histone modifications and their effects on chromatin state and accessibility play a key role in the regulation of gene expression in eukaryotes. To gain insights into their functions during plant growth and development, the distribution of histone modifications can be analyzed at a genome-wide scale through chromatin immunoprecipitation assays followed by sequencing of the isolated genomic DNA. Here, we present a protocol for systematic analysis of the distribution and dynamic changes of selected histone modifications, during flower development in the model plant Arabidopsis thaliana. This protocol utilizes a previously established floral induction system to synchronize flower development, which allows the collection of sufficient plant material for analysis by genomic technologies. In this chapter, we describe how to use this system to study, from the same set of samples, chromatin and transcriptome dynamics during early stages of flower formation.

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.