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.

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.

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.

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.

Single-cell transcriptomics allows gene expression heterogeneity to be assessed at cellular resolution but the original location of each cell is unknown. Here the authors combine single nuclei RNA-seq with 3D spatial reconstruction of floral meristems to link gene activities with morphology.

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.

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.

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.

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.

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.

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 CCW₆GG), 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.

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.

Strigolactone biosynthesis in Medicago truncatula and rice requires the symbiotic GRAS-type transcription factors NSP1 and NSP2

Legume GRAS (GAI, RGA, SCR)-type transcription factors NODULATION SIGNALING PATHWAY1 (NSP1) and NSP2 are essential for rhizobium Nod factor-induced nodulation. Both proteins are considered to be Nod factor response factors regulating gene expression after symbiotic signaling. However, legume NSP1 and NSP2 can be functionally replaced by nonlegume orthologs, including rice (Oryza sativa) NSP1 and NSP2, indicating that both proteins are functionally conserved in higher plants. Here, we show that NSP1 and NSP2 are indispensable for strigolactone (SL) biosynthesis in the legume Medicago truncatula and in rice. Mutant nsp1 plants do not produce SLs, whereas in M. truncatula, NSP2 is essential for conversion of orobanchol into didehydro-orobanchol, which is the main SL produced by this species. The disturbed SL biosynthesis in nsp1 nsp2 mutant backgrounds correlates with reduced expression of DWARF27, a gene essential for SL biosynthesis. Rice and M. truncatula represent distinct phylogenetic lineages that split approximately 150 million years ago. Therefore, we conclude that regulation of SL biosynthesis by NSP1 and NSP2 is an ancestral function conserved in higher plants. NSP1 and NSP2 are single-copy genes in legumes, which implies that both proteins fulfill dual regulatory functions to control downstream targets after rhizobium-induced signaling as well as SL biosynthesis in nonsymbiotic conditions.

Arabidopsis thaliana transcription factors bZIP19 and bZIP23 regulate the adaptation to zinc deficiency

Zinc is an essential micronutrient for all living organisms. When facing a shortage in zinc supply, plants adapt by enhancing the zinc uptake capacity. The molecular regulators controlling this adaptation are not known. We present the identification of two closely related members of the Arabidopsis thaliana basic-region leucine-zipper (bZIP) transcription factor gene family, bZIP19 and bZIP23, that regulate the adaptation to low zinc supply. They were identified, in a yeast-one-hybrid screening, to associate to promoter regions of the zinc deficiency-induced ZIP4 gene of the Zrt- and Irt-related protein (ZIP) family of metal transporters. Although mutation of only one of the bZIP genes hardly affects plants, we show that the bzip19 bzip23 double mutant is hypersensitive to zinc deficiency. Unlike the wild type, the bzip19 bzip23 mutant is unable to induce the expression of a small set of genes that constitutes the primary response to zinc deficiency, comprising additional ZIP metal transporter genes. This set of target genes is characterized by the presence of one or more copies of a 10-bp imperfect palindrome in their promoter region, to which both bZIP proteins can bind. The bZIP19 and bZIP23 transcription factors, their target genes, and the characteristic cis zinc deficiency response elements they can bind to are conserved in higher plants. These findings are a significant step forward to unravel the molecular mechanism of zinc homeostasis in plants, allowing the improvement of zinc bio-fortification to alleviate human nutrition problems and phytoremediation strategies to clean contaminated soils.

Target genes of the MADS transcription factor SEPALLATA3: integration of developmental and hormonal pathways in the Arabidopsis flower

The molecular mechanisms by which floral homeotic genes act as major developmental switches to specify the identity of floral organs are still largely unknown. Floral homeotic genes encode transcription factors of the MADS-box family, which are supposed to assemble in a combinatorial fashion into organ-specific multimeric protein complexes. Major mediators of protein interactions are MADS-domain proteins of the SEPALLATA subfamily, which play a crucial role in the development of all types of floral organs. In order to characterize the roles of the SEPALLATA3 transcription factor complexes at the molecular level, we analyzed genome-wide the direct targets of SEPALLATA3. We used chromatin immunoprecipitation followed by ultrahigh-throughput sequencing or hybridization to whole-genome tiling arrays to obtain genome-wide DNA-binding patterns of SEPALLATA3. The results demonstrate that SEPALLATA3 binds to thousands of sites in the genome. Most potential target sites that were strongly bound in wild-type inflorescences are also bound in the floral homeotic agamous mutant, which displays only the perianth organs, sepals, and petals. Characterization of the target genes shows that SEPALLATA3 integrates and modulates different growth-related and hormonal pathways in a combinatorial fashion with other MADS-box proteins and possibly with non-MADS transcription factors. In particular, the results suggest multiple links between SEPALLATA3 and auxin signaling pathways. Our gene expression analyses link the genomic binding site data with the phenotype of plants expressing a dominant repressor version of SEPALLATA3, suggesting that it modulates auxin response to facilitate floral organ outgrowth and morphogenesis. Furthermore, the binding of the SEPALLATA3 protein to cis-regulatory elements of other MADS-box genes and expression analyses reveal that this protein is a key component in the regulatory transcriptional network underlying the formation of floral organs.

Most regulatory genes encode transcription factors, which modulate gene expression by binding to regulatory sequences of their target genes. In plants in particular, which genes are directly controlled by these transcription factors, and the molecular mechanisms of target gene recognition in vivo, are still largely unexplored. One of the best-understood developmental processes in plants is flower development. In different combinations, transcription factors of the MADS-box family control the identities of the different types of floral organs: sepals, petals, stamens, and carpels. Here, we present the first genome-wide analysis of binding sites of a MADS-box transcription factor in plants. We show that the MADS-domain protein SEPALLATA3 (SEP3) binds to the regulatory regions of thousands of potential target genes, many of which are also transcription factors. We provide insight into mechanisms of DNA recognition by SEP3, and suggest roles for other transcription factor families in SEP3 target gene regulation. In addition to effects on genes involved in floral organ identity, our data suggest that SEP3 binds to, and modulates, the transcription of target genes involved in hormonal signaling pathways.

The key floral regulator SEPALLATA3 binds to the promoters of a large number of potential direct target genes to integrate different growth-related and hormonal pathways in flower development.