MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

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.

Hybrid allele-specific ChIP-seq analysis identifies variation in brassinosteroid-responsive transcription factor binding linked to traits in maize

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.

Conserved and variable heat stress responses of the Heat Shock Factor transcription factor family in maize and Setaria viridis

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.

Mapping responsive genomic elements to heat stress in a maize diversity panel

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.

Functional dissection of the DORNRÖSCHEN-LIKE enhancer 2 during embryonic and phyllotactic patterning

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.

Specific chromatin changes mark lateral organ founder cells in the Arabidopsis inflorescence meristem

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.

Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis

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.

Enhancer elements can control spatial and temporal patterns of gene expression. Here the authors profile chromatin accessibility, histone modifications and gene expression during Arabidopsis flower development providing evidence for sets of distal enhancers acting in a highly stage-specific manner.

Dynamic and spatial restriction of Polycomb activity by plant histone demethylases

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.

Interactions between transcription factors and chromatin regulators in the control of flower development

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.

Profiling Histone Modifications in Synchronized Floral Tissues for Quantitative Resolution of Chromatin and Transcriptome Dynamics

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.

ULTRAPETALA1 and LEAFY pathways function independently in specifying identity and determinacy at the Arabidopsis floral meristem

BACKGROUND AND AIMS

The morphological variability of the flower in angiosperms, combined with its relatively simple structure, makes it an excellent model to study cell specification and the establishment of morphogenetic patterns. Flowers are the products of floral meristems, which are determinate structures that generate four different types of floral organs before terminating. The precise organization of the flower in whorls, each defined by the identity and number of organs it contains, is controlled by a multi-layered network involving numerous transcriptional regulators. In particular, the AGAMOUS (AG) MADS domain-containing transcription factor plays a major role in controlling floral determinacy in Arabidopsis thaliana in addition to specifying reproductive organ identity. This study aims to characterize the genetic interactions between the ULTRAPETALA1 (ULT1) and LEAFY (LFY) transcriptional regulators during flower morphogenesis, with a focus on AG regulation.

METHODS

Genetic and molecular approaches were used to address the question of redundancy and reciprocal interdependency for the establishment of flower meristem initiation, identity and termination. In particular, the effects of loss of both ULT1 and LFY function were determined by analysing flower developmental phenotypes of double-mutant plants. The dependency of each factor on the other for activating developmental genes was also investigated in gain-of-function experiments.

KEY RESULTS

The ULT1 and LFY pathways, while both activating AG expression in the centre of the flower meristem, functioned independently in floral meristem determinacy. Ectopic transcriptional activation by ULT1 of AG and AP3, another gene encoding a MADS domain-containing flower architect, did not depend on LFY function. Similarly, LFY did not require ULT1 function to ectopically determine floral fate.

CONCLUSIONS

The results indicate that the ULT1 and LFY pathways act separately in regulating identity and determinacy at the floral meristem. In particular, they independently induce AG expression in the centre of the flower to terminate meristem activity. A model is proposed whereby these independent contributions bring about a switch at the AG locus from an inactive to an active transcriptional state at the correct time and place during flower development.

Gene activation and cell fate control in plants: a chromatin perspective

In plants, environment-adaptable organogenesis extends throughout the lifespan, and iterative development requires repetitive rounds of activation and repression of several sets of genes. Eukaryotic genome compaction into chromatin forms a physical barrier for transcription; therefore, induction of gene expression requires alteration in chromatin structure. One of the present great challenges in molecular and developmental biology is to understand how chromatin is brought from a repressive to permissive state on specific loci and in a very specific cluster of cells, as well as how this state is further maintained and propagated through time and cell division in a cell lineage. In this review, we report recent discoveries implementing our knowledge on chromatin dynamics that modulate developmental gene expression. We also discuss how new data sets highlight plant specificities, likely reflecting requirement for a highly dynamic chromatin.

Natural variation of H3K27me3 distribution between two Arabidopsis accessions and its association with flanking transposable elements

Background

Histone H3 lysine 27 tri-methylation and lysine 9 di-methylation are independent repressive chromatin modifications in Arabidopsis thaliana. H3K27me3 is established and maintained by Polycomb repressive complexes whereas H3K9me2 is catalyzed by SUVH histone methyltransferases. Both modifications can spread to flanking regions after initialization and were shown to be mutually exclusive in Arabidopsis.

Results

We analyzed the extent of natural variation of H3K27me3 in the two accessions Landsberg erecta (Ler) and Columbia (Col) and their F1 hybrids. The majority of H3K27me3 target genes in Col were unchanged in Ler and F1 hybrids. A small number of Ler-specific targets were detected and confirmed. Consistent with a cis-regulatory mechanism for establishing H3K27me3, differential targets showed allele-specific H3K27me3 in hybrids. Five Ler-specific targets showed the active mark H3K4me3 in Col and for this group, differential H3K27me3 enrichment accorded to expression variation. On the other hand, the majority of Ler-specific targets were not expressed in Col, Ler or 17 other accessions. Instead of H3K27me3, the antagonistic mark H3K9me2 and other heterochromatic features were observed at these loci in Col. These loci were frequently flanked by transposable elements, which were often missing in the Ler genome assembly.

Conclusion

There is little variation in H3K27me3 occupancy within the species, although H3K27me3 targets were previously shown as overrepresented among differentially expressed genes. The existing variation in H3K27me3 seems mostly explained by flanking polymorphic transposable elements. These could nucleate heterochromatin, which then spreads into neighboring H3K27me3 genes, thus converting them to H3K9me2 targets.

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.

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.

Metaanalysis of ChIP-chip data

Genome-wide analysis of histone modifications via ChIP-chip (chromatin immunoprecipitation followed by whole genome tiling array hybridization) may generate lists of up to several thousand potential target genes. In the case of the model organism Arabidopsis thaliana, several databases are available to alleviate further characterization and classification of genomic data sets. The term metaanalysis has been coined for this type of multidatabase comparison. In this chapter, we describe open source software and web tools that perform transcriptional and functional analysis of target genes. Sources of transcription data and clustering tools to subdivide genes according to their expression pattern are described. The user is guided through all necessary steps, including data download and formatting. In addition, the Gene Ontology (GO) vocabulary and methods to uncover over- or underrepresented functions among target genes are introduced. Genomic targets of the histone H3K27me3 modification are presented as a case study to demonstrate that metaanalysis can uncover novel functions that were hidden in genomic data sets.

High-yield expression and purification of isotopically labeled norcoclaurine synthase, a Bet v 1-homologous enzyme, from Thalictrum flavum for NMR studies

The enzyme norcoclaurine synthase (NCS) found in the common meadow rue, Thalictrum flavum, and other plants shows sequence homology to members of the class 10 of pathogenesis related (PR 10) proteins that contains allergens such as the major birch pollen allergen Bet v 1, the major cherry allergen Pru av 1, and the major apple allergen Mal d 1. The enzyme is involved in the plant's secondary metabolism and is required for the production of bioactive secondary metabolites like morphine. Whereas the physiological function of PR 10 class allergens is still unknown, NCS activity has been studied in detail. Investigation of the structural properties of NCS by NMR spectroscopy can thus not only provide new information concerning the reaction mechanism of the enzyme, but is also expected to help clarify the long standing and heavily debated question on the physiological function as well as the reasons for the allergenic potential of members of this protein family. As the first important step towards the three-dimensional solution structure, we optimized expression of recombinant NCS in Escherichia coli and established an efficient purification protocol yielding high amounts of pure isotopically labeled active enzyme. The identity of NCS was confirmed by electrospray ionization mass spectrometry, and activity of the purified enzyme was determined by an assay detecting the radiolabeled reaction product. Spectroscopic analysis by NMR spectroscopy showed that the protein was properly folded with well defined tertiary structure.