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

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.

MADS8 is indispensable for female reproductive development at high ambient temperatures in cereal crops

Temperature is a major factor that regulates plant growth and phenotypic diversity. To ensure reproductive success at a range of temperatures, plants must maintain developmental stability of their sexual organs when exposed to temperature fluctuations. However, the mechanisms integrating plant floral organ development and temperature responses are largely unknown. Here, we generated barley and rice loss-of-function mutants in the SEPALLATA-like MADS-box gene MADS8. The mutants in both species form multiple carpels that lack ovules at high ambient temperatures. Tissue-specific markers revealed that HvMADS8 is required to maintain floral meristem determinacy and ovule initiation at high temperatures, and transcriptome analyses confirmed that temperature-dependent differentially expressed genes in Hvmads8 mutants predominantly associate with floral organ and meristem regulation. HvMADS8 temperature-responsive activity relies on increased binding to promoters of downstream targets, as revealed by a cleavage under targets and tagmentation (CUT&Tag) analysis. We also demonstrate that HvMADS8 directly binds to 2 orthologs of D-class floral homeotic genes to activate their expression. Overall, our findings revealed a new, conserved role for MADS8 in maintaining pistil number and ovule initiation in cereal crops, extending the known function of plant MADS-box proteins in floral organ regulation.

Antagonistic regulation of target genes by the SISTER OF TM3-JOINTLESS2 complex in tomato inflorescence branching

Inflorescence branch number is a yield-related trait controlled by cell fate determination in meristems. Two MADS-box transcription factors (TFs)-SISTER OF TM3 (STM3) and JOINTLESS 2 (J2)-have opposing regulatory roles in inflorescence branching. However, the mechanisms underlying their regulatory functions in inflorescence determinacy remain unclear. Here, we characterized the functions of these TFs in tomato (Solanum lycopersicum) floral meristem and inflorescence meristem (IM) through chromatin immunoprecipitation and sequencing analysis of their genome-wide occupancy. STM3 and J2 activate or repress the transcription of a set of common putative target genes, respectively, through recognition and binding to CArG box motifs. FRUITFULL1 (FUL1) is a shared putative target of STM3 and J2 and these TFs antagonistically regulate FUL1 in inflorescence branching. Moreover, STM3 physically interacts with J2 to mediate its cytosolic redistribution and restricts J2 repressor activity by reducing its binding to target genes. Conversely, J2 limits STM3 regulation of target genes by transcriptional repression of the STM3 promoter and reducing STM3-binding activity. Our study thus reveals an antagonistic regulatory relationship in which STM3 and J2 control tomato IM determinacy and branch number.

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.

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.

Excess nitrogen responsive HvMADS27 transcription factor controls barley root architecture by regulating abscisic acid level

Nitrogen (N) is an important element for plant growth and development. Although several studies have examined plants' response to N deficiency, studies on plants' response to excess N, which is common in fertilizer-based agrosystems, are limited. Therefore, the aim of this study was to examine the response of barley to excess N conditions, specifically the root response. Additionally, genomic mechanism of excess N response in barley was elucidated using transcriptomic technologies. The results of the study showed that barley MADS27 transcription factor was mainly expressed in the roots and its gene contained N-responsive cis-regulatory elements in the promoter region. Additionally, there was a significant decrease in HvMADS27 expression under excess N condition; however, its expression was not significantly affected under low N condition. Phenotypic analysis of the root system of HvMADS27 knockdown and overexpressing barley plants revealed that HvMADS27 regulates barley root architecture under excess N stress. Further analysis of wild-type (WT) and transgenic barley plants (hvmads27 kd and hvmads27 c-Myc OE) revealed that HvMADS27 regulates the expression of HvBG1 β-glucosidase, which in turn regulates abscisic acid (ABA) level in roots. Overall, the findings of this study showed that HvMADS27 expression is downregulated in barley roots under excess N stress, which induces HvBG1 expression, leading to the release of ABA from ABA-glucose conjugate, and consequent shortening of the roots.

MADS1 maintains barley spike morphology at high ambient temperatures

Temperature stresses affect plant phenotypic diversity. The developmental stability of the inflorescence, required for reproductive success, is tightly regulated by the interplay of genetic and environmental factors. However, the mechanisms underpinning how plant inflorescence architecture responds to temperature are largely unknown. We demonstrate that the barley SEPALLATA MADS-box protein HvMADS1 is responsible for maintaining an unbranched spike architecture at high temperatures, while the loss-of-function mutant forms a branched inflorescence-like structure. HvMADS1 exhibits increased binding to target promoters via A-tract CArG-box motifs, which change conformation with temperature. Target genes for high-temperature-dependent HvMADS1 activation are predominantly associated with inflorescence differentiation and phytohormone signalling. HvMADS1 directly regulates the cytokinin-degrading enzyme HvCKX3 to integrate temperature response and cytokinin homeostasis, which is required to repress meristem cell cycle/division. Our findings reveal a mechanism by which genetic factors direct plant thermomorphogenesis, extending the recognized role of plant MADS-box proteins in floral development.

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.

Multiple levels of crosstalk in hormone networks regulating plant defense

Plant hormones are essential for regulating the interactions between plants and their complex biotic and abiotic environments. Each hormone initiates a specific molecular pathway and these different hormone pathways are integrated in a complex network of synergistic, antagonistic and additive interactions. This inter-pathway communication is called hormone crosstalk. By influencing the immune network topology, hormone crosstalk is essential for tailoring plant responses to diverse microbes and insects in diverse environmental and internal contexts. Crosstalk provides robustness to the immune system but also drives specificity of induced defense responses against the plethora of biotic interactors. Recent advances in dry-lab and wet-lab techniques have greatly enhanced our understanding of the broad-scale effects of hormone crosstalk on immune network functioning and have revealed underlying principles of crosstalk mechanisms. Molecular studies have demonstrated that hormone crosstalk is modulated at multiple levels of regulation, such as by affecting protein stability, gene transcription and hormone homeostasis. These new insights into hormone crosstalk regulation of plant defense are reviewed here, with a focus on crosstalk acting on the jasmonic acid pathway in Arabidopsis thaliana, highlighting the transcription factors MYC2 and ORA59 as major targets for modulation by other hormones.

How Hormones and MADS-Box Transcription Factors Are Involved in Controlling Fruit Set and Parthenocarpy in Tomato

Fruit set is the earliest phase of fruit growth and represents the onset of ovary growth after successful fertilization. In parthenocarpy, fruit formation is less affected by environmental factors because it occurs in the absence of pollination and fertilization, making parthenocarpy a highly desired agronomic trait. Elucidating the genetic program controlling parthenocarpy, and more generally fruit set, may have important implications in agriculture, considering the need for crops to be adaptable to climate changes. Several phytohormones play an important role in the transition from flower to fruit. Further complexity emerges from functional analysis of floral homeotic genes. Some homeotic MADS-box genes are implicated in fruit growth and development, displaying an expression pattern commonly observed for ovary growth repressors. Here, we provide an overview of recent discoveries on the molecular regulatory gene network underlying fruit set in tomato, the model organism for fleshy fruit development due to the many genetic and genomic resources available. We describe how the genetic modification of components of this network can cause parthenocarpy, discussing the contribution of hormonal signals and MADS-box transcription factors.

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.

Landscape of DNA binding signatures of myocyte enhancer factor-2B reveals a unique interplay of base and shape readout

Myocyte enhancer factor-2B (MEF2B) has the unique capability of binding to its DNA target sites with a degenerate motif, while still functioning as a gene-specific transcriptional regulator. Identifying its DNA targets is crucial given regulatory roles exerted by members of the MEF2 family and MEF2B's involvement in B-cell lymphoma. Analyzing structural data and SELEX-seq experimental results, we deduced the DNA sequence and shape determinants of MEF2B target sites on a high-throughput basis in vitro for wild-type and mutant proteins. Quantitative modeling of MEF2B binding affinities and computational simulations exposed the DNA readout mechanisms of MEF2B. The resulting binding signature of MEF2B revealed distinct intricacies of DNA recognition compared to other transcription factors. MEF2B uses base readout at its half-sites combined with shape readout at the center of its degenerate motif, where A-tract polarity dictates nuances of binding. The predominant role of shape readout at the center of the core motif, with most contacts formed in the minor groove, differs from previously observed protein-DNA readout modes. MEF2B, therefore, represents a unique protein for studies of the role of DNA shape in achieving binding specificity. MEF2B-DNA recognition mechanisms are likely representative for other members of the MEF2 family.

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.

OsMADS25 regulates root system development via auxin signalling in rice

The phytohormone auxin is essential for root development in plants. OsMADS25 is a homologue of the AGL17-clade MADS-box genes in rice. Despite recent progress, the molecular mechanisms underlying the regulation of root development by OsMADS25 are not well known. It is unclear whether OsMADS25 regulates root development via auxin signalling. In this study, we examined the role of OsMADS25 in root development and characterized the signalling pathway through which OsMADS25 regulates root system development in rice. OsMADS25 overexpression significantly increased, but RNAi gene silencing repressed primary root (PR) length and lateral root (LR) density. Moreover, OsMADS25 promoted LR development in response to NO₃^- . Further study showed that OsMADS25 increased auxin accumulation in the root system by enhancing auxin biosynthesis and transport, while also reducing auxin degradation, therefore stimulating root development. More importantly, OsMADS25 was found to regulate OsIAA14 expression directly by binding to the CArG-box in the promoter region of OsIAA14, which encodes an Aux/indole acetic acid (IAA) transcriptional repressor of auxin signalling. Elevated auxin levels and decreased OsIAA14 expression might lead to reduced OsIAA14 protein accumulation, as a mechanism to regulate auxin signalling. Therefore, our findings reveal a molecular mechanism by which OsMADS25 modulates root system growth and development in rice, at least partilly, via Aux/IAA-based auxin signalling.

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.

Comparative analysis of binding patterns of MADS-domain proteins in Arabidopsis thaliana

BACKGROUND

Correct flower formation requires highly specific temporal and spatial regulation of gene expression. In Arabidopsis thaliana the majority of the master regulators that determine flower organ identity belong to the MADS-domain transcription factor family. The canonical DNA binding motif for this transcription factor family is the CArG-box, which has the consensus CC(A/T)₆GG. However, so far, a comprehensive analysis of MADS-domain binding patterns has not yet been performed.

RESULTS

Eight publicly available ChIP-seq datasets of MADS-domain proteins that regulate the floral transition and flower formation were analyzed. Surprisingly, the preferred DNA binding motif of each protein was a CArG-box with an NAA extension. Furthermore, motifs of other transcription factors were found in the vicinity of binding sites of MADS-domain transcription factors, suggesting that interaction of MADS-domain proteins with other transcription factors is important for target gene regulation. Finally, conservation of CArG-boxes between Arabidopsis ecotypes was assessed to obtain information about their evolutionary importance. CArG-boxes that fully matched the consensus were more conserved than other CArG-boxes, suggesting that the perfect CArG-box is evolutionary more important than other CArG-box variants.

CONCLUSION

Our analysis provides detailed insight into MADS-domain protein binding patterns. The results underline the importance of an extended version of the CArG-box and provide a first view on evolutionary conservation of MADS-domain protein binding sites in Arabidopsis ecotypes.

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.

Divergence of regulatory networks governed by the orthologous transcription factors FLC and PEP1 in Brassicaceae species

Genome-wide landscapes of transcription factor (TF) binding sites (BSs) diverge during evolution, conferring species-specific transcriptional patterns. The rate of divergence varies in different metazoan lineages but has not been widely studied in plants. We identified the BSs and assessed the effects on transcription of FLOWERING LOCUS C (FLC) and PERPETUAL FLOWERING 1 (PEP1), two orthologous MADS-box TFs that repress flowering and confer vernalization requirement in the Brassicaceae species Arabidopsis thaliana and Arabis alpina, respectively. We found that only 14% of their BSs were conserved in both species and that these contained a CArG-box that is recognized by MADS-box TFs. The CArG-box consensus at conserved BSs was extended compared with the core motif. By contrast, species-specific BSs usually lacked the CArG-box in the other species. Flowering-time genes were highly overrepresented among conserved targets, and their CArG-boxes were widely conserved among Brassicaceae species. Cold-regulated (COR) genes were also overrepresented among targets, but the cognate BSs and the identity of the regulated genes were usually different in each species. In cold, COR gene transcript levels were increased in flc and pep1-1 mutants compared with WT, and this correlated with reduced growth in pep1-1 Therefore, FLC orthologs regulate a set of conserved target genes mainly involved in reproductive development and were later independently recruited to modulate stress responses in different Brassicaceae lineages. Analysis of TF BSs in these lineages thus distinguishes widely conserved targets representing the core function of the TF from those that were recruited later in evolution.

MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution

The floral quartet model of floral organ specification poses that different tetramers of MIKC-type MADS-domain transcription factors control gene expression and hence the identity of floral organs during development. Here, we provide a brief history of the floral quartet model and review several lines of recent evidence that support the model. We also describe how the model has been used in contemporary developmental and evolutionary biology to shed light on enigmatic topics such as the origin of land and flowering plants. Finally, we suggest a novel hypothesis describing how floral quartet-like complexes may interact with chromatin during target gene activation and repression.

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.

Floral Organogenesis: When Knowing Your ABCs Is Not Enough

The use of new experimental approaches enhances the understanding of floral organogenesis.

DNA Shape Features Improve Transcription Factor Binding Site Predictions In Vivo

Interactions of transcription factors (TFs) with DNA comprise a complex interplay between base-specific amino acid contacts and readout of DNA structure. Recent studies have highlighted the complementarity of DNA sequence and shape in modeling TF binding in vitro. Here, we have provided a comprehensive evaluation of in vivo datasets to assess the predictive power obtained by augmenting various DNA sequence-based models of TF binding sites (TFBSs) with DNA shape features (helix twist, minor groove width, propeller twist, and roll). Results from 400 human ChIP-seq datasets for 76 TFs show that combining DNA shape features with position-specific scoring matrix (PSSM) scores improves TFBS predictions. Improvement has also been observed using TF flexible models and a machine-learning approach using a binary encoding of nucleotides in lieu of PSSMs. Incorporating DNA shape information is most beneficial for E2F and MADS-domain TF families. Our findings indicate that incorporating DNA sequence and shape information benefits the modeling of TF binding under complex in vivo conditions.

Genome-Wide Targets Regulated by the OsMADS1 Transcription Factor Reveals Its DNA Recognition Properties

OsMADS1 controls rice (Oryza sativa) floral fate and organ development. Yet, its genome-wide targets and the mechanisms underlying its role as a transcription regulator controlling developmental gene expression are unknown. We identify 3112 gene-associated OsMADS1-bound sites in the floret genome. These occur in the vicinity of transcription start sites, within gene bodies, and in intergenic regions. Majority of the bound DNA contained CArG motif variants or, in several cases, only A-tracts. Sequences flanking the binding peak had a higher AT nucleotide content, implying that broader DNA structural features may define in planta binding. Sequences for binding by other transcription factor families like MYC, AP2/ERF, bZIP, etc. are enriched in OsMADS1-bound DNAs. Target genes implicated in transcription, chromatin remodeling, cellular processes, and hormone metabolism were enriched. Combining expression data from OsMADS1 knockdown florets with these DNA binding data, a snapshot of a gene regulatory network was deduced where targets, such as AP2/ERF and bHLH transcription factors and chromatin remodelers form nodes. We show that the expression status of these nodal factors can be altered by inducing the OsMADS1-GR fusion protein and present a model for a regulatory cascade where the direct targets of OsMADS1, OsbHLH108/SPT, OsERF034, and OsHSF24, in turn control genes such as OsMADS32 and OsYABBY5 This cascade, with other similar relationships, cumulatively contributes to floral organ development. Overall, OsMADS1 binds to several regulatory genes and, probably in combination with other factors, controls a gene regulatory network that ensures rice floret development.

Molecular mechanisms of floral organ specification by MADS domain proteins

Flower development is a model system to understand organ specification in plants. The identities of different types of floral organs are specified by homeotic MADS transcription factors that interact in a combinatorial fashion. Systematic identification of DNA-binding sites and target genes of these key regulators show that they have shared and unique sets of target genes. DNA binding by MADS proteins is not based on 'simple' recognition of a specific DNA sequence, but depends on DNA structure and combinatorial interactions. Homeotic MADS proteins regulate gene expression via alternative mechanisms, one of which may be to modulate chromatin structure and accessibility in their target gene promoters.

Novel alleles of the VERNALIZATION1 genes in wheat are associated with modulation of DNA curvature and flexibility in the promoter region

BACKGROUND

In wheat, the vernalization requirement is mainly controlled by the VRN genes. Different species of hexaploid and tetraploid wheat are widely used as genetic source for new mutant variants and alleles for fundamental investigations and practical breeding programs. In this study, VRN-A1 and VRN-B1 were analysed for 178 accessions representing six tetraploid wheat species (Triticum dicoccoides, T. dicoccum, T. turgidum, T. polonicum, T. carthlicum, T. durum) and five hexaploid species (T. compactum, T. sphaerococcum, T. spelta, T. macha, T. vavilovii).

RESULTS

Novel allelic variants in the promoter region of VRN-A1 and VRN-B1 were identified based on the change in curvature and flexibility of the DNA molecules. The new variants of VRN-A1 (designated as Vrn-A1a.2, Vrn-A1b.2 - Vrn-A1b.6 and Vrn-A1i) were found to be widely distributed in hexaploid and tetraploid wheat, and in fact were predominant over the known VRN-A1 alleles. The greatest diversity of the new variants of VRN-B1 (designated as VRN-B1.f, VRN-B1.s and VRN-B1.m) was found in the tetraploid and some hexaploid wheat species. For the first time, minor differences within the sequence motif known as the VRN-box of VRN1 were correlated with wheat growth habit. Thus, vrn-A1b.3 and vrn-A1b.4 were revealed in winter wheat in contrast to Vrn-A1b.2, Vrn-A1b.5, Vrn-A1b.6 and Vrn-A1i. It was found that single nucleotide mutation in the VRN-box can influence the vernalization requirement and growth habit of wheat. Our data suggest that both the A-tract and C-rich segment within the VRN-box contribute to its functionality, and provide a new view of the hypothesised role of the VRN-box in regulating transcription of the VRN1 genes. Specifically, it is proposed that combination of mutations in this region can modulate vernalization sensitivity and flowering time of wheat.

CONCLUSIONS

New allelic variants of the VRN-A1 and VRN-B1 genes were identified in hexaploid and tetraploid wheat. Mutations in A-tract and C-rich segments within the VRN-box of VRN-A1 are associated with modulation of the vernalization requirement and flowering time. New allelic variants will be useful in fundamental investigations into the regulation of VRN1 expression, and provide a valuable genetic resource for practical breeding of wheat.

Evolution of DNA-Binding Sites of a Floral Master Regulatory Transcription Factor

Flower development is controlled by the action of key regulatory transcription factors of the MADS-domain family. The function of these factors appears to be highly conserved among species based on mutant phenotypes. However, the conservation of their downstream processes is much less well understood, mostly because the evolutionary turnover and variation of their DNA-binding sites (BSs) among plant species have not yet been experimentally determined. Here, we performed comparative ChIP (chromatin immunoprecipitation)-seq experiments of the MADS-domain transcription factor SEPALLATA3 (SEP3) in two closely related Arabidopsis species: Arabidopsis thaliana and A. lyrata which have very similar floral organ morphology. We found that BS conservation is associated with DNA sequence conservation, the presence of the CArG-box BS motif and on the relative position of the BS to its potential target gene. Differences in genome size and structure can explain that SEP3 BSs in A. lyrata can be located more distantly to their potential target genes than their counterparts in A. thaliana. In A. lyrata, we identified transposition as a mechanism to generate novel SEP3 binding locations in the genome. Comparative gene expression analysis shows that the loss/gain of BSs is associated with a change in gene expression. In summary, this study investigates the evolutionary dynamics of DNA BSs of a floral key-regulatory transcription factor and explores factors affecting this phenomenon.

Combinatorial activities of SHORT VEGETATIVE PHASE and FLOWERING LOCUS C define distinct modes of flowering regulation in Arabidopsis

Background

The initiation of flowering is an important developmental transition as it marks the beginning of the reproductive phase in plants. The MADS-box transcription factors (TFs) FLOWERING LOCUS C (FLC) and SHORT VEGETATIVE PHASE (SVP) form a complex to repress the expression of genes that initiate flowering in Arabidopsis. Both TFs play a central role in the regulatory network by conferring seasonal patterns of flowering. However, their interdependence and biological relevance when acting as a complex have not been extensively studied.

Results

We characterized the effects of both TFs individually and as a complex on flowering initiation using transcriptome profiling and DNA-binding occupancy. We find four major clusters regulating transcriptional responses, and that DNA binding scenarios are highly affected by the presence of the cognate partner. Remarkably, we identify genes whose regulation depends exclusively on simultaneous action of both proteins, thus distinguishing between the specificity of the SVP:FLC complex and that of each TF acting individually. The downstream targets of the SVP:FLC complex include a higher proportion of genes regulating floral induction, whereas those bound by either TF independently are biased towards floral development. Many genes involved in gibberellin-related processes are bound by the SVP:FLC complex, suggesting that direct regulation of gibberellin metabolism by FLC and SVP contributes to their effects on flowering.

Conclusions

The regulatory codes controlled by SVP and FLC were deciphered at the genome-wide level revealing substantial flexibility based on dependent and independent DNA binding that may contribute to variation and robustness in the regulation of flowering.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0597-1) contains supplementary material, which is available to authorized users.

How is FLC repression initiated by cold?

Vernalization is the promotion of flowering in response to prolonged exposure to low temperatures. In Arabidopsis, FLOWERING LOCUS C (FLC), a suppressor of flowering, is repressed by low temperatures but the mechanism leading to the initial decrease in FLC transcription remains a mystery. No mutants that block the repression of FLC at low temperatures have been identified to date. If the failure to identify such a mutant is assumed to imply that no such mutant exists, then it follows that the first response to the drop in temperature is physical, not genetic. In this Opinion article we propose that the drop in temperature first causes a simple change in the topology of the chromatin polymer, which in turn initiates the repression of FLC transcription.

Arabidopsis SEPALLATA proteins differ in cooperative DNA-binding during the formation of floral quartet-like complexes

The SEPALLATA (SEP) genes of Arabidopsis thaliana encode MADS-domain transcription factors that specify the identity of all floral organs. The four Arabidopsis SEP genes function in a largely yet not completely redundant manner. Here, we analysed interactions of the SEP proteins with DNA. All of the proteins were capable of forming tetrameric quartet-like complexes on DNA fragments carrying two sequence elements termed CArG-boxes. Distances between the CArG-boxes for strong cooperative DNA-binding were in the range of 4-6 helical turns. However, SEP1 also bound strongly to CArG-box pairs separated by smaller or larger distances, whereas SEP2 preferred large and SEP4 preferred small inter-site distances for binding. Cooperative binding of SEP3 was comparatively weak for most of the inter-site distances tested. All SEP proteins constituted floral quartet-like complexes together with the floral homeotic proteins APETALA3 (AP3) and PISTILLATA (PI) on the target genes AP3 and SEP3. Our results suggest an important part of an explanation for why the different SEP proteins have largely, but not completely redundant functions in determining floral organ identity: they may bind to largely overlapping, but not identical sets of target genes that differ in the arrangement and spacing of the CArG-boxes in their cis-regulatory regions.

Structural basis for the oligomerization of the MADS domain transcription factor SEPALLATA3 in Arabidopsis

In plants, MADS domain transcription factors act as central regulators of diverse developmental pathways. In Arabidopsis thaliana, one of the most central members of this family is SEPALLATA3 (SEP3), which is involved in many aspects of plant reproduction, including floral meristem and floral organ development. SEP3 has been shown to form homo and heterooligomeric complexes with other MADS domain transcription factors through its intervening (I) and keratin-like (K) domains. SEP3 function depends on its ability to form specific protein-protein complexes; however, the atomic level determinants of oligomerization are poorly understood. Here, we report the 2.5-Å crystal structure of a small portion of the intervening and the complete keratin-like domain of SEP3. The domains form two amphipathic alpha helices separated by a rigid kink, which prevents intramolecular association and presents separate dimerization and tetramerization interfaces comprising predominantly hydrophobic patches. Mutations to the tetramerization interface demonstrate the importance of highly conserved hydrophobic residues for tetramer stability. Atomic force microscopy was used to show SEP3-DNA interactions and the role of oligomerization in DNA binding and conformation. Based on these data, the oligomerization patterns of the larger family of MADS domain transcription factors can be predicted and manipulated based on the primary sequence.