Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene

Evolution of gene regulatory networks in insects

Changes in gene regulatory networks (GRNs) underlying the evolution of traits have been intensively studied, with insects providing excellent model cases. In studies using Drosophila, butterflies, and other insects, several well-known cases have shown that changes in the cis-regulatory region of a gene controlling a trait can result in the co-option of the gene for a role different from that in its original developmental context. When the expression of a regulatory gene that controls the expression of multiple downstream genes is altered, the expression of these downstream genes changes accordingly, representing the simplest form of GRN co-option. Many studies have explored the applicability of this model to the acquisition of new traits, yielding substantial insights. However, no study has yet comprehensively elucidated the co-option of a GRN or the evolution of a network architecture, including associated genes and their regulatory relationships. In the near future, the use of single-cell multiomics and machine learning will allow for larger-scale data analysis, leading to a better understanding of the evolution of traits through the evolution of GRNs.

Genome-wide identification of chitinase gene family in Hordeum vulgare: insights into stress response mechanisms and evolutionary dynamics

BACKGROUND

Chitinase, a key enzyme family within the pathogenesis-related (PR) protein, plays a crucial role in plant defense by degrading chitin, a major component of fungal cell walls. The HvCHT genes in barley are involved in responses to biotic and abiotic stresses, although their full range of functions is not yet fully understood.

RESULTS

In this study, we identified 24 potential HvCHT genes through a genome-wide analysis. The comparative synteny analysis showed conserved relationships between HvCHT genes and their homologs in Sorghum bicolor, Oryza sativa, and Arabidopsis thaliana. Chromosomal mapping, gene structure, characterization, protein motif analysis, and miRNA regulation were performed to gain insight into the genetic structures of these genes. Segmental duplication events observed in the HvCHT family suggest an important role in the evolutionary development of these genes. Additionally, cis-regulatory element analysis revealed the presence of light-responsive elements, and regulators for Abscisic acid, methyl jasmonate (MeJA), salicylic acid, and gibberellins, indicating potential involvement in stress responses. Transcriptomic data showed differential expression of HvCHT genes in response to salt stress, with distinct patterns observed in leaf and root tissues. Furthermore, the genes defensive responses to drought stress and Fusarium infection were characterized across multiple time points. Notably, qRT-PCR analysis confirmed the upregulation of HvCHT1, HvCHT4, and HvCHT17, highlighting their potential involvement in stress-related pathways.

CONCLUSION

These findings provide a comprehensive overview of the HvCHT genes role in barley defense mechanisms, underlining their regulatory functions in biotic and abiotic stressors. The results lay the groundwork for future functional studies on HvCHT genes, with the potential to enhance stress tolerance in crops.

CLINICAL TRIAL NUMBER

Not applicable.

Molecular signature of primate astrocytes reveals pathways and regulatory changes contributing to human brain evolution

Astrocytes contribute to the development and regulation of the higher-level functions of the brain, the critical targets of evolution. However, how astrocytes evolve in primates is unsettled. Here, we obtain human, chimpanzee, and macaque induced pluripotent stem-cell-derived astrocytes (iAstrocytes). Human iAstrocytes are bigger and more complex than the non-human primate iAstrocytes. We identify new loci contributing to the increased human astrocyte. We show that genes and pathways implicated in long-range intercellular signaling are activated in the human iAstrocytes and partake in controlling iAstrocyte complexity. Genes downregulated in human iAstrocytes frequently relate to neurological disorders and were decreased in adult brain samples. Through regulome analysis and machine learning, we uncover that functional activation of enhancers coincides with a previously unappreciated, pervasive gain of "stripe" transcription factor binding sites. Altogether, we reveal the transcriptomic signature of primate astrocyte evolution and a mechanism driving the acquisition of the regulatory potential of enhancers.

Novelty versus innovation of gene regulatory elements in human evolution and disease

It is not currently understood how much of human evolution is due to modifying existing functional elements in the genome versus forging novel elements from nonfunctional DNA. Many early experiments that aimed to assign genetic changes on the human lineage to their resulting phenotypic change have focused on mutations that modify existing elements. However, a number of recent studies have highlighted the potential ease and importance of forging novel gene regulatory elements from nonfunctional sequences on the human lineage. In this review, we distinguish gene regulatory element novelty from innovation. We propose definitions for these terms and emphasize their importance in studying the genetic basis of human uniqueness. We discuss why the forging of novel regulatory elements may have been less emphasized during the previous decades, and why novel regulatory elements are likely to play a significant role in both human adaptation and disease.

Wing spot in a tropical and a temperate drosophilid: C = C enrichment and conserved thermal response

Wings are primarily used in flight but also play a role in mating behaviour in many insects. Drosophila species exhibit a variety of pigmentation patterns on their wings. In some sexually dimorphic Drosophilids, a pigmented spot pattern is found at the top-right edge of the male wings. Our understanding of wing spot thermal plasticity in sexually dimorphic species is limited with wing spots being primarily associated with sexual selection. Here, we investigated the wing pigmentation response of two species with wing spots: D. biarmipes and D. suzukii species to thermal variation. We exposed freshly hatched larvae of both the species to three different growth temperatures and checked for wing pigmentation in adult males. Our results indicate wing pigmentation is a plastic trait in the species studied and that wing pigmentation is negatively correlated with higher temperature. In both species, wings were darker at lower temperature compared to higher temperature. Further, D. suzukii exhibits darker wing pigmentation compared to D. biarmipes. Variation in wing pigmentation in both D. suzukii and D. biarmipes could reflect habitat level differences; indicating a strong G*E interaction. Raman spectral analysis indicated a shift in chemical profiles of pigmented vs. non-pigmented areas of the wing. The wing spot was found enriched with carbon-carbon double-bond compared to the non-pigmented wing area. We report that C = C formation in spotted area is thermally controlled and conserved in two members of the suzukii subgroup i.e. D. biarmipes and D. suzukii. Our study indicated a conserved mechanism of the spot formation in two Drosophila species coming from contrasting distribution ranges.

Disentangling variational bias: the roles of development, mutation, and selection

The extraordinary diversity and adaptive fit of organisms to their environment depends fundamentally on the availability of variation. While most population genetic frameworks assume that random mutations produce isotropic phenotypic variation, the distribution of variation available to natural selection is more restricted, as the distribution of phenotypic variation is affected by a range of factors in developmental systems. Here, we revisit the concept of developmental bias - the observation that the generation of phenotypic variation is biased due to the structure, character, composition, or dynamics of the developmental system - and argue that a more rigorous investigation into the role of developmental bias in the genotype-to-phenotype map will produce fundamental insights into evolutionary processes, with potentially important consequences on the relation between micro- and macro-evolution. We discuss the hierarchical relationships between different types of variational biases, including mutation bias and developmental bias, and their roles in shaping the realized phenotypic space. Furthermore, we highlight the challenges in studying variational bias and propose potential approaches to identify developmental bias using modern tools.

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.

Evolution of the regulation of developmental gene expression in blind Mexican cavefish

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

A chromosome-level genome assembly of Drosophila madeirensis, a fruit fly species endemic to the island of Madeira

Drosophila subobscura is distributed across Europe, the Near East, and the Americas, while its sister species, Drosophila madeirensis, is endemic to the island of Madeira in the Atlantic Ocean. D. subobscura is known for its strict light-dependence in mating and its unique courtship displays, including nuptial gift-giving. D. subobscura has also attracted the interest of researchers because of its abundant variations in chromosomal polymorphisms correlated to the latitude and season, which have been used as a tool to track global climate warming. Although D. madeirensis can be an important resource for understanding the evolutionary underpinning of these genetic characteristics of D. subobscura, little work has been done on the biology of this species. Here, we used a HiFi long-read sequencing data set to produce a de novo genome assembly for D. madeirensis. This assembly comprises a total of 111 contigs spanning 135.5 Mb and has an N50 of 24.2 Mb and a BUSCO completeness score of 98.6%. Each of the 6 chromosomes of D. madeirensis consisted of a single contig except for some centromeric regions. Breakpoints of the chromosomal inversions between D. subobscura and D. madeirensis were characterized using this genome assembly, updating some of the previously identified locations.

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The convergently evolved fish antifreeze proteins provide excellent opportunities to investigate evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner and 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage. Each independently evolved a coding region for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms are consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

Regulatory evolution tuning pigmentation intensity quantitatively in Drosophila

Quantitative variation in attributes such as color, texture, or stiffness dominates morphological diversification. It results from combinations of alleles at many Mendelian loci. Here, we identify an additional source of quantitative variation among species, continuous evolution in a gene regulatory region. Specifically, we examined the modulation of wing pigmentation in a group of fly species and showed that inter-species variation correlated with the quantitative expression of the pigmentation gene yellow. This variation results from an enhancer of yellow determining darkness through species-specific activity. We mapped the divergent activities between two sister species and found the changes to be broadly distributed along the enhancer. Our results demonstrate that enhancers can act as dials fueling quantitative morphological diversification by modulating trait properties.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Changes in locus wide repression underlie the evolution of Drosophila abdominal pigmentation

Changes in gene regulation represent an important path to generate developmental differences affecting anatomical traits. Interspecific divergence in gene expression often results from changes in transcription-stimulating enhancer elements. While gene repression is crucial for precise spatiotemporal expression patterns, the relative contribution of repressive transcriptional silencers to regulatory evolution remains to be addressed. Here, we show that the Drosophila pigmentation gene ebony has mainly evolved through changes in the spatial domains of silencers patterning its abdominal expression. By precisely editing the endogenous ebony locus of D. melanogaster, we demonstrate the requirement of two redundant abdominal enhancers and three silencers that repress the redundant enhancers in a patterned manner. We observe a role for changes in these silencers in every case of ebony evolution observed to date. Our findings suggest that negative regulation by silencers likely has an under-appreciated role in gene regulatory evolution.

Evolution of modular and pleiotropic enhancers

Cis-regulatory elements (CREs), or enhancers, are segments of noncoding DNA that regulate the spatial and temporal expression of nearby genes. Sometimes, genes are expressed in more than one tissue, and this can be driven by two main types of CREs: tissue-specific "modular" CREs, where different CREs drive expression of the gene in the different tissues, or by "pleiotropic" CREs, where the same CRE drives expression in the different tissues. In this perspective, we will discuss some of the ways (i) modular and pleiotropic CREs might originate; (ii) propose that modular CREs might derive from pleiotropic CREs via a process of duplication, degeneration, and complementation (the CRE-DDC model); and (iii) propose that hotspot loci of evolution are associated with the origin of modular CREs belonging to any gene in a regulatory network.

Widespread cis- and trans-regulatory evolution underlies the origin, diversification, and loss of a sexually dimorphic fruit fly pigmentation trait

Changes in gene expression are a prominent feature of morphological evolution. These changes occur to hierarchical gene regulatory networks (GRNs) of transcription factor genes that regulate the expression of trait-building differentiation genes. While changes in the expression of differentiation genes are essential to phenotypic evolution, they can be caused by mutations within cis-regulatory elements (CREs) that drive their expression (cis-evolution) or within genes for CRE-interacting transcription factors (trans-evolution). Locating these mutations remains a challenge, especially when experiments are limited to one species that possesses the ancestral or derived phenotype. We investigated CREs that control the expression of the differentiation genes tan and yellow, the expression of which evolved during the gain, modification, and loss of dimorphic pigmentation among Sophophora fruit flies. We show these CREs to be necessary components of a pigmentation GRN, as deletion from Drosophila melanogaster (derived dimorphic phenotype) resulted in lost expression and lost male-specific pigmentation. We evaluated the ability of orthologous CRE sequences to drive reporter gene expression in species with modified (Drosophila auraria), secondarily lost (Drosophila ananassae), and ancestrally absent (Drosophila willistoni) pigmentation. We show that the transgene host frequently determines CRE activity, implicating trans-evolution as a significant factor for this trait's diversity. We validated the gain of dimorphic Bab transcription factor expression as a trans-change contributing to the dimorphic trait. Our findings suggest an amenability to change for the landscape of trans-regulators and begs for an explanation as to why this is so common compared to the evolution of differentiation gene CREs.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

Increased chromatin accessibility promotes the evolution of a transcriptional silencer in Drosophila

The loss of discrete morphological traits, the most common evolutionary transition, is typically driven by changes in developmental gene expression. Mutations accumulating in regulatory elements of these genes can disrupt DNA binding sites for transcription factors patterning their spatial expression, or delete entire enhancers. Regulatory elements, however, may be silenced through changes in chromatin accessibility or the emergence of repressive elements. Here, we show that increased chromatin accessibility at the gene yellow, combined with the gain of a repressor site, underlies the loss of a wing spot pigmentation pattern in a Drosophila species. The gain of accessibility of this repressive element is regulated by E93, a transcription factor governing the progress of metamorphosis. This convoluted evolutionary scenario contrasts with the parsimonious mutational paths generally envisioned and often documented for morphological losses. It illustrates how evolutionary changes in chromatin accessibility may directly contribute to morphological diversification.

The function and evolution of a genetic switch controlling sexually dimorphic eye differentiation in honeybees

Animals develop sex-specific morphological structures that are diverse between organisms. However, understanding the developmental and evolutionary mechanisms governing these traits is still limited and largely restricted to DM domain genes, which are conserved, sex-specific developmental regulators identified in genetic models. Here, we report a sex-specific developmental regulator gene, glubschauge (glu) that selectively regulates sexually dimorphic eye differentiation in honeybees. We found that the sex determination gene feminizer (fem) controls sex-specific splicing of glu transcripts, establishing a genetic switch in which Glu proteins with a zinc finger (ZnF) domain are only expressed in females. We showed that female coding sequence was essential and sufficient for partial feminization. Comparative sequence and functional studies revealed that the evolutionary origination of the genetic switch was followed by the mutational origin of the essential ZnF domain. Our results demonstrate that glu is a newly evolved sex-specific genetic switch for region-specific regulation of a dimorphic character.

It is time to move on from homology in comparative biology

The continued use of the idea of homology is questionable on philosophical and scientific grounds. It is based on the widespread idea that a "homologue" in extant taxa can be "traced back" to a feature in common ancestor. In contrast, Richard Owen, who first used the term in 1846, saw homology (homologue) differently, as "sameness": "the same organ in different animals under every variety of form and function." At that point in time, he was not influenced by evolutionary thinking, and more focused on the details and approaches to biological comparison and description. His was a perceptive approach to comparison. This paper argues that the concept of homology no longer plays a useful role in comparative biology. It is a conceptual idea with little or no empirical implications for modern comparisons among phenotypes. Comparative biology now uses formal phylogenetic analysis in which similar features in individuals of two or more taxa are treated as characters on a tree and tested for historical "sameness" in terms of the concept of synapomorphy. If we are to understand the complexities of phenotypic evolution, applying this method to detailed comparative data will be essential. At the same time, a deep understanding of the phenotype and its history will emerge only through the use of multidisciplinary approaches that address historical changes at different hierarchical levels.

The regulation of a pigmentation gene in the formation of complex color patterns in Drosophila abdomens

Changes in the control of developmental gene expression patterns have been implicated in the evolution of animal morphology. However, the genetic mechanisms underlying complex morphological traits remain largely unknown. Here we investigated the molecular mechanisms that induce the pigmentation gene yellow in a complex color pattern on the abdomen of Drosophila guttifera. We show that at least five developmental genes may collectively activate one cis-regulatory module of yellow in distinct spot rows and a dark shade to assemble the complete abdominal pigment pattern of Drosophila guttifera. One of these genes, wingless, may play a conserved role in the early phase of spot pattern development in several species of the quinaria group. Our findings shed light on the evolution of complex animal color patterns through modular changes of gene expression patterns.

Can changes in 3D genome architecture create new regulatory landscapes that contribute to phenotypic evolution?

Animal genomes are compartmentalized into insulated regulatory units named topology-associated domains (TADs). TADs insulate gene promoters from enhancers that occupy neighboring TADs. Chromosomal rearrangements that disrupt TAD structure can generate new regulatory interactions between enhancers and promoters that were once separated into different TADs, which might lead to new gene expression patterns. On the one hand, TAD rearrangements are known to cause deleterious phenotypes, but, on the other hand, rearrangements can also create novel expression patterns that may be selected during evolution because they generate advantageous phenotypes. Here, we review recent studies that explore the effects of chromosomal rearrangements and genetic perturbations on TAD structure and gene regulation in the context of development and evolution. We discuss the possible contribution of evolutionary breakpoints (EBRs) that affect TAD structure to the evolution of gene regulation and the phenotype.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Relocation to avoid costs: A hypothesis on red carotenoid-based signals based on recent CYP2J19 gene expression data

In many vertebrates, the enzymatic oxidation of dietary yellow carotenoids generates red keto-carotenoids giving color to ornaments. The oxidase CYP2J19 is here a key effector. Its purported intracellular location suggests a shared biochemical pathway between trait expression and cell functioning. This might guarantee the reliability of red colorations as individual quality signals independent of production costs. We hypothesize that the ornament type (feathers vs. bare parts) and production costs (probably CYP2J19 activity compromising vital functions) could have promoted tissue-specific gene relocation. We review current avian tissue-specific CYP2J19 expression data. Among the ten red-billed species showing CYP2J19 bill expression, only one showed strong hepatic expression. Moreover, a phylogenetically-controlled analysis of 25 red-colored species shows that those producing red bare parts are less likely to have strong hepatic CYP2J19 expression than species with only red plumages. Thus, both production costs and shared pathways might have contributed to the evolution of red signals.

Divergence and correlated evolution of male wing spot and courtship display between Drosophila nepalensis and D. trilutea

Male-specific wing spots are usually associated with wing displays in the courtship behavior of Drosophila and may play important roles in sexual selection. Two closely related species, D. nepalensis and D. trilutea, differ in wing spots and scissoring behavior. Here, we compare male morphological characters, pigmentation intensity of male wing spots, wing-scissoring behavior, courtship songs, and reproductive isolation between 2 species. F1 fertile females and sterile males result from the cross between females of D. nepalensis and males of D. trilutea. The pigmentation of wing spots is significantly weaker in D. trilutea than in D. nepalensis and the F1 hybrid. Males scissor both wings in front of the female during courtship, with a posture spreading wings more widely, and at a faster frequency in D. nepalensis than in D. trilutea and the F1s. Males of D. trilutea vibrate wings to produce 2 types (A and B) of pulse songs, whereas D. nepalensis and the F1s sing only type B songs. The incidence of wing vibration and scissoring during courtship suggests that wing vibration is essential but scissoring is a facultative courtship element for successful mating in both species. The association between the darker wing spots with more elaborate scissoring might be the consequence of correlated evolution of these traits in D. nepalensis; however, D. trilutea retains wing scissoring during courtship despite having weaker pigmentation of wing spots. The genetic architecture of 2 traits differs in the F1s, consistent with maternal or sex-linked effects for spots but nonadditive effects for scissoring.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Biology and ecology of the Oriental flower-breeding Drosophila elegans and related species

Animals adapt to their environments in the course of evolution. One effective approach to elucidate mechanisms of adaptive evolution is to compare closely related species with model organisms in which knowledge of the molecular and physiological bases of various traits has been accumulated. Drosophila elegans and its close relatives, belonging to the same species group as the model organism D. melanogaster, exhibit various unique characteristics such as flower-breeding habit, courtship display, territoriality, sexual dimorphism, and colour polymorphism. Their ease of culturing and availability of genomic information makes them a useful model for understanding mechanisms of adaptive evolution. Here, we review the morphology, distribution, and phylogenetic relationships of D. elegans and related species, as well as their characteristic flower-dependent biology, food habits, and life-history traits. We also describe their unique mating and territorial behaviours and note their distinctive karyotype and the genetic mechanisms of morphological diversity that have recently been revealed.

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.

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.

Evaluation of candidate point mutation of Kisspeptin 1 gene associated with litter size in Indian Goat breeds and its effect on transcription factor binding sites

Kisspeptin gene (Kiss1) has a significant role in reproductive processes in mammals. However, only little information is available about the association of Kiss1 gene with litter size in Indian goat breeds. Thus, blood samples from 285 randomly selected animals were collected for DNA isolation and SNP profiling. The PCR product of 242 bp size harboring g.2540C>T mutation of Kiss1 gene was digested with the restriction enzyme Sac1. Least squares analysis revealed that Barbari goats showed significantly higher average litter size (2.86±0.08) compared to Beetal, Sirohi and Sojat breeds (P < 0.01). SNP locus g.2540C>T of Kiss1 gene also showed significant effect on litter size (P < 0.01). Goats with Genotype CT (2.66 ± 0.07) and TT (2.67 ± 0.26) had significantly higher (P < 0.01) litter size than CC (1.50 ± 0.05). From the transcription factor binding site analysis, it was predicted that due to g.2540C>T SNP, both native and mutant variant forms coded for putative binding sites for different transcription Factor. Allele T had putative binding sites for the androgen receptor which plays a significant role in the signaling pathway involved in increase in ovulation rate; which consequently can have a tremendous effect on average litter size.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Fossil dragonfly-type larva with lateral abdominal protrusions and implications on the early evolution of Pterygota

Aquatic larvae are known in three early branches of Pterygota: Ephemeroptera (mayflies), Plecoptera (stoneflies), and Odonata (dragonflies, damselflies). A common origin of these larvae has been suggested, yet also counterarguments have been put forward, for example, the different position of larval gills: laterally on the abdomen in Ephemeroptera, terminally in Odonata, variably in Plecoptera. We discuss recent fossil findings and report a new dragonfly-type larva from Kachin amber (Myanmar), which possesses ancestral characters such as a terminal filum, maintained in ephemeropterans, but lost in modern odonatan larvae. The new larva possesses lateral protrusions on the abdominal segments where in other lineages gills occur. Together with other fossils, such as a plecopteran retaining lateral gills on the abdomen, this indicates that lateral protrusions on the abdomen might have well been an ancestral feature, removing one important argument against the idea of an aquatic larva in the ground pattern of Pterygota.

•

A new dragonfly-type larva was found in Kachin amber (Myanmar, 99 million years old).

•

The larva possesses a terminal filum, which is not known in modern dragonfly larvae

•

It also exhibits lateral abdominal protrusions where in other lineages gills occur

•

This find makes an aquatic larva in the ground pattern of Pterygota more likely

Local tissue interactions govern pLL patterning in medaka

Vertebrate organs are arranged in a stereotypic, species-specific position along the animal body plan. Substantial morphological variation exists between related species, especially so in the vastly diversified teleost clade. It is still unclear how tissues, organs and systems can accommodate such diverse scaffolds. Here, we use the distinctive arrangement of neuromasts in the posterior lateral line (pLL) system of medaka fish to address the tissue-interactions defining a pattern. We show that patterning in this peripheral nervous system is established by autonomous organ precursors independent of neuronal wiring. In addition, we target the keratin 15 gene to generate stuck-in-the-midline (siml) mutants, which display epithelial lesions and a disrupted pLL patterning. By using siml/wt chimeras, we determine that the aberrant siml pLL pattern depends on the mutant epithelium, since a wild type epithelium can rescue the siml phenotype. Inducing epithelial lesions by 2-photon laser ablation during pLL morphogenesis phenocopies siml genetic mutants and reveals that epithelial integrity defines the final position of the embryonic pLL neuromasts. Our results using the medaka pLL disentangle intrinsic from extrinsic properties during the establishment of a sensory system. We speculate that intrinsic programs guarantee proper organ morphogenesis, while instructive interactions from surrounding tissues facilitates the accommodation of sensory organs to the diverse body plans found among teleosts.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Modeling the evolutionary architectures of transcribed human enhancer sequences reveals distinct origins, functions, and associations with human-trait variation

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits.

To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers.

Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.

An Optimized Transformation System and Functional Test of CYC-Like TCP Gene CpCYC in Chirita pumila (Gesneriaceae)

The development of an ideal model plant located at a key phylogenetic node is critically important to advance functional and regulatory studies of key regulatory genes in the evolutionary developmental (evo-devo) biology field. In this study, we selected Chirita pumila in the family Gesneriaceae, a basal group in Lamiales, as a model plant to optimize its genetic transformation system established previously by us through investigating a series of factors and further conduct functional test of the CYC-like floral symmetry gene CpCYC. By transforming a RNAi:CpCYC vector, we successfully achieved the desired phenotypes of upright actinomorphic flowers, which suggest that CpCYC actually determines the establishment of floral zygomorphy and the horizontal orientation of flowers in C. pumila. We also confirmed the activities of CpCYC promoter in dorsal petals, dorsal/lateral staminodes, as well as the pedicel by transferring a CpCYC promoter:GUS vector into C. pumila. Furthermore, we testified the availability of a transient gene expression system using C. pumila mesophyll protoplasts. The improved transformation system together with the inherent biological features would make C. pumila an attractive new model in functional and regulatory studies for a broad range of evo-devo issues.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Decoding the Reproductive System of the Olive Fruit Fly, Bactrocera oleae

In most diploid organisms, mating is a prerequisite for reproduction and, thus, critical to the maintenance of their population and the perpetuation of the species. Besides the importance of understanding the fundamentals of reproduction, targeting the reproductive success of a pest insect is also a promising method for its control, as a possible manipulation of the reproductive system could affect its destructive activity. Here, we used an integrated approach for the elucidation of the reproductive system and mating procedures of the olive fruit fly, Bactrocera oleae. Initially, we performed a RNAseq analysis in reproductive tissues of virgin and mated insects. A comparison of the transcriptomes resulted in the identification of genes that are differentially expressed after mating. Functional annotation of the genes showed an alteration in the metabolic, catalytic, and cellular processes after mating. Moreover, a functional analysis through RNAi silencing of two differentially expressed genes, yellow-g and troponin C, resulted in a significantly reduced oviposition rate. This study provided a foundation for future investigations into the olive fruit fly's reproductive biology to the development of new exploitable tools for its control.

Remarkable sexually dimorphic features ofConiceromyia(Diptera: Phoridae): evolution in the light of phylogeny and comparative evidence about their function

Coniceromyia is a genus of 100 species of phorid flies mostly distributed in the Neotropical region. The genus is distinguishable based on several male-exclusive features in different parts of the body, many of which are unique among the Phoridae. In addition, many species of Coniceromyia have nearly identical morphology of their male copulatory apparatus (i.e. hypopygium). The co-occurrence of these unusual characteristics suggests an evolutionary correlation between them. To investigate this possible correlation and to understand other aspects of the evolution of these puzzling male-exclusive characters, we performed the first phylogenetic analysis of Coniceromyia, based on morphological and molecular data. Ancestral state reconstructions and comparative analyses then allowed us to infer the evolution of these characters and search for general evolutionary patterns and correlated histories. We demonstrate that these male-exclusive features varied from highly homoplastic to uniquely derived on the phylogenetic history of Coniceromyia. For some characters, we found evidence of a biased evolution favouring gains over losses of the feature, but no male characteristics were significantly correlated with hypopygium morphology. The evolutionary patterns of the male-exclusive features and comparative evidence with other better known groups suggest possible functions for these features related to sexual selection.

Fitness maps to a large-effect locus in introduced stickleback populations

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Individual differences in honey bee behavior enabled by plasticity in brain gene regulatory networks

Understanding the regulatory architecture of phenotypic variation is a fundamental goal in biology, but connections between gene regulatory network (GRN) activity and individual differences in behavior are poorly understood. We characterized the molecular basis of behavioral plasticity in queenless honey bee (Apis mellifera) colonies, where individuals engage in both reproductive and non-reproductive behaviors. Using high-throughput behavioral tracking, we discovered these colonies contain a continuum of phenotypes, with some individuals specialized for either egg-laying or foraging and 'generalists' that perform both. Brain gene expression and chromatin accessibility profiles were correlated with behavioral variation, with generalists intermediate in behavior and molecular profiles. Models of brain GRNs constructed for individuals revealed that transcription factor (TF) activity was highly predictive of behavior, and behavior-associated regulatory regions had more TF motifs. These results provide new insights into the important role played by brain GRN plasticity in the regulation of behavior, with implications for social evolution.

Evolved Differences in cis and trans Regulation Between the Maternal and Zygotic mRNA Complements in the Drosophila Embryo

How gene expression can evolve depends on the mechanisms driving gene expression. Gene expression is controlled in different ways in different developmental stages; here we ask whether different developmental stages show different patterns of regulatory evolution. To explore the mode of regulatory evolution, we used the early stages of embryonic development controlled by two different genomes, that of the mother and that of the zygote. During embryogenesis in all animals, initial developmental processes are driven entirely by maternally provided gene products deposited into the oocyte. The zygotic genome is activated later, when developmental control is handed off from maternal gene products to the zygote during the maternal-to-zygotic transition. Using hybrid crosses between sister species of Drosophila (Dsimulans, D. sechellia, and D. mauritiana) and transcriptomics, we find that the regulation of maternal transcript deposition and zygotic transcription evolve through different mechanisms. We find that patterns of transcript level inheritance in hybrids, relative to parental species, differ between maternal and zygotic transcripts, and maternal transcript levels are more likely to be conserved. Changes in transcript levels occur predominantly through differences in trans regulation for maternal genes, while changes in zygotic transcription occur through a combination of both cis and trans regulatory changes. Differences in the underlying regulatory landscape in the mother and the zygote are likely the primary determinants for how maternal and zygotic transcripts evolve.

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.

Enhancer evolutionary co-option through shared chromatin accessibility input

Form diversity is fueled by changes in the expression of genes that build organisms. New expression often results from the emergence of new DNA switches, known as transcriptional enhancers. Many enhancers are thought to appear through the recycling of older enhancers, a process called evolutionary co-option. Enhancer co-option is difficult to assess, and the molecular mechanisms explaining its prevalence are elusive. Using state-of-the-art quantification and analyses, we reveal that the sequences of an ancestral and a derived enhancer overlap extensively. They contain specific binding sites for regulators imparting spatial activities. We found that the two enhancers also share a site facilitating access to chromatin in a region where they overlap.

The diversity of forms in multicellular organisms originates largely from the spatial redeployment of developmental genes [S. B. Carroll, Cell 134, 25–36 (2008)]. Several scenarios can explain the emergence of cis-regulatory elements that govern novel aspects of a gene expression pattern [M. Rebeiz, M. Tsiantis, Curr. Opin. Genet. Dev. 45, 115–123 (2017)]. One scenario, enhancer co-option, holds that a DNA sequence producing an ancestral regulatory activity also becomes the template for a new regulatory activity, sharing regulatory information. While enhancer co-option might fuel morphological diversification, it has rarely been documented [W. J. Glassford et al., Dev. Cell 34, 520–531 (2015)]. Moreover, if two regulatory activities are borne from the same sequence, their modularity, considered a defining feature of enhancers [J. Banerji, L. Olson, W. Schaffner, Cell 33, 729–740 (1983)], might be affected by pleiotropy. Sequence overlap may thereby play a determinant role in enhancer function and evolution. Here, we investigated this problem with two regulatory activities of the Drosophila gene yellow, the novel spot enhancer and the ancestral wing blade enhancer. We used precise and comprehensive quantification of each activity in Drosophila wings to systematically map their sequences along the locus. We show that the spot enhancer has co-opted the sequences of the wing blade enhancer. We also identified a pleiotropic site necessary for DNA accessibility of a shared regulatory region. While the evolutionary steps leading to the derived activity are still unknown, such pleiotropy suggests that enhancer accessibility could be one of the molecular mechanisms seeding evolutionary co-option.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.