The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

Quantification provides a conceptual basis for convergent evolution

While much of evolutionary biology attempts to explain the processes of diversification, there is an important place for the study of phenotypic similarity across life forms. When similar phenotypes evolve independently in different lineages this is referred to as convergent evolution. Although long recognised, evolutionary convergence is receiving a resurgence of interest. This is in part because new genomic data sets allow detailed and tractable analysis of the genetic underpinnings of convergent phenotypes, and in part because of renewed recognition that convergence may reflect limitations in the diversification of life. In this review we propose that although convergent evolution itself does not require a new evolutionary framework, none the less there is room to generate a more systematic approach which will enable evaluation of the importance of convergent phenotypes in limiting the diversity of life's forms. We therefore propose that quantification of the frequency and strength of convergence, rather than simply identifying cases of convergence, should be considered central to its systematic comprehension. We provide a non-technical review of existing methods that could be used to measure evolutionary convergence, bringing together a wide range of methods. We then argue that quantification also requires clear specification of the level at which the phenotype is being considered, and argue that the most constrained examples of convergence show similarity both in function and in several layers of underlying form. Finally, we argue that the most important and impressive examples of convergence are those that pertain, in form and function, across a wide diversity of selective contexts as these persist in the likely presence of different selection pressures within the environment.

Convergent evolution of SWS2 opsin facilitates adaptive radiation of threespine stickleback into different light environments

Repeated adaptation to a new environment often leads to convergent phenotypic changes whose underlying genetic mechanisms are rarely known. Here, we study adaptation of color vision in threespine stickleback during the repeated postglacial colonization of clearwater and blackwater lakes in the Haida Gwaii archipelago. We use whole genomes from 16 clearwater and 12 blackwater populations, and a selection experiment, in which stickleback were transplanted from a blackwater lake into an uninhabited clearwater pond and resampled after 19 y to test for selection on cone opsin genes. Patterns of haplotype homozygosity, genetic diversity, site frequency spectra, and allele-frequency change support a selective sweep centered on the adjacent blue- and red-light sensitive opsins SWS2 and LWS. The haplotype under selection carries seven amino acid changes in SWS2, including two changes known to cause a red-shift in light absorption, and is favored in blackwater lakes but disfavored in the clearwater habitat of the transplant population. Remarkably, the same red-shifting amino acid changes occurred after the duplication of SWS2 198 million years ago, in the ancestor of most spiny-rayed fish. Two distantly related fish species, bluefin killifish and black bream, express these old paralogs divergently in black- and clearwater habitats, while sticklebacks lost one paralog. Our study thus shows that convergent adaptation to the same environment can involve the same genetic changes on very different evolutionary time scales by reevolving lost mutations and reusing them repeatedly from standing genetic variation.

Genomic Scans across Three Eucalypts Suggest that Adaptation to Aridity is a Genome-Wide Phenomenon

Widespread species spanning strong environmental (e.g., climatic) gradients frequently display morphological and physiological adaptations to local conditions. Some adaptations are common to different species that occupy similar environments. However, the genomic architecture underlying such convergent traits may not be the same between species. Using genomic data from previous studies of three widespread eucalypt species that grow along rainfall gradients in southern Australia, our probabilistic approach provides evidence that adaptation to aridity is a genome-wide phenomenon, likely to involve multiple and diverse genes, gene families and regulatory regions that affect a multitude of complex genetic and biochemical processes.

The ultimate and proximate mechanisms driving the evolution of long tails in forest deer mice

Understanding both the role of selection in driving phenotypic change and its underlying genetic basis remain major challenges in evolutionary biology. Here, we use modern tools to revisit a classic system of local adaptation in the North American deer mouse, Peromyscus maniculatus, which occupies two main habitat types: prairie and forest. Using historical collections, we find that forest‐dwelling mice have longer tails than those from nonforested habitat, even when we account for individual and population relatedness. Using genome‐wide SNP data, we show that mice from forested habitats in the eastern and western parts of their range form separate clades, suggesting that increased tail length evolved independently. We find that forest mice in the east and west have both more and longer caudal vertebrae, but not trunk vertebrae, than nearby prairie forms. By intercrossing prairie and forest mice, we show that the number and length of caudal vertebrae are not correlated in this recombinant population, indicating that variation in these traits is controlled by separate genetic loci. Together, these results demonstrate convergent evolution of the long‐tailed forest phenotype through two distinct genetic mechanisms, affecting number and length of vertebrae, and suggest that these morphological changes—either independently or together—are adaptive.

Genomic Trajectories to Desiccation Resistance: Convergence and Divergence Among Replicate Selected Drosophila Lines

Adaptation to environmental stress is critical for long-term species persistence. With climate change and other anthropogenic stressors compounding natural selective pressures, understanding the nature of adaptation is as important as ever in evolutionary biology. In particular, the number of alternative molecular trajectories available for an organism to reach the same adaptive phenotype remains poorly understood. Here, we investigate this issue in a set of replicated Drosophila melanogaster lines selected for increased desiccation resistance-a classical physiological trait that has been closely linked to Drosophila species distributions. We used pooled whole-genome sequencing (Pool-Seq) to compare the genetic basis of their selection responses, using a matching set of replicated control lines for characterizing laboratory (lab-)adaptation, as well as the original base population. The ratio of effective population size to census size was high over the 21 generations of the experiment at 0.52-0.88 for all selected and control lines. While selected SNPs in replicates of the same treatment (desiccation-selection or lab-adaptation) tended to change frequency in the same direction, suggesting some commonality in the selection response, candidate SNP and gene lists often differed among replicates. Three of the five desiccation-selection replicates showed significant overlap at the gene and network level. All five replicates showed enrichment for ovary-expressed genes, suggesting maternal effects on the selected trait. Divergence between pairs of replicate lines for desiccation-candidate SNPs was greater than between pairs of control lines. This difference also far exceeded the divergence between pairs of replicate lines for neutral SNPs. Overall, while there was overlap in the direction of allele frequency changes and the network and functional categories affected by desiccation selection, replicates showed unique responses at all levels, likely reflecting hitchhiking effects, and highlighting the challenges in identifying candidate genes from these types of experiments when traits are likely to be polygenic.

Evolutionary convergence in experimental Pseudomonas populations

Model microbial systems provide opportunity to understand the genetic bases of ecological traits, their evolution, regulation and fitness contributions. Experimental populations of Pseudomonas fluorescens rapidly diverge in spatially structured microcosms producing a range of surface-colonising forms. Despite divergent molecular routes, wrinkly spreader (WS) niche specialist types overproduce a cellulosic polymer allowing mat formation at the air-liquid interface and access to oxygen. Given the range of ways by which cells can form mats, such phenotypic parallelism is unexpected. We deleted the cellulose-encoding genes from the ancestral genotype and asked whether this mutant could converge on an alternate phenotypic solution. Two new traits were discovered. The first involved an exopolysaccharide encoded by pgaABCD that functions as cell-cell glue similar to cellulose. The second involved an activator of an amidase (nlpD) that when defective causes cell chaining. Both types form mats, but were less fit in competition with cellulose-based WS types. Surprisingly, diguanylate cyclases linked to cellulose overexpression underpinned evolution of poly-beta-1,6-N-acetyl-d-glucosamine (PGA)-based mats. This prompted genetic analyses of the relationships between the diguanylate cyclases WspR, AwsR and MwsR, and both cellulose and PGA. Our results suggest that c-di-GMP regulatory networks may have been shaped by evolution to accommodate loss and gain of exopolysaccharide modules facilitating adaptation to new environments.

Pigeonetics takes flight: Evolution, development, and genetics of intraspecific variation

Intensive artificial selection over thousands of years has produced hundreds of varieties of domestic pigeon. As Charles Darwin observed, the morphological differences among breeds can rise to the magnitude of variation typically observed among different species. Nevertheless, different pigeon varieties are interfertile, thereby enabling forward genetic and genomic approaches to identify genes that underlie derived traits. Building on classical genetic studies of pigeon variation, recent molecular investigations find a spectrum of coding and regulatory alleles controlling derived traits, including plumage color, feather growth polarity, and limb identity. Developmental and genetic analyses of pigeons are revealing the molecular basis of variation in a classic example of extreme intraspecific diversity, and have the potential to nominate genes that control variation among other birds and vertebrates in general.

Potential Direct Regulators of the Drosophila yellow Gene Identified by Yeast One-Hybrid and RNAi Screens

The regulation of gene expression controls development, and changes in this regulation often contribute to phenotypic evolution. Drosophila pigmentation is a model system for studying evolutionary changes in gene regulation, with differences in expression of pigmentation genes such as yellow that correlate with divergent pigment patterns among species shown to be caused by changes in cis- and trans-regulation. Currently, much more is known about the cis-regulatory component of divergent yellow expression than the trans-regulatory component, in part because very few trans-acting regulators of yellow expression have been identified. This study aims to improve our understanding of the trans-acting control of yellow expression by combining yeast-one-hybrid and RNAi screens for transcription factors binding to yellow cis-regulatory sequences and affecting abdominal pigmentation in adults, respectively. Of the 670 transcription factors included in the yeast-one-hybrid screen, 45 showed evidence of binding to one or more sequence fragments tested from the 5′ intergenic and intronic yellow sequences from D. melanogaster, D. pseudoobscura, and D. willistoni, suggesting that they might be direct regulators of yellow expression. Of the 670 transcription factors included in the yeast-one-hybrid screen, plus another TF previously shown to be genetically upstream of yellow, 125 were also tested using RNAi, and 32 showed altered abdominal pigmentation. Nine transcription factors were identified in both screens, including four nuclear receptors related to ecdysone signaling (Hr78, Hr38, Hr46, and Eip78C). This finding suggests that yellow expression might be directly controlled by nuclear receptors influenced by ecdysone during early pupal development when adult pigmentation is forming.

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.

Dynamic Convergent Evolution Drives the Passage Adaptation across 48 Years' History of H3N2 Influenza Evolution

Influenza viruses are often propagated in a diverse set of culturing media and additional substitutions known as passage adaptation can cause extra evolution in the target strain, leading to ineffective vaccines. Using 25,482 H3N2 HA1 sequences curated from Global Initiative on Sharing All Influenza Data and National Center for Biotechnology Information databases, we found that passage adaptation is a very dynamic process that changes over time and evolves in a seesaw like pattern. After crossing the species boundary from bird to human in 1968, the influenza H3N2 virus evolves to be better adapted to the human environment and passaging them in embryonated eggs (i.e., an avian environment) leads to increasingly stronger positive selection. On the contrary, passage adaptation to the mammalian cell lines changes from positive selection to negative selection. Using two statistical tests, we identified 19 codon positions around the receptor binding domain strongly contributing to passage adaptation in the embryonated egg. These sites show strong convergent evolution and overlap extensively with positively selected sites identified in humans, suggesting that passage adaptation can confound many of the earlier studies on influenza evolution. Interestingly, passage adaptation in recent years seems to target a few codon positions in antigenic surface epitopes, which makes it difficult to produce antigenically unaltered vaccines using embryonic eggs. Our study outlines another interesting scenario whereby both convergent and adaptive evolution are working in synchrony driving viral adaptation. Future studies from sequence analysis to vaccine production need to take careful consideration of passage adaptation.

Phenotypic Plasticity through Transcriptional Regulation of the Evolutionary Hotspot Gene tan in Drosophila melanogaster

Phenotypic plasticity is the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. Phenotypic plasticity can be adaptive. Furthermore, it is thought to facilitate evolution. Although phenotypic plasticity is a widespread phenomenon, its molecular mechanisms are only beginning to be unravelled. Environmental conditions can affect gene expression through modification of chromatin structure, mainly via histone modifications, nucleosome remodelling or DNA methylation, suggesting that phenotypic plasticity might partly be due to chromatin plasticity. As a model of phenotypic plasticity, we study abdominal pigmentation of Drosophila melanogaster females, which is temperature sensitive. Abdominal pigmentation is indeed darker in females grown at 18°C than at 29°C. This phenomenon is thought to be adaptive as the dark pigmentation produced at lower temperature increases body temperature. We show here that temperature modulates the expression of tan (t), a pigmentation gene involved in melanin production. t is expressed 7 times more at 18°C than at 29°C in female abdominal epidermis. Genetic experiments show that modulation of t expression by temperature is essential for female abdominal pigmentation plasticity. Temperature modulates the activity of an enhancer of t without modifying compaction of its chromatin or level of the active histone mark H3K27ac. By contrast, the active mark H3K4me3 on the t promoter is strongly modulated by temperature. The H3K4 methyl-transferase involved in this process is likely Trithorax, as we show that it regulates t expression and the H3K4me3 level on the t promoter and also participates in female pigmentation and its plasticity. Interestingly, t was previously shown to be involved in inter-individual variation of female abdominal pigmentation in Drosophila melanogaster, and in abdominal pigmentation divergence between Drosophila species. Sensitivity of t expression to environmental conditions might therefore give more substrate for selection, explaining why this gene has frequently been involved in evolution of pigmentation.

Environmental conditions can strongly modulate the phenotype produced by a particular genotype. This process, called phenotypic plasticity, has major implications in medicine and agricultural sciences, and is thought to facilitate evolution. Phenotypic plasticity is observed in many animals and plants but its mechanisms are only partially understood. As a model of phenotypic plasticity, we study the effect of temperature on female abdominal pigmentation in the fruit fly Drosophila melanogaster. Here we show that temperature affects female abdominal pigmentation by modulating the expression of tan (t), a gene involved in melanin production, in female abdominal epidermis. This effect is mediated at least partly by a particular regulatory sequence of t, the t_MSE enhancer. However we detected no modulation of chromatin structure of t_MSE by temperature. By contrast, the level of the active chromatin mark H3K4me3 on the t promoter is strongly increased at lower temperature. We show that the H3K4 methyl-transferase Trithorax is involved in female abdominal pigmentation and its plasticity and regulates t expression and H3K4me3 level on the t promoter. Several studies have linked t to pigmentation evolution within and between Drosophila species. Our results suggest that sensitivity of t expression to temperature might facilitate its role in pigmentation evolution.

Ectodysplasin signalling genes and phenotypic evolution in sculpins (Cottus)

Despite their deeply conserved function among vertebrates, ectodysplasin (Eda) signalling genes are involved in microevolutionary change in humans and sticklebacks. If such a dual role is common, Eda signalling genes constitute hotspots for morphological evolution. Variation in sculpin (Cottus) skin prickling and body shape resembles patterns caused by variation in Eda signalling in sticklebacks. We mapped Eda signalling genes and performed quantitative trait locus mapping in crosses between Cottus rhenanus and Cottus perifretum. A genomic region containing the Eda receptor (Edar) was strongly associated with prickling and contributed to shape. The expression of Edar in developing prickles and skeletal elements in Cottus was confirmed by in situ hybridization. Coding sequence changes between Edar alleles in C. rhenanus and C. perifretum exceeded sequence differentiation in other vertebrates. However, it is likely that additional genetic elements besides coding changes affect the phenotypic variation. Although the phenotype in a natural hybrid lineage between C. rhenanus and C. perifretum resembles C. perifretum, the respective coding Edar alleles are not fully fixed (88.6%). Hence, our results support an involvement of Eda signalling in microevolutionary changes, but imply that the Edar gene is affected by multiple evolutionary processes that vary among freshwater sculpins.

What, if anything, are hybrids: enduring truths and challenges associated with population structure and gene flow

Hybridization is a potent evolutionary process that can affect the origin, maintenance, and loss of biodiversity. Because of its ecological and evolutionary consequences, an understanding of hybridization is important for basic and applied sciences, including conservation biology and agriculture. Herein, we review and discuss ideas that are relevant to the recognition of hybrids and hybridization. We supplement this discussion with simulations. The ideas we present have a long history, particularly in botany, and clarifying them should have practical consequences for managing hybridization and gene flow in plants. One of our primary goals is to illustrate what we can and cannot infer about hybrids and hybridization from molecular data; in other words, we ask when genetic analyses commonly used to study hybridization might mislead us about the history or nature of gene flow and selection. We focus on patterns of variation when hybridization is recent and populations are polymorphic, which are particularly informative for applied issues, such as contemporary hybridization following recent ecological change. We show that hybridization is not a singular process, but instead a collection of related processes with variable outcomes and consequences. Thus, it will often be inappropriate to generalize about the threats or benefits of hybridization from individual studies, and at minimum, it will be important to avoid categorical thinking about what hybridization and hybrids are. We recommend potential sampling and analytical approaches that should help us confront these complexities of hybridization.

The industrial melanism mutation in British peppered moths is a transposable element

Discovering the mutational events that fuel adaptation to environmental change remains an important challenge for evolutionary biology. The classroom example of a visible evolutionary response is industrial melanism in the peppered moth (Biston betularia): the replacement, during the Industrial Revolution, of the common pale typica form by a previously unknown black (carbonaria) form, driven by the interaction between bird predation and coal pollution. The carbonaria locus has been coarsely localized to a 200-kilobase region, but the specific identity and nature of the sequence difference controlling the carbonaria-typica polymorphism, and the gene it influences, are unknown. Here we show that the mutation event giving rise to industrial melanism in Britain was the insertion of a large, tandemly repeated, transposable element into the first intron of the gene cortex. Statistical inference based on the distribution of recombined carbonaria haplotypes indicates that this transposition event occurred around 1819, consistent with the historical record. We have begun to dissect the mode of action of the carbonaria transposable element by showing that it increases the abundance of a cortex transcript, the protein product of which plays an important role in cell-cycle regulation, during early wing disc development. Our findings fill a substantial knowledge gap in the iconic example of microevolutionary change, adding a further layer of insight into the mechanism of adaptation in response to natural selection. The discovery that the mutation itself is a transposable element will stimulate further debate about the importance of 'jumping genes' as a source of major phenotypic novelty.

Characterizing the adult and larval transcriptome of the multicolored Asian lady beetle, Harmonia axyridis

The reasons for the evolution and maintenance of striking visual phenotypes are as widespread as the species that display these phenotypes. While study systems such as Heliconius and Dendrobatidae have been well characterized and provide critical information about the evolution of these traits, a breadth of new study systems, in which the phenotype of interest can be easily manipulated and quantified, are essential for gaining a more general understanding of these specific evolutionary processes. One such model is the multicolored Asian lady beetle, Harmonia axyridis, which displays significant elytral spot and color polymorphism. Using transcriptome data from two life stages, adult and larva, we characterize the transcriptome, thereby laying a foundation for further analysis and identification of the genes responsible for the continual maintenance of spot variation in H. axyridis.

Craniofacial shape transition across the house mouse hybrid zone: implications for the genetic architecture and evolution of between-species differences

Craniofacial shape differences between taxa have often been linked to environmental adaptation, e.g., new food sources, or have been studied in the context of domestication. Evidence for the genetic basis of such phenotypic differences to date suggests that between-species as well as between-population variation has an oligogenic basis, i.e., few loci of large effect explain most of the variation. In mice, it has been shown that within-population craniofacial variation has a highly polygenic basis, but there are no data regarding the genetic basis of between-species differences in natural populations. Here, we address this question using a phenotype-focused approach. Using 3D geometric morphometrics, we phenotyped a panel of mice derived from a natural hybrid zone between Mus musculus domesticus and Mus mus musculus and quantify the transition of craniofacial shape along the hybridization gradient. We find a continuous shape transition along the hybridization gradient and unaltered developmental stability associated with hybridization. This suggests that the morphospace between the two subspecies is continuous despite reproductive isolation and strong barriers to gene flow. We show that quantitative changes in overall genome composition generate quantitative changes in craniofacial shape; this supports a highly polygenic basis for between-species craniofacial differences in the house mouse. We discuss our findings in the context of oligogenic versus polygenic models of the genetic architecture of morphological traits.

Molecular evolution of anthocyanin pigmentation genes following losses of flower color

Background

Phenotypic transitions, such as trait gain or loss, are predicted to carry evolutionary consequences for the genes that control their development. For example, trait losses can result in molecular decay of the pathways underlying the trait. Focusing on the Iochrominae clade (Solanaceae), we examine how repeated losses of floral anthocyanin pigmentation associated with flower color transitions have affected the molecular evolution of three anthocyanin pathway genes (Chi, F3h, and Dfr).

Results

We recovered intact coding regions for the three genes in all of the lineages that have lost floral pigmentation, suggesting that molecular decay is not associated with these flower color transitions. However, two of the three genes (Chi, F3h) show significantly elevated dN/dS ratios in lineages without floral pigmentation. Maximum likelihood analyses suggest that this increase is due to relaxed constraint on anthocyanin genes in the unpigmented lineages as opposed to positive selection. Despite the increase, the values for dN/dS in both pigmented and unpigmented lineages were consistent overall with purifying selection acting on these loci.

Conclusions

The broad conservation of anthocyanin pathway genes across lineages with and without floral anthocyanins is consistent with the growing consensus that losses of pigmentation are largely achieved by changes in gene expression as opposed to structural mutations. Moreover, this conservation maintains the potential for regain of flower color, and indicates that evolutionary losses of floral pigmentation may be readily reversible.

Electronic supplementary material

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

Evolution of Schooling Behavior in Threespine Sticklebacks Is Shaped by the Eda Gene

Despite longstanding interest in the genetic mechanisms that underlie behavioral evolution, very few genes that underlie naturally occurring variation in behavior between individuals or species are known, particularly in vertebrates. Here, we build on our previous forward genetic mapping experiments and use transgenic approaches to identify Ectodysplasin as a gene that causes differences in schooling behavior between wild populations of threespine stickleback (Gasterosteus aculeatus) fish. This work provides rare insight into the proximate mechanisms that have shaped the evolution of vertebrate behavior.

Recurrent specialization on a toxic fruit in an island Drosophila population

Recurrent specialization on similar host plants offers a unique opportunity to unravel the evolutionary and genetic mechanisms underlying dietary shifts. Recent studies have focused on ecological races belonging to the same species, but it is hard in many cases to untangle the role of adaptive introgression versus distinct mutations in facilitating recurrent evolution. We discovered on the island of Mayotte a population of the generalist fly Drosophila yakuba that is strictly associated with noni (Morinda citrifolia). This case strongly resembles Drosophila sechellia, a genetically isolated insular relative of D. yakuba whose intensely studied specialization on toxic noni fruits has always been considered a unique event in insect evolution. Experiments revealed that unlike mainland D. yakuba strains, Mayotte flies showed strong olfactory attraction and significant toxin tolerance to noni. Island females strongly discriminated against mainland males, suggesting that dietary adaptation has been accompanied by partial reproductive isolation. Population genomic analysis indicated a recent colonization (∼29 kya), at a time when year-round noni fruits may have presented a predictable resource on the small island, with ongoing migration after colonization. This relatively recent time scale allowed us to search for putatively adaptive loci based on genetic variation. Strong signals of genetic differentiation were found for several detoxification genes, including a major toxin tolerance locus in D. sechellia Our results suggest that recurrent evolution on a toxic resource can involve similar historical events and common genetic bases, and they establish an important genetic system for the study of early stages of ecological specialization and speciation.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

Competition between anthocyanin and flavonol biosynthesis produces spatial pattern variation of floral pigments between Mimulus species

Flower color patterns have long served as a model for developmental genetics because pigment phenotypes are visually striking, yet generally not required for plant viability, facilitating the genetic analysis of color and pattern mutants. The evolution of novel flower colors and patterns has played a key role in the adaptive radiation of flowering plants via their specialized interactions with different pollinator guilds (e.g., bees, butterflies, birds), motivating the search for allelic differences affecting flower color pattern in closely related plant species with different pollinators. We have identified LIGHT AREAS1 (LAR1), encoding an R2R3-MYB transcription factor, as the causal gene underlying the spatial pattern variation of floral anthocyanin pigmentation between two sister species of monkeyflower: the bumblebee-pollinated Mimulus lewisii and the hummingbird-pollinated Mimulus cardinalis. We demonstrated that LAR1 positively regulates FLAVONOL SYNTHASE (FLS), essentially eliminating anthocyanin biosynthesis in the white region (i.e., light areas) around the corolla throat of M. lewisii flowers by diverting dihydroflavonol into flavonol biosynthesis from the anthocyanin pigment pathway. FLS is preferentially expressed in the light areas of the M. lewisii flower, thus prepatterning the corolla. LAR1 expression in M. cardinalis flowers is much lower than in M. lewisii, explaining the unpatterned phenotype and recessive inheritance of the M. cardinalis allele. Furthermore, our gene-expression analysis and genetic mapping results suggest that cis-regulatory change at the LAR1 gene played a critical role in the evolution of different pigmentation patterns between the two species.

Divergent MLS1 Promoters Lie on a Fitness Plateau for Gene Expression

Qualitative patterns of gene activation and repression are often conserved despite an abundance of quantitative variation in expression levels within and between species. A major challenge to interpreting patterns of expression divergence is knowing which changes in gene expression affect fitness. To characterize the fitness effects of gene expression divergence, we placed orthologous promoters from eight yeast species upstream of malate synthase (MLS1) in Saccharomyces cerevisiae. As expected, we found these promoters varied in their expression level under activated and repressed conditions as well as in their dynamic response following loss of glucose repression. Despite these differences, only a single promoter driving near basal levels of expression caused a detectable loss of fitness. We conclude that the MLS1 promoter lies on a fitness plateau whereby even large changes in gene expression can be tolerated without a substantial loss of fitness.

Can the experimental evolution programme help us elucidate the genetic basis of adaptation in nature?

There have been a variety of approaches taken to try to characterize and identify the genetic basis of adaptation in nature, spanning theoretical models, experimental evolution studies and direct tests of natural populations. Theoretical models can provide formalized and detailed hypotheses regarding evolutionary processes and patterns, from which experimental evolution studies can then provide important proofs of concepts and characterize what is biologically reasonable. Genetic and genomic data from natural populations then allow for the identification of the particular factors that have and continue to play an important role in shaping adaptive evolution in the natural world. Further to this, experimental evolution studies allow for tests of theories that may be difficult or impossible to test in natural populations for logistical and methodological reasons and can even generate new insights, suggesting further refinement of existing theories. However, as experimental evolution studies often take place in a very particular set of controlled conditions--that is simple environments, a small range of usually asexual species, relatively short timescales--the question remains as to how applicable these experimental results are to natural populations. In this review, we discuss important insights coming from experimental evolution, focusing on four key topics tied to the evolutionary genetics of adaptation, and within those topics, we discuss the extent to which the experimental work compliments and informs natural population studies. We finish by making suggestions for future work in particular a need for natural population genomic time series data, as well as the necessity for studies that combine both experimental evolution and natural population approaches.

Widespread flower color convergence in Solanaceae via alternate biochemical pathways

Phenotypic convergence is rampant throughout the tree of life. While recent studies have made significant progress in ascertaining the proximate mechanisms underlying convergent phenotypes, less is known about the frequency and predictability with which convergent phenotypes arise via the same or multiple pathways at the macroevolutionary scale. We investigated the proximate causes and evolutionary patterns of red flower color in the tomato family, Solanaceae, using large-scale data mining and new sequence data to reconstruct a megaphylogeny of 1341 species. We then combined spectral and anatomical data to assess how many times red flowers have evolved, the relative contribution of different pathways to independent origins of red, and whether the underlying pathway is predicted by phylogenetic relatedness. We estimated at least 30 relatively recent origins of red flowers using anthocyanins, carotenoids, or a dual production of both pigments, with significant phylogenetic signal in the use of anthocyanins and dual production, indicating that closely related red-flowered species tend to employ the same mechanism for coloration. Our study is the first to test whether developmental pathways exhibit phylogenetic signal and implies that historical contingency strongly influences the evolution of new phenotypes.

Is there more than one way to skin a newt? Convergent toxin resistance in snakes is not due to a common genetic mechanism

Convergent evolution of tetrodotoxin (TTX) resistance, at both the phenotypic and genetic levels, characterizes coevolutionary arms races between amphibians and their snake predators around the world, and reveals remarkable predictability in the process of adaptation. Here we examine the repeatability of the evolution of TTX resistance in an undescribed predator-prey relationship between TTX-bearing Eastern Newts (Notophthalmus viridescens) and Eastern Hog-nosed Snakes (Heterodon platirhinos). We found that that local newts contain levels of TTX dangerous enough to dissuade most predators, and that Eastern Hog-nosed Snakes within newt range are highly resistant to TTX. In fact, these populations of Eastern Hog-nosed Snakes are so resistant to TTX that the potential for current reciprocal selection might be limited. Unlike all other cases of TTX resistance in vertebrates, H. platirhinos lacks the adaptive amino acid substitutions in the skeletal muscle sodium channel that reduce TTX binding, suggesting that physiological resistance in Eastern Hog-nosed Snakes is conferred by an alternate genetic mechanism. Thus, phenotypic convergence in this case is not due to parallel molecular evolution, indicating that there may be more than one way for this adaptation to arise, even among closely related species.

Natural bone fragmentation in the blind cave-dwelling fish, Astyanax mexicanus: candidate gene identification through integrative comparative genomics

Animals that colonize dark and nutrient-poor subterranean environments evolve numerous extreme phenotypes. These include dramatic changes to the craniofacial complex, many of which are under genetic control. These phenotypes can demonstrate asymmetric genetic signals wherein a QTL is detected on one side of the face but not the other. The causative gene(s) underlying QTL are difficult to identify with limited genomic resources. We approached this task by searching for candidate genes mediating fragmentation of the third suborbital bone (SO3) directly inferior to the orbit of the eye. We integrated positional genomic information using emerging Astyanax resources, and linked these intervals to homologous (syntenic) regions of the Danio rerio genome. We identified a discrete, approximately 6 Mb, conserved region wherein the gene causing SO3 fragmentation likely resides. We interrogated this interval for genes demonstrating significant differential expression using mRNA-seq analysis of cave and surface morphs across life history. We then assessed genes with known roles in craniofacial evolution and development based on GO term annotation. Finally, we screened coding sequence alterations in this region, identifying two key genes: transforming growth factor β3 (tgfb3) and bone morphogenetic protein 4 (bmp4). Of these candidates, tgfb3 is most promising as it demonstrates significant differential expression across multiple stages of development, maps close (<1 Mb) to the fragmentation critical locus, and is implicated in a variety of other animal systems (including humans) in non-syndromic clefting and malformations of the cranial sutures. Both abnormalities are analogous to the failure-to-fuse phenotype that we observe in SO3 fragmentation. This integrative approach will enable discovery of the causative genetic lesions leading to complex craniofacial features analogous to human craniofacial disorders. This work underscores the value of cave-dwelling fish as a powerful evolutionary model of craniofacial disease, and demonstrates the power of integrative system-level studies for informing the genetic basis of craniofacial aberrations in nature.

Replaying the tape of life in the twenty-first century

Should the tape of life be replayed, would it produce similar living beings? A classical answer has long been ‘no’, but accumulating data are now challenging this view. Repeatability in experimental evolution, in phenotypic evolution of diverse species and in the genes underlying phenotypic evolution indicates that despite unpredictability at the level of basic evolutionary processes (such as apparition of mutations), a certain kind of predictability can emerge at higher levels over long time periods. For instance, a survey of the alleles described in the literature that cause non-deleterious phenotypic differences among animals, plants and yeasts indicates that similar phenotypes have often evolved in distinct taxa through independent mutations in the same genes. Does this mean that the range of possibilities for evolution is limited? Does this mean that we can predict the outcomes of a replayed tape of life? Imagining other possible paths for evolution runs into four important issues: (i) resolving the influence of contingency, (ii) imagining living organisms that are different from the ones we know, (iii) finding the relevant concepts for predicting evolution, and (iv) estimating the probability of occurrence for complex evolutionary events that occurred only once during the evolution of life on earth.

Convergent Evolution During Local Adaptation to Patchy Landscapes

Species often encounter, and adapt to, many patches of similar environmental conditions across their range. Such adaptation can occur through convergent evolution if different alleles arise in different patches, or through the spread of shared alleles by migration acting to synchronize adaptation across the species. The tension between the two reflects the constraint imposed on evolution by the underlying genetic architecture versus how effectively selection and geographic isolation act to inhibit the geographic spread of locally adapted alleles. This paper studies the balance between these two routes to adaptation in a model of continuous environments with patchy selection pressures. We address the following questions: How long does it take for a novel allele to appear in a patch where it is locally adapted through mutation? Or, through migration from another, already adapted patch? Which is more likely to occur, as a function of distance between the patches? What population genetic signal is left by the spread of migrant alleles? To answer these questions we examine the family structure underlying migration-selection equilibrium surrounding an already adapted patch, treating those rare families that reach new patches as spatial branching processes. A main result is that patches further apart than a critical distance will likely evolve independent locally adapted alleles; this distance is proportional to the spatial scale of selection ([Formula: see text], where σ is the dispersal distance and sm is the selective disadvantage of these alleles between patches), and depends linearly on log(sm/μ), where μ is the mutation rate. This provides a way to understand the role of geographic separation between patches in promoting convergent adaptation and the genomic signals it leaves behind. We illustrate these ideas using the convergent evolution of cryptic coloration in the rock pocket mouse, Chaetodipus intermedius, as an empirical example.

The genomics of ecological vicariance in threespine stickleback fish

Populations occurring in similar habitats and displaying similar phenotypes are increasingly used to explore parallel evolution at the molecular level. This generally ignores the possibility that parallel evolution can be mimicked by the fragmentation of an ancestral population followed by genetic exchange with ecologically different populations. Here we demonstrate such an ecological vicariance scenario in multiple stream populations of threespine stickleback fish divergent from a single adjacent lake population. On the basis of demographic and population genomic analyses, we infer the initial spread of a stream-adapted ancestor followed by the emergence of a lake-adapted population, that selective sweeps have occurred mainly in the lake population, that adaptive lake–stream divergence is maintained in the face of gene flow from the lake into the streams, and that this divergence involves major inversion polymorphisms also important to marine-freshwater stickleback divergence. Overall, our study highlights the need for a robust understanding of the demographic and selective history in evolutionary investigations.

Threespine stickleback fish are adapted to lake and stream habitats in Central Europe. Here, the authors show colonization of a lake basin by a stream-adapted ancestor, followed by the emergence of a lake-adapted population in the face of gene flow across lake–stream boundaries.

Extent of QTL Reuse During Repeated Phenotypic Divergence of Sympatric Threespine Stickleback

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution.

The Role of Standing Variation in Geographic Convergent Adaptation

The extent to which populations experiencing shared selective pressures adapt through a shared genetic response is relevant to many questions in evolutionary biology. In this article, we explore how standing genetic variation contributes to convergent genetic responses in a geographically spread population. Geographically limited dispersal slows the spread of each selected allele, hence allowing other alleles to spread before any one comes to dominate the population. When selectively equivalent alleles meet, their progress is substantially slowed, dividing the species range into a random tessellation, which can be well understood by analogy to a Poisson process model of crystallization. In this framework, we derive the geographic scale over which an allele dominates and the proportion of adaptive alleles that arise from standing variation. Finally, we explore how negative pleiotropic effects of alleles can bias the subset of alleles that contribute to the species' adaptive response. We apply the results to the malaria-resistance glucose-6-phosphate dehydrogenase-deficiency alleles, where the large mutational target size makes it a likely candidate for adaptation from deleterious standing variation. Our results suggest that convergent adaptation may be common. Therefore, caution must be exercised when arguing that strongly geographically restricted alleles are the outcome of local adaptation. We close by discussing the implications of these results for ideas of species coherence and the nature of divergence between species.

Elucidating the molecular architecture of adaptation via evolve and resequence experiments

Evolve and resequence (E&R) experiments use experimental evolution to adapt populations to a novel environment, then next-generation sequencing to analyse genetic changes. They enable molecular evolution to be monitored in real time on a genome-wide scale. Here, we review the field of E&R experiments across diverse systems, ranging from simple non-living RNA to bacteria, yeast and the complex multicellular organism Drosophila melanogaster. We explore how different evolutionary outcomes in these systems are largely consistent with common population genetics principles. Differences in outcomes across systems are largely explained by different starting population sizes, levels of pre-existing genetic variation, recombination rates and adaptive landscapes. We highlight emerging themes and inconsistencies that future experiments must address.

Divergence and gene flow among Darwin's finches: A genome-wide view of adaptive radiation driven by interspecies allele sharing

A recent analysis of the genomes of Darwin's finches revealed extensive interspecies allele sharing throughout the history of the radiation and identified a key locus responsible for morphological evolution in this group. The radiation of Darwin's finches on the Galápagos archipelago has long been regarded as an iconic study system for field ecology and evolutionary biology. Coupled with an extensive history of field work, these latest findings affirm the increasing acceptance of introgressive hybridization, or gene flow between species, as a significant contributor to adaptive evolution. Here, we review and discuss these findings in relation to both classical work on Darwin's finches and contemporary work showing similar evolutionary signatures in other biological systems. The continued unification of genomic data with field biology promises to further elucidate the molecular basis of adaptation in Darwin's finches and well beyond.

Heterochrony underpins natural variation in Cardamine hirsuta leaf form

A key problem in biology is whether the same processes underlie morphological variation between and within species. Here, by using plant leaves as an example, we show that the causes of diversity at these two evolutionary scales can be divergent. Some species like the model plant Arabidopsis thaliana have simple leaves, whereas others like the A. thaliana relative Cardamine hirsuta bear complex leaves comprising leaflets. Previous work has shown that these interspecific differences result mostly from variation in local tissue growth and patterning. Now, by cloning and characterizing a quantitative trait locus (QTL) for C. hirsuta leaf shape, we find that a different process, age-dependent progression of leaf form, underlies variation in this trait within species. This QTL effect is caused by cis-regulatory variation in the floral repressor ChFLC, such that genotypes with low-expressing ChFLC alleles show both early flowering and accelerated age-dependent changes in leaf form, including faster leaflet production. We provide evidence that this mechanism coordinates leaf development with reproductive timing and may help to optimize resource allocation to the next generation.

The genetic architecture of sexually selected traits in two natural populations of Drosophila montana

We investigated the genetic architecture of courtship song and cuticular hydrocarbon traits in two phygenetically distinct populations of Drosophila montana. To study natural variation in these two important traits, we analysed within-population crosses among individuals sampled from the wild. Hence, the genetic variation analysed should represent that available for natural and sexual selection to act upon. In contrast to previous between-population crosses in this species, no major quantitative trait loci (QTLs) were detected, perhaps because the between-population QTLs were due to fixed differences between the populations. Partitioning the trait variation to chromosomes suggested a broadly polygenic genetic architecture of within-population variation, although some chromosomes explained more variation in one population compared with the other. Studies of natural variation provide an important contrast to crosses between species or divergent lines, but our analysis highlights recent concerns that segregating variation within populations for important quantitative ecological traits may largely consist of small effect alleles, difficult to detect with studies of moderate power.

Genetic Changes to a Transcriptional Silencer Element Confers Phenotypic Diversity within and between Drosophila Species

The modification of transcriptional regulation has become increasingly appreciated as a major contributor to morphological evolution. However, the role of negative-acting control elements (e.g. silencers) in generating morphological diversity has been generally overlooked relative to positive-acting “enhancer” elements. The highly variable body coloration patterns among Drosophilid insects represents a powerful model system in which the molecular alterations that underlie phenotypic diversity can be defined. In a survey of pigment phenotypes among geographically disparate Japanese populations of Drosophila auraria, we discovered a remarkable degree of variation in male-specific abdominal coloration. In testing the expression patterns of the major pigment-producing enzymes, we found that phenotypes uniquely correlated with differences in the expression of ebony, a gene required for yellow-colored cuticle. Assays of ebony’s transcriptional control region indicated that a lightly pigmented strain harbored cis-regulatory mutations that caused correlated changes in its expression. Through a series of chimeric reporter constructs between light and dark strain alleles, we localized function-altering mutations to a conserved silencer that mediates a male-specific pattern of ebony repression. This suggests that the light allele was derived through the loss of this silencer’s activity. Furthermore, examination of the ebony gene of D. serrata, a close relative of D. auraria which secondarily lost male-specific pigmentation revealed the parallel loss of this silencer element. These results demonstrate how loss-of-function mutations in a silencer element resulted in increased gene expression. We propose that the mutational inactivation of silencer elements may represent a favored path to evolve gene expression, impacting morphological traits.

One of the greatest challenges in understanding the relationship between genotype and phenotype is to discern how changes in DNA affect the normal functioning of genes. Mutations may generate a new function for a gene, yet it is frequently observed that they inactivate some aspect of a gene’s normal capacity. Investigations focused on understanding the developmental basis for the evolution of anatomical structures has found a prevalent role for mutations that alter developmental gene regulation. In animals, genes are transcriptionally activated in specific tissues during development by regulatory sequences distributed across their expansive non-protein coding regions. Regulatory elements known as silencers act to prevent genes from being expressed in certain tissues, providing a mechanism for precise control. Here, we show how a silencer that prevents expression of a pigment-producing enzyme in certain Drosophila species has repeatedly been subject to inactivating mutations that increased this gene’s expression. This example illustrates how such negative-acting regulatory sequences can represent a convenient target for increasing gene expression through the loss of a genetic element.

Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes

The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.

The differential view of genotype-phenotype relationships

An integrative view of diversity and singularity in the living world requires a better understanding of the intricate link between genotypes and phenotypes. Here we re-emphasize the old standpoint that the genotype-phenotype (GP) relationship is best viewed as a connection between two differences, one at the genetic level and one at the phenotypic level. As of today, predominant thinking in biology research is that multiple genes interact with multiple environmental variables (such as abiotic factors, culture, or symbionts) to produce the phenotype. Often, the problem of linking genotypes and phenotypes is framed in terms of genotype and phenotype maps, and such graphical representations implicitly bring us away from the differential view of GP relationships. Here we show that the differential view of GP relationships is a useful explanatory framework in the context of pervasive pleiotropy, epistasis, and environmental effects. In such cases, it is relevant to view GP relationships as differences embedded into differences. Thinking in terms of differences clarifies the comparison between environmental and genetic effects on phenotypes and helps to further understand the connection between genotypes and phenotypes.

Wax, sex and the origin of species: Dual roles of insect cuticular hydrocarbons in adaptation and mating

Evolutionary changes in traits that affect both ecological divergence and mating signals could lead to reproductive isolation and the formation of new species. Insect cuticular hydrocarbons (CHCs) are potential examples of such dual traits. They form a waxy layer on the cuticle of the insect to maintain water balance and prevent desiccation, while also acting as signaling molecules in mate recognition and chemical communication. Because the synthesis of these hydrocarbons in insect oenocytes occurs through a common biochemical pathway, natural or sexual selection on one role may affect the other. In this review, we explore how ecological divergence in insect CHCs can lead to divergence in mating signals and reproductive isolation. We suggest that the evolution of insect CHCs may be ripe models for understanding ecological speciation.

What to compare and how: Comparative transcriptomics for Evo-Devo

Evolutionary developmental biology has grown historically from the capacity to relate patterns of evolution in anatomy to patterns of evolution of expression of specific genes, whether between very distantly related species, or very closely related species or populations. Scaling up such studies by taking advantage of modern transcriptomics brings promising improvements, allowing us to estimate the overall impact and molecular mechanisms of convergence, constraint or innovation in anatomy and development. But it also presents major challenges, including the computational definitions of anatomical homology and of organ function, the criteria for the comparison of developmental stages, the annotation of transcriptomics data to proper anatomical and developmental terms, and the statistical methods to compare transcriptomic data between species to highlight significant conservation or changes. In this article, we review these challenges, and the ongoing efforts to address them, which are emerging from bioinformatics work on ontologies, evolutionary statistics, and data curation, with a focus on their implementation in the context of the development of our database Bgee (http://bgee.org). J. Exp. Zool. (Mol. Dev. Evol.) 324B: 372–382, 2015. © 2015 The Authors. J. Exp. Zool. (Mol. Dev. Evol.) published by Wiley Periodicals, Inc.

The evolutionary origination and diversification of a dimorphic gene regulatory network through parallel innovations in cis and trans

The origination and diversification of morphological characteristics represents a key problem in understanding the evolution of development. Morphological traits result from gene regulatory networks (GRNs) that form a web of transcription factors, which regulate multiple cis-regulatory element (CRE) sequences to control the coordinated expression of differentiation genes. The formation and modification of GRNs must ultimately be understood at the level of individual regulatory linkages (i.e., transcription factor binding sites within CREs) that constitute the network. Here, we investigate how elements within a network originated and diversified to generate a broad range of abdominal pigmentation phenotypes among Sophophora fruit flies. Our data indicates that the coordinated expression of two melanin synthesis enzymes, Yellow and Tan, recently evolved through novel CRE activities that respond to the spatial patterning inputs of Hox proteins and the sex-specific input of Bric-à-brac transcription factors. Once established, it seems that these newly evolved activities were repeatedly modified by evolutionary changes in the network’s trans-regulators to generate large-scale changes in pigment pattern. By elucidating how yellow and tan are connected to the web of abdominal trans-regulators, we discovered that the yellow and tan abdominal CREs are composed of distinct regulatory inputs that exhibit contrasting responses to the same Hox proteins and Hox cofactors. These results provide an example in which CRE origination underlies a recently evolved novel trait, and highlights how coordinated expression patterns can evolve in parallel through the generation of unique regulatory linkages.

The genomic content of regulatory genes such as transcription factors is surprisingly conserved between diverse animal species, raising the paradox of how new traits emerge, and are subsequently modified and lost. In this study we make a connection between the developmental basis for the formation of a fruit fly trait and the evolutionary basis for that trait’s origin, diversification, and loss. We show how the origin of a novel pigmentation trait is associated with the evolution of two regulatory sequences that control the co-expression of two key pigmentation genes. These sequences interact in unique ways with evolutionarily conserved Hox transcription factors to drive gene co-expression. Once these unique connections evolved, the alteration of this trait appears to have proceeded through changes to regulatory genes rather than regulatory sequences of the pigmentation genes. Thus, our findings support a scenario where regulatory sequence evolution provided new functions to old transcription factors, how co-expression can emerge from different utilizations of the same transcription factors, and that trait diversity was surprisingly shaped by changes in some manner to the deeply conserved regulatory genes.

Conservatism and novelty in the genetic architecture of adaptation in Heliconius butterflies

Understanding the genetic architecture of adaptive traits has been at the centre of modern evolutionary biology since Fisher; however, evaluating how the genetic architecture of ecologically important traits influences their diversification has been hampered by the scarcity of empirical data. Now, high-throughput genomics facilitates the detailed exploration of variation in the genome-to-phenotype map among closely related taxa. Here, we investigate the evolution of wing pattern diversity in Heliconius, a clade of neotropical butterflies that have undergone an adaptive radiation for wing-pattern mimicry and are influenced by distinct selection regimes. Using crosses between natural wing-pattern variants, we used genome-wide restriction site-associated DNA (RAD) genotyping, traditional linkage mapping and multivariate image analysis to study the evolution of the architecture of adaptive variation in two closely related species: Heliconius hecale and H. ismenius. We implemented a new morphometric procedure for the analysis of whole-wing pattern variation, which allows visualising spatial heatmaps of genotype-to-phenotype association for each quantitative trait locus separately. We used the H. melpomene reference genome to fine-map variation for each major wing-patterning region uncovered, evaluated the role of candidate genes and compared genetic architectures across the genus. Our results show that, although the loci responding to mimicry selection are highly conserved between species, their effect size and phenotypic action vary throughout the clade. Multilocus architecture is ancestral and maintained across species under directional selection, whereas the single-locus (supergene) inheritance controlling polymorphism in H. numata appears to have evolved only once. Nevertheless, the conservatism in the wing-patterning toolkit found throughout the genus does not appear to constrain phenotypic evolution towards local adaptive optima.

How consistent are the transcriptome changes associated with cold acclimation in two species of the Drosophila virilis group?

For many organisms the ability to cold acclimate with the onset of seasonal cold has major implications for their fitness. In insects, where this ability is widespread, the physiological changes associated with increased cold tolerance have been well studied. Despite this, little work has been done to trace changes in gene expression during cold acclimation that lead to an increase in cold tolerance. We used an RNA-Seq approach to investigate this in two species of the Drosophila virilis group. We found that the majority of genes that are differentially expressed during cold acclimation differ between the two species. Despite this, the biological processes associated with the differentially expressed genes were broadly similar in the two species. These included: metabolism, cell membrane composition, and circadian rhythms, which are largely consistent with previous work on cold acclimation/cold tolerance. In addition, we also found evidence of the involvement of the rhodopsin pathway in cold acclimation, a pathway that has been recently linked to thermotaxis. Interestingly, we found no evidence of differential expression of stress genes implying that long-term cold acclimation and short-term stress response may have a different physiological basis.

Two developmentally temporal quantitative trait loci underlie convergent evolution of increased branchial bone length in sticklebacks

In convergent evolution, similar phenotypes evolve repeatedly in independent populations, often reflecting adaptation to similar environments. Understanding whether convergent evolution proceeds via similar or different genetic and developmental mechanisms offers insight towards the repeatability and predictability of evolution. Oceanic populations of threespine stickleback fish, Gasterosteus aculeatus, have repeatedly colonized countless freshwater lakes and streams, where new diets lead to morphological adaptations related to feeding. Here, we show that heritable increases in branchial bone length have convergently evolved in two independently derived freshwater stickleback populations. In both populations, an increased bone growth rate in juveniles underlies the convergent adult phenotype, and one population also has a longer cartilage template. Using F2 crosses from these two freshwater populations, we show that two quantitative trait loci (QTL) control branchial bone length at distinct points in development. In both populations, a QTL on chromosome 21 controls bone length throughout juvenile development, and a QTL on chromosome 4 controls bone length only in adults. In addition to these similar developmental profiles, these QTL show similar chromosomal locations in both populations. Our results suggest that sticklebacks have convergently evolved longer branchial bones using similar genetic and developmental programmes in two independently derived populations.

A recurrent regulatory change underlying altered expression and Wnt response of the stickleback armor plates gene EDA

Armor plate changes in sticklebacks are a classic example of repeated adaptive evolution. Previous studies identified ectodysplasin (EDA) gene as the major locus controlling recurrent plate loss in freshwater fish, though the causative DNA alterations were not known. Here we show that freshwater EDA alleles have cis-acting regulatory changes that reduce expression in developing plates and spines. An identical T → G base pair change is found in EDA enhancers of divergent low-plated fish. Recreation of the T → G change in a marine enhancer strongly reduces expression in posterior armor plates. Bead implantation and cell culture experiments show that Wnt signaling strongly activates the marine EDA enhancer, and the freshwater T → G change reduces Wnt responsiveness. Thus parallel evolution of low-plated sticklebacks has occurred through a shared DNA regulatory change, which reduces the sensitivity of an EDA enhancer to Wnt signaling, and alters expression in developing armor plates while preserving expression in other tissues.

DOI:http://dx.doi.org/10.7554/eLife.05290.001

Stickleback fish develop bony plates on their surface to protect themselves from predators. The extent and pattern of their bony armor depends on their habitat: marine sticklebacks are typically covered from head to tail with bony plates, but freshwater sticklebacks retain only a few plates on their sides.

One gene that promotes the formation of the bony plates is called ectodysplasin (EDA). This encodes a signaling protein that is important for the development of the skeleton, skin and many other tissues. Variations in the sequence of this gene are shared among different stickleback populations worldwide. However, it has not been clear which genetic changes can explain how lightly armored freshwater sticklebacks could have evolved from their well-armored marine ancestors on several separate occasions.

Here, O'Brown et al. studied EDA in marine and groups of freshwater sticklebacks that have evolved in different locations around the world. The experiments show that the level of expression of EDA in the developing plates and spines is lower in the freshwater fish. O'Brown et al. thought this could be due to genetic changes in regions of EDA that lie outside the region that encodes the protein, so called ‘regulatory elements’.

Indeed, further experiments found that all freshwater fish have a small change in the DNA of a regulatory element that switches on the gene in plate-forming regions of the body. When this change was introduced into marine sticklebacks, the fish had lower levels of gene expression in these plate-forming regions.

These findings demonstrate that lightly armored sticklebacks have evolved multiple times from their well-armored marine ancestors through the same small change in their DNA that alters the expression of the EDA gene. The next challenge will be to understand why this particular small change in DNA appears to be favored over all the other changes that could occur in the regulatory element, and to see if factors that act through this regulatory switch also modify armor structures in natural populations.

DOI:http://dx.doi.org/10.7554/eLife.05290.002