Evolutionary origin of a novel gene expression pattern through co-option of the latent activities of existing regulatory sequences

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

Unified bursting strategies in ectopic and endogenous even-skipped expression patterns

Transcription often occurs in bursts as gene promoters switch stochastically between active and inactive states. Enhancers can dictate transcriptional activity in animal development through the modulation of burst frequency, duration, or amplitude. Previous studies observed that different enhancers can achieve a wide range of transcriptional outputs through the same strategies of bursting control. For example, despite responding to different transcription factors, all even-skipped enhancers increase transcription by upregulating burst frequency and amplitude while burst duration remains largely constant. These shared bursting strategies suggest that a unified molecular mechanism constraints how enhancers modulate transcriptional output. Alternatively, different enhancers could have converged on the same bursting control strategy because of natural selection favoring one of these particular strategies. To distinguish between these two scenarios, we compared transcriptional bursting between endogenous and ectopic gene expression patterns. Because enhancers act under different regulatory inputs in ectopic patterns, dissimilar bursting control strategies between endogenous and ectopic patterns would suggest that enhancers adapted their bursting strategies to their trans-regulatory environment. Here, we generated ectopic even-skipped transcription patterns in fruit fly embryos and discovered that bursting strategies remain consistent in endogenous and ectopic even-skipped expression. These results provide evidence for a unified molecular mechanism shaping even-skipped bursting strategies and serve as a starting point to uncover the realm of strategies employed by other enhancers.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Cysteine Enrichment Mediates Co-Option of Uricase in Reptilian Skin and Transition to Uricotelism

Uric acid is the main means of nitrogen excretion in uricotelic vertebrates (birds and reptiles) and the end product of purine catabolism in humans and a few other mammals. While uricase is inactivated in mammals unable to degrade urate, the presence of orthologous genes without inactivating mutations in avian and reptilian genomes is unexplained. Here we show that the Gallus gallus gene we name cysteine-rich urate oxidase (CRUOX) encodes a functional protein representing a unique case of cysteine enrichment in the evolution of vertebrate orthologous genes. CRUOX retains the ability to catalyze urate oxidation to hydrogen peroxide and 5-hydroxyisourate (HIU), albeit with a 100-fold reduced efficiency. However, differently from all uricases hitherto characterized, it can also facilitate urate regeneration from HIU, a catalytic property that we propose depends on its enrichment in cysteine residues. X-ray structural analysis highlights differences in the active site compared to known orthologs and suggests a mechanism for cysteine-mediated self-aggregation under H2O2-oxidative conditions. Cysteine enrichment was concurrent with the transition to uricotelism and a shift in gene expression from the liver to the skin where CRUOX is co-expressed with β-keratins. Therefore, the loss of urate degradation in amniotes has followed opposite evolutionary trajectories: while uricase has been eliminated by pseudogenization in some mammals, it has been repurposed as a redox-sensitive enzyme in the reptilian skin.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Evolutionary mechanisms modulating the mammalian skull development

Mammals possess impressive craniofacial variation that mirrors their adaptation to diverse ecological niches, feeding behaviour, physiology and overall lifestyle. The spectrum of craniofacial geometries is established mainly during embryonic development. The formation of the head represents a sequence of events regulated on genomic, molecular, cellular and tissue level, with each step taking place under tight spatio-temporal control. Even minor variations in timing, position or concentration of the molecular drivers and the resulting events can affect the final shape, size and position of the skeletal elements and the geometry of the head. Our knowledge of craniofacial development increased substantially in the last decades, mainly due to research using conventional vertebrate model organisms. However, how developmental differences in head formation arise specifically within mammals remains largely unexplored. This review highlights three evolutionary mechanisms acknowledged to modify ontogenesis: heterochrony, heterotopy and heterometry. We present recent research that links changes in developmental timing, spatial organization or gene expression levels to the acquisition of species-specific skull morphologies. We highlight how these evolutionary modifications occur on the level of the genes, molecules and cellular processes, and alter conserved developmental programmes to generate a broad spectrum of skull shapes characteristic of the class Mammalia. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

Differentially Active and Conserved Neural Enhancers Define Two Forms of Adaptive Noncoding Evolution in Humans

The human and chimpanzee genomes are strikingly similar, but our neural phenotypes are very different. Many of these differences are likely driven by changes in gene expression, and some of those changes may have been adaptive during human evolution. Yet, the relative contributions of positive selection on regulatory regions or other functional regulatory changes are unclear. Where are these changes located throughout the human genome? Are functional regulatory changes near genes or are they in distal enhancer regions? In this study, we experimentally combined both human and chimpanzee cis-regulatory elements (CREs) that showed either (1) signs of accelerated evolution in humans or (2) that have been shown to be active in the human brain. Using a massively parallel reporter assay, we tested the ability of orthologous human and chimpanzee CREs to activate transcription in induced pluripotent stem-cell-derived neural progenitor cells and neurons. With this assay, we identified 179 CREs with differential activity between human and chimpanzee; in contrast, we found 722 CREs with signs of positive selection in humans. Selection and differentially expressed CREs strikingly differ in level of expression, size, and genomic location. We found a subset of 69 CREs in loci with genetic variants associated with neuropsychiatric diseases, which underscores the consequence of regulatory activity in these loci for proper neural development and function. By combining CREs that either experienced recent selection in humans or CREs that are functional brain enhancers, presents a novel way of studying the evolution of noncoding elements that contribute to human neural phenotypes.

Developmental biology is poised to discover altogether new principles in biology in the 21st century

In the 20th century, developmental biology spearheaded a revolution in our understanding of complex biological problems. Its success rests in great part on a truly unique approach that has recruited a diversity of systems and research organisms rather than focusing on isolated cells or molecules, while also employing a wide variety of technological and intellectual approaches. But what will developmental biology contribute to this century? Advances in technology and instrumentation are presently moving at neck-breaking speed and herald the advent of an age of technological wonders in which previously inaccessible biology is now tangibly within our grasps. For instance, single-cell RNAseq has revealed novel, transient cell states in both stem and differentiated cells that are specified by defined changes in gene expression frequency during regeneration. Additionally, genome-wide epigenetic analyses combined with gene editing and transgenic methodologies have identified the existence of regeneration responsive enhancers in adult vertebrate tissues. These circumstances combined with our discipline’s diversity of experimental and intellectual approaches offer unimaginable opportunities for developmental biologists not only to discover new biology but also to reveal entirely new principles of biology.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Repression precedes independent evolutionary gains of a highly specific gene expression pattern

Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns.

Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.

Hidden paths to endless forms most wonderful: parasite-blind diversification of host quality

Evolutionary diversification can occur in allopatry or sympatry, can be driven by selection or unselected, and can be phenotypically manifested immediately or remain latent until manifested in a newly encountered environment. Diversification of host-parasite interactions is frequently studied in the context of intrinsically selective coevolution, but the potential for host-parasite interaction phenotypes to diversify latently during parasite-blind host evolution is rarely considered. Here, we use a social bacterium experimentally adapted to several environments in the absence of phage to analyse allopatric diversification of host quality-the degree to which a host population supports a viral epidemic. Phage-blind evolution reduced host quality overall, with some bacteria becoming completely resistant to growth suppression by phage. Selective-environment differences generated only mild divergence in host quality. However, selective environments nonetheless played a major role in shaping evolution by determining the degree of stochastic diversification among replicate populations within treatments. Ancestral motility genotype was also found to strongly shape patterns of latent host-quality evolution and diversification. These outcomes show that (i) adaptive landscapes can differ in how they constrain stochastic diversification of a latent phenotype and (ii) major effects of selection on biological diversification can be missed by focusing on trait means. Collectively, our findings suggest that latent-phenotype evolution should inform host-parasite evolution theory and that diversification should be conceived broadly to include latent phenotypes.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Changes in regeneration-responsive enhancers shape regenerative capacities in vertebrates

Vertebrates vary in their ability to regenerate, and the genetic mechanisms underlying such disparity remain elusive. Comparative epigenomic profiling and single-cell sequencing of two related teleost fish uncovered species-specific and evolutionarily conserved genomic responses to regeneration. The conserved response revealed several regeneration-responsive enhancers (RREs), including an element upstream to inhibin beta A (inhba), a known effector of vertebrate regeneration. This element activated expression in regenerating transgenic fish, and its genomic deletion perturbed caudal fin regeneration and abrogated cardiac regeneration altogether. The enhancer is present in mammals, shares functionally essential activator protein 1 (AP-1)-binding motifs, and responds to injury, but it cannot rescue regeneration in fish. This work suggests that changes in AP-1-enriched RREs are likely a crucial source of loss of regenerative capacities in vertebrates.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

On the specificity of gene regulatory networks: How does network co-option affect subsequent evolution?

The process of multicellular organismal development hinges upon the specificity of developmental programs: for different parts of the organism to form unique features, processes must exist to specify each part. This specificity is thought to be hardwired into gene regulatory networks, which activate cohorts of genes in particular tissues at particular times during development. However, the evolution of gene regulatory networks sometimes occurs by mechanisms that sacrifice specificity. One such mechanism is network co-option, in which existing gene networks are redeployed in new developmental contexts. While network co-option may offer an efficient mechanism for generating novel phenotypes, losses of tissue specificity at redeployed network genes could restrict the ability of the affected traits to evolve independently. At present, there has not been a detailed discussion regarding how tissue specificity of network genes might be altered due to gene network co-option at its initiation, as well as how trait independence can be retained or restored after network co-option. A lack of clarity about network co-option makes it more difficult to speculate on the long-term evolutionary implications of this mechanism. In this review, we will discuss the possible initial outcomes of network co-option, outline the mechanisms by which networks may retain or subsequently regain specificity after network co-option, and comment on some of the possible evolutionary consequences of network co-option. We place special emphasis on the need to consider selectively-neutral outcomes of network co-option to improve our understanding of the role of this mechanism in trait evolution.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Red Light/Green Light, a Dual Fluorescent Protein Reporter System To Study Enhancer-Promoter Specificity in Drosophila

Enhancers activate gene transcription in spatial and temporal patterns by interactions with gene promoters. These elements typically reside distal to their target promoter, with which they must interact selectively. Additional elements may contribute to enhancer-promoter specificity, including remote control element sequences within enhancers, tethering elements near promoters, and insulator/boundary elements that disrupt off-target interactions. However, few of these elements have been mapped, and as a result, the mechanisms by which these elements interact remain poorly understood. One impediment is their method of study, namely reporter transgenes in which enhancers are placed adjacent to a heterologous promoter, which may circumvent mechanisms controlling enhancer-promoter specificity and long-range interactions. Here, we report an optimized dual reporter transgene system in Drosophila melanogaster that allows the simultaneous comparison of an enhancer’s ability to activate proximal and distal fluorescent reporter genes. Testing a panel of fluorescent transgenes in vivo, we found a two-protein combination that allows simultaneous measurement with minimal detection interference. We note differences among four tested enhancers in their ability to regulate a distally placed reporter transgene. These results suggest that enhancers differ in their requirements for promoter interaction and raise important practical considerations when studying enhancer function.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Evolution: How Many Phenotypes Do Regulatory Mutations Affect?

Mutations in gene regulatory regions are thought to play an important role in the evolution of morphological structures. This is largely due to their minimal pleiotropic effects, limiting their impact to one particular body part. A recent study finds that one such regulatory mutation may affect two particular morphological structures.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

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.

Chromatin Evolution-Key Innovations Underpinning Morphological Complexity

The history of life consists of a series of major evolutionary transitions, including emergence and radiation of complex multicellular eukaryotes from unicellular ancestors. The cells of multicellular organisms, with few exceptions, contain the same genome, however, their organs are composed of a variety of cell types that differ in both structure and function. This variation is largely due to the transcriptional activity of different sets of genes in different cell types. This indicates that complex transcriptional regulation played a key role in the evolution of complexity in eukaryotes. In this review, we summarize how gene duplication and subsequent evolutionary innovations, including the structural evolution of nucleosomes and chromatin-related factors, contributed to the complexity of the transcriptional system and provided a basis for morphological diversity.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Evolutionary Origins of Pseudogenes and Their Association with Regulatory Sequences in Plants

Pseudogenes (Ψs), nonfunctional relatives of functional genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Evolutionary analysis provides clues to Ψ origins and effects on gene regulation. However, few systematic studies of plant Ψs have been conducted, hampering comparative analyses. Here, we examined the origin, evolution, and expression patterns of Ψs and their relationships with noncoding sequences in seven angiosperm plants. We identified ∼250,000 Ψs, most of which are more lineage specific than protein-coding genes. The distribution of Ψs on the chromosome indicates that genome recombination may contribute to Ψ elimination. Most Ψs evolve rapidly in terms of sequence and expression levels, showing tissue- or stage-specific expression patterns. We found that a surprisingly large fraction of nontransposable element regulatory noncoding RNAs (microRNAs and long noncoding RNAs) originate from transcription of Ψ proximal upstream regions. We also found that transcription factor binding sites preferentially occur in putative Ψ proximal upstream regions compared with random intergenic regions, suggesting that Ψs have conditioned genome evolution by providing transcription factor binding sites that serve as promoters and enhancers. We therefore propose that rapid rewiring of Ψ transcriptional regulatory regions is a major mechanism driving the origin of novel regulatory modules.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Enhancer reprogramming in mammalian genomes

Background

Transcription factor binding site (TFBS) loss, gain, and reshuffling within the sequence of a regulatory element could alter the function of that regulatory element. Some of the changes will be detrimental to the fitness of the species and will result in gradual removal from a population, while other changes would be either beneficial or just a part of genetic drift and end up being fixed in a population. This “reprogramming” of regulatory elements results in modification of the gene regulatory landscape during evolution.

Results

We identified reprogrammed enhancers (RPEs) by comparing the distribution of tissue-specific enhancers in the human and mouse genomes. We observed that around 30% of mammalian enhancers have been reprogrammed after the human-mouse speciation. In 79% of cases, the reprogramming of an enhancer resulted in a quantifiably different expression of a flanking gene. In the case of the Thy-1 cell surface antigen gene, for example, enhancer reprogramming is associated with cortex to thymus change in gene expression. To understand the mechanisms of enhancer reprogramming, we profiled the evolutionary changes in the TFBS enhancer content and found that enhancer reprogramming took place through the acquisition of new TFBSs in 72% of reprogramming events.

Conclusions

Our results suggest that enhancer reprogramming takes place within well-established regulatory loci with RPEs contributing additively to fine-tuning of the gene regulatory program in mammals. We also found evidence for acquisition of novel gene function through enhancer reprogramming, which allows expansion of gene regulatory landscapes into new regulatory domains.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2343-7) contains supplementary material, which is available to authorized users.

Complexin in ivermectin resistance in body lice

Ivermectin has emerged as very promising pediculicide, particularly in cases of resistance to commonly used pediculicides. Recently, however, the first field-evolved ivermectin-resistance in lice was reported. To gain insight into the mechanisms underlying ivermectin-resistance, we both looked for mutations in the ivermectin-target site (GluCl) and searched the entire proteome for potential new loci involved in resistance from laboratory susceptible and ivermectin-selected resistant body lice. Polymorphism analysis of cDNA GluCl showed no non-silent mutations. Proteomic analysis identified 22 differentially regulated proteins, of which 13 were upregulated and 9 were downregulated in the resistant strain. We evaluated the correlation between mRNA and protein levels by qRT-PCR and found that the trend in transcriptional variation was consistent with the proteomic changes. Among differentially expressed proteins, a complexin i.e. a neuronal protein which plays a key role in regulating neurotransmitter release, was shown to be the most significantly down-expressed in the ivermectin-resistant lice. Moreover, DNA-mutation analysis revealed that some complexin transcripts from resistant lice gained a premature stop codon, suggesting that this down-expression might be due, in part, to secondary effects of a nonsense mutation inside the gene. We further confirmed the association between complexin and ivermectin-resistance by RNA-interfering and found that knocking down the complexin expression induces resistance to ivermectin in susceptible lice. Our results provide evidence that complexin plays a significant role in regulating ivermectin resistance in body lice and represents the first evidence that links complexin to insecticide resistance.

Diverse Cis-Regulatory Mechanisms Contribute to Expression Evolution of Tandem Gene Duplicates

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression.

The Evolution of Gene Expression in cis and trans

There is abundant variation in gene expression between individuals, populations, and species. The evolution of gene regulation and expression within and between species is thought to frequently contribute to adaptation. Yet considerable evidence suggests that the primary evolutionary force acting on variation in gene expression is stabilizing selection. We review here the results of recent studies characterizing the evolution of gene expression occurring in cis (via linked polymorphisms) or in trans (through diffusible products of other genes) and their contribution to adaptation and response to the environment. We review the evidence for buffering of variation in gene expression at the level of both transcription and translation, and the possible mechanisms for this buffering. Lastly, we summarize unresolved questions about the evolution of gene regulation.

Evolutionary emergence of the rac3b/rfng/sgca regulatory cluster refined mechanisms for hindbrain boundaries formation

Developmental programs often rely on parallel morphogenetic mechanisms that guarantee precise tissue architecture. While redundancy constitutes an obvious selective advantage, little is known on how novel morphogenetic mechanisms emerge during evolution. In zebrafish, rhombomeric boundaries behave as an elastic barrier, preventing cell intermingling between adjacent compartments. Here, we identify the fundamental role of the small-GTPase Rac3b in actomyosin cable assembly at hindbrain boundaries. We show that the novel rac3b/rfng/sgca regulatory cluster, which is specifically expressed at the boundaries, emerged in the Ostariophysi superorder by chromosomal rearrangement that generated new cis-regulatory interactions. By combining 4C-seq, ATAC-seq, transgenesis, and CRISPR-induced deletions, we characterized this regulatory domain, identifying hindbrain boundary-specific cis-regulatory elements. Our results suggest that the capacity of boundaries to act as an elastic mesh for segregating rhombomeric cells evolved by cooption of critical genes to a novel regulatory block, refining the mechanisms for hindbrain segmentation.

Enhancer adoption caused by genomic insertion elicits interdigital Shh expression and syndactyly in mouse

In this study, we reexamined an old mouse mutant named Hammer toe (Hm), which arose spontaneously almost a half century ago and exhibits a limb phenotype with webbing. We revealed that a 150-kb noncoding genomic fragment that was originally located in chromosome 14 has been inserted into a genomic region proximal to Sonic hedgehog (Shh), located in chromosome 5. This inserted fragment possesses enhancer activity to induce Shh expression in the interdigital regions in Hm, which in turn down-regulates bone morphogenetic protein signaling and eventually results in syndactyly and web formation. Since the donor fragment residing in chromosome 14 has enhancer activity to induce interdigital gene expression, the Hm mutation appears to be an archetypal case of enhancer adoption.

Acquisition of new cis-regulatory elements (CREs) can cause alteration of developmental gene regulation and may introduce morphological novelty in evolution. Although structural variation in the genome generated by chromosomal rearrangement is one possible source of new CREs, only a few examples are known, except for cases of retrotransposition. In this study, we show the acquisition of novel regulatory sequences as a result of large genomic insertion in the spontaneous mouse mutation Hammer toe (Hm). Hm mice exhibit syndactyly with webbing, due to suppression of interdigital cell death in limb development. We reveal that, in the Hm genome, a 150-kb noncoding DNA fragment from chromosome 14 is inserted into the region upstream of the Sonic hedgehog (Shh) promoter in chromosome 5. Phenotyping of mouse embryos with a series of CRISPR/Cas9-aided partial deletion of the 150-kb insert clearly indicated that two different regions are necessary for the syndactyly phenotype of Hm. We found that each of the two regions contains at least one enhancer for interdigital regulation. These results show that a set of enhancers brought by the large genomic insertion elicits the interdigital Shh expression and the Hm phenotype. Transcriptome analysis indicates that ectopic expression of Shh up-regulates Chordin (Chrd) that antagonizes bone morphogenetic protein signaling in the interdigital region. Indeed, Chrd-overexpressing transgenic mice recapitulated syndactyly with webbing. Thus, the Hm mutation provides an insight into enhancer acquisition as a source of creation of novel gene regulation.

Enhancer evolution and the origins of morphological novelty

A central goal of evolutionary biology is to understand the genetic origin of morphological novelties-i.e. anatomical structures unique to a taxonomic group. Elaboration of morphology during development depends on networks of regulatory genes that activate patterned gene expression through transcriptional enhancer regions. We summarize recent case studies and genome-wide investigations that have uncovered diverse mechanisms though which new enhancers arise. We also discuss how these enhancer-originating mechanisms have clarified the history of genetic networks underlying diversification of genital structures in flies, limbs and neural crest in chordates, and plant leaves. These studies have identified enhancers that were pivotal for morphological divergence and highlighted how novel genetic networks shaping form emerged from pre-existing ones.

Developmental constraints shape the evolution of the nematode mid-developmental transition

Evolutionary theory assumes that genetic variation is uniform and gradual in nature, yet morphological and gene expression studies have revealed that different life-stages exhibit distinct levels of cross-species conservation. In particular, a stage in mid-embryogenesis is highly conserved across species of the same phylum, suggesting that this stage is subject to developmental constraints, either by increased purifying selection or by a strong mutational bias. An alternative explanation, however, holds that the same 'hourglass' pattern of variation may result from increased positive selection at the earlier and later stages of development. To distinguish between these scenarios, we examined gene expression variation in a population of the nematode Caenorhabditis elegans using an experimental design that eliminated the influence of positive selection. By measuring gene expression for all genes throughout development in 20 strains, we found that variations were highly uneven throughout development, with a significant depletion during mid-embryogenesis. In particular, the family of homeodomain transcription factors, whose expression generally coincides with mid-embryogenesis, evolved under high constraint. Our data further show that genes responsible for the integration of germ layers during morphogenesis are the most constrained class of genes. Together, these results provide strong evidence for developmental constraints as the mechanism underlying the hourglass model of animal evolution. Understanding the pattern and mechanism of developmental constraints provides a framework to understand how evolutionary processes have interacted with embryogenesis and led to the diversity of animal life on Earth.

Cis-regulatory PLETHORA promoter elements directing root and nodule expression are conserved between Arabidopsis thaliana and Medicago truncatula

Nodules are unique organs formed on roots of legumes by soil-borne bacteria, collectively known as rhizobium. Recently, we have shown that orthologs of the AINTEGUMENTA-like (AIL) AP2 transcription factors PLETHORA (PLT) 1 to 4, that redundantly regulate Arabidopsis thaliana root development are involved in root and nodule growth in Medicago truncatula. Hence, it is conceivable that rhizobium has co-opted these genes for nodule development. Whether this co-option requires the presence of specific cis-elements in the promoters and/or specialization of PLT protein function is not clear. Here, we analyzed the qualitative expression patterns of the Arabidopsis PLT1 to 4 promoters in Medicago roots and nodules and compared these with the described expression patterns of the Medicago PLT genes. Our studies reveal that the expression patterns of the investigated promoters and their Medicago orthologs are very similar, indicating that at least all cis-elements regulating spatial PLT expression are conserved among the Arabidopsis and Medicago PLT1 to 4 promoters.

The evolutionary origination of a novel expression pattern through an extreme heterochronic shift

The evolutionary origins of morphological structures are thought to often depend upon the redeployment of old genes into new developmental settings. Although many examples of cis-regulatory divergence have shown how pre-existing patterns of gene expression have been altered, only a small number of case studies have traced the origins of cis-regulatory elements that drive new expression domains. Here, we elucidate the evolutionary history of a novel expression pattern of the yellow gene within the Zaprionus genus of fruit flies. We observed a unique pattern of yellow transcript accumulation in the wing disc during the third larval instar, a stage that precedes its typical expression pattern associated with cuticular melanization by about a week. The region of the Zaprionus wing disc that expresses yellow subsequently develops into a portion of the thorax, a tissue for which yellow expression has been reported for several fruit fly species. Tests of GFP reporter transgenes containing the Zaprionus yellow regulatory region revealed that the wing disc pattern arose by changes in the cis-regulatory region of yellow. Moreover, the wing disc enhancer activity of yellow depends upon a short conserved sequence with ancestral thoracic functions, suggesting that the pupal thorax regulatory sequence was genetically reprogrammed to drive expression that commences much earlier during development. These results highlight how novel domains of gene expression may arise by extreme shifts in timing during the origins of novel traits.

Enhancer modularity and the evolution of new traits

Animals have modular cis-regulatory regions in their genomes, and expression of a single gene is often regulated by multiple enhancers residing in such a region. In the laboratory, and also in natural populations, loss of an enhancer can result in a loss of gene expression. Although only a few examples have been well characterized to date, some studies have suggested that an evolutionary gain of a new enhancer function can establish a new gene expression domain. Our recent study showed that Drosophila guttifera has more enhancers and additional expression domains of the wingless gene during the pupal stage, compared to D. melanogaster, and that these new features appear to have evolved in the ancestral lineage leading to D. guttifera. (1) Gain of a new expression domain of a developmental regulatory gene (toolkit gene), such as wingless, can cause co-option of the expression of its downstream genes to the new domain, resulting in duplication of a preexisting structure at this new body position. Recently, with the advancement of evo-devo studies, we have learned that the developmental regulatory systems are strikingly similar across various animal taxa, in spite of the great diversity of the animals' morphology. Even behind "new" traits, co-options of essential developmental genes from known systems are very common. We previously provided concrete evidence of gains of enhancer activities of a developmental regulatory gene underlying gains of new traits. (1) Broad occurrence of this scenario is testable and should be validated in the future.

The White-Knight Hypothesis, or Does the Environment Limit Innovations?

Organisms often harbor latent traits that are byproducts of other adaptations. Such latent traits are not themselves adaptive but can become adaptive in the right environment. Here I discuss several examples of such traits. Their abundance suggests that environmental change rather than new mutations might often limit the origin of evolutionary adaptations and innovations. This is important, because environments can change much faster than new mutations arise. I introduce a conceptual model that distinguishes between mutation-limited and environment-limited trait origins and suggest how experiments could help discriminate between them. Wherever latent traits are plentiful, ecology rather than genetics might determine how fast new adaptations originate and thus how fast adaptive Darwinian evolution proceeds.

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

Coupled enhancer and coding sequence evolution of a homeobox gene shaped leaf diversity

In this study, Vuolo et al. investigate the mechanisms underlying the genetic basis for morphological diversity in leaf shape. They show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base.

Here we investigate mechanisms underlying the diversification of biological forms using crucifer leaf shape as an example. We show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base. A single amino acid substitution evolved together with this regulatory change, which reduced RCO protein stability, preventing pleiotropic effects caused by its altered gene expression. We detected hallmarks of positive selection in these evolved regulatory and coding sequence variants and showed that modulating RCO activity can improve plant physiological performance. Therefore, interplay between enhancer and coding sequence evolution created a potentially adaptive path for morphological evolution.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Sensitivity of Allelic Divergence to Genomic Position: Lessons from the Drosophila tan Gene

To identify genetic variants underlying changes in phenotypes within and between species, researchers often utilize transgenic animals to compare the function of alleles in different genetic backgrounds. In Drosophila, targeted integration mediated by the ΦC31 integrase allows activity of alternative alleles to be compared at the same genomic location. By using the same insertion site for each transgene, position effects are generally assumed to be controlled for because both alleles are surrounded by the same genomic context. Here, we test this assumption by comparing the activity of tan alleles from two Drosophila species, D. americana and D. novamexicana, at five different genomic locations in D. melanogaster. We found that the relative effects of these alleles varied among insertion sites, with no difference in activity observed between them at two sites. One of these sites simply silenced both transgenes, but the other allowed expression of both alleles that was sufficient to rescue a mutant phenotype yet failed to reveal the functional differences between the two alleles. These results suggest that more than one insertion site should be used when comparing the activity of transgenes because failing to do so could cause functional differences between alleles to go undetected.

Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.

A de novo silencer causes elimination of MITF-M expression and profound hearing loss in pigs

Background

Genesis of novel gene regulatory modules is largely responsible for morphological and functional evolution. De novo generation of novel cis-regulatory elements (CREs) is much rarer than genomic events that alter existing CREs such as transposition, promoter switching or co-option. Only one case of de novo generation has been reported to date, in fish and without involvement of phenotype alteration. Yet, this event likely occurs in other animals and helps drive genetic/phenotypic variation.

Results

Using a porcine model of spontaneous hearing loss not previously characterized we performed gene mapping and mutation screening to determine the genetic foundation of the phenotype. We identified a mutation in the non-regulatory region of the melanocyte-specific promoter of microphthalmia-associated transcription factor (MITF) gene that generated a novel silencer. The consequent elimination of expression of the MITF-M isoform led to early degeneration of the intermediate cells of the cochlear stria vascularis and profound hearing loss, as well as depigmentation, all of which resemble the typical phenotype of Waardenburg syndrome in humans. The mutation exclusively affected MITF-M and no other isoforms. The essential function of Mitf-m in hearing development was further validated using a knock-out mouse model.

Conclusions

Elimination of the MITF-M isoform alone is sufficient to cause deafness and depigmentation. To our knowledge, this study provides the first evidence of a de novo CRE in mammals that produces a systemic functional effect.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0273-2) contains supplementary material, which is available to authorized users.

Simulating evolution of protein complexes through gene duplication and co-option

We present a model of the evolution of protein complexes with novel functions through gene duplication, mutation, and co-option. Under a wide variety of input parameters, digital organisms evolve complexes of 2-5 bound proteins which have novel functions but whose component proteins are not independently functional. Evolution of complexes with novel functions happens more quickly as gene duplication rates increase, point mutation rates increase, protein complex functional probability increases, protein complex functional strength increases, and protein family size decreases. Evolution of complexity is inhibited when the metabolic costs of making proteins exceeds the fitness gain of having functional proteins, or when point mutation rates get so large the functional proteins undergo deleterious mutations faster than new functional complexes can evolve.