Resurrecting the role of transcription factor change in developmental evolution

Terminal regions of a protein are a hotspot for low complexity regions and selection

Volatile low complexity regions (LCRs) are a novel source of adaptive variation, functional diversification and evolutionary novelty. An interplay of selection and mutation governs the composition and length of low complexity regions. High %GC and mutations provide length variability because of mechanisms like replication slippage. Owing to the complex dynamics between selection and mutation, we need a better understanding of their coexistence. Our findings underscore that positively selected sites (PSS) and low complexity regions prefer the terminal regions of genes, co-occurring in most Tetrapoda clades. We observed that positively selected sites within a gene have position-specific roles. Central-positively selected site genes primarily participate in defence responses, whereas terminal-positively selected site genes exhibit non-specific functions. Low complexity region-containing genes in the Tetrapoda clade exhibit a significantly higher %GC and lower ω (dN/dS: non-synonymous substitution rate/synonymous substitution rate) compared with genes without low complexity regions. This lower ω implies that despite providing rapid functional diversity, low complexity region-containing genes are subjected to intense purifying selection. Furthermore, we observe that low complexity regions consistently display ubiquitous prevalence at lower purity levels, but exhibit a preference for specific positions within a gene as the purity of the low complexity region stretch increases, implying a composition-dependent evolutionary role. Our findings collectively contribute to the understanding of how genetic diversity and adaptation are shaped by the interplay of selection and low complexity regions in the Tetrapoda clade.

Divergence of Grainy head affects chromatin accessibility, gene expression, and embryonic viability in Drosophila melanogaster

Pioneer factors are critical for gene regulation and development because they bind chromatin and make DNA more accessible for binding by other transcription factors. The pioneer factor Grainy head (Grh) is present across metazoans and has been shown to retain a role in epithelium development in fruit flies, nematodes, and mice despite extensive divergence in both amino acid sequence and length. Here, we investigate the evolution of Grh function by comparing the effects of the fly ( Drosophila melanogaster ) and worm ( Caenorhabditis elegans ) Grh orthologs on chromatin accessibility, gene expression, embryonic development, and viability in transgenic D. melanogaster . We found that the Caenorhabditis elegans ortholog rescued cuticle development but not full embryonic viability in Drosophila melanogaster grh null mutants. At the molecular level, the C. elegans ortholog only partially rescued chromatin accessibility and gene expression. Divergence in the disordered N-terminus of the Grh protein contributes to these differences in embryonic viability and molecular phenotypes. These data show how pioneer factors can diverge in sequence and function at the molecular level while retaining conserved developmental functions at the organismal level.

SUMMARY STATEMENT

Despite divergence in a disordered region that affects function at both molecular and organismal levels, the Caenorhabditis elegans Grainy head (Grh) protein rescued cuticle morphology in D. melanogaster embryos.

Rapid evolution of genes with anti-cancer functions during the origins of large bodies and cancer resistance in elephants

Elephants have emerged as a model system to study the evolution of body size and cancer resistance because, despite their immense size, they have a very low prevalence of cancer. Previous studies have found that duplication of tumor suppressors at least partly contributes to the evolution of anti-cancer cellular phenotypes in elephants. Still, many other mechanisms must have contributed to their augmented cancer resistance. Here, we use a suite of codon-based maximum-likelihood methods and a dataset of 13,310 protein-coding gene alignments from 261 Eutherian mammals to identify positively selected and rapidly evolving elephant genes. We found 496 genes (3.73% of alignments tested) with statistically significant evidence for positive selection and 660 genes (4.96% of alignments tested) that likely evolved rapidly in elephants. Positively selected and rapidly evolving genes are statistically enriched in gene ontology terms and biological pathways related to regulated cell death mechanisms, DNA damage repair, cell cycle regulation, epidermal growth factor receptor (EGFR) signaling, and immune functions, particularly neutrophil granules and degranulation. All of these biological factors are plausibly related to the evolution of cancer resistance. Thus, these positively selected and rapidly evolving genes are promising candidates for genes contributing to elephant-specific traits, including the evolution of molecular and cellular characteristics that enhance cancer resistance.

Cooption of polyalanine tract into a repressor domain in the mammalian transcription factor HoxA11

An enduring problem in biology is explaining how novel functions of genes originated and how those functions diverge between species. Despite detailed studies on the functional evolution of a few proteins, the molecular mechanisms by which protein functions have evolved are almost entirely unknown. Here, we show that a polyalanine tract in the homeodomain transcription factor HoxA11 arose in the stem-lineage of mammals and functions as an autonomous repressor module by physically interacting with the PAH domains of SIN3 proteins. These results suggest that long polyalanine tracts, which are common in transcription factors and often associated with disease, may tend to function as repressor domains and can contribute to the diversification of transcription factor functions despite the deleterious consequences of polyalanine tract expansion.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Updated Phylogeny and Protein Structure Predictions Revise the Hypothesis on the Origin of MADS-box Transcription Factors in Land Plants

MADS-box transcription factors (TFs), among the first TFs extensively studied, exhibit a wide distribution across eukaryotes and play diverse functional roles. Varying by domain architecture, MADS-box TFs in land plants are categorized into Type I (M-type) and Type II (MIKC-type). Type I and II genes have been considered orthologous to the SRF and MEF2 genes in animals, respectively, presumably originating from a duplication before the divergence of eukaryotes. Here, we exploited the increasing availability of eukaryotic MADS-box sequences and reassessed their evolution. While supporting the ancient duplication giving rise to SRF- and MEF2-types, we found that Type I and II genes originated from the MEF2-type genes through another duplication in the most recent common ancestor (MRCA) of land plants. Protein structures predicted by AlphaFold2 and OmegaFold support our phylogenetic analyses, with plant Type I and II TFs resembling the MEF2-type structure, rather than SRFs. We hypothesize that the ancestral SRF-type TFs were lost in the MRCA of Archaeplastida (the kingdom Plantae sensu lato). The retained MEF2-type TFs acquired a Keratin-like domain and became MIKC-type before the divergence of Streptophyta. Subsequently in the MRCA of land plants, M-type TFs evolved from a duplicated MIKC-type precursor through loss of the Keratin-like domain, leading to the Type I clade. Both Type I and II TFs expanded and functionally differentiated in concert with the increasing complexity of land plant body architecture. The recruitment of these originally stress-responsive TFs into developmental programs, including those underlying reproduction, may have facilitated the adaptation to the terrestrial environment.

Accelerated Evolution Analysis Uncovers PKNOX2 as a Key Transcription Factor in the Mammalian Cochlea

The genetic bases underlying the evolution of morphological and functional innovations of the mammalian inner ear are poorly understood. Gene regulatory regions are thought to play an important role in the evolution of form and function. To uncover crucial hearing genes whose regulatory machinery evolved specifically in mammalian lineages, we mapped accelerated noncoding elements (ANCEs) in inner ear transcription factor (TF) genes and found that PKNOX2 harbors the largest number of ANCEs within its transcriptional unit. Using reporter gene expression assays in transgenic zebrafish, we determined that four PKNOX2-ANCEs drive differential expression patterns when compared with ortholog sequences from close outgroup species. Because the functional role of PKNOX2 in cochlear hair cells has not been previously investigated, we decided to study Pknox2 null mice generated by CRISPR/Cas9 technology. We found that Pknox2-/- mice exhibit reduced distortion product otoacoustic emissions (DPOAEs) and auditory brainstem response (ABR) thresholds at high frequencies together with an increase in peak 1 amplitude, consistent with a higher number of inner hair cells (IHCs)-auditory nerve synapsis observed at the cochlear basal region. A comparative cochlear transcriptomic analysis of Pknox2-/- and Pknox2+/+ mice revealed that key auditory genes are under Pknox2 control. Hence, we report that PKNOX2 plays a critical role in cochlear sensitivity at higher frequencies and that its transcriptional regulation underwent lineage-specific evolution in mammals. Our results provide novel insights about the contribution of PKNOX2 to normal auditory function and to the evolution of high-frequency hearing in mammals.

The START domain potentiates HD-ZIPIII transcriptional activity

The CLASS III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIPIII) transcription factors (TFs) were repeatedly deployed over 725 million years of evolution to regulate central developmental innovations. The START domain of this pivotal class of developmental regulators was recognized over 20 years ago, but its putative ligands and functional contributions remain unknown. Here, we demonstrate that the START domain promotes HD-ZIPIII TF homodimerization and increases transcriptional potency. Effects on transcriptional output can be ported onto heterologous TFs, consistent with principles of evolution via domain capture. We also show the START domain binds several species of phospholipids, and that mutations in conserved residues perturbing ligand binding and/or its downstream conformational readout abolish HD-ZIPIII DNA-binding competence. Our data present a model in which the START domain potentiates transcriptional activity and uses ligand-induced conformational change to render HD-ZIPIII dimers competent to bind DNA. These findings resolve a long-standing mystery in plant development and highlight the flexible and diverse regulatory potential coded within this widely distributed evolutionary module.

Evolution of a new form of haploid-specific gene regulation appearing in a limited clade of ascomycete yeast species

Over evolutionary timescales, the logic and pattern of cell-type specific gene expression can remain constant, yet the molecular mechanisms underlying such regulation can drift between alternative forms. Here, we document a new example of this principle in the regulation of the haploid-specific genes in a small clade of fungal species. For most ascomycete fungal species, transcription of these genes is repressed in the a/α cell type by a heterodimer of two homeodomain proteins, Mata1 and Matα2. We show that in the species Lachancea kluyveri, most of the haploid-specific genes are regulated in this way, but repression of one haploid-specific gene (GPA1) requires, in addition to Mata1 and Matα2, a third regulatory protein, Mcm1. Model building, based on x-ray crystal structures of the three proteins, rationalizes the requirement for all three proteins: no single pair of the proteins is optimally arranged, and we show that no single pair can bring about repression. This case study exemplifies the idea that the energy of DNA binding can be "shared out" in different ways and can result in different DNA-binding solutions across different genes-while maintaining the same overall pattern of gene expression.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

Spatiotemporal Regulation of a Single Adaptively Evolving Trans-Regulatory Element Contributes to Spermatogenetic Expression Divergence in Drosophila

Due to extensive pleiotropy, trans-acting elements are often thought to be evolutionarily constrained. While the impact of trans-acting elements on gene expression evolution has been extensively studied, relatively little is understood about the contribution of a single trans regulator to interspecific expression and phenotypic divergence. Here, we disentangle the effects of genomic context and miR-983, an adaptively evolving young microRNA, on expression divergence between Drosophila melanogaster and D. simulans. We show miR-983 effects promote interspecific expression divergence in testis despite its antagonism with the often-predominant context effects. Single-cyst RNA-seq reveals that distinct sets of genes gain and lose miR-983 influence under disruptive or diversifying selection at different stages of spermatogenesis, potentially helping minimize antagonistic pleiotropy. At the round spermatid stage, the effects of miR-983 are weak and distributed, coincident with the transcriptome undergoing drastic expression changes. Knocking out miR-983 causes reduced sperm length with increased within-individual variation in D. melanogaster but not in D. simulans, and the D. melanogaster knockout also exhibits compromised sperm defense ability. Our results provide empirical evidence for the resolution of antagonistic pleiotropy and also have broad implications for the function and evolution of new trans regulators.

Correlation Between Improved Mating Efficiency and Weakened Scaffold-Kinase Interaction in the Mating Pheromone Response Pathway Revealed by Interspecies Complementation

Scaffold protein Ste5 and associated kinases, including Ste11, Ste7, and Fus3, are core components of the mating pheromone pathway, which is required to induce a mating response. Orthologs of these proteins are widely present in fungi, but to which extent one protein can be replaced by its ortholog is less well understood. Here, interspecies complementation was carried out to evaluate the functional homology of Ste5 and associated kinases in Kluyveromyces lactis, K. marxianus, and Saccharomyces cerevisiae. These three species occupy important positions in the evolution of hemiascomycetes. Results indicated that Ste5 and associated kinases in K. lactis and K. marxianus could be functionally replaced by their orthologs to different extents. However, the extent of sequence identity, either between full-length proteins or between domains, did not necessarily indicate the extent of functional replaceability. For example, Ste5, the most unconserved protein in sequence, achieved the highest average functional replaceability. Notably, swapping Ste5 between K. lactis and K. marxianus significantly promoted mating in both species and the weakened interaction between the Ste5 and Ste7 might contribute to this phenotype. Consistently, chimeric Ste5 displaying a higher affinity for Ste7 decreased the mating efficiency, while chimeric Ste5 displaying a lower affinity for Ste7 improved the mating efficiency. Furthermore, the length of a negatively charged segment in the Ste7-binding domain of Ste5 was negatively correlated with the mating efficiency in K. lactis and K. marxianus. Extending the length of the segment in KlSte5 improved its interaction with Ste7 and that might contribute to the reduced mating efficiency. Our study suggested a novel role of Ste5-Ste7 interaction in the negative regulation of the pheromone pathway. Meanwhile, Ste5 mutants displaying improved mating efficiency facilitated the breeding and selection of Kluyveromyces strains for industrial applications.

Lessons from a transcription factor: Alx1 provides insights into gene regulatory networks, cellular reprogramming, and cell type evolution

The skeleton-forming cells of sea urchins and other echinoderms have been studied by developmental biologists as models of cell specification and morphogenesis for many decades. The gene regulatory network (GRN) deployed in the embryonic skeletogenic cells of euechinoid sea urchins is one of the best understood in any developing animal. Recent comparative studies have leveraged the information contained in this GRN, bringing renewed attention to the diverse patterns of skeletogenesis within the phylum and the evolutionary basis for this diversity. The homeodomain-containing transcription factor, Alx1, was originally shown to be a core component of the skeletogenic GRN of the sea urchin embryo. Alx1 has since been found to be key regulator of skeletal cell identity throughout the phylum. As such, Alx1 is currently serving as a lens through which multiple developmental processes are being investigated. These include not only GRN organization and evolution, but also cell reprogramming, cell type evolution, and the gene regulatory control of morphogenesis. This review summarizes our current state of knowledge concerning Alx1 and highlights the insights it is yielding into these important developmental and evolutionary processes.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

Evolution and expansion of the RUNX2 QA repeat corresponds with the emergence of vertebrate complexity

Runt-related transcription factor 2 (RUNX2) is critical for the development of the vertebrate bony skeleton. Unlike other RUNX family members, RUNX2 possesses a variable poly-glutamine, poly-alanine (QA) repeat domain. Natural variation within this repeat is able to alter the transactivation potential of RUNX2, acting as an evolutionary 'tuning knob' suggested to influence mammalian skull shape. However, the broader role of the RUNX2 QA repeat throughout vertebrate evolution is unknown. In this perspective, we examine the role of the RUNX2 QA repeat during skeletal development and discuss how its emergence and expansion may have facilitated the evolution of morphological novelty in vertebrates.

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.

A tale of two leeches: Toward the understanding of the evolution and development of behavioral neural circuits

In the animal kingdom, behavioral traits encompass a broad spectrum of biological phenotypes that have critical roles in adaptive evolution, but an EvoDevo approach has not been broadly used to study behavior evolution. Here, we propose that, by integrating two leech model systems, each of which has already attained some success in its respective field, it is possible to take on behavioral traits with an EvoDevo approach. We first identify the developmental changes that may theoretically lead to behavioral evolution and explain why an EvoDevo study of behavior is challenging. Next, we discuss the pros and cons of the two leech model species, Hirudo, a classic model for invertebrate neurobiology, and Helobdella, an emerging model for clitellate developmental biology, as models for behavioral EvoDevo research. Given the limitations of each leech system, neither is particularly strong for behavioral EvoDevo. However, the two leech systems are complementary in their technical accessibilities, and they do exhibit some behavioral similarities and differences. By studying them in parallel and together with additional leech species such as Haementeria, it is possible to explore the different levels of behavioral development and 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.

CHD9 upregulates RUNX2 and has a potential role in skeletal evolution

BACKGROUND

Changes in gene regulation are widely recognized as an important driver of adaptive phenotypic evolution. However, the specific molecular mechanisms that underpin such changes are still poorly understood. Chromatin state plays an essential role in gene regulation, by influencing the accessibility of coding loci to the transcriptional machinery. Changes in the function of chromatin remodellers are therefore strong candidates to drive changes in gene expression associated with phenotypic adaptation. Here, we identify amino acid homoplasies in the chromatin remodeller CHD9, shared between the extinct marsupial thylacine and eutherian wolf which show remarkable skull convergence. CHD9 is involved in osteogenesis, though its role in the process is still poorly understood. We examine whether CHD9 is able to regulate the expression of osteogenic target genes and examine the function of a key substitution in the CHD9 DNA binding domain.

RESULTS

We examined whether CHD9 was able to upregulate its osteogenic target genes, RUNX2, Osteocalcin (OC) and ALP in HEK293T cells. We found that overexpression of CHD9 upregulated RUNX2, the master regulator of osteoblast cell fate, but not the downstream genes OC or ALP, supporting the idea that CHD9 regulates osteogenic progenitors rather than terminal osteoblasts. We also found that the evolutionary substitution in the CHD9 DNA binding domain does not alter protein secondary structure, but was able to drive a small but insignificant increase in RUNX2 activation. Finally, CHD9 was unable to activate an episomal RUNX2 promoter-reporter construct, suggesting that CHD9 requires the full chromatin complement for its function.

CONCLUSIONS

We provide new evidence to the role of CHD9 in osteogenic differentiation through its newly observed ability to upregulate the expression of RUNX2. Though we were unable to identify significant functional consequences of the evolutionary substitution in HEK293T cells, our study provides important steps forward in the functional investigation of protein homoplasy and its role in developmental processes. Mutations in coding genes may be a mechanism for driving adaptive changes in gene expression, and their validation is essential towards determining the functional consequences of evolutionary homoplasy.

Compound Dynamics and Combinatorial Patterns of Amino Acid Repeats Encode a System of Evolutionary and Developmental Markers

Homopolymeric amino acid repeats (AARs) like polyalanine (polyA) and polyglutamine (polyQ) in some developmental proteins (DPs) regulate certain aspects of organismal morphology and behavior, suggesting an evolutionary role for AARs as developmental "tuning knobs." It is still unclear, however, whether these are occasional protein-specific phenomena or hints at the existence of a whole AAR-based regulatory system in DPs. Using novel approaches to trace their functional and evolutionary history, we find quantitative evidence supporting a generalized, combinatorial role of AARs in developmental processes with evolutionary implications. We observe nonrandom AAR distributions and combinations in HOX and other DPs, as well as in their interactomes, defining elements of a proteome-wide combinatorial functional code whereby different AARs and their combinations appear preferentially in proteins involved in the development of specific organs/systems. Such functional associations can be either static or display detectable evolutionary dynamics. These findings suggest that progressive changes in AAR occurrence/combination, by altering embryonic development, may have contributed to taxonomic divergence, leaving detectable traces in the evolutionary history of proteomes. Consistent with this hypothesis, we find that the evolutionary trajectories of the 20 AARs in eukaryotic proteomes are highly interrelated and their individual or compound dynamics can sharply mark taxonomic boundaries, or display clock-like trends, carrying overall a strong phylogenetic signal. These findings provide quantitative evidence and an interpretive framework outlining a combinatorial system of AARs whose compound dynamics mark at the same time DP functions and evolutionary transitions.

Looking back to look forward: protein-protein interactions and the evolution of development

The evolutionary modification of development was fundamental in generating extant plant diversity. Similarly, the modification of development is a path forward to engineering the plants of the future, provided we know enough about what to modify. Understanding how extant diversity was generated will reveal productive pathways forward for modifying development. Here, I discuss four examples of developmental pathways that have been remodeled by changes to protein-protein interactions. These are cases where changes to developmental pathways have been paralleled by recent changes, selected for or engineered by humans. Extant plant diversity represents a vast treasure trove of molecular solutions to ecological problems. Mining this treasure trove will allow for the intentional modification of plant development for solving future problems.

Physiological and genomic evidence that selection on the transcription factor Epas1 has altered cardiovascular function in high-altitude deer mice

Evolutionary adaptation to extreme environments often requires coordinated changes in multiple intersecting physiological pathways, but how such multi-trait adaptation occurs remains unresolved. Transcription factors, which regulate the expression of many genes and can simultaneously alter multiple phenotypes, may be common targets of selection if the benefits of induced changes outweigh the costs of negative pleiotropic effects. We combined complimentary population genetic analyses and physiological experiments in North American deer mice (Peromyscus maniculatus) to examine links between genetic variation in transcription factors that coordinate physiological responses to hypoxia (hypoxia-inducible factors, HIFs) and multiple physiological traits that potentially contribute to high-altitude adaptation. First, we sequenced the exomes of 100 mice sampled from different elevations and discovered that several SNPs in the gene Epas1, which encodes the oxygen sensitive subunit of HIF-2α, exhibited extreme allele frequency differences between highland and lowland populations. Broader geographic sampling confirmed that Epas1 genotype varied predictably with altitude throughout the western US. We then discovered that Epas1 genotype influences heart rate in hypoxia, and the transcriptomic responses to hypoxia (including HIF targets and genes involved in catecholamine signaling) in the heart and adrenal gland. Finally, we used a demographically-informed selection scan to show that Epas1 variants have experienced a history of spatially varying selection, suggesting that differences in cardiovascular function and gene regulation contribute to high-altitude adaptation. Our results suggest a mechanism by which Epas1 may aid long-term survival of high-altitude deer mice and provide general insights into the role that highly pleiotropic transcription factors may play in the process of environmental adaptation.

Variation across Species and Levels: Implications for Model Species Research

The selection of model species tends to involve two typically unstated assumptions, namely: (1) that the similarity between species decreases steadily with phylogenetic distance, and (2) that similarities are greater at lower levels of biological organization. The first assumption holds on average, but species similarities tend to decrease with the square root of divergence time, rather than linearly, and lineages with short generation times (which includes most model species) tend to diverge faster than average, making the decrease in similarity non-monotonic. The second assumption is more difficult to test. Comparative molecular research has traditionally emphasized species similarities over differences, whereas comparative research at higher levels of organization frequently highlights the species differences. However, advances in comparative genomics have brought to light a great variety of species differences, not just in gene regulation but also in protein coding genes. Particularly relevant are cases in which homologous high-level characters are based on non-homologous genes. This phenomenon of non-orthologous gene displacement, or "deep non-homology," indicates that species differences at the molecular level can be surprisingly large. Given these observations, it is not surprising that some findings obtained in model species do not generalize across species as well as researchers had hoped, even if the research is molecular.

Successive Domain Rearrangements Underlie the Evolution of a Regulatory Module Controlled by a Small Interfering Peptide

The establishment of new interactions between transcriptional regulators increases the regulatory diversity that drives phenotypic novelty. To understand how such interactions evolve, we have studied a regulatory module (DDR) composed by three MYB-like proteins: DIVARICATA (DIV), RADIALIS (RAD), and DIV-and-RAD-Interacting Factor (DRIF). The DIV and DRIF proteins form a transcriptional complex that is disrupted in the presence of RAD, a small interfering peptide, due to the formation of RAD–DRIF dimers. This dynamic interaction result in a molecular switch mechanism responsible for the control of distinct developmental processes in plants. Here, we have determined how the DDR regulatory module was established by analyzing the origin and evolution of the DIV, DRIF, and RAD protein families and the evolutionary history of their interactions. We show that duplications of a pre-existing MYB domain originated the DIV and DRIF protein families in the ancestral lineage of green algae, and, later, the RAD family in seed plants. Intraspecies interactions between the MYB domains of DIV and DRIF proteins are detected in green algae, whereas the earliest evidence of an interaction between DRIF and RAD proteins occurs in the gymnosperms, coincident with the establishment of the RAD family. Therefore, the DDR module evolved in a stepwise progression with the DIV–DRIF transcription complex evolving prior to the antagonistic RAD–DRIF interaction that established the molecular switch mechanism. Our results suggest that the successive rearrangement and divergence of a single protein domain can be an effective evolutionary mechanism driving new protein interactions and the establishment of novel regulatory modules.

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.

Genetic architecture and sex-specific selection govern modular, male-biased evolution of doublesex

doublesex regulates early embryonic sex differentiation in holometabolous insects, along with the development of species-, sex-, and morph-specific adaptations during pupal stages. How does a highly conserved gene with a critical developmental role also remain functionally dynamic enough to gain ecologically important adaptations that are divergent in sister species? We analyzed patterns of exon-level molecular evolution and protein structural homology of doublesex from 145 species of four insect orders representing 350 million years of divergence. This analysis revealed that evolution of doublesex was governed by a modular architecture: Functional domains and female-specific regions were highly conserved, whereas male-specific sequences and protein structures evolved up to thousand-fold faster, with sites under pervasive and/or episodic positive selection. This pattern of sex bias was reversed in Hymenoptera. Thus, highly conserved yet dynamic master regulators such as doublesex may partition specific conserved and novel functions in different genic modules at deep evolutionary time scales.

The Ciona myogenic regulatory factor functions as a typical MRF but possesses a novel N-terminus that is essential for activity

Electroporation-based assays were used to test whether the myogenic regulatory factor (MRF) of Ciona intestinalis (CiMRF) interferes with endogenous developmental programs, and to evaluate the importance of its unusual N-terminus for muscle development. We found that CiMRF suppresses both notochord and endoderm development when it is expressed in these tissues by a mechanism that may involve activation of muscle-specific microRNAs. Because these results add to a large body of evidence demonstrating the exceptionally high degree of functional conservation among MRFs, we were surprised to discover that non-ascidian MRFs were not myogenic in Ciona unless they formed part of a chimeric protein containing the CiMRF N-terminus. Equally surprising, we found that despite their widely differing primary sequences, the N-termini of MRFs of other ascidian species could form chimeric MRFs that were also myogenic in Ciona. This domain did not rescue the activity of a Brachyury protein whose transcriptional activation domain had been deleted, and so does not appear to constitute such a domain. Our results indicate that ascidians have previously unrecognized and potentially novel requirements for MRF-directed myogenesis. Moreover, they provide the first example of a domain that is essential to the core function of an important family of gene regulatory proteins, one that, to date, has been found in only a single branch of the family.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

UNVEILing connections between genotype, phenotype, and fitness in natural populations

Understanding the links between genetic variation and fitness in natural populations is a central goal of evolutionary genetics. This monumental task spans the fields of classical and molecular genetics, population genetics, biochemistry, physiology, developmental biology, and ecology. Advances to our molecular and developmental toolkits are facilitating integrative approaches across these traditionally separate fields, providing a more complete picture of the genotype-phenotype map in natural and non-model systems. Here, we summarize research presented at the first annual symposium of the UNVEIL Network, an NSF-funded collaboration between the University of Montana and the University of Nebraska, Lincoln, which took place from the 1st to the 3rd of June, 2018. We discuss how this body of work advances basic evolutionary science, what it implies for our ability to predict evolutionary change, and how it might inform novel conservation strategies.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

Borders of Cis-Regulatory DNA Sequences Preferentially Harbor the Divergent Transcription Factor Binding Motifs in the Human Genome

Changes in cis-regulatory DNA sequences and transcription factor (TF) repertoires provide major sources of phenotypic diversity that shape the evolution of gene regulation in eukaryotes. The DNA-binding specificities of TFs may be diversified or produce new variants in different eukaryotic species. However, it is currently unclear how various levels of divergence in TF DNA-binding specificities or motifs became introduced into the cis-regulatory DNA regions of the genome over evolutionary time. Here, we first estimated the evolutionary divergence levels of TF binding motifs and quantified their occurrence at DNase I-hypersensitive sites. Results from our in silico motif scan and experimentally derived chromatin immunoprecipitation (TF-ChIP) show that the divergent motifs tend to be introduced in the edges of cis-regulatory regions, which is probably accompanied by the expansion of the accessible core of promoter-associated regulatory elements during evolution. We also find that the genes neighboring the expanded cis-regulatory regions with the most divergent motifs are associated with functions like development and morphogenesis. Accordingly, we propose that the accumulation of divergent motifs in the edges of cis-regulatory regions provides a functional mechanism for the evolution of divergent regulatory circuits.

Heterologous Hsp90 promotes phenotypic diversity through network evolution

Biological processes in living cells are often carried out by gene networks in which signals and reactions are integrated through network hubs. Despite their functional importance, it remains unclear to what extent network hubs are evolvable and how alterations impact long-term evolution. We investigated these issues using heat shock protein 90 (Hsp90), a central hub of proteostasis networks. When native Hsp90 in Saccharomyces cerevisiae cells was replaced by the ortholog from hypersaline-tolerant Yarrowia lipolytica that diverged from S. cerevisiae about 270 million years ago, the cells exhibited improved growth in hypersaline environments but compromised growth in others, indicating functional divergence in Hsp90 between the two yeasts. Laboratory evolution shows that evolved Y. lipolytica-HSP90–carrying S. cerevisiae cells exhibit a wider range of phenotypic variation than cells carrying native Hsp90. Identified beneficial mutations are involved in multiple pathways and are often pleiotropic. Our results show that cells adapt to a heterologous Hsp90 by modifying different subnetworks, facilitating the evolution of phenotypic diversity inaccessible to wild-type cells.

Biological processes in living cells are often carried out by gene networks. Hubs are highly connected network components important for integrating signal inputs and generating responsive functional outputs. Heat shock protein 90 (Hsp90), a versatile hub in the protein homeostasis network, is a molecular chaperone essential for cell viability in all tested eukaryotic cells. In yeast, about a quarter of the expressed proteins are profoundly influenced when Hsp90 activity is reduced. Despite its pivotal role, we found that the function of Hsp90 has diverged between two yeast species, Yarrowia lipolytica and Saccharomyces cerevisiae, which split about 270 million years ago. To understand the impacts and adaptive strategies in cells with an altered network hub, we conducted laboratory evolution experiments using a S. cerevisiae strain in which native Hsp90 is replaced by its counterpart in Y. lipolytica. We observed different fitness gain or loss under various stress conditions in individual evolved clones, suggesting that cells adapted via different evolutionary paths. Genome sequencing and mutation reconstitution experiments show that beneficial mutations occurred in multiple Hsp90-related pathways that interact with each other. Our results show that a perturbed network allows cells to evolve a broader range of phenotypic diversity unavailable to wild-type cells.

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Rethinking Living Fossils

Biologists would be mistaken if they relegated living fossils to paleontological inquiry or assumed that the concept is dead. It is now used to describe entities ranging from viruses to higher taxa, despite recent warnings of misleading inferences. Current work on character evolution illustrates how analyzing living fossils and stasis in terms of parts (characters) and wholes (e.g., organisms and lineages) advances our understanding of prolonged stasis at many hierarchical levels. Instead of viewing the concept's task as categorizing living fossils, we show how its primary role is to mark out what is in need of explanation, accounting for the persistence of both molecular and morphological traits. Rethinking different conceptions of living fossils as specific hypotheses reveals novel avenues for research that integrate phylogenetics, ecological and evolutionary modeling, and evo-devo to produce a more unified theoretical outlook.

Tissue-Specific Transcriptome for Poeciliopsis prolifica Reveals Evidence for Genetic Adaptation Related to the Evolution of a Placental Fish

The evolution of the placenta is an excellent model to examine the evolutionary processes underlying adaptive complexity due to the recent, independent derivation of placentation in divergent animal lineages. In fishes, the family Poeciliidae offers the opportunity to study placental evolution with respect to variation in degree of post-fertilization maternal provisioning among closely related sister species. In this study, we present a detailed examination of a new reference transcriptome sequence for the live-bearing, matrotrophic fish, Poeciliopsis prolifica, from multiple-tissue RNA-seq data. We describe the genetic components active in liver, brain, late-stage embryo, and the maternal placental/ovarian complex, as well as associated patterns of positive selection in a suite of orthologous genes found in fishes. Results indicate the expression of many signaling transcripts, "non-coding" sequences and repetitive elements in the maternal placental/ovarian complex. Moreover, patterns of positive selection in protein sequence evolution were found associated with live-bearing fishes, generally, and the placental P. prolifica, specifically, that appear independent of the general live-bearer lifestyle. Much of the observed patterns of gene expression and positive selection are congruent with the evolution of placentation in fish functionally converging with mammalian placental evolution and with the patterns of rapid evolution facilitated by the teleost-specific whole genome duplication event.

Changing MADS-Box Transcription Factor Protein-Protein Interactions as a Mechanism for Generating Floral Morphological Diversity

Flowers display fantastic morphological diversity. Despite extreme variability in form, floral organ identity is specified by a core set of deeply conserved proteins-the floral MADS-box transcription factors. This indicates that while core gene function has been maintained, MADS-box transcription factors have evolved to regulate different downstream genes. Thus, the evolution of gene regulation downstream of the MADS-box transcription factors is likely central to the evolution of floral form. Gene regulation is determined by the combination of transcriptional regulators present at a particular cis-regulatory element at a particular time. Therefore, the interactions between transcription factors can be of profound importance in determining patterns of gene regulation. Here, after a short primer on flowers and floral morphology, I discuss the centrality of protein-protein interactions to MADS-box transcription factor function, and review the evidence that the evolution of MADS-box protein-protein interactions is a key driver in the evolution of gene regulation downstream of the MADS-box genes.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

Taxon-restricted genes at the origin of a novel trait allowing access to a new environment

Taxon-restricted genes make up a considerable proportion of genomes, yet their contribution to phenotypic evolution is poorly understood. We combined gene expression with functional and behavioral assays to study the origin and adaptive value of an evolutionary innovation exclusive to the water strider genus Rhagovelia: the propelling fan. We discovered that two taxon-restricted genes, which we named geisha and mother-of-geisha, specifically control fan development. geisha originated through a duplication event at the base of the Rhagovelia lineage, and both duplicates acquired a novel expression in a specific cell population prefiguring fan development. These gene duplicates played a central role in Rhagovelia's adaptation to a new physical environment, demonstrating that the evolution of taxon-restricted genes can contribute directly to evolutionary novelties that allow access to unexploited ecological niches.

Transforming a transcription factor

A transcription factor that regulates skeleton formation in sea urchin embryos has evolved a new domain that is essential for this process.

Functional divergence of paralogous transcription factors supported the evolution of biomineralization in echinoderms

Alx1 is a pivotal transcription factor in a gene regulatory network that controls skeletogenesis throughout the echinoderm phylum. We performed a structure-function analysis of sea urchin Alx1 using a rescue assay and identified a novel, conserved motif (Domain 2) essential for skeletogenic function. The paralogue of Alx1, Alx4, was not functionally interchangeable with Alx1, but insertion of Domain 2 conferred robust skeletogenic function on Alx4. We used cross-species expression experiments to show that Alx1 proteins from distantly related echinoderms are not interchangeable, although the sequence and function of Domain 2 are highly conserved. We also found that Domain 2 is subject to alternative splicing and provide evidence that this domain was originally gained through exonization. Our findings show that a gene duplication event permitted the functional specialization of a transcription factor through changes in exon-intron organization and thereby supported the evolution of a major morphological novelty.

No Excess of Cis-Regulatory Variation Associated with Intraspecific Selection in Wild Pearl Millet (Cenchrus americanus)

Several studies suggest that cis-regulatory mutations are the favorite target of evolutionary changes, one reason being that cis-regulatory mutations might have fewer deleterious pleiotropic effects than protein-coding mutations. A review of the process also suggests that this bias towards adaptive cis-regulatory variation might be less pronounced at the intraspecific level compared with the interspecific level. In this study, we assessed the contribution of cis-regulatory variation to adaptation at the intraspecific level using populations of wild pearl millet (Cenchrus americanus ssp. monodii) sampled along an environmental gradient in Niger. From RNA sequencing of hybrids to assess allele-specific expression, we identified genes with cis-regulatory divergence between two parental accessions collected in contrasted environmental conditions. This revealed that ∼15% of transcribed genes showed cis-regulatory variation. Intersecting the gene set exhibiting cis-regulatory variation with the gene set identified as targets of selection revealed no excess of cis-acting mutations among the selected genes. We additionally found no excess of cis-regulatory variation among genes associated with adaptive traits. As our approach relied on methods identifying mainly genes submitted to strong selection pressure or with high phenotypic effect, the contribution of cis-regulatory changes to soft selection or polygenic adaptive traits remains to be tested. However our results favor the hypothesis that enrichment of adaptive cis-regulatory divergence builds up over time. For short evolutionary time-scales, cis-acting mutations are not predominantly involved in adaptive evolution associated with strong selective signal.

Developmental plasticity: re-conceiving the genotype

In recent decades, the phenotype of an organism (i.e. its traits and behaviour) has been studied as the outcome of a developmental 'programme' coded in its genotype. This deterministic view is implicit in the Modern Synthesis approach to adaptive evolution as a sorting process among genetic variants. Studies of developmental pathways have revealed that genotypes are in fact differently expressed depending on environmental conditions. Accordingly, the genotype can be understood as a repertoire of potential developmental outcomes or norm of reaction. Reconceiving the genotype as an environmental response repertoire rather than a fixed developmental programme leads to three critical evolutionary insights. First, plastic responses to specific conditions often comprise functionally appropriate trait adjustments, resulting in an individual-level, developmental mode of adaptive variation. Second, because genotypes are differently expressed depending on the environment, the genetic diversity available to natural selection is itself environmentally contingent. Finally, environmental influences on development can extend across multiple generations via cytoplasmic and epigenetic factors transmitted to progeny individuals, altering their responses to their own, immediate environmental conditions and, in some cases, leading to inherited but non-genetic adaptations. Together, these insights suggest a more nuanced understanding of the genotype and its evolutionary role, as well as a shift in research focus to investigating the complex developmental interactions among genotypes, environments and previous environments.