SEPALLATA-driven MADS transcription factor tetramerization is required for inner whorl floral organ development.

MADS transcription factors are master regulators of plant reproduction and flower development. The SEPALLATA (SEP) subfamily of MADS transcription factors is required for the development of floral organs and plays roles in inflorescence architecture and development of the floral meristem. SEPALLATAs act as organizers of MADS complexes, forming both heterodimers and heterotetramers in vitro. To date, the MADS complexes characterized in angiosperm floral organ development contain at least one SEPALLATA protein. Whether DNA-binding by SEPALLATA-containing dimeric MADS complexes is sufficient for launching floral organ identity programs, however, is not clear as only defects in floral meristem determinacy were observed in tetramerization--impaired SEPALLATA mutant proteins. Here, we used a combination of genome-wide binding studies, high resolution structural studies of the SEP3/AGAMOUS (AG) tetramerization domain, structure-based mutagenesis and complementation experiments in Arabidopsis (Arabidopsis thaliana) sep1 sep2 sep3 and sep1 sep2 sep3 ag-4 plants transformed with versions of SEP3 encoding tetramerization mutants. We demonstrate that while SEP3 heterodimers can bind DNA both in vitro and in vivo and recognize the majority of SEP3 wild-type binding sites genome-wide, tetramerization is required not only for floral meristem determinacy, but also for floral organ identity in the second, third and fourth whorls.

The ALOG domain defines a family of plant-specific transcription factors acting during Arabidopsis flower development.

The ALOG (Arabidopsis LIGHT-DEPENDENT SHORT HYPOCOTYLS 1 (LSH1) and Oryza G1) proteins are conserved plant-specific Transcription Factors (TFs). They play critical roles in the development of various plant organs (meristems, inflorescences, floral organs, and nodules) from bryophytes to higher flowering plants. Despite the fact that the first members of this family were originally discovered in Arabidopsis, their role in this model plant has remained poorly characterized. Moreover, how these transcriptional regulators work at the molecular level is unknown. Here, we study the redundant function of the ALOG proteins LSH1,3,4 from Arabidopsis. We uncover their role in the repression of bract development and position them within a gene regulatory network controlling this process and involving the floral regulators LEAFY, BLADE-ON-PETIOLE, and PUCHI. Next, using in vitro genome-wide studies, we identified the conserved DNA motif bound by ALOG proteins from evolutionarily distant species (the liverwort Marchantia polymorpha and the flowering plants Arabidopsis, tomato, and rice). Resolution of the crystallographic structure of the ALOG DNA-binding domain in complex with DNA revealed the domain is a four-helix bundle with a disordered NLS and a zinc ribbon insertion between helices 2 and 3. The majority of DNA interactions are mediated by specific contacts made by the third alpha helix and the NLS. Taken together, this work provides the biochemical and structural basis for DNA-binding specificity of an evolutionarily conserved TF family and reveals its role as a key player in Arabidopsis flower development.

Plant-TFClass: a structural classification for plant transcription factors.

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TF binding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

The F-box protein UFO controls flower development by redirecting the master transcription factor LEAFY to new cis-elements.

In angiosperms, flower development requires the combined action of the transcription factor LEAFY (LFY) and the ubiquitin ligase adaptor F-box protein, UNUSUAL FLORAL ORGANS (UFO), but the molecular mechanism underlying this synergy has remained unknown. Here we show in transient assays and stable transgenic plants that the connection to ubiquitination pathways suggested by the UFO F-box domain is mostly dispensable. On the basis of biochemical and genome-wide studies, we establish that UFO instead acts by forming an active transcriptional complex with LFY at newly discovered regulatory elements. Structural characterization of the LFY-UFO-DNA complex by cryo-electron microscopy further demonstrates that UFO performs this function by directly interacting with both LFY and DNA. Finally, we propose that this complex might have a deep evolutionary origin, largely predating flowering plants. This work reveals a unique mechanism of an F-box protein directly modulating the DNA binding specificity of a master transcription factor.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

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.

Cauliflower fractal forms arise from perturbations of floral gene networks.

Throughout development, plant meristems regularly produce organs in defined spiral, opposite, or whorl patterns. Cauliflowers present an unusual organ arrangement with a multitude of spirals nested over a wide range of scales. How such a fractal, self-similar organization emerges from developmental mechanisms has remained elusive. Combining experimental analyses in an Arabidopsis thaliana cauliflower-like mutant with modeling, we found that curd self-similarity arises because the meristems fail to form flowers but keep the "memory" of their transient passage in a floral state. Additional mutations affecting meristem growth can induce the production of conical structures reminiscent of the conspicuous fractal Romanesco shape. This study reveals how fractal-like forms may emerge from the combination of key, defined perturbations of floral developmental programs and growth dynamics.

The LEAFY floral regulator displays pioneer transcription factor properties.

Pioneer transcription factors (TFs) are a special category of TFs with the capacity to bind to closed chromatin regions in which DNA is wrapped around histones and may be highly methylated. Subsequently, pioneer TFs are able to modify the chromatin state to initiate gene expression. In plants, LEAFY (LFY) is a master floral regulator and has been suggested to act as a pioneer TF in Arabidopsis. Here, we demonstrate that LFY is able to bind both methylated and non-methylated DNA using a combination of in vitro genome-wide binding experiments and structural modeling. Comparisons between regions bound by LFY in vivo and chromatin accessibility data suggest that a subset of LFY bound regions is occupied by nucleosomes. We confirm that LFY is able to bind nucleosomal DNA in vitro using reconstituted nucleosomes. Finally, we show that constitutive LFY expression in seedling tissues is sufficient to induce chromatin accessibility in the LFY direct target genes APETALA1 and AGAMOUS. Taken together, our study suggests that LFY possesses key pioneer TF features that contribute to launching the floral gene expression program.

OsFD4 promotes the rice floral transition via florigen activation complex formation in the shoot apical meristem.

In rice, the florigens Heading Date 3a (Hd3a) and Rice Flowering Locus T 1 (RFT1), OsFD-like basic leucine zipper (bZIP) transcription factors, and Gf14 proteins assemble into florigen activation/repressor complexes (FACs/FRCs), which regulate transition to flowering in leaves and apical meristem. Only OsFD1 has been described as part of complexes promoting flowering at the meristem, and little is known about the role of other bZIP transcription factors, the combinatorial complexity of FAC formation, and their DNA-binding properties. Here, we used mutant analysis, protein-protein interaction assays and DNA affinity purification (DAP) sequencing coupled to in silico prediction of binding syntaxes to study several bZIP proteins that assemble into FACs or FRCs. We identified OsFD4 as a component of a FAC promoting flowering at the shoot apical meristem, downstream of OsFD1. The osfd4 mutants are late flowering and delay expression of genes promoting inflorescence development. Protein-protein interactions indicate an extensive network of contacts between several bZIPs and Gf14 proteins. Finally, we identified genomic regions bound by bZIPs with promotive and repressive effects on flowering. We conclude that distinct bZIPs orchestrate floral induction at the meristem and that FAC formation is largely combinatorial. While binding to the same consensus motif, their DNA-binding syntax is different, suggesting discriminatory functions.

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.

Impact of whole genome triplication on the evolutionary history and the functional dynamics of regulatory genes involved in Brassica self-incompatibility signalling pathway.

Polyploidy or whole genome duplication is a frequent and recurrent phenomenon in flowering plants that has played a major role in their diversification, adaptation and speciation. The adaptive success of polyploids relates to the different evolutionary fates of duplicated genes. In this study, we explored the impact of the whole genome triplication (WGT) event in the Brassiceae tribe on the genes involved in the self-incompatibility (SI) signalling pathway, a mechanism allowing recognition and rejection of self-pollen in hermaphrodite plants. By taking advantage of the knowledge acquired on this pathway as well as of several reference genomes in Brassicaceae species, we determined copy number of the different genes involved in this pathway and investigated their structural and functional evolutionary dynamics. We could infer that whereas most genes involved in the SI signalling returned to single copies after the WGT event (i.e. ARC1, JDP1, THL1, THL2, Exo70A01) in diploid Brassica species, a few were retained in duplicated (GLO1 and PLDα) or triplicated copies (MLPK). We also carefully studied the gene structure of these latter duplicated genes (including the conservation of functional domains and active sites) and tested their transcription in the stigma to identify which copies seem to be involved in the SI signalling pathway. By taking advantage of these analyses, we then explored the putative origin of a contrasted SI phenotype between two Brassica rapa varieties that have been fully sequenced and shared the same S-allele (S60).

Cytonuclear interactions remain stable during allopolyploid evolution despite repeated whole-genome duplications in Brassica.

Several plastid macromolecular protein complexes are encoded by both nuclear and plastid genes. Therefore, cytonuclear interactions are held in place to prevent genomic conflicts that may lead to incompatibilities. Allopolyploidy resulting from hybridization and genome doubling of two divergent species can disrupt these fine-tuned interactions, as newly formed allopolyploid species confront biparental nuclear chromosomes with a uniparentally inherited plastid genome. To avoid any deleterious effects of unequal genome inheritance, preferential transcription of the plastid donor over the other donor has been hypothesized to occur in allopolyploids. We used Brassica as a model to study the effects of paleopolyploidy in diploid parental species, as well as the effects of recent and ancient allopolyploidy in Brassica napus, on genes implicated in plastid protein complexes. We first identified redundant nuclear copies involved in those complexes. Compared with cytosolic protein complexes and with genome-wide retention rates, genes involved in plastid protein complexes show a higher retention of genes in duplicated and triplicated copies. Those redundant copies are functional and are undergoing strong purifying selection. We then compared transcription patterns and sequences of those redundant gene copies between resynthesized allopolyploids and their diploid parents. The neopolyploids showed no biased subgenome expression or maternal homogenization via gene conversion, despite the presence of some non-synonymous substitutions between plastid genomes of parental progenitors. Instead, subgenome dominance was observed regardless of the maternal progenitor. Our results provide new insights on the evolution of plastid protein complexes that could be tested and generalized in other allopolyploid species.

Capturing Auxin Response Factors Syntax Using DNA Binding Models.

Auxin is a key hormone performing a wealth of functions throughout the life cycle of plants. It acts largely by regulating genes at the transcriptional level through a family of transcription factors called auxin response factors (ARFs). Even though all ARF monomers analyzed so far bind a similar DNA sequence, there is evidence that ARFs differ in their target genomic regions and regulated genes. Here, we report the use of position weight matrices (PWMs) to model ARF DNA binding specificity based on published DNA affinity purification sequencing (DAP-seq) data. We found that the genome binding of two ARFs (ARF2 and ARF5/Monopteros [MP]) differ largely because these two factors have different preferred ARF binding site (ARFbs) arrangements (orientation and spacing). We illustrated why PWMs are more versatile to reliably identify ARFbs than the widely used consensus sequences and demonstrated their power with biochemical experiments in the identification of the regulatory regions of IAA19, an well-characterized auxin-responsive gene. Finally, we combined gene regulation by auxin with ARF-bound regions and identified specific ARFbs configurations that are over-represented in auxin-upregulated genes, thus deciphering the ARFbs syntax functional for regulation. Our study provides a general method to exploit the potential of genome-wide DNA binding assays and to decode gene regulation.