A promiscuous intermediate underlies the evolution of LEAFY DNA binding specificity

A conserved ARF-DNA interface underlies auxin-triggered transcriptional response

Auxin Response Factor (ARF) plant transcription factors are the key effectors in auxin signaling. Their DNA-Binding Domain (DBD) contains a B3 domain that allows base-specific interactions with Auxin Response Elements (AuxREs) in DNA target sites. Land plants encode three phylogenetically distinct ARF classes: the closely related A- and B-classes have overlapping DNA binding properties, contrasting with the different DNA-binding properties of the divergent C-class ARFs. ARF DNA-binding divergence likely occurred early in the evolution of the gene family, but the molecular determinants underlying it remain unclear. Here, we show that the B3 DNA-binding residues are deeply conserved in ARFs, and variability within these is only present in tracheophytes, correlating with greatly expanded ARF families. Using the liverwort Marchantia polymorpha, we confirm the essential role of conserved DNA-contacting residues for ARF function. We further show that ARF B3-AuxRE interfaces are not mutation-tolerant, suggesting low evolvability that has led to the conservation of the B3-DNA interface between ARF classes. Our data support the almost complete interchangeability between A/B-class ARF B3 by performing interspecies domain swaps in M. polymorpha, even between ARF lineages that diverged over half a billion years ago. Our analysis further suggests that C-class ARF DNA-binding specificity diverged early during ARF evolution in a common streptophyte ancestor, followed by strong selection in A and B-class ARFs as part of a competition-based auxin response system.

Specificity landscapes of 40 R2R3-MYBs reveal how paralogs target different cis-elements by homodimeric binding

Paralogous transcription factors (TFs) frequently recognize highly similar DNA motifs. Homodimerization can help distinguish them according to their different dimeric configurations. Here, by studying R2R3-MYB TFs, we show that homodimerization can also directly change the recognized DNA motifs to distinguish between similar TFs. By high-throughput SELEX, we profiled the specificity landscape for 40 R2R3-MYBs of subfamily VIII and curated 833 motif models. The dimeric models show that homodimeric binding has evoked specificity changes for AtMYBs. Focusing on AtMYB2 as an example, we show that homodimerization has modified its specificity and allowed it to recognize additional cis-regulatory sequences that are different from the closely related CCWAA-box AtMYBs and are unique among all AtMYBs. Genomic sites described by the modified dimeric specificities of AtMYB2 are conserved in evolution and involved in AtMYB2-specific transcriptional activation. Collectively, this study provides rich data on sequence preferences of VIII R2R3-MYBs and suggests an alternative mechanism that guides closely related TFs to respective cis-regulatory sites.

Evolutionary trajectory of transcription factors and selection of targets for metabolic engineering

Transcription factors (TFs) provide potentially powerful tools for plant metabolic engineering as they often control multiple genes in a metabolic pathway. However, selecting the best TF for a particular pathway has been challenging, and the selection often relies significantly on phylogenetic relationships. Here, we offer examples where evolutionary relationships have facilitated the selection of the suitable TFs, alongside situations where such relationships are misleading from the perspective of metabolic engineering. We argue that the evolutionary trajectory of a particular TF might be a better indicator than protein sequence homology alone in helping decide the best targets for plant metabolic engineering efforts. This article is part of the theme issue 'The evolution of plant metabolism'.

Following the Evolutionary Paths of Dscam1 Proteins toward Highly Specific Homophilic Interactions

Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Differential gene expression during floral transition in pineapple

Pineapple (Ananas comosus var. comosus) and ornamental bromeliads are commercially induced to flower by treatment with ethylene or its analogs. The apex is transformed from a vegetative to a floral meristem and shows morphological changes in 8 to 10 days, with flowers developing 8 to 10 weeks later. During eight sampling stages ranging from 6 h to 8 days after treatment, 7961 genes were found to exhibit differential expression (DE) after the application of ethylene. In the first 3 days after treatment, there was little change in ethylene synthesis or in the early stages of the ethylene response. Subsequently, three ethylene response transcription factors (ERTF) were up-regulated and the potential gene targets were predicted to be the positive flowering regulator CONSTANS-like 3 (CO), a WUSCHEL gene, two APETALA1/FRUITFULL (AP1/FUL) genes, an epidermal patterning gene, and a jasmonic acid synthesis gene. We confirm that pineapple has lost the flowering repressor FLOWERING LOCUS C. At the initial stages, the SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was not significantly involved in this transition. Another WUSCHEL gene and a PHD homeobox transcription factor, though not apparent direct targets of ERTF, were up-regulated within a day of treatment, their predicted targets being the up-regulated CO, auxin response factors, SQUAMOSA, and histone H3 genes with suppression of abscisic acid response genes. The FLOWERING LOCUS T (FT), TERMINAL FLOWER (TFL), AGAMOUS-like APETELAR (AP2), and SEPETALA (SEP) increased rapidly within 2 to 3 days after ethylene treatment. Two FT genes were up-regulated at the apex and not at the leaf bases after treatment, suggesting that transport did not occur. These results indicated that the ethylene response in pineapple and possibly most bromeliads act directly to promote the vegetative to flower transition via APETALA1/FRUITFULL (AP1/FUL) and its interaction with SPL, FT, TFL, SEP, and AP2. A model based on AP2/ERTF DE and predicted DE target genes was developed to give focus to future research. The identified candidate genes are potential targets for genetic manipulation to determine their molecular role in flower transition.

Roles and evolution of four LEAFY homologs in floral patterning and leaf development in woodland strawberry

The plant-specific transcription factor LEAFY (LFY), generally maintained as a single-copy gene in most angiosperm species, plays critical roles in flower development. The woodland strawberry (Fragaria vesca) possesses four LFY homologs in the genome; however, their respective functions and evolution remain unknown. Here, we identified and validated that mutations in one of the four LFY homologs, FveLFYa, cause homeotic conversion of floral organs and reiterative outgrowth of ectopic flowers. In contrast to FveLFYa, FveLFYb/c/d appear dispensable under normal growth conditions, as fvelfyc mutants are indistinguishable from wild type and FveLFYb and FveLFYd are barely expressed. Transgenic analysis and yeast one-hybrid assay showed that FveLFYa and FveLFYb, but not FveLFYc and FveLFYd, are functionally conserved with AtLFY in Arabidopsis (Arabidopsis thaliana). Unexpectedly, LFY-binding site prediction and yeast one-hybrid assay revealed that the transcriptional links between LFY and the APETALA1 (AP1) promoter/the large AGAMOUS (AG) intron are missing in F. vesca, which is due to the loss of LFY-binding sites. The data indicate that mutations in cis-regulatory elements could contribute to LFY evolution. Moreover, we showed that FveLFYa is involved in leaf development, as approximately 30% of mature leaves have smaller or fewer leaflets in fvelfya. Phylogenetic analysis indicated that LFY homologs in Fragaria species may arise from recent duplication events in their common ancestor and are undergoing convergent gene loss. Together, these results provide insight into the role of LFY in flower and leaf development in strawberry and have important implications for the evolution of LFY.

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.

Interspecies transfer of RAMOSA1 orthologs and promoter cis sequences impacts maize inflorescence architecture

Grass inflorescences support floral structures that each bear a single grain, where variation in branch architecture directly impacts yield. The maize (Zea mays) RAMOSA1 (ZmRA1) transcription factor acts as a key regulator of inflorescence development by imposing branch meristem determinacy. Here, we show RA1 transcripts accumulate in boundary domains adjacent to spikelet meristems in sorghum (Sorghum bicolor, Sb) and green millet (Setaria viridis, Sv) inflorescences similar as in the developing maize tassel and ear. To evaluate the functional conservation of syntenic RA1 orthologs and promoter cis sequences in maize, sorghum, and setaria, we utilized interspecies gene transfer and assayed genetic complementation in a common inbred background by quantifying recovery of normal branching in highly ramified ra1-R mutants. A ZmRA1 transgene that includes endogenous upstream and downstream flanking sequences recovered normal tassel and ear branching in ra1-R. Interspecies expression of two transgene variants of the SbRA1 locus, modeled as the entire endogenous tandem duplication or just the nonframeshifted downstream copy, complemented ra1-R branching defects and induced unusual fasciation and branch patterns. The SvRA1 locus lacks conserved, upstream noncoding cis sequences found in maize and sorghum; interspecies expression of a SvRA1 transgene did not or only partially recovered normal inflorescence forms. Driving expression of the SvRA1 coding region by the ZmRA1 upstream region, however, recovered normal inflorescence morphology in ra1-R. These data leveraging interspecies gene transfer suggest that cis-encoded temporal regulation of RA1 expression is a key factor in modulating branch meristem determinacy that ultimately impacts grass inflorescence architecture.

Navigating the CoGe Online Software Suite for Polyploidy Research

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.

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.

How was apical growth regulated in the ancestral land plant? Insights from the development of non-seed plants

Land plant life cycles are separated into distinct haploid gametophyte and diploid sporophyte stages. Indeterminate apical growth evolved independently in bryophyte (moss, liverwort, and hornwort) and fern gametophytes, and tracheophyte (vascular plant) sporophytes. The extent to which apical growth in tracheophytes co-opted conserved gametophytic gene networks, or exploited ancestral sporophytic networks, is a long-standing question in plant evolution. The recent phylogenetic confirmation of bryophytes and tracheophytes as sister groups has led to a reassessment of the nature of the ancestral land plant. Here, we review developmental genetic studies of apical regulators and speculate on their likely evolutionary history.

ABERRANT PANICLE ORGANIZATION2 controls multiple steps in panicle formation through common direct-target genes

At the transition from vegetative to reproductive growth in rice (Oryza sativa), a developmental program change occurs, resulting in panicle (rice inflorescence) formation. The initial event of the transition is the change of the shoot apical meristem to an inflorescence meristem (IM), accompanied by a rapid increase in the meristem size. Suppression of leaf growth also occurs, resulting in the formation of bracts. The IM generates branch meristems (BMs), indeterminate meristems that reiteratively generate next-order meristems. All meristems eventually acquire a determinate spikelet meristem identity and terminate after producing a floret. ABERRANT PANICLE ORGANIZATION2 (APO2) is the rice ortholog of Arabidopsis (Arabidopsis thaliana) LEAFY (LFY), a plant-specific transcription factor (TF). APO2 is a positive regulator of panicle branch formation. Here, we show that APO2 is also required to increase the meristem size of the IM and suppress bract outgrowth. We identified genes directly and indirectly regulated by APO2 and identified APO2-binding sites. These analyses showed that APO2 directly controls known regulators of panicle development, including SQUAMOSA PROMOTER BINDING PROTEIN LIKE14 and NECK LEAF1. Furthermore, we revealed that a set of genes act as downstream regulators of APO2 in controlling meristem cell proliferation during reproductive transition, bract suppression, and panicle branch formation. Our findings indicate that APO2 acts as a master regulator of rice panicle development by regulating multiple steps in the reproductive transition through directly controlling a set of genes.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Temperature-dependent genetics of thermotolerance between yeast species

Many traits of industrial and basic biological interest arose long ago, and manifest now as fixed differences between a focal species and its reproductively isolated relatives. In these systems, extant individuals can hold clues to the mechanisms by which phenotypes evolved in their ancestors. We harnessed yeast thermotolerance as a test case for such molecular-genetic inferences. In viability experiments, we showed that extant Saccharomyces cerevisiae survived at temperatures where cultures of its sister species S. paradoxus died out. Then, focusing on loci that contribute to this difference, we found that the genetic mechanisms of high-temperature growth changed with temperature. We also uncovered an enrichment of low-frequency variants at thermotolerance loci in S. cerevisiae population sequences, suggestive of a history of non-neutral selective forces acting at these genes. We interpret these results in light of models of the evolutionary mechanisms by which the thermotolerance trait arose in the S. cerevisiae lineage. Together, our results and interpretation underscore the power of genetic approaches to explore how an ancient trait came to be.

Plant transcription factors - being in the right place with the right company

Transcriptional regulation underlies many of the growth and developmental processes that shape plants as well as their adaptation to their environment. Key to transcriptional control are transcription factors, DNA-binding proteins that serve two essential functions: to find the appropriate DNA contact sites in their target genes; and to recruit other proteins to execute transcriptional transactions. In recent years, protein structural, genomic, bioinformatic, and proteomic analyses have led to new insights into how these central functions are regulated. Here, we review new findings relating to plant transcription factor function and to their role in shaping transcription in the context of chromatin.

The Genetic Control of the Compound Leaf Patterning in Medicago truncatula

Simple and compound which are the two basic types of leaves are distinguished by the pattern of the distribution of blades on the petiole. Compared to simple leaves comprising a single blade, compound leaves have multiple blade units and exhibit more complex and diverse patterns of organ organization, and the molecular mechanisms underlying their pattern formation are receiving more and more attention in recent years. Studies in model legume Medicago truncatula have led to an improved understanding of the genetic control of the compound leaf patterning. This review is an attempt to summarize the current knowledge about the compound leaf morphogenesis of M. truncatula, with a focus on the molecular mechanisms involved in pattern formation. It also includes some comparisons of the molecular mechanisms between leaf morphogenesis of different model species and offers useful information for the molecular design of legume crops.

Secreted Phospholipases A2 - not just Enzymes: Revisited

Secreted phospholipases A₂ (sPLA₂s) participate in a very broad spectrum of biological processes through their enzymatic activity and as ligands for membrane and soluble receptors. The physiological roles of sPLA₂s as enzymes have been very well described, while their functions as ligands are still poorly known. Since the last overview of sPLA₂-binding proteins (sPLA₂-BPs) 10 years ago, several important discoveries have occurred in this area. New and more sensitive analytical tools have enabled the discovery of additional sPLA₂-BPs, which are presented and critically discussed here. The structural diversity of sPLA₂-BPs reveals sPLA₂s as very promiscuous proteins, and we offer some structural explanations for this nature that makes these proteins evolutionarily highly advantageous. Three areas of physiological engagement of sPLA₂-BPs have appeared most clearly: cellular transport and signalling, and regulation of the enzymatic activity of sPLA₂s. Due to the multifunctionality of sPLA₂s, they appear to be exceptional pharmacological targets. We reveal the potential to exploit interactions of sPLA₂s with other proteins in medical terms, for the development of original diagnostic and therapeutic procedures. We conclude this survey by suggesting the priority questions that need to be answered.

Evolution and expression of LEAFY genes in ferns and lycophytes

BACKGROUND

The LEAFY (LFY) transcription factors are present in algae and across land plants. The available expression and functional data of these genes in embryophytes suggest that LFY genes control a plethora of processes including the first zygotic cell division in bryophytes, shoot cell divisions of the gametophyte and sporophyte in ferns, cone differentiation in gymnosperms and floral meristem identity in flowering plants. However, their putative plesiomorphic role in plant reproductive transition in vascular plants remains untested.

RESULTS

We perform Maximum Likelihood (ML) phylogenetic analyses for the LFY gene lineage in embryophytes with expanded sampling in lycophytes and ferns. We recover the previously identified seed plant duplication that results in LEAFY and NEEDLY paralogs. In addition, we recover multiple species-specific duplications in ferns and lycophytes and large-scale duplications possibly correlated with the occurrence of whole genome duplication (WGD) events in Equisetales and Salviniales. To test putative roles in diverse ferns and lycophytes we perform LFY expression analyses in Adiantum raddianum, Equisetum giganteum and Selaginella moellendorffii. Our results show that LFY genes are active in vegetative and reproductive tissues, with higher expression in early fertile developmental stages and during sporangia differentiation.

CONCLUSIONS

Our data point to previously unrecognized roles of LFY genes in sporangia differentiation in lycophytes and ferns and suggests that functions linked to reproductive structure development are not exclusive to seed plant LFY homologs.

Transitions in development – an interview with Marie Monniaux

Marie Monniaux is a permanent CNRS researcher in the ‘Evo-devo of the flower’ group at the Laboratory for Plant Reproduction and Development (RDP) at the École normale supérieure (ENS) in Lyon, France. Marie uses Petunia to understand the development and evolution of the flower petal. We met Marie over Teams for a virtual chat about her career path, finding a permanent position and her ideas for the future.

HvLEAFY controls the early stages of floral organ specification and inhibits the formation of multiple ovaries in barley

Transition to the reproductive phase, inflorescence formation and flower development are crucial elements that ensure maximum reproductive success in a plant's life cycle. To understand the regulatory mechanisms underlying correct flower development in barley (Hordeum vulgare), we characterized the multiovary 5 (mov5.o) mutant. This mutant develops abnormal flowers that exhibit mosaic floral organs typified by multiple carpels at the total or partial expense of stamens. Genetic mapping positioned mov5 on the long arm of chromosome 2H, incorporating a region that encodes HvLFY, the barley orthologue of LEAFY from Arabidopsis. Sequencing revealed that, in mov5.o plants, HvLFY contains a single amino acid substitution in a highly conserved proline residue. CRISPR-mediated knockout of HvLFY replicated the mov5.o phenotype, suggesting that HvLFY_mov5 represents a loss of function allele. In heterologous assays, the HvLFY_mov5 polymorphism influenced protein-protein interactions and affinity for a putative binding site in the promoter of HvMADS58, a C-class MADS-box gene. Moreover, molecular analysis indicated that HvLFY interacts with HvUFO and regulates the expression of floral homeotic genes including HvMADS2, HvMADS4 and HvMADS16. Other distinct changes in expression differ from those reported in the rice LFY mutants apo2/rfl, suggesting that LFY function in the grasses is modulated in a species-specific manner. This pathway provides a key entry point for the study of LFY function and multiple ovary formation in barley, as well as cereal species in general.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Evolving origin-of-transfer sequences on staphylococcal conjugative and mobilizable plasmids-who's mimicking whom?

In Staphylococcus aureus, most multiresistance plasmids lack conjugation or mobilization genes for horizontal transfer. However, most are mobilizable due to carriage of origin-of-transfer (oriT) sequences mimicking those of conjugative plasmids related to pWBG749. pWBG749-family plasmids have diverged to carry five distinct oriT subtypes and non-conjugative plasmids have been identified that contain mimics of each. The relaxasome accessory factor SmpO, encoded by each conjugative plasmid, determines specificity for its cognate oriT. Here we characterized the binding of SmpO proteins to each oriT. SmpO proteins predominantly formed tetramers in solution and bound 5′-GNNNNC-3′ sites within each oriT. Four of the five SmpO proteins specifically bound their cognate oriT. An F7K substitution in pWBG749 SmpO switched oriT-binding specificity in vitro. In vivo, the F7K substitution reduced but did not abolish self-transfer of pWBG749. Notably, the substitution broadened the oriT subtypes that were mobilized. Thus, this substitution represents a potential evolutionary intermediate with promiscuous DNA-binding specificity that could facilitate a switch between oriT specificities. Phylogenetic analysis suggests pWBG749-family plasmids have switched oriT specificity more than once during evolution. We hypothesize the convergent evolution of oriT specificity in distinct branches of the pWBG749-family phylogeny reflects indirect selection pressure to mobilize plasmids carrying non-cognate oriT-mimics.

HAG1 and SWI3A/B control of male germ line development in P. patens suggests conservation of epigenetic reproductive control across land plants

KEY MESSAGE

Bryophytes as models to study the male germ line: loss-of-function mutants of epigenetic regulators HAG1 and SWI3a/b demonstrate conserved function in sexual reproduction. With the water-to-land transition, land plants evolved a peculiar haplodiplontic life cycle in which both the haploid gametophyte and the diploid sporophyte are multicellular. The switch between these phases was coined alternation of generations. Several key regulators that control the bauplan of either generation are already known. Analyses of such regulators in flowering plants are difficult due to the highly reduced gametophytic generation, and the fact that loss of function of such genes often is embryo lethal in homozygous plants. Here we set out to determine gene function and conservation via studies in bryophytes. Bryophytes are sister to vascular plants and hence allow evolutionary inferences. Moreover, embryo lethal mutants can be grown and vegetatively propagated due to the dominance of the bryophyte gametophytic generation. We determined candidates by selecting single copy orthologs that are involved in transcriptional control, and of which flowering plant mutants show defects during sexual reproduction, with a focus on the under-studied male germ line. We selected two orthologs, SWI3a/b and HAG1, and analyzed loss-of-function mutants in the moss P. patens. In both mutants, due to lack of fertile spermatozoids, fertilization and hence the switch to the diploid generation do not occur. Pphag1 additionally shows arrested male and impaired female gametangia development. We analyzed HAG1 in the dioecious liverwort M. polymorpha and found that in Mphag1 the development of gametangiophores is impaired. Taken together, we find that involvement of both regulators in sexual reproduction is conserved since the earliest divergence of land plants.

Plant design gets its details: Modulating plant architecture by phase transitions

Plants evolved different strategies to better adapt to the environmental conditions in which they live: the control of their body architecture and the timing of phase change are two important processes that can improve their fitness. As they age, plants undergo two major phase changes (juvenile to adult and adult to reproductive) that are a response to environmental and endogenous signals. These phase transitions are accompanied by alterations in plant morphology and also by changes in physiology and the behavior of gene regulatory networks. Six main pathways involving environmental and endogenous cues that crosstalk with each other have been described as responsible for the control of plant phase transitions: the photoperiod pathway, the autonomous pathway, the vernalization pathway, the temperature pathway, the GA pathway, and the age pathway. However, studies have revealed that sugar is also involved in phase change and the control of branching behavior. In this review, we discuss recent advances in plant biology concerning the genetic and molecular mechanisms that allow plants to regulate phase transitions in response to the environment. We also propose connections between phase transition and plant architecture control.

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.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

What can lycophytes teach us about plant evolution and development? Modern perspectives on an ancient lineage

All Evo-Devo studies rely on representative sampling across the tree of interest to elucidate evolutionary trajectories through time. In land plants, genetic resources are well established in model species representing lineages including bryophytes (mosses, liverworts, and hornworts), monilophytes (ferns and allies), and seed plants (gymnosperms and flowering plants), but few resources are available for lycophytes (club mosses, spike mosses, and quillworts). Living lycophytes are a sister group to the euphyllophytes (the fern and seed plant clade), and have retained several ancestral morphological traits despite divergence from a common ancestor of vascular plants around 420 million years ago. This sister relationship offers a unique opportunity to study the conservation of traits such as sporophyte branching, vasculature, and indeterminacy, as well as the convergent evolution of traits such as leaves and roots which have evolved independently in each vascular plant lineage. To elucidate the evolution of vascular development and leaf formation, molecular studies using RNA Seq, quantitative reverse transcription polymerase chain reaction, in situ hybridisation and phylogenetics have revealed the diversification and expression patterns of KNOX, ARP, HD-ZIP, KANADI, and WOX gene families in lycophytes. However, the molecular basis of further trait evolution is not known. Here we describe morphological traits of living lycophytes and their extinct relatives, consider the molecular underpinnings of trait evolution and discuss future research required in lycophytes to understand the key evolutionary innovations enabling the growth and development of all vascular plants.

Reconstructing trait evolution in plant evo-devo studies

Our planet is teeming with an astounding diversity of plants. In a mere single group of closely related species, tremendous diversity can be observed in their form and function - the colour of petals in flowering plants, the shape of the fronds in ferns, and the branching pattern of the gametophyte in mosses. Diversity can also be found in subtler traits, such as the resistance to pathogens or the ability to recruit symbiotic microbes from the environment. Plant traits can also be highly conserved - at the cellular and metabolic levels, entire biosynthetic pathways are present in all plant groups, and morphological characteristics such as vascular tissues have been conserved for hundreds of millions of years. The research community that seeks to understand these traits - both the diverse and the conserved - by taking an evolutionary point-of-view on plant biology is growing. Here, we summarize a subset of the different aspects of plant evolutionary biology, provide a guide for structuring comparative biology approaches and discuss the pitfalls that (plant) researchers should avoid when embarking on such studies.

A Bacterial Tower of Babel: Quorum-Sensing Signaling Diversity and Its Evolution

Quorum sensing is a process in which bacteria secrete and sense a diffusible molecule, thereby enabling bacterial groups to coordinate their behavior in a density-dependent manner. Quorum sensing has evolved multiple times independently, utilizing different molecular pathways and signaling molecules. A common theme among many quorum-sensing families is their wide range of signaling diversity-different variants within a family code for different signal molecules with a cognate receptor specific to each variant. This pattern of vast allelic polymorphism raises several questions-How do different signaling variants interact with one another? How is this diversity maintained? And how did it come to exist in the first place? Here we argue that social interactions between signaling variants can explain the emergence and persistence of signaling diversity throughout evolution. Finally, we extend the discussion to include cases where multiple diverse systems work in concert in a single bacterium.

Evolution of Plant Metabolism: A (Bio)synthesis

Studying the evolution of metabolism is technically challenging. A new study combining in silico metabolic maps and phylogenomics allows reconstructing the diversification of plant metabolism across one billion years of evolution.

Sequencing and Analyzing the Transcriptomes of a Thousand Species Across the Tree of Life for Green Plants

The 1,000 Plants (1KP) initiative was the first large-scale effort to collect next-generation sequencing (NGS) data across a phylogenetically representative sampling of species for a major clade of life, in this case theViridiplantae, or green plants. As an international multidisciplinary consortium, we focused on plant evolution and its practical implications. Among the major outcomes were the inference of a reference species tree for green plants by phylotranscriptomic analysis of low-copy genes, a survey of paleopolyploidy (whole-genome duplications) across the Viridiplantae, the inferred evolutionary histories for many gene families and biological processes, the discovery of novel light-sensitive proteins for optogenetic studies in mammalian neuroscience, and elucidation of the genetic network for a complex trait (C₄ photosynthesis). Altogether, 1KP demonstrated how value can be extracted from a phylodiverse sequencing data set, providing a template for future projects that aim to generate even more data, including complete de novo genomes, across the tree of life.

Evolution of Plant Hormone Response Pathways

This review focuses on the evolution of plant hormone signaling pathways. Like the chemical nature of the hormones themselves, the signaling pathways are diverse. Therefore, we focus on a group of hormones whose primary perception mechanism involves an Skp1/Cullin/F-box-type ubiquitin ligase: auxin, jasmonic acid, gibberellic acid, and strigolactone. We begin with a comparison of the core signaling pathways of these four hormones, which have been established through studies conducted in model organisms in the Angiosperms. With the advent of next-generation sequencing and advanced tools for genetic manipulation, the door to understanding the origins of hormone signaling mechanisms in plants beyond these few model systems has opened. For example, in-depth phylogenetic analyses of hormone signaling components are now being complemented by genetic studies in early diverging land plants. Here we discuss recent investigations of how basal land plants make and sense hormones. Finally, we propose connections between the emergence of hormone signaling complexity and major developmental transitions in plant evolution.

Protein-coding changes preceded cis-regulatory gains in a newly evolved transcription circuit

Changes in both the coding sequence of transcriptional regulators and in the cis-regulatory sequences recognized by these regulators have been implicated in the evolution of transcriptional circuits. However, little is known about how they evolved in concert. We describe an evolutionary pathway in fungi where a new transcriptional circuit (a-specific gene repression by the homeodomain protein Matα2) evolved by coding changes in this ancient regulator, followed millions of years later by cis-regulatory sequence changes in the genes of its future regulon. By analyzing a group of species that has acquired the coding changes but not the cis-regulatory sites, we show that the coding changes became necessary for the regulator's deeply conserved function, thereby poising the regulator to jump-start formation of the new circuit.

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.

High-Throughput Protein Production Combined with High- Throughput SELEX Identifies an Extensive Atlas of Ciona robusta Transcription Factor DNA-Binding Specificities

Transcription factors (TFs) control gene transcription, binding to specific DNA motifs located in cis-regulatory elements across the genome. The identification of TF-binding motifs is thus an important aspect to understand the role of TFs in gene regulation. SELEX, Systematic Evolution of Ligands by EXponential enrichment, is an efficient in vitro method, which can be used to determine the DNA-binding specificity of TFs. Thanks to the development of high-throughput (HT) DNA cloning system and protein production technology, the classical SELEX assay has be extended to high-throughput scale (HT-SELEX).We report here the detailed protocol for the cloning, production, and purification of 420 Ciona robusta DNA BD. 263 Ciona robusta TF DNA-binding domain proteins were purified in milligram quantities and analyzed by HT-SELEX. The identification of 139 recognition sequences generates an atlas of protein-DNA-binding specificities that is crucial for the understanding of the gene regulatory network (GRN) of Ciona robusta. Overall, our analysis suggests that the Ciona robusta repertoire of sequence-specific transcription factors comprises less than 500 genes. The protocols for high-throughput protein production and HT-SELEX described in this article for the study of Ciona robusta TF DNA-binding specificity are generic and have been successfully applied to a wide range of TFs from other species, including human, mouse, and Drosophila.

Anthoceros genomes illuminate the origin of land plants and the unique biology of hornworts

Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbon-concentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Unwinding BRAHMA Functions in Plants

The ATP-dependent Switch/Sucrose non-fermenting (SWI/SNF) chromatin remodeling complex (CRC) regulates the transcription of many genes by destabilizing interactions between DNA and histones. In plants, BRAHMA (BRM), one of the two catalytic ATPase subunits of the complex, is the closest homolog of the yeast and animal SWI2/SNF2 ATPases. We summarize here the advances describing the roles of BRM in plant development as well as its recently reported chromatin-independent role in pri-miRNA processing in vitro and in vivo. We also enlighten the roles of plant-specific partners that physically interact with BRM. Three main types of partners can be distinguished: (i) DNA-binding proteins such as transcription factors which mostly cooperate with BRM in developmental processes, (ii) enzymes such as kinases or proteasome-related proteins that use BRM as substrate and are often involved in response to abiotic stress, and (iii) an RNA-binding protein which is involved with BRM in chromatin-independent pri-miRNA processing. This overview contributes to the understanding of the central position occupied by BRM within regulatory networks controlling fundamental biological processes in plants.

Extreme divergence between one-to-one orthologs: the structure of N15 Cro bound to operator DNA and its relationship to the λ Cro complex

The gene cro promotes lytic growth of phages through binding of Cro protein dimers to regulatory DNA sites. Most Cro proteins are one-to-one orthologs, yet their sequence, structure and binding site sequences are quite divergent across lambdoid phages. We report the cocrystal structure of bacteriophage N15 Cro with a symmetric consensus site. We contrast this complex with an orthologous structure from phage λ, which has a dissimilar binding site sequence and a Cro protein that is highly divergent in sequence, dimerization interface and protein fold. The N15 Cro complex has less DNA bending and smaller DNA-induced changes in protein structure. N15 Cro makes fewer direct contacts and hydrogen bonds to bases, relying mostly on water-mediated and Van der Waals contacts to recognize the sequence. The recognition helices of N15 Cro and λ Cro make mostly nonhomologous and nonanalogous contacts. Interface alignment scores show that half-site binding geometries of N15 Cro and λ Cro are less similar to each other than to distantly related CI repressors. Despite this divergence, the Cro family shows several code-like protein–DNA sequence covariations. In some cases, orthologous genes can achieve a similar biological function using very different specific molecular interactions.

Asymmetrical diversification of the receptor-ligand interaction controlling self-incompatibility in Arabidopsis

How two-component genetic systems accumulate evolutionary novelty and diversify in the course of evolution is a fundamental problem in evolutionary systems biology. In the Brassicaceae, self-incompatibility (SI) is a spectacular example of a diversified allelic series in which numerous highly diverged receptor-ligand combinations are segregating in natural populations. However, the evolutionary mechanisms by which new SI specificities arise have remained elusive. Using in planta ancestral protein reconstruction, we demonstrate that two allelic variants segregating as distinct receptor-ligand combinations diverged through an asymmetrical process whereby one variant has retained the same recognition specificity as their (now extinct) putative ancestor, while the other has functionally diverged and now represents a novel specificity no longer recognized by the ancestor. Examination of the structural determinants of the shift in binding specificity suggests that qualitative rather than quantitative changes of the interaction are an important source of evolutionary novelty in this highly diversified receptor-ligand system.

Phylogenetic Conflicts, Combinability, and Deep Phylogenomics in Plants

Studies have demonstrated that pervasive gene tree conflict underlies several important phylogenetic relationships where different species tree methods produce conflicting results. Here, we present a means of dissecting the phylogenetic signal for alternative resolutions within a data set in order to resolve recalcitrant relationships and, importantly, identify what the data set is unable to resolve. These procedures extend upon methods for isolating conflict and concordance involving specific candidate relationships and can be used to identify systematic error and disambiguate sources of conflict among species tree inference methods. We demonstrate these on a large phylogenomic plant data set. Our results support the placement of Amborella as sister to the remaining extant angiosperms, Gnetales as sister to pines, and the monophyly of extant gymnosperms. Several other contentious relationships, including the resolution of relationships within the bryophytes and the eudicots, remain uncertain given the low number of supporting gene trees. To address whether concatenation of filtered genes amplified phylogenetic signal for relationships, we implemented a combinatorial heuristic to test combinability of genes. We found that nested conflicts limited the ability of data filtering methods to fully ameliorate conflicting signal amongst gene trees. These analyses confirmed that the underlying conflicting signal does not support broad concatenation of genes. Our approach provides a means of dissecting a specific data set to address deep phylogenetic relationships while also identifying the inferential boundaries of the data set. [Angiosperms; coalescent; gene-tree conflict; genomics; phylogenetics; phylogenomics.]

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

Ancient duplications and grass-specific transposition influenced the evolution of LEAFY transcription factor genes

The LFY transcription factor gene family are important in the promotion of cell proliferation and floral development. Understanding their evolution offers an insight into floral development in plant evolution. Though a promiscuous transition intermediate and a gene duplication event within the LFY family had been identified previously, the early evolutionary path of this family remained elusive. Here, we reconstructed the LFY family phylogeny using maximum-likelihood and Bayesian inference methods incorporating LFY genes from all major lineages of streptophytes. The well-resolved phylogeny unveiled a high-confidence duplication event before the functional divergence of types I and II LFY genes in the ancestry of liverworts, mosses and tracheophytes, supporting sub-functionalization of an ancestral promiscuous gene. The identification of promiscuous genes in Osmunda suggested promiscuous LFY genes experienced an ancient transient duplication. Genomic synteny comparisons demonstrated a deep genomic positional conservation of LFY genes and an ancestral lineage-specific transposition activity in grasses.

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.

Similarity regression predicts evolution of transcription factor sequence specificity

Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.