MADS transcription factors cooperate: complexities of complex formation

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.

Genome-wide identification and characterization of TCP gene family in Dendrobium nobile and their role in perianth development

TCP is a widely distributed, essential plant transcription factor that regulates plant growth and development. An in-depth study of TCP genes in Dendrobium nobile, a crucial parent in genetic breeding and an excellent model material to explore perianth development in Dendrobium, has not been conducted. We identified 23 DnTCP genes unevenly distributed across 19 chromosomes and classified them as Class I PCF (12 members), Class II: CIN (10 members), and CYC/TB1 (1 member) based on the conserved domain and phylogenetic analysis. Most DnTCPs in the same subclade had similar gene and motif structures. Segmental duplication was the predominant duplication event for TCP genes, and no tandem duplication was observed. Seven genes in the CIN subclade had potential miR319 and -159 target sites. Cis-acting element analysis showed that most DnTCP genes contained many developmental stress-, light-, and phytohormone-responsive elements in their promoter regions. Distinct expression patterns were observed among the 23 DnTCP genes, suggesting that these genes have diverse regulatory roles at different stages of perianth development or in different organs. For instance, DnTCP4 and DnTCP18 play a role in early perianth development, and DnTCP5 and DnTCP10 are significantly expressed during late perianth development. DnTCP17, 20, 21, and 22 are the most likely to be involved in perianth and leaf development. DnTCP11 was significantly expressed in the gynandrium. Specially, MADS-specific binding sites were present in most DnTCP genes putative promoters, and two Class I DnTCPs were in the nucleus and interacted with each other or with the MADS-box. The interactions between TCP and the MADS-box have been described for the first time in orchids, which broadens our understanding of the regulatory network of TCP involved in perianth development in orchids.

Functional conservation and divergence of SEPALLATA-like genes in floral development in Cymbidium sinense

Cymbidium sinense is one of the most important traditional Chinese Orchids due to its unique and highly ornamental floral organs. Although the ABCDE model for flower development is well-established in model plant species, the precise roles of these genes in C. sinense are not yet fully understood. In this study, four SEPALLATA-like genes were isolated and identified from C. sinense. CsSEP1 and CsSEP3 were grouped into the AGL9 clade, while CsSEP2 and CsSEP4 were included in the AGL2/3/4 clade. The expression pattern of CsSEP genes showed that they were significantly accumulated in reproductive tissues and expressed during flower bud development but only mildly detected or even undetected in vegetative organs. Subcellular localization revealed that CsSEP1 and CsSEP4 were localized to the nucleus, while CsSEP2 and CsSEP3 were located at the nuclear membrane. Promoter sequence analysis predicted that CsSEP genes contained a number of hormone response elements (HREs) and MADS-box binding sites. The early flowering phenotype observed in transgenic Arabidopsis plants expressing four CsSEP genes, along with the expression profiles of endogenous genes, such as SOC1, LFY, AG, FT, SEP3 and TCPs, in both transgenic Arabidopsis and C. sinense protoplasts, suggested that the CsSEP genes played a regulatory role in the flowering transition by influencing downstream genes related to flowering. However, only transgenic plants overexpressing CsSEP3 and CsSEP4 caused abnormal phenotypes of floral organs, while CsSEP1 and CsSEP2 had no effect on floral organs. Protein-protein interaction assays indicated that CsSEPs formed a protein complex with B-class CsAP3-2 and CsSOC1 proteins, affecting downstream genes to regulate floral organs and flowering time. Our findings highlighted both the functional conservation and divergence of SEPALLATA-like genes in C. sinense floral development. These results provided a valuable foundation for future studies of the molecular network underlying floral development in C. sinense.

Double DAP-seq uncovered synergistic DNA binding of interacting bZIP transcription factors

Many eukaryotic transcription factors (TF) form homodimer or heterodimer complexes to regulate gene expression. Dimerization of BASIC LEUCINE ZIPPER (bZIP) TFs are critical for their functions, but the molecular mechanism underlying the DNA binding and functional specificity of homo- versus heterodimers remains elusive. To address this gap, we present the double DNA Affinity Purification-sequencing (dDAP-seq) technique that maps heterodimer binding sites on endogenous genomic DNA. Using dDAP-seq we profile twenty pairs of C/S1 bZIP heterodimers and S1 homodimers in Arabidopsis and show that heterodimerization significantly expands the DNA binding preferences of these TFs. Analysis of dDAP-seq binding sites reveals the function of bZIP9 in abscisic acid response and the role of bZIP53 heterodimer-specific binding in seed maturation. The C/S1 heterodimers show distinct preferences for the ACGT elements recognized by plant bZIPs and motifs resembling the yeast GCN4 cis-elements. This study demonstrates the potential of dDAP-seq in deciphering the DNA binding specificities of interacting TFs that are key for combinatorial gene regulation.

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.

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.

Utilizing MIKC-type MADS-box protein SOC1 for yield potential enhancement in maize

Overexpression of Zea mays SOC gene promotes flowering, reduces plant height, and leads to no reduction in grain production per plant, suggesting enhanced yield potential, at least, through increasing planting density. MIKC-type MADS-box gene SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) is an integrator conserved in the plant flowering pathway. In this study, the maize SOC1 (ZmSOC1) gene was cloned and overexpressed in transgenic maize Hi-II genotype. The T₀ plants were backcrossed with nontransgenic inbred B73 to produce first generation backcross (BC₁) seeds. Phenotyping of both transgenic and null segregant (NT) BC₁ plants was conducted in three independent experiments. The BC₁ transgenic plants showed new attributes such as increased vegetative growth, accelerated flowering time, reduced overall plant height, and increased grain weight. Second generation backcross (BC₂) plants were evaluated in the field using two planting densities. Compared to BC₂ NT plants, BC₂ transgenic plants, were 12-18% shorter, flowered 5 days earlier, and showed no reduction in grain production per plant and an increase in fat, starch, and simple sugars in the grain. Transcriptome comparison in young leaves of 56-day-old BC₁ plants revealed that the overexpressed ZmSOC1 resulted in 107 differentially expressed genes. The upregulated transcription factor DNA BINDING WITH ONE FINGER 5.4 (DOF5.4) was among the genes responsible for the reduced plant height. Modulating expression of SOC1 opens a new and effective approach to promote flowering and reduce plant height, which may have potential to enhance crop yield and improve grain quality.

Functional conservation and divergence of SEPALLATA-like genes in the development of two-type florets in marigold

Marigold (Tagetes erecta), as one member of Asteraceae family, bears a typical capitulum with two morphologically distinct florets. The SEPALLATA genes are involved in regulating the floral meristem determinacy, organ identity, fruit maturation, seed formation, and plant architecture. Here, five SEP-like genes were cloned and identified from marigold. Sequence alignment and phylogenetic analysis demonstrated that TeSEP3-1, TeSEP3-2, and TeSEP3-3 proteins were grouped into SEP3 clade, and TeSEP1 and TeSEP4 proteins were clustered into SEP1/2/4 clade. Quantitative real-time PCR analysis revealed that TeSEP1 and TeSEP3-3 were broadly expressed in floral organs, and that TeSEP3-2 and TeSEP4 were mainly expressed in pappus and corollas, while TeSEP3-1 was mainly expressed in two inner whorls. Ectopic expression of TeSEP1, TeSEP3-2, TeSEP3-3, and TeSEP4 in arabidopsis and tobacco resulted in early flowering. However, overexpression of TeSEP3-1 in arabidopsis and tobacco caused no visible phenotypic changes. Notably, overexpression of TeSEP4 in tobacco decreased the number of petals and stamens. Overexpression of TeSEP1 in tobacco led to longer sepals and simpler inflorescence architecture. The comprehensive pairwise interaction analysis suggested that TeSEP proteins had a broad interaction with class A, C, D, E proteins to form dimers. The yeast three-hybrid analysis suggested that in ternary complexes, class B proteins interacted with TeSEP3 by forming heterodimer TePI-TeAP3-2. The regulatory network analysis of MADS-box genes in marigold further indicated that TeSEP proteins played a "glue" role in regulating floral organ development, implying functional conservation and divergence of MADS box genes in regulating two-type floret developments. This study provides an insight into the formation mechanism of floral organs of two-type florets, thus broadening our knowledge of the genetic basis of flower evolution.

Meeting at the DNA: Specifying Cytokinin Responses through Transcription Factor Complex Formation

Cytokinin is a plant hormone regulating numerous biological processes. Its diverse functions are realized through the expression control of specific target genes. The transcription of the immediate early cytokinin target genes is regulated by type-B response regulator proteins (RRBs), which are transcription factors (TFs) of the Myb family. RRB activity is controlled by phosphorylation and protein degradation. Here, we focus on another step of regulation, the interaction of RRBs among each other or with other TFs to form active or repressive TF complexes. Several examples in Arabidopsis thaliana illustrate that RRBs form homodimers or complexes with other TFs to specify the cytokinin response. This increases the variability of the output response and provides opportunities of crosstalk between the cytokinin signaling pathway and other cellular signaling pathways. We propose that a targeted approach is required to uncover the full extent and impact of RRB interaction with other TFs.

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.