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

Reflections on the ABC model of flower development

The formulation of the ABC model by a handful of pioneer plant developmental geneticists was a seminal event in the quest to answer a seemingly simple question: how are flowers formed? Fast forward 30 years and this elegant model has generated a vibrant and diverse community, capturing the imagination of developmental and evolutionary biologists, structuralists, biochemists and molecular biologists alike. Together they have managed to solve many floral mysteries, uncovering the regulatory processes that generate the characteristic spatio-temporal expression patterns of floral homeotic genes, elucidating some of the mechanisms allowing ABC genes to specify distinct organ identities, revealing how evolution tinkers with the ABC to generate morphological diversity, and even shining a light on the origins of the floral gene regulatory network itself. Here we retrace the history of the ABC model, from its genesis to its current form, highlighting specific milestones along the way before drawing attention to some of the unsolved riddles still hidden in the floral alphabet.

Banana MabHLH28 positively regulates the expression of softening-related genes to mediate fruit ripening independently or via cooperating with MaWRKY49/111

Texture softening is a physiological indicator of fruit ripening, which eventually contributes to fruit quality and the consumer's acceptance. Despite great progress having been made in identification of the genes related to fruit softening, the upstream transcriptional regulatory pathways of these softening-related genes are not fully elucidated. Here, a novel bHLH gene, designated as MabHLH28, was identified because of its significant upregulation in banana fruit ripening. DAP-Seq analysis revealed that MabHLH28 bound to the core sequence of 'CAYGTG' presented in promoter regions of fruit softening-associated genes, such as the genes related to cell wall modification (MaPG3, MaPE1, MaPL5, MaPL8, MaEXP1, MaEXP2, MaEXPA2, and MaEXPA15) and starch degradation (MaGWD1 and MaLSF2), and these bindings were validated by EMSA and DLR assays. Transient overexpression and knockdown of MabHLH28 in banana fruit resulted in up- and down-regulation of softening-related genes, thereby hastening and postponing fruit ripening. Furthermore, overexpression of MabHLH28 in tomato accelerated the ripening process by elevating the accumulation of softening-associated genes. In addition, MabHLH28 showed interaction withMaWRKY49/111 and itself to form protein complexes, which could combinatorically strengthen the transcription of softening-associated genes. Taken together, our findings suggest that MabHLH28 mediates fruit softening by upregulating the expression of softening-related genes either alone or in combination with MaWRKY49/111.

Research advances on the gene regulation of floral development and color in orchids

Orchidaceae is one of the largest monocotyledon families and contributes significantly to worldwide biodiversity, with value in the fields of landscaping, medicine, and ecology. The diverse phenotypes and vibrant colors of orchid floral organs make them excellent research objects for investigating flower development and pigmentation. In recent years, a number of orchid genomes have been published, laying the molecular foundation for revealing flower development and color presentation. In this article, we review transcription factors, the structural genes responsible for the floral pigment synthesis pathways, the molecular mechanisms of flower morphogenesis, and the potential relationship between flower type and flower color. This study provides a theoretical reference for the research on molecular mechanisms related to flower morphogenesis and color presentation, genetic improvement, and new variety creation in orchids.

An ecological transcriptome approach to capture the molecular and physiological mechanisms of mass flowering in Shorea curtisii

Climatic factors have commonly been attributed as the trigger of general flowering, a unique community-level mass flowering phenomenon involving most dipterocarp species that forms the foundation of Southeast Asian tropical rainforests. This intriguing flowering event is often succeeded by mast fruiting, which provides a temporary yet substantial burst of food resources for animals, particularly frugivores. However, the physiological mechanism that triggers general flowering, particularly in dipterocarp species, is not well understood largely due to its irregular and unpredictable occurrences in the tall and dense forests. To shed light on this mechanism, we employed ecological transcriptomic analyses on an RNA-seq dataset of a general flowering species, Shorea curtisii (Dipterocarpaceae), sequenced from leaves and buds collected at multiple vegetative and flowering phenological stages. We assembled 64,219 unigenes from the transcriptome of which 1,730 and 3,559 were differentially expressed in the leaf and the bud, respectively. Differentially expressed unigene clusters were found to be enriched with homologs of Arabidopsis thaliana genes associated with response to biotic and abiotic stresses, nutrient level, and hormonal treatments. When combined with rainfall data, our transcriptome data reveals that the trees were responding to a brief period of drought prior to the elevated expression of key floral promoters and followed by differential expression of unigenes that indicates physiological changes associated with the transition from vegetative to reproductive stages. Our study is timely for a representative general flowering dipterocarp species that occurs in forests that are under the constant threat of deforestation and climate change as it pinpoints important climate sensitive and flowering-related homologs and offers a glimpse into the cascade of gene expression before and after the onset of floral initiation.

Sea-ATI unravels novel vocabularies of plant active cistrome

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants.

Evolutionary studies of the bHLH transcription factors belonging to MBW complex: their role in seed development

BACKGROUND AND AIMS

The MBW complex consist of proteins belonging to three major families (MYB, bHLH and WDR) involved in various processes throughout plant development: epidermal cell development, mucilage secretory cells and flavonoid biosynthesis. Recently, it has been reported that TT8, encoding a bHLH transcription factor, is involved in the biosynthesis of flavonoids in the seed coat and it also plays a role in bypassing the postzygotic barrier resulting from an unbalance in genetic loads of the parental lines. Here, we focus on the functional evolution, in seed development, of the bHLH proteins that are part of the MBW complex, complemented with a literature review.

METHODS

Phylogenetic analyses performed across seed plants and expression analyses in the reproductive tissues of four selected angiosperms (Arabidopsis thaliana, Brassica napus, Capsella rubella and Solanum lycopersicum) allow us to hypothesize on the evolution of its functions.

KEY RESULTS

TT8 expression in the innermost layer of the seed coat is conserved in the selected angiosperms. However, except for Arabidopsis, TT8 is also expressed in ovules, carpels and fruits. The homologues belonging to the sister clade of TT8, EGL3/GL3, involved in trichome development, are expressed in the outermost layer of the seed coat, suggesting potential roles in mucilage.

CONCLUSIONS

The ancestral function of these genes appears to be flavonoid biosynthesis, and the conservation of TT8 expression patterns in the innermost layer of the seed coat in angiosperms suggests that their function in postzygotic barriers might also be conserved. Moreover, the literature review and the results of the present study suggest a sophisticated association, linking the mechanisms of action of these genes to the cross-communication activity between the different tissues of the seed. Thus, it provides avenues to study the mechanisms of action of TT8 in the postzygotic triploid block, which is crucial because it impacts seed development in unbalanced crosses.

Floral Development on Vitis vinifera Is Associated with MADS-Box Transcription Factors through the Transcriptional Regulation of VviZIP3

Several grapevine (Vitis vinifera L.) cultivars show a tendency to develop parthenocarpic seedless grapes, affecting fruit yield and quality. This reproductive disorder originates in defective ovule fertilization due to a failure in pollen tube growth. Zinc (Zn) is a crucial trace element, playing a vital role in various physiological and metabolic processes. It is particularly essential for the healthy growth of flowers and fruits. Insufficient zinc has been suggested as a potential reason for issues in this development process. This microelement is taken up through a mechanism that involves transporters, including the ZRT-IRT-like protein (ZIP) gene family, associated with the influx of Zn into the cell. In grapevines, 20 genes for ZIP-type transporters have been described. In this study, we analyzed the expression pattern of VviZIP3 during flower development and employ transgenic methods to assess its transcriptional regulation. Furthermore, through computational examination of the promoter region, we identified two CArG boxes, recognized as responsive elements to MADS transcription factors. These factors play a key role in shaping various components of a flower, such as pollen. Our investigation of the VviZIP3 promoter confirms the functionality of these CArG boxes. Overall, our results suggest that the increased expression of VviZIP3 during flowering is likely under the influence of MADS transcription factors.

TM3 and STM3 Promote Flowering Together with FUL2 and MBP20, but Act Antagonistically in Inflorescence Branching in Tomato

The moment at which a plant transitions to reproductive development is paramount to its life cycle and is strictly controlled by many genes. The transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) plays a central role in this process in Arabidopsis. However, the role of SOC1 in tomato (Solanum lycopersicum) has been sparsely studied. Here, we investigated the function of four tomato SOC1 homologs in the floral transition and inflorescence development. We thoroughly characterized the SOC1-like clade throughout the Solanaceae and selected four tomato homologs that are dynamically expressed upon the floral transition. We show that of these homologs, TOMATO MADS 3 (TM3) and SISTER OF TM3 (STM3) promote the primary and sympodial transition to flowering, while MADS-BOX PROTEIN 23 (MBP23) and MBP18 hardly contribute to flowering initiation in the indeterminate cultivar Moneyberg. Protein-protein interaction assays and whole-transcriptome analysis during reproductive meristem development revealed that TM3 and STM3 interact and share many targets with FRUITFULL (FUL) homologs, including cytokinin regulators. Furthermore, we observed that mutating TM3/STM3 affects inflorescence development, but counteracts the inflorescence-branching phenotype of ful2 mbp20. Collectively, this indicates that TM3/STM3 promote the floral transition together with FUL2/MBP20, while these transcription factors have opposite functions in inflorescence development.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)₆GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKC^C-type MTFs (M^C-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. M^C-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which M^C-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of M^C-MTFs). A special focus of our review is on the structural prerequisites of M^C-MTFs that enable the specific recognition of target genes.

Floral Homeotic Factors: A Question of Specificity

MADS-domain transcription factors are involved in the control of a multitude of processes in eukaryotes, and in plants, they play particularly important roles during reproductive development. Among the members of this large family of regulatory proteins are the floral organ identity factors, which specify the identities of the different types of floral organs in a combinatorial manner. Much has been learned over the past three decades about the function of these master regulators. For example, it has been shown that they have similar DNA-binding activities and that their genome-wide binding patterns exhibit large overlaps. At the same time, it appears that only a minority of binding events lead to changes in gene expression and that the different floral organ identity factors have distinct sets of target genes. Thus, binding of these transcription factors to the promoters of target genes alone may not be sufficient for their regulation. How these master regulators achieve specificity in a developmental context is currently not well understood. Here, we review what is known about their activities and highlight open questions that need to be addressed to gain more detailed insights into the molecular mechanisms underlying their functions. We discuss evidence for the involvement of cofactors as well as the results from studies on transcription factors in animals that may be instructive for a better understanding of how the floral organ identity factors achieve regulatory specificity.

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.

Designing artificial synthetic promoters for accurate, smart, and versatile gene expression in plants

With the development of high-throughput biology techniques and artificial intelligence, it has become increasingly feasible to design and construct artificial biological parts, modules, circuits, and even whole systems. To overcome the limitations of native promoters in controlling gene expression, artificial promoter design aims to synthesize short, inducible, and conditionally controlled promoters to coordinate the expression of multiple genes in diverse plant metabolic and signaling pathways. Synthetic promoters are versatile and can drive gene expression accurately with smart responses; they show potential for enhancing desirable traits in crops, thereby improving crop yield, nutritional quality, and food security. This review first illustrates the importance of synthetic promoters, then introduces promoter architecture and thoroughly summarizes advances in synthetic promoter construction. Restrictions to the development of synthetic promoters and future applications of such promoters in synthetic plant biology and crop improvement are also discussed.

Flower Development in Arabidopsis

Like in other angiosperms, the development of flowers in Arabidopsis starts right after the floral transition, when the shoot apical meristem (SAM) stops producing leaves and makes flowers instead. On the flanks of the SAM emerge the flower meristems (FM) that will soon differentiate into the four main floral organs, sepals, petals, stamens, and pistil, stereotypically arranged in concentric whorls. Each phase of flower development-floral transition, floral bud initiation, and floral organ development-is under the control of specific gene networks. In this chapter, we describe these different phases and the gene regulatory networks involved, from the floral transition to the floral termination.

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.

Context-specific functions of transcription factors controlling plant development: From leaves to flowers

Plant development is regulated by transcription factors that often act in more than one process and stage of development. Yet the molecular mechanisms that govern the functional diversity and specificity of these proteins remains far from understood. Flower development provides an ideal context to study these mechanisms since the development of distinct floral organs depends on similar but distinct combinations of transcriptional regulators. Recent work also highlights the importance of leaf polarity regulators as additional key factors in flower initiation, floral organ morphogenesis, and possibly floral organ positioning. A detailed understanding of how these factors work in combination will enable us to address outstanding questions in flower development including how distinct shapes and positions of floral organs are generated. Experimental approaches and computer-based modeling will be required to characterize gene-regulatory networks at the level of single cells.

Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Transcription factors (TFs) play a critical role in determining cell fate decisions by integrating developmental and environmental signals through binding to specific cis-regulatory modules and regulating spatio-temporal specificity of gene expression patterns. Precise identification of functional TF binding sites in time and space not only will revolutionize our understanding of regulatory networks governing cell fate decisions but is also instrumental to uncover how genetic variations cause morphological diversity or disease. In this review, we discuss recent advances in mapping TF binding sites and characterizing the various parameters underlying the complexity of binding site recognition by TFs.

Unraveling the role of MADS transcription factor complexes in apple tree dormancy

A group of MADS transcription factors (TFs) are believed to control temperature-mediated bud dormancy. These TFs, called DORMANCY-ASSOCIATED MADS-BOX (DAM), are encoded by genes similar to SHORT VEGETATIVE PHASE (SVP) from Arabidopsis. MADS proteins form transcriptional complexes whose combinatory composition defines their molecular function. However, how MADS multimeric complexes control the dormancy cycle in trees is unclear. Apple MdDAM and other dormancy-related MADS proteins form complexes with MdSVPa, which is essential for the ability of transcriptional complexes to bind to DNA. Sequential DNA-affinity purification sequencing (seq-DAP-seq) was performed to identify the genome-wide binding sites of apple MADS TF complexes. Target genes associated with the binding sites were identified by combining seq-DAP-seq data with transcriptomics datasets obtained using a glucocorticoid receptor fusion system, and RNA-seq data related to apple dormancy. We describe a gene regulatory network (GRN) formed by MdSVPa-containing complexes, which regulate the dormancy cycle in response to environmental cues and hormonal signaling pathways. Additionally, novel molecular evidence regarding the evolutionary functional segregation between DAM and SVP proteins in the Rosaceae is presented. MdSVPa sequentially forms complexes with the MADS TFs that predominate at each dormancy phase, altering its DNA-binding specificity and, therefore, the transcriptional regulation of its target genes.

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.

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.

MADS transcription factors regulate multiple aspects of plant development. Here the authors show that the intervening I domain is conserved in both type I and type II plant MADS lineages and contributes to the functional identity of the protein by influencing both DNA binding activity and dimerisation specificity.

ZmOrphan94 Transcription Factor Downregulates ZmPEPC1 Gene Expression in Maize Bundle Sheath Cells

Spatial separation of the photosynthetic reactions is a key feature of C₄ metabolism. In most C₄ plants, this separation requires compartmentation of photosynthetic enzymes between mesophyll (M) and bundle sheath (BS) cells. The upstream region of the gene encoding the maize PHOSPHOENOLPYRUVATE CARBOXYLASE 1 (ZmPEPC1) has been shown sufficient to drive M-specific ZmPEPC1 gene expression. Although this region has been well characterized, to date, only few trans-factors involved in the ZmPEPC1 gene regulation were identified. Here, using a yeast one-hybrid approach, we have identified three novel maize transcription factors ZmHB87, ZmCPP8, and ZmOrphan94 as binding to the ZmPEPC1 upstream region. Bimolecular fluorescence complementation assays in maize M protoplasts unveiled that ZmOrphan94 forms homodimers and interacts with ZmCPP8 and with two other ZmPEPC1 regulators previously reported, ZmbHLH80 and ZmbHLH90. Trans-activation assays in maize M protoplasts unveiled that ZmHB87 does not have a clear transcriptional activity, whereas ZmCPP8 and ZmOrphan94 act as activator and repressor, respectively. Moreover, we observed that ZmOrphan94 reduces the trans-activation activity of both activators ZmCPP8 and ZmbHLH90. Using the electromobility shift assay, we showed that ZmOrphan94 binds to several cis-elements present in the ZmPEPC1 upstream region and one of these cis-elements overlaps with the ZmbHLH90 binding site. Gene expression analysis revealed that ZmOrphan94 is preferentially expressed in the BS cells, suggesting that ZmOrphan94 is part of a transcriptional regulatory network downregulating ZmPEPC1 transcript level in the BS cells. Based on both this and our previous work, we propose a model underpinning the importance of a regulatory mechanism within BS cells that contributes to the M-specific ZmPEPC1 gene expression.

Transcriptome analysis of gynoecium morphogenesis uncovers the chronology of gene regulatory network activity

The gynoecium is the most complex organ formed by the flowering plants. It encloses the ovules, provides a surface for pollen contact and self-incompatibility reactions, allows pollen tube growth, and, post fertilization, develops into the fruit. Consequently, the regulation of gynoecium morphogenesis is complex and appropriate timing of this process in part determines reproductive success. However, little is known about the global control of gynoecium development, even though many regulatory genes have been characterized. Here, we characterized dynamic gene expression changes using laser-microdissected gynoecium tissue from four developmental stages in Arabidopsis. We provide a high-resolution map of global expression dynamics during gynoecium morphogenesis and link these to the gynoecium interactome. We reveal groups of genes acting together early and others acting late in morphogenesis. Clustering of co-expressed genes enables comparisons between the leaf, shoot apex, and gynoecium transcriptomes, allowing the dissection of common and distinct regulators. Furthermore, our results lead to the discovery of genes with putative transcription factor activity (B3LF1, -2, DOFLF1), which, when mutated, lead to impaired gynoecium expansion, illustrating that global transcriptome analyses reveal yet unknown developmental regulators. Our data show that genes encoding highly interacting proteins, such as SEPALLATA3, AGAMOUS, and TOPLESS, are expressed evenly during development but switch interactors over time, whereas stage-specific proteins tend to have fewer interactors. Our analysis connects specific transcriptional regulator activities, protein interactions, and underlying metabolic processes, contributing toward a dynamic network model for gynoecium development.

Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

Arabidopsis Abscisic Acid Repressor 1 is a susceptibility hub that interacts with multiple Pseudomonas syringae effectors

Pathogens secrete effector proteins into host cells to suppress host immunity and promote pathogen virulence, although many features at the molecular interface of host-pathogen interactions remain to be characterized. In a yeast two-hybrid assay, we found that the Pseudomonas syringae effector HopZ1a interacts with the Arabidopsis transcriptional regulator Abscisic Acid Repressor 1 (ABR1). Further analysis revealed that ABR1 interacts with multiple P. syringae effectors, suggesting that it may be targeted as a susceptibility hub. Indeed, loss-of-function abr1 mutants exhibit reduced susceptibility to a number of P. syringae strains. The ABR1 protein comprises a conserved APETALA2 (AP2) domain flanked by long regions of predicted structural disorder. We verified the DNA-binding activity of the AP2 domain and demonstrated that the disordered domains act redundantly to enhance DNA binding and to facilitate transcriptional activation by ABR1. Finally, we compared gene expression profiles from wild-type and abr1 plants following inoculation with P. syringae, which suggested that the reduced susceptibility of abr1 mutants is due to the loss of a virulence target rather than an enhanced immune response. These data highlight ABR1 as a functionally important component at the host-pathogen interface.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq

We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

Genome-wide binding of SEPALLATA3 and AGAMOUS complexes determined by sequential DNA-affinity purification sequencing

The MADS transcription factors (TF), SEPALLATA3 (SEP3) and AGAMOUS (AG) are required for floral organ identity and floral meristem determinacy. While dimerization is obligatory for DNA binding, SEP3 and SEP3–AG also form tetrameric complexes. How homo and hetero-dimerization and tetramerization of MADS TFs affect genome-wide DNA-binding and gene regulation is not known. Using sequential DNA affinity purification sequencing (seq-DAP-seq), we determined genome-wide binding of SEP3 homomeric and SEP3–AG heteromeric complexes, including SEP3^Δtet-AG, a complex with a SEP3 splice variant, SEP3^Δtet, which is largely dimeric and SEP3–AG tetramer. SEP3 and SEP3–AG share numerous bound regions, however each complex bound unique sites, demonstrating that protein identity plays a role in DNA-binding. SEP3–AG and SEP3^Δtet-AG share a similar genome-wide binding pattern; however the tetrameric form could access new sites and demonstrated a global increase in DNA-binding affinity. Tetramerization exhibited significant cooperative binding with preferential distances between two sites, allowing efficient binding to regions that are poorly recognized by dimeric SEP3^Δtet-AG. By intersecting seq-DAP-seq with ChIP-seq and expression data, we identified unique target genes bound either in SEP3–AG seq-DAP-seq or in SEP3/AG ChIP-seq. Seq-DAP-seq is a versatile genome-wide technique and complements in vivo methods to identify putative direct regulatory targets.

GABA signaling in plants: targeting the missing pieces of the puzzle

The adaptation of plants to unstable environments relies on their ability to sense their surroundings and to generate and transmit corresponding signals to different parts of the plant to evoke changes necessary for optimizing growth and defense. Plants, like animals, contain a huge repertoire of intra- and intercellular signals, including organic and inorganic molecules. The occurrence of neurotransmitter-like signaling molecules in plants has been an intriguing field of research. Among these, γ-aminobutyric acid (GABA) was discovered in plants over half a century ago, and studies of its roles as a primary metabolite have been well documented, particularly in the context of stress responses. In contrast, evidence of the potential mechanism by which GABA acts as a signaling molecule in plants has only recently been reported. In spite of this breakthrough, the roles of GABA as a signaling molecule in plants have yet to be established and several aspects of the complexity of the GABA signaling system remain obscure. This review summarizes the uncertainties in GABA signaling in plants and suggests research directions and technologies that would help in answering unsolved questions.

Contrasted evolutionary trajectories of plant transcription factors

Because of their prominent roles in plant development, transcription factors (TF) play central roles as drivers of innovation in the evolution of the green lineage (viridiplantae). The advent of massive sequencing combined with comparative genetics/genomics allows a rigorous investigation of how TF families have contributed to plant diversification from charophyte algae to bryophytes to angiosperms. Here, we review recent progress on TF family reconstruction and the identification of distantly related TFs present throughout the evolutionary timeline from algae to angiosperms. These data provide examples of contrasting evolutionary trajectories of TF families and illustrate how conserved TFs adopt diverse roles over the course of evolution.

The MADS-box transcription factor PHERES1 controls imprinting in the endosperm by binding to domesticated transposons

MADS-box transcription factors (TFs) are ubiquitous in eukaryotic organisms and play major roles during plant development. Nevertheless, their function in seed development remains largely unknown. Here, we show that the imprinted Arabidopsis thaliana MADS-box TF PHERES1 (PHE1) is a master regulator of paternally expressed imprinted genes, as well as of non-imprinted key regulators of endosperm development. PHE1 binding sites show distinct epigenetic modifications on maternal and paternal alleles, correlating with parental-specific transcriptional activity. Importantly, we show that the CArG-box-like DNA-binding motifs that are bound by PHE1 have been distributed by RC/Helitron transposable elements. Our data provide an example of the molecular domestication of these elements which, by distributing PHE1 binding sites throughout the genome, have facilitated the recruitment of crucial endosperm regulators into a single transcriptional network.

Inference of plant gene regulatory networks using data-driven methods: A practical overview

Transcriptional regulation is a complex and dynamic process that plays a vital role in plant growth and development. A key component in the regulation of genes is transcription factors (TFs), which coordinate the transcriptional control of gene activity. A gene regulatory network (GRN) is a collection of regulatory interactions between TFs and their target genes. The accurate delineation of GRNs offers a significant contribution to our understanding about how plant cells are organized and function, and how individual genes are regulated in various conditions, organs or cell types. During the past decade, important progress has been made in the identification of GRNs using experimental and computational approaches. However, a detailed overview of available platforms supporting the analysis of GRNs in plants is missing. Here, we review current databases, platforms and tools that perform data-driven analyses of gene regulation in Arabidopsis. The platforms are categorized into two sections, 1) promoter motif analysis tools that use motif mapping approaches to find TF motifs in the regulatory sequences of genes of interest and 2) network analysis tools that identify potential regulators for a set of input genes using a range of data types in order to generate GRNs. We discuss the diverse datasets integrated and highlight the strengths and caveats of different platforms. Finally, we shed light on the limitations of the above approaches and discuss future perspectives, including the need for integrative approaches to unravel complex GRNs in plants.

Structural Basis for Plant MADS Transcription Factor Oligomerization

MADS transcription factors (TFs) are DNA binding proteins found in almost all eukaryotes that play essential roles in diverse biological processes. While present in animals and fungi as a small TF family, the family has dramatically expanded in plants over the course of evolution, with the model flowering plant, Arabidopsis thaliana, possessing over 100 type I and type II MADS TFs. All MADS TFs contain a core and highly conserved DNA binding domain called the MADS or M domain. Plant MADS TFs have diversified this domain with plant-specific auxiliary domains. Plant type I MADS TFs have a highly diverse and largely unstructured Carboxy-terminal (C domain), whereas type II MADS have added oligomerization domains, called Intervening (I domain) and Keratin-like (K domain), in addition to the C domain. In this mini review, we describe the overall structure of the type II "MIKC" type MADS TFs in plants, with a focus on the K domain, a critical oligomerization module. We summarize the determining factors for oligomerization and provide mechanistic insights on how secondary structural elements are required for oligomerization capability and specificity. Using MADS TFs that are involved in flower organ specification as an example, we provide case studies and homology modeling of MADS TFs complex formation. Finally, we highlight outstanding questions in the field.

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.

Gynoecium development: networks in Arabidopsis and beyond

Life has always found a way to preserve itself. One strategy that has been developed for this purpose is sexual reproduction. In land plants, the gynoecium is considered to be at the top of evolutionary innovation, since it has been a key factor in the success of the angiosperms. The gynoecium is composed of carpels with different tissues that need to develop and differentiate in the correct way. In order to control and guide gynoecium development, plants have adapted elements of pre-existing gene regulatory networks (GRNs) but new ones have also evolved. The GRNs can interact with internal factors (e.g. hormones and other metabolites) and external factors (e.g. mechanical signals and temperature) at different levels, giving robustness and flexibility to gynoecium development. Here, we review recent findings regarding the role of cytokinin-auxin crosstalk and the genes that connect these hormonal pathways during early gynoecium development. We also discuss some examples of internal and external factors that can modify GRNs. Finally, we make a journey through the flowering plant lineage to determine how conserved are these GRNs that regulate gynoecium and fruit development.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity

BACKGROUND

Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system.

RESULTS

The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae.

CONCLUSIONS

We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

Architecture of gene regulatory networks controlling flower development in Arabidopsis thaliana

Floral homeotic transcription factors (TFs) act in a combinatorial manner to specify the organ identities in the flower. However, the architecture and the function of the gene regulatory network (GRN) controlling floral organ specification is still poorly understood. In particular, the interconnections of homeotic TFs, microRNAs (miRNAs) and other factors controlling organ initiation and growth have not been studied systematically so far. Here, using a combination of genome-wide TF binding, mRNA and miRNA expression data, we reconstruct the dynamic GRN controlling floral meristem development and organ differentiation. We identify prevalent feed-forward loops (FFLs) mediated by floral homeotic TFs and miRNAs that regulate common targets. Experimental validation of a coherent FFL shows that petal size is controlled by the SEPALLATA3-regulated miR319/TCP4 module. We further show that combinatorial DNA-binding of homeotic factors and selected other TFs is predictive of organ-specific patterns of gene expression. Our results provide a valuable resource for studying molecular regulatory processes underlying floral organ specification in plants.

Homeotic transcription factors and miRNAs promote floral organ specification. Here Chen et al. reconstruct gene regulatory networks in Arabidopsis flowers and find evidence for feed forward loops between transcription factors, miRNAs and their targets that determine organ-specific gene expression.

The HTPmod Shiny application enables modeling and visualization of large-scale biological data

The wave of high-throughput technologies in genomics and phenomics are enabling data to be generated on an unprecedented scale and at a reasonable cost. Exploring the large-scale data sets generated by these technologies to derive biological insights requires efficient bioinformatic tools. Here we introduce an interactive, open-source web application (HTPmod) for high-throughput biological data modeling and visualization. HTPmod is implemented with the Shiny framework by integrating the computational power and professional visualization of R and including various machine-learning approaches. We demonstrate that HTPmod can be used for modeling and visualizing large-scale, high-dimensional data sets (such as multiple omics data) under a broad context. By reinvestigating example data sets from recent studies, we find not only that HTPmod can reproduce results from the original studies in a straightforward fashion and within a reasonable time, but also that novel insights may be gained from fast reinvestigation of existing data by HTPmod.

Dijun Chen et al. present HTPmod, a Shiny web application for modeling and visualization of large-scale genomic and phenomic datasets. The authors show that HTPmod can quickly reproduce analyses of high-throughput biological datasets and produce publication-quality figures.

The floral homeotic protein SEPALLATA3 recognizes target DNA sequences by shape readout involving a conserved arginine residue in the MADS-domain

SEPALLATA3 of Arabidopsis thaliana is a MADS-domain transcription factor (TF) and a key regulator of flower development. MADS-domain proteins bind to sequences termed 'CArG-boxes' [consensus 5'-CC(A/T)₆ GG-3']. Because only a fraction of the CArG-boxes in the Arabidopsis genome are bound by SEPALLATA3, more elaborate principles have to be discovered to better understand which features turn CArG-boxes into genuine recognition sites. Here, we investigate to what extent the shape of the DNA is involved in a 'shape readout' that contributes to the binding of SEPALLATA3. We determined in vitro binding affinities of SEPALLATA3 to DNA probes that all contain the CArG-box motif, but differ in their predicted DNA shape. We found that binding affinity correlates well with a narrow minor groove of the DNA. Substitution of canonical bases with non-standard bases supports the hypothesis of minor groove shape readout by SEPALLATA3. Analysis of mutant SEPALLATA3 proteins further revealed that a highly conserved arginine residue, which is expected to contact the DNA minor groove, contributes significantly to the shape readout. Our studies show that the specific recognition of cis-regulatory elements by a plant MADS-domain TF, and by inference probably also of other TFs of this type, heavily depends on shape readout mechanisms.

Synergistic Binding of bHLH Transcription Factors to the Promoter of the Maize NADP-ME Gene Used in C4 Photosynthesis Is Based on an Ancient Code Found in the Ancestral C3 State

C4 photosynthesis has evolved repeatedly from the ancestral C3 state to generate a carbon concentrating mechanism that increases photosynthetic efficiency. This specialized form of photosynthesis is particularly common in the PACMAD clade of grasses, and is used by many of the world's most productive crops. The C4 cycle is accomplished through cell-type-specific accumulation of enzymes but cis-elements and transcription factors controlling C4 photosynthesis remain largely unknown. Using the NADP-Malic Enzyme (NADP-ME) gene as a model we tested whether mechanisms impacting on transcription in C4 plants evolved from ancestral components found in C3 species. Two basic Helix-Loop-Helix (bHLH) transcription factors, ZmbHLH128 and ZmbHLH129, were shown to bind the C4NADP-ME promoter from maize. These proteins form heterodimers and ZmbHLH129 impairs trans-activation by ZmbHLH128. Electrophoretic mobility shift assays indicate that a pair of cis-elements separated by a seven base pair spacer synergistically bind either ZmbHLH128 or ZmbHLH129. This pair of cis-elements is found in both C3 and C4 Panicoid grass species of the PACMAD clade. Our analysis is consistent with this cis-element pair originating from a single motif present in the ancestral C3 state. We conclude that C4 photosynthesis has co-opted an ancient C3 regulatory code built on G-box recognition by bHLH to regulate the NADP-ME gene. More broadly, our findings also contribute to the understanding of gene regulatory networks controlling C4 photosynthesis.

MADS transcription factors cooperate: complexities of complex formation

This article comments on: Rümpler F, Theißen G, Melzer R. 2018. A conserved leucine zipper-like motif accounts for strong tetramerization capabilities of SEPALLATA-like MADS-domain transcription factors. Journal of Experimental Botany 69, 1943–1954.

MADS-box genes and crop domestication: the jack of all traits

MADS-box genes are key regulators of virtually every aspect of plant reproductive development. They play especially prominent roles in flowering time control, inflorescence architecture, floral organ identity determination, and seed development. The developmental and evolutionary importance of MADS-box genes is widely acknowledged. However, their role during flowering plant domestication is less well recognized. Here, we provide an overview illustrating that MADS-box genes have been important targets of selection during crop domestication and improvement. Numerous examples from a diversity of crop plants show that various developmental processes have been shaped by allelic variations in MADS-box genes. We propose that new genomic and genome editing resources provide an excellent starting point for further harnessing the potential of MADS-box genes to improve a variety of reproductive traits in crops. We also suggest that the biophysics of MADS-domain protein-protein and protein-DNA interactions, which is becoming increasingly well characterized, makes them especially suited to exploit coding sequence variations for targeted breeding approaches.

Master Regulatory Transcription Factors in Plant Development: A Blooming Perspective

Transcription factors that trigger major developmental decisions in plants and animals are termed "master regulators". Such master regulators are classically seen as acting on the top of a regulatory hierarchy that determines a complete developmental program, and they usually encode transcription factors. Here, we introduce master regulators of flowering time and flower development as examples to show how analysis of molecular interactions and gene-regulatory networks in plants has changed our view on the molecular mechanisms by which these factors control developmental processes. A picture has emerged that emphasizes a complex combinatorial interplay in determining cell-type transcriptional programs, and a high level of feedback control. The expression of master regulators themselves is usually regulated by multiple factors integrating environmental and endogenous spatiotemporal cues. Master regulatory transcription factors regulate gene expression by different mechanisms, including modifications in chromatin status in the bound regions. A poorly understood phenomenon is how developmental master regulators exert functions in different cell- and organ types. This is especially relevant for those factors that have important functions in several developmental processes.