Heat stress causes chromatin accessibility and related gene expression changes in crown tissues of barley (Hordeum vulgare).

Plant responses to stress caused by high temperatures involve changes occurring at the molecular, metabolic, and physiological levels. Understanding the mechanisms by which plants recognize signals to activate this response is a prerequisite for identifying key genes and signaling pathways and for obtaining heat-tolerant plants. We demonstrated the first implementation of an assay for transposase-accessible chromatin to identify open chromatin regions (OCRs) in crown tissues of barley using three genotypes carrying different allelic forms of the sdw1 gene encoding gibberellin 20-oxidase subjected to elevated temperatures. In parallel, we performed gene expression analysis, which allowed us to relate changes in chromatin state to changes in transcriptional activity. The obtained data revealed that the hypersensitive chromatin regions within the genes were more repeatable than those outside the gene intervals. We observed that prolonged exposure to high temperatures increased chromatin accessibility. Genes with OCRs in their regulatory regions were involved in stress signaling and tolerance, including calcium-dependent protein kinase, mitogen-activated protein kinase (MAPK3), receptor-like cytoplasmic kinase (RLK), TIFY domain-containing transcriptional regulator, bZIP transcription factor, and regulatory protein NPR1. The effect of genotype on gene expression was not as pronounced as that of temperature. By combining results from the differential analysis of chromatin accessibility and expression profiles, we identified genes with high temperature-induced changes in chromatin accessibility associated with expression alterations. Importantly, our data revealed a relationship between the loss of chromatin accessibility in response to heat and the downregulation of genes related to gibberellin signaling.

Estimating DNA-Binding Specificities of Transcription Factors Using SELEX-Seq.

Here we provide an updated protocol for the Systematic Evolution of Ligands followed by massively parallel sequencing (SELEX-seq) method to study protein-DNA interaction specificities. This in vitro method is used to characterize DNA-binding specificities of transcription factors (TFs). The procedure is based on cycles of immunoprecipitation of protein-DNA complexes, starting with a randomized DNA library of defined fragment length, followed by massively parallel sequencing. The updated protocol includes aspects of experimental design and procedure as well as basic instructions on data analysis.

Immunoprecipitation-Mass Spectrometry (IP-MS) of Protein-Protein Interactions of Nuclear-Localized Plant Proteins.

Transcription factors that act within a gene regulatory network (GRN) often interact with other proteins such as chromatin remodeling factors, histone modifiers, and other co-regulators. Characterizing these interactions is crucial for understanding the function and mechanism of action of a transcription factor. Here, a method for the identification of protein-protein interactions of nuclear-localized, transcription-associated factors is described. The method is based on the immunoprecipitation (IP) of a fluorophore-tagged target, followed by mass spectrometry (MS), peptide identification, and quantification of interacting proteins. By applying label-free quantification to IPs and their input protein extracts, statistically controlled protein enrichment ratios uncover high-confidence interaction partners of the target. A complete step-by-step procedure, including sample preparation, MS settings, data analysis, and visualization is provided.

Dual specificity and target gene selection by the MADS-domain protein FRUITFULL.

How transcription factors attain their target gene specificity and how this specificity may be modulated, acquiring different regulatory functions through the development of plant tissues, is an open question. Here we characterized different regulatory roles of the MADS-domain transcription factor FRUITFULL (FUL) in flower development and mechanisms modulating its activity. We found that the dual role of FUL in regulating floral transition and pistil development is associated with its different in vivo patterns of DNA binding in both tissues. Characterization of FUL protein complexes by liquid chromatography-tandem mass spectrometry and SELEX-seq experiments shows that aspects of tissue-specific target site selection can be predicted by tissue-specific variation in the composition of FUL protein complexes with different DNA binding specificities, without considering the chromatin status of the target region. This suggests a role for dynamic changes in FUL TF complex composition in reshaping the regulatory functions of FUL during flower development.

A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data.

Cellular heterogeneity in growth and differentiation results in organ patterning. Single-cell transcriptomics allows characterization of gene expression heterogeneity in developing organs at unprecedented resolution. However, the original physical location of the cell is lost during this methodology. To recover the original location of cells in the developing organ is essential to link gene activity with cellular identity and function in plants. Here, we propose a method to reconstruct genome-wide gene expression patterns of individual cells in a 3D flower meristem by combining single-nuclei RNA-seq with microcopy-based 3D spatial reconstruction. By this, gene expression differences among meristematic domains giving rise to different tissue and organ types can be determined. As a proof of principle, the method is used to trace the initiation of vascular identity within the floral meristem. Our work demonstrates the power of spatially reconstructed single cell transcriptome atlases to understand plant morphogenesis. The floral meristem 3D gene expression atlas can be accessed at http://threed-flower-meristem.herokuapp.com .

The intervening domain is required for DNA-binding and functional identity of plant MADS transcription factors.

The MADS transcription factors (TF) are an ancient eukaryotic protein family. In plants, the family is divided into two main lineages. Here, we demonstrate that DNA binding in both lineages absolutely requires a short amino acid sequence C-terminal to the MADS domain (M domain) called the Intervening domain (I domain) that was previously defined only in type II lineage MADS. Structural elucidation of the MI domains from the floral regulator, SEPALLATA3 (SEP3), shows a conserved fold with the I domain acting to stabilise the M domain. Using the floral organ identity MADS TFs, SEP3, APETALA1 (AP1) and AGAMOUS (AG), domain swapping demonstrate that the I domain alters genome-wide DNA-binding specificity and dimerisation specificity. Introducing AG carrying the I domain of AP1 in the Arabidopsis ap1 mutant resulted in strong complementation and restoration of first and second whorl organs. Taken together, these data demonstrate that the I domain acts as an integral part of the DNA-binding domain and significantly contributes to the functional identity of the MADS TF.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system.

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Cell identity specification in plants: lessons from flower development.

Multicellular organisms display a fascinating complexity of cellular identities and patterns of diversification. The concept of 'cell type' aims to describe and categorize this complexity. In this review, we discuss the traditional concept of cell types and highlight the impact of single-cell technologies and spatial omics on the understanding of cellular differentiation in plants. We summarize and compare position-based and lineage-based mechanisms of cell identity specification using flower development as a model system. More than understanding ontogenetic origins of differentiated cells, an important question in plant science is to understand their position- and developmental stage-specific heterogeneity. Combinatorial action and crosstalk of external and internal signals is the key to cellular heterogeneity, often converging on transcription factors that orchestrate gene expression programs.

The Chromatin-Associated Protein PWO1 Interacts with Plant Nuclear Lamin-like Components to Regulate Nuclear Size.

Spatial organization of chromatin contributes to gene regulation of many cellular processes and includes a connection of chromatin with the nuclear lamina (NL). The NL is a protein mesh that resides underneath the inner nuclear membrane and consists of lamins and lamina-associated proteins. Chromatin regions associated with lamins in animals are characterized mostly by constitutive heterochromatin, but association with facultative heterochromatin mediated by Polycomb-group (PcG) proteins has been reported as well. In contrast with animals, plant NL components are largely not conserved and NL association with chromatin is poorly explored. Here, we present the connection between the lamin-like protein, CROWDED NUCLEI1 (CRWN1), and the chromatin- and PcG-associated component, PROLINE-TRYPTOPHANE-TRYPTOPHANE-PROLINE INTERACTOR OF POLYCOMBS1, in Arabidopsis (Arabidopsis thaliana). We show that PWO1 and CRWN1 proteins associate physically with each other, act in the same pathway to maintain nuclear morphology, and control expression of a similar set of target genes. Moreover, we demonstrate that transiently expressed PWO1 proteins form foci located partially at the subnuclear periphery. Ultimately, as CRWN1 and PWO1 are plant-specific, our results argue that plants might have developed an equivalent, rather than homologous, mechanism of linking chromatin repression and NL.

Building Transcription Factor Binding Site Models to Understand Gene Regulation in Plants.

Transcription factors (TFs) are key cellular components that control gene expression. They recognize specific DNA sequences, the TF binding sites (TFBSs), and thus are targeted to specific regions of the genome where they can recruit transcriptional co-factors and/or chromatin regulators to fine-tune spatiotemporal gene regulation. Therefore, the identification of TFBSs in genomic sequences and their subsequent quantitative modeling is of crucial importance for understanding and predicting gene expression. Here, we review how TFBSs can be determined experimentally, how the TFBS models can be constructed in silico, and how they can be optimized by taking into account features such as position interdependence within TFBSs, DNA shape, and/or by introducing state-of-the-art computational algorithms such as deep learning methods. In addition, we discuss the integration of context variables into the TFBS modeling, including nucleosome positioning, chromatin states, methylation patterns, 3D genome architectures, and TF cooperative binding, in order to better predict TF binding under cellular contexts. Finally, we explore the possibilities of combining the optimized TFBS model with technological advances, such as targeted TFBS perturbation by CRISPR, to better understand gene regulation, evolution, and plant diversity.

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 manner1-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.

Duplication of an upstream silencer of FZP increases grain yield in rice.

Transcriptional silencer and copy number variants (CNVs) are associated with gene expression. However, their roles in generating phenotypes have not been well studied. Here we identified a rice quantitative trait locus, SGDP7 (Small Grain and Dense Panicle 7). SGDP7 is identical to FZP (FRIZZY PANICLE), which represses the formation of axillary meristems. The causal mutation of SGDP7 is an 18-bp fragment, named CNV-18bp, which was inserted ~5.3 kb upstream of FZP and resulted in a tandem duplication in the cultivar Chuan 7. The CNV-18bp duplication repressed FZP expression, prolonged the panicle branching period and increased grain yield by more than 15% through substantially increasing the number of spikelets per panicle (SPP) and slightly decreasing the 1,000-grain weight (TGW). The transcription repressor OsBZR1 binds the CGTG motifs in CNV-18bp and thereby represses FZP expression, indicating that CNV-18bp is the upstream silencer of FZP. These findings showed that the silencer CNVs coordinate a trade-off between SPP and TGW by fine-tuning FZP expression, and balancing the trade-off could enhance yield potential.

Differences in DNA Binding Specificity of Floral Homeotic Protein Complexes Predict Organ-Specific Target Genes.

Floral organ identities in plants are specified by the combinatorial action of homeotic master regulatory transcription factors. However, how these factors achieve their regulatory specificities is still largely unclear. Genome-wide in vivo DNA binding data show that homeotic MADS domain proteins recognize partly distinct genomic regions, suggesting that DNA binding specificity contributes to functional differences of homeotic protein complexes. We used in vitro systematic evolution of ligands by exponential enrichment followed by high-throughput DNA sequencing (SELEX-seq) on several floral MADS domain protein homo- and heterodimers to measure their DNA binding specificities. We show that specification of reproductive organs is associated with distinct binding preferences of a complex formed by SEPALLATA3 and AGAMOUS. Binding specificity is further modulated by different binding site spacing preferences. Combination of SELEX-seq and genome-wide DNA binding data allows differentiation between targets in specification of reproductive versus perianth organs in the flower. We validate the importance of DNA binding specificity for organ-specific gene regulation by modulating promoter activity through targeted mutagenesis. Our study shows that intrafamily protein interactions affect DNA binding specificity of floral MADS domain proteins. Differential DNA binding of MADS domain protein complexes plays a role in the specificity of target gene regulation.

SELEX-Seq: A Method to Determine DNA Binding Specificities of Plant Transcription Factors.

Systematic evolution of ligands by exponential enrichment (SELEX) is a method that allows isolating specific nucleotide sequences that interact with a DNA binding protein of choice. By using a transcription factor (TF) and a randomized pool of double-stranded DNA, this technique can be used to characterize TF DNA binding specificities and affinities. The method is based on protein-DNA complex immunoprecipitation with protein-specific antibodies and subsequent DNA selection and amplification. Application of massively parallel sequencing (-seq) at each cycle of SELEX allows determining the relative affinities to any DNA sequence for any transcription factor or TF complex. The resulting TF DNA binding motifs can be used to predict potential DNA binding sites in genomes and thereby direct target genes of TFs.

Post-transcriptional control of GRF transcription factors by microRNA miR396 and GIF co-activator affects leaf size and longevity.

The growth-regulating factors (GRFs) are plant-specific transcription factors. They form complexes with GRF-interacting factors (GIFs), a small family of transcriptional co-activators. In Arabidopsis thaliana, seven out of the nine GRFs are controlled by microRNA miR396. Analysis of Arabidopsis plants carrying a GRF3 allele insensitive to miR396 revealed a strong boost in the number of cells in leaves, which was further enhanced synergistically by an additional increase of GIF1 levels. Genetic experiments revealed that GRF3 can still increase cell number in gif1 mutants, albeit to a much lesser extent. Genome-wide transcript profiling indicated that the simultaneous increase of GRF3 and GIF1 levels causes additional effects in gene expression compared to either of the transgenes alone. We observed that GIF1 interacts in vivo with GRF3, as well as with chromatin-remodeling complexes, providing a mechanistic explanation for the synergistic activities of a GRF3-GIF1 complex. Interestingly, we found that, in addition to the leaf size, the GRF system also affects the organ longevity. Genetic and molecular analysis revealed that the functions of GRFs in leaf growth and senescence can be uncoupled, demonstrating that the miR396-GRF-GIF network impinges on different stages of leaf development. Our results integrate the post-transcriptional control of the GRF transcription factors with the progression of leaf development.

Structural determinants of DNA recognition by plant MADS-domain transcription factors.

Plant MADS-domain transcription factors act as key regulators of many developmental processes. Despite the wealth of information that exists about these factors, the mechanisms by which they recognize their cognate DNA-binding site, called CArG-box (consensus CCW6GG), and how different MADS-domain proteins achieve DNA-binding specificity, are still largely unknown. We used information from in vivo ChIP-seq experiments, in vitro DNA-binding data and evolutionary conservation to address these important questions. We found that structural characteristics of the DNA play an important role in the DNA binding of plant MADS-domain proteins. The central region of the CArG-box largely resembles a structural motif called 'A-tract', which is characterized by a narrow minor groove and may assist bending of the DNA by MADS-domain proteins. Periodically spaced A-tracts outside the CArG-box suggest additional roles for this structure in the process of DNA binding of these transcription factors. Structural characteristics of the CArG-box not only play an important role in DNA-binding site recognition of MADS-domain proteins, but also partly explain differences in DNA-binding specificity of different members of this transcription factor family and their heteromeric complexes.

Naturally occurring allele diversity allows potato cultivation in northern latitudes.

Potato (Solanum tuberosum L.) originates from the Andes and evolved short-day-dependent tuber formation as a vegetative propagation strategy. Here we describe the identification of a central regulator underlying a major-effect quantitative trait locus for plant maturity and initiation of tuber development. We show that this gene belongs to the family of DOF (DNA-binding with one finger) transcription factors and regulates tuberization and plant life cycle length, by acting as a mediator between the circadian clock and the StSP6A mobile tuberization signal. We also show that natural allelic variants evade post-translational light regulation, allowing cultivation outside the geographical centre of origin of potato. Potato is a member of the Solanaceae family and is one of the world's most important food crops. This annual plant originates from the Andean regions of South America. Potato develops tubers from underground stems called stolons. Its equatorial origin makes potato essentially short-day dependent for tuberization and potato will not make tubers in the long-day conditions of spring and summer in the northern latitudes. When introduced in temperate zones, wild material will form tubers in the course of the autumnal shortening of day-length. Thus, one of the first selected traits in potato leading to a European potato type is likely to have been long-day acclimation for tuberization. Potato breeders can exploit the naturally occurring variation in tuberization onset and life cycle length, allowing varietal breeding for different latitudes, harvest times and markets.

Proteomics-based identification of low-abundance signaling and regulatory protein complexes in native plant tissues.

Owing to the low abundance of signaling proteins and transcription factors, their protein complexes are not easily identified by classical proteomics. The isolation of these protein complexes from endogenous plant tissues (rather than plant cell cultures) is therefore an important technical challenge. Here, we describe a sensitive, quantitative proteomics-based procedure to determine the composition of plant protein complexes. The method makes use of fluorophore-tagged protein immunoprecipitation (IP) and label-free mass spectrometry (MS)-based quantification to correct for nonspecifically precipitated proteins. We provide procedures for the isolation of membrane-bound receptor complexes and transcriptional regulators from nuclei. The protocol consists of an IP step (~6 h) and sample preparation for liquid chromatography-tandem MS (LC-MS/MS; 2 d). We also provide a guide for data analysis. Our single-step affinity purification protocol is a good alternative to two-step tandem affinity purification (TAP), as it is shorter and relatively easy to perform. The data analysis by label-free quantification (LFQ) requires a cheaper and less challenging experimental setup compared with known labeling techniques in plants.

Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies.

Members of the MADS-box transcription factor family play essential roles in almost every developmental process in plants. Many MADS-box genes have conserved functions across the flowering plants, but some have acquired novel functions in specific species during evolution. The analyses of MADS-domain protein interactions and target genes have provided new insights into their molecular functions. Here, we review recent findings on MADS-box gene functions in Arabidopsis and discuss the evolutionary history and functional diversification of this gene family in plants. We also discuss possible mechanisms of action of MADS-domain proteins based on their interactions with chromatin-associated factors and other transcriptional regulators.

Characterization of MADS-domain transcription factor complexes in Arabidopsis flower development.

Floral organs are specified by the combinatorial action of MADS-domain transcription factors, yet the mechanisms by which MADS-domain proteins activate or repress the expression of their target genes and the nature of their cofactors are still largely unknown. Here, we show using affinity purification and mass spectrometry that five major floral homeotic MADS-domain proteins (AP1, AP3, PI, AG, and SEP3) interact in floral tissues as proposed in the "floral quartet" model. In vitro studies confirmed a flexible composition of MADS-domain protein complexes depending on relative protein concentrations and DNA sequence. In situ bimolecular fluorescent complementation assays demonstrate that MADS-domain proteins interact during meristematic stages of flower development. By applying a targeted proteomics approach we were able to establish a MADS-domain protein interactome that strongly supports a mechanistic link between MADS-domain proteins and chromatin remodeling factors. Furthermore, members of other transcription factor families were identified as interaction partners of floral MADS-domain proteins suggesting various specific combinatorial modes of action.