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

A long noncoding RNA at the cortex locus controls adaptive coloration in butterflies

Evolutionary variation in the wing pigmentation of butterflies and moths offers striking examples of adaptation by crypsis and mimicry. The cortex locus has been independently mapped as the locus controlling color polymorphisms in 15 lepidopteran species, suggesting that it acts as a genomic hotspot for the diversification of wing patterns, but functional validation through protein-coding knockouts has proven difficult to obtain. Our study unveils the role of a long noncoding RNA (lncRNA) which we name ivory, transcribed from the cortex locus, in modulating color patterning in butterflies. Strikingly, ivory expression prefigures most melanic patterns during pupal development, suggesting an early developmental role in specifying scale identity. To test this, we generated CRISPR mosaic knock-outs in five nymphalid butterfly species and show that ivory mutagenesis yields transformations of dark pigmented scales into white or light-colored scales. Genotyping of Vanessa cardui germline mutants associates these phenotypes to small on-target deletions at the conserved first exon of ivory. In contrast, cortex germline mutant butterflies with confirmed null alleles lack any wing phenotype and exclude a color patterning role for this adjacent gene. Overall, these results show that a lncRNA gene acts as a master switch of color pattern specification and played key roles in the adaptive diversification of wing patterns in butterflies.

Optix and cortex/ivory/mir-193 again: the repeated use of two mimicry hotspot loci

The extent to which evolution is repeatable has been a debated topic among evolutionary biologists. Although rewinding the tape of life perhaps would not lead to the same outcome every time, repeated evolution of analogous genes for similar functions has been extensively reported. Wing phenotypes of butterflies and moths have provided a wealth of examples of gene re-use, with certain 'hotspot loci' controlling wing patterns across diverse taxa. Here, we present an example of convergent evolution in the molecular genetic basis of Batesian wing mimicry in two Hypolimnas butterfly species. We show that mimicry is controlled by variation near cortex/ivory/mir-193, a known butterfly hotspot locus. By dissecting the genetic architecture of mimicry in Hypolimnas misippus and Hypolimnas bolina, we present evidence that distinct non-coding regions control the development of white pattern elements in the forewing and hindwing of the two species, suggesting independent evolution, and that no structural variation is found at the locus. Finally, we also show that orange coloration in H. bolina is associated with optix, a well-known patterning gene. Overall, our study once again implicates variation near the hotspot loci cortex/ivory/mir-193 and optix in butterfly wing mimicry and thereby highlights the repeatability of adaptive evolution.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Stable inversion clines in a grasshopper species group despite complex geographical history

Chromosomal inversions are known to play roles in adaptation and differentiation in many species. They involve clusters of correlated genes (i.e., loci in linkage disequilibrium, LD) possibly associated with environmental variables. The grasshopper "species complex" Trimerotropis pallidipennis comprises several genetic lineages distributed from North to South America in arid and semi-arid high-altitude environments. The southernmost lineage, Trimerotropis sp., segregates for four to seven putative inversions that display clinal variation, possibly through adaptation to temperate environments. We analysed chromosomal, mitochondrial and genome-wide single nucleotide polymorphism data in 19 Trimerotropis sp. populations mainly distributed along two altitudinal gradients (MS and Ju). Populations across Argentina comprise two main chromosomally and genetically differentiated lineages: one distributed across the southernmost border of the "Andes Centrales," adding evidence for a differentiation hotspot in this area; and the other widely distributed in Argentina. Within the latter, network analytical approaches to LD found three clusters of correlated loci (LD-clusters), with inversion karyotypes explaining >79% of the genetic variation. Outlier loci associated with environmental variables mapped to two of these LD-clusters. Furthermore, despite the complex geographical history indicated by population genetic analyses, the clines in inversion karyotypes have remained stable for more than 20 generations, implicating their role in adaptation and differentiation within this lineage. We hypothesize that these clines could be the consequence of a coupling between extrinsic postzygotic barriers and spatially varying selection along environmental gradients resulting in a hybrid zone. These results provide a framework for future investigations about candidate genes implicated in rapid adaptation to new environments.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

Inversion breakpoints and the evolution of supergenes

The coexistence of discrete morphs that differ in multiple traits is common within natural populations of many taxa. Such morphs are often associated with chromosomal inversions, presumably because the recombination suppressing effects of inversions help maintain alternate adaptive combinations of alleles across the multiple loci affecting these traits. However, inversions can also harbour selected mutations at their breakpoints, leading to their rise in frequency in addition to (or independent from) their role in recombination suppression. In this review, we first describe the different ways that breakpoints can create mutations. We then critically examine the evidence for the breakpoint-mutation and recombination suppression hypotheses for explaining the existence of discrete morphs associated with chromosomal inversions. We find that the evidence that inversions are favoured due to recombination suppression is often indirect. The evidence that breakpoints harbour mutations that are adaptive is also largely indirect, with the characterization of inversion breakpoints at the sequence level being incomplete in most systems. Direct tests of the role of suppressed recombination and breakpoint mutations in inversion evolution are thus needed. Finally, we emphasize how the two hypotheses of recombination suppression and breakpoint mutation can act in conjunction, with implications for understanding the emergence of supergenes and their evolutionary dynamics. We conclude by discussing how breakpoint characterization could improve our understanding of complex, discrete phenotypic forms in nature.

Porcupine/Wntless-dependent trafficking of the conserved WntA ligand in butterflies

Wnt ligands are key signaling molecules in animals, but little is known about the evolutionary dynamics and mode of action of the WntA orthologs, which are not present in the vertebrates or in Drosophila. Here we show that the WntA subfamily evolved at the base of the Bilateria + Cnidaria clade, and conserved the thumb region and Ser209 acylation site present in most other Wnts, suggesting WntA requires the core Wnt secretory pathway. WntA proteins are distinguishable from other Wnts by a synapomorphic Iso/Val/Ala216 amino-acid residue that replaces the otherwise ubiquitous Thr216 position. WntA embryonic expression is conserved between beetles and butterflies, suggesting functionality, but the WntA gene was lost three times within arthropods, in podoplean copepods, in the cyclorrhaphan fly radiation, and in ensiferan crickets and katydids. Finally, CRISPR mosaic knockouts (KOs) of porcupine and wntless phenocopied the pattern-specific effects of WntA KOs in the wings of Vanessa cardui butterflies. These results highlight the molecular conservation of the WntA protein across invertebrates, and imply it functions as a typical Wnt ligand that is acylated and secreted through the Porcupine/Wntless secretory pathway.

Genomic variation, population history and within-archipelago adaptation between island bird populations

Oceanic island archipelagos provide excellent models to understand evolutionary processes. Colonization events and gene flow can interact with selection to shape genetic variation at different spatial scales. Landscape-scale variation in biotic and abiotic factors may drive fine-scale selection within islands, while long-term evolutionary processes may drive divergence between distantly related populations. Here, we examine patterns of population history and selection between recently diverged populations of the Berthelot's pipit (Anthus berthelotii), a passerine endemic to three North Atlantic archipelagos. First, we use demographic trees and f3 statistics to show that genome-wide divergence across the species range is largely shaped by colonization and bottlenecks, with evidence of very weak gene flow between populations. Then, using a genome scan approach, we identify signatures of divergent selection within archipelagos at single nucleotide polymorphisms (SNPs) in genes potentially associated with craniofacial development and DNA repair. We did not detect within-archipelago selection at the same SNPs as were detected previously at broader spatial scales between archipelagos, but did identify signatures of selection at loci associated with similar biological functions. These findings suggest that similar ecological factors may repeatedly drive selection between recently separated populations, as well as at broad spatial scales across varied landscapes.

Convergence in sympatry: Evolution of blue-banded wing pattern in Morpho butterflies

Species interactions such as mimicry can promote trait convergence but disentangling this effect from those of shared ecology, evolutionary history, and niche conservatism is often challenging. Here by focusing on wing colour pattern variation within and between three butterfly species living in sympatry in a large proportion of their range, we tested the effect of species interactions on trait diversification. These butterflies display a conspicuous iridescent blue coloration on the dorsal side of their wings and a cryptic brownish colour on the ventral side. Combined with an erratic and fast flight, these colour patterns increase the difficulty of capture by predators and contribute to the high escape abilities of these butterflies. We hypothesize that, beyond their direct contribution to predator escape, these wing patterns can be used as signals of escape abilities by predators, resulting in positive frequency-dependent selection favouring convergence in wing pattern in sympatry. To test this hypothesis, we quantified dorsal wing pattern variations of 723 butterflies from the three species sampled throughout their distribution, including sympatric and allopatric situations and compared the phenotypic distances between species, sex and localities. We detected a significant effect of localities on colour pattern, and higher inter-specific resemblance in sympatry as compared to allopatry, consistent with the hypothesis of local convergence of wing patterns. Our results provide support to the existence of escape mimicry in the wild and stress the importance of estimating trait variation within species to understand trait variation between species, and to a larger extent, trait diversification at the macro-evolutionary scale.

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.

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first- and second-generation hybrids produced by controlled crosses between either two co-mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major-effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.

The combined use of raw and phylogenetically independent methods of outlier detection uncovers genome-wide dynamics of local adaptation in a lizard

Local adaptation is a dynamic process by which different allele combinations are selected in different populations at different times, and whose genetic signature can be inferred by genome-wide outlier analyses. We combined gene flow estimates with two methods of outlier detection, one of them independent of population coancestry (CIOA) and the other one not (ROA), to identify genetic variants favored when ecology promotes phenotypic convergence. We analyzed genotyping-by-sequencing data from five populations of a lizard distributed over an environmentally heterogeneous range that has been changing since the split of eastern and western lineages ca. 3 mya. Overall, western lizards inhabit forest habitat and are unstriped, whereas eastern ones inhabit shrublands and are striped. However, one population (Lerma) has unstriped phenotype despite its eastern ancestry. The analysis of 73,291 SNPs confirmed the east-west division and identified nonoverlapping sets of outliers (12 identified by ROA and 9 by CIOA). ROA revealed ancestral adaptive variation in the uncovered outliers that were subject to divergent selection and differently fixed for eastern and western populations at the extremes of the environmental gradient. Interestingly, such variation was maintained in Lerma, where we found high levels of heterozygosity for ROA outliers, whereas CIOA uncovered innovative variants that were selected only there. Overall, it seems that both the maintenance of ancestral variation and asymmetric migration have counterbalanced adaptive lineage splitting in our model species. This scenario, which is likely promoted by a changing and heterogeneous environment, could hamper ecological speciation of locally adapted populations despite strong genetic structure between lineages.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency-dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator-prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once-paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

Discrete or indiscrete? Redefining the colour polymorphism of the land snail Cepaea nemoralis

Biologists have long tried to describe and name the different phenotypes that make up the shell polymorphism of the land snail Cepaea nemoralis. Traditionally, the view is that the ground colour of the shell is one of a few major colour classes, either yellow, pink or brown, but in practise it is frequently difficult to distinguish the colours, and define different shades of the same colour. To understand whether colour variation is in reality continuous, and to investigate how the variation may be perceived by an avian predator, we applied psychophysical models of colour vision to shell reflectance measures. We found that both achromatic and chromatic variation are indiscrete in Cepaea nemoralis, being continuously distributed over many perceptual units. Nonetheless, clustering analysis based on the density of the distribution did reveal three groups, roughly corresponding to human-perceived yellow, pink and brown shells. We also found large-scale geographic variation in the frequency of these groups across Europe, and some covariance between shell colour and banding patterns. Although further studies are necessary, the observation of continuous variation in colour is intriguing because the traditional theory is that the underlying supergene that determines colour has evolved to prevent phenotypes from “dissolving” into continuous trait distributions. The findings thus have significance for understanding the Cepaea polymorphism, and the nature of the selection that acts upon it, as well as more generally highlighting the need to measure colour objectively in other systems.

Subtle variation in size and shape of the whole forewing and the red band among co-mimics revealed by geometric morphometric analysis in Heliconius butterflies

Heliconius are unpalatable butterflies that exhibit remarkable intra- and interspecific variation in wing color pattern, specifically warning coloration. Species that have converged on the same pattern are often clustered in Müllerian mimicry rings. Overall, wing color patterns are nearly identical among co-mimics. However, fine-scale differences exist, indicating that factors in addition to natural selection may underlie wing phenotype. Here, we investigate differences in shape and size of the forewing and the red band in the Heliconius postman mimicry ring (H. erato phyllis and the co-mimics H. besckei, H. melpomene burchelli, and H. melpomene nanna) using a landmark-based approach. If phenotypic evolution is driven entirely by predation pressure, we expect nonsignificant differences among co-mimics in terms of wing shape. Also, a reinforcement of wing pattern (i.e., greater similarity) could occur when co-mimics are in sympatry. We also examined variation in the red forewing band because this trait is critical for both mimicry and sexual communication. Morphometric results revealed significant but small differences among species, particularly in the shape of the forewing of co-mimics. Although we did not observe greater similarity when co-mimics were in sympatry, nearly identical patterns provided evidence of convergence for mimicry. In contrast, mimetic pairs could be distinguished based on the shape (but not the size) of the red band, suggesting an "advergence" process. In addition, sexual dimorphism in the red band shape (but not size) was found for all lineages. Thus, we infer that natural selection due to predation by birds might not be the only mechanism responsible for variation in color patterns, and sexual selection could be an important driver of wing phenotypic evolution in this mimicry ring.

Waiting in the wings: what can we learn about gene co-option from the diversification of butterfly wing patterns?

A major challenge is to understand how conserved gene regulatory networks control the wonderful diversity of form that we see among animals and plants. Butterfly wing patterns are an excellent example of this diversity. Butterfly wings form as imaginal discs in the caterpillar and are constructed by a gene regulatory network, much of which is conserved across the holometabolous insects. Recent work in Heliconius butterflies takes advantage of genomic approaches and offers insights into how the diversification of wing patterns is overlaid onto this conserved network. WntA is a patterning morphogen that alters spatial information in the wing. Optix is a transcription factor that acts later in development to paint specific wing regions red. Both of these loci fit the paradigm of conserved protein-coding loci with diverse regulatory elements and developmental roles that have taken on novel derived functions in patterning wings. These discoveries offer insights into the 'Nymphalid Ground Plan', which offers a unifying hypothesis for pattern formation across nymphalid butterflies. These loci also represent 'hotspots' for morphological change that have been targeted repeatedly during evolution. Both convergent and divergent evolution of a great diversity of patterns is controlled by complex alleles at just a few genes. We suggest that evolutionary change has become focused on one or a few genetic loci for two reasons. First, pre-existing complex cis-regulatory loci that already interact with potentially relevant transcription factors are more likely to acquire novel functions in wing patterning. Second, the shape of wing regulatory networks may constrain evolutionary change to one or a few loci. Overall, genomic approaches that have identified wing patterning loci in these butterflies offer broad insight into how gene regulatory networks evolve to produce diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Predator Perspective Drives Geographic Variation in Frequency-Dependent Polymorphism

Color polymorphism in natural populations can manifest as a striking patchwork of phenotypes in space, with neighboring populations characterized by dramatic differences in morph composition. These geographic mosaics can be challenging to explain in the absence of localized selection because they are unlikely to result from simple isolation-by-distance or clinal variation in selective regimes. To identify processes that can lead to the formation of geographic mosaics, we developed a simulation-based model to explore the influence of predator perspective, selection, migration, and genetic linkage of color loci on allele frequencies in polymorphic populations over space and time. Using simulated populations inspired by the biology of Heliconius longwing butterflies, Cepaea land snails, Oophaga poison frogs, and Sonora ground snakes, we found that the relative sizes of predator and prey home ranges can produce large differences in morph composition between neighboring populations under both positive and negative frequency-dependent selection. We also demonstrated the importance of the interaction of predator perspective with the type of frequency dependence and localized directional selection across migration and selection intensities. Our results show that regional-scale predation can promote the formation of phenotypic mosaics in prey species, without the need to invoke spatial variation in selective regimes. We suggest that predator behavior can play an important and underappreciated role in the formation and maintenance of geographic mosaics in polymorphic species.

Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures

Explaining the maintenance of adaptive diversity within populations is a long-standing goal in evolutionary biology, with important implications for conservation, medicine, and agriculture. Adaptation often leads to the fixation of beneficial alleles, and therefore it erodes local diversity so that understanding the coexistence of multiple adaptive phenotypes requires deciphering the ecological mechanisms that determine their respective benefits. Here, we show how antagonistic frequency-dependent selection (FDS), generated by natural and sexual selection acting on the same trait, maintains mimicry polymorphism in the toxic butterfly Heliconius numata Positive FDS imposed by predators on mimetic signals favors the fixation of the most abundant and best-protected wing-pattern morph, thereby limiting polymorphism. However, by using mate-choice experiments, we reveal disassortative mate preferences of the different wing-pattern morphs. The resulting negative FDS on wing-pattern alleles is consistent with the excess of heterozygote genotypes at the supergene locus controlling wing-pattern variation in natural populations of H. numata The combined effect of positive and negative FDS on visual signals is sufficient to maintain a diversity of morphs displaying accurate mimicry with other local prey, although some of the forms only provide moderate protection against predators. Our findings help understand how alternative adaptive phenotypes can be maintained within populations and emphasize the need to investigate interactions between selective pressures in other cases of puzzling adaptive polymorphism.

Young inversion with multiple linked QTLs under selection in a hybrid zone

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favoured alleles at multiple loci. However, it is unknown whether favoured mutations slowly accumulate on older inversions or if young inversions spread because they capture pre-existing adaptive quantitative trait loci (QTLs). By genetic mapping, chromosome painting and genome sequencing, we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Crossing fitness valleys: empirical estimation of a fitness landscape associated with polymorphic mimicry

Characterizing fitness landscapes associated with polymorphic adaptive traits enables investigation of mechanisms allowing transitions between fitness peaks. Here, we explore how natural selection can promote genetic mechanisms preventing heterozygous phenotypes from falling into non-adaptive valleys. Polymorphic mimicry is an ideal system to investigate such fitness landscapes, because the direction of selection acting on complex mimetic colour patterns can be predicted by the local mimetic community composition. Using more than 5000 artificial butterflies displaying colour patterns exhibited by the polymorphic Müllerian mimic Heliconius numata, we directly tested the role of wild predators in shaping fitness landscapes. We compared predation rates on mimetic phenotypes (homozygotes at the supergene controlling colour pattern), intermediate phenotypes (heterozygotes), exotic morphs (absent from the local community) and palatable cryptic phenotypes. Exotic morphs were significantly more attacked than local morphs, highlighting predators' discriminatory capacities. Overall, intermediates were attacked twice as much as local homozygotes, suggesting the existence of deep fitness valleys promoting strict dominance and reduced recombination between supergene alleles. By including information on predators' colour perception, we also showed that protection on intermediates strongly depends on their phenotypic similarity to homozygous phenotypes and that ridges exist between similar phenotypes, which may facilitate divergence in colour patterns.

Complex modular architecture around a simple toolkit of wing pattern genes

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

Genes controlling mimetic colour pattern variation in butterflies

Butterfly wing patterns are made up of arrays of coloured scales. There are two genera in which within-species variation in wing patterning is common and has been investigated at the molecular level, Heliconius and Papilio. Both of these species have mimetic relationships with other butterfly species that increase their protection from predators. Heliconius have a 'tool-kit' of five genetic loci that control colour pattern, three of which have been identified at the gene level, and which have been repeatedly used to modify colour pattern by different species in the genus. By contrast, the three Papilio species that have been investigated each have different genetic mechanisms controlling their polymorphic wing patterns.

The genetic basis of discrete and quantitative colour variation in the polymorphic lizard, Ctenophorus decresii

Background

Colour polymorphic species provide invaluable insight into processes that generate and maintain intra-specific variation. Despite an increasing understanding of the genetic basis of discrete morphs, sources of colour variation within morphs remain poorly understood. Here we use the polymorphic tawny dragon lizard Ctenophorus decresii to test simple Mendelian models for the inheritance of discrete morphs, and to investigate the genetic basis of continuous variation among individuals across morphs. Males of this species express either orange, yellow, orange surrounded by yellow, or grey throats. Although four discrete morphs are recognised, the extent of orange and yellow varies greatly. We artificially elevated testosterone in F0 females and F1 juveniles to induce them to express the male throat colour polymorphism, and quantified colour variation across the pedigree.

Results

Inheritance of discrete morphs in C. decresii best fit a model whereby two autosomal loci with complete dominance respectively determine the presence of orange and yellow. However, a single locus model with three co-dominant alleles for orange, yellow and grey could not be definitively rejected. Additionally, quantitative expression of the proportion of orange and yellow on the throat was strongly heritable (orange: h² = 0.84 ± 0.14; yellow: h² = 0.67 ± 0.19), with some evidence for covariance between the two.

Conclusions

Our study supports the theoretical prediction that polymorphism should be governed by few genes of major effect, but implies broader genetic influence on variation in constituent morph traits.

Electronic supplementary material

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

Does 3D Phenotyping Yield Substantial Insights in the Genetics of the Mouse Mandible Shape?

We describe the application of high-resolution 3D microcomputed tomography, together with 3D landmarks and geometric morphometrics, to validate and further improve previous quantitative genetic studies that reported QTL responsible for variation in the mandible shape of laboratory mice using a new backcross between C57BL/6J and A/J inbred strains. Despite the increasing availability of 3D imaging techniques, artificial flattening of the mandible by 2D imaging techniques seems at first an acceptable compromise for large-scale phenotyping protocols, thanks to an abundance of low-cost digital imaging systems such as microscopes or digital cameras. We evaluated the gain of information from considering explicitly this additional third dimension, and also from capturing variation on the bone surface where no precise anatomical landmark can be marked. Multivariate QTL mapping conducted with different landmark configurations (2D vs. 3D; manual vs. semilandmarks) broadly agreed with the findings of previous studies. Significantly more QTL (23) were identified and more precisely mapped when the mandible shape was captured with a large set of semilandmarks coupled with manual landmarks. It appears that finer phenotypic characterization of the mandibular shape with 3D landmarks, along with higher density genotyping, yields better insights into the genetic architecture of mandibular development. Most of the main variation is, nonetheless, preferentially embedded in the natural 2D plane of the hemi-mandible, reinforcing the results of earlier influential investigations.

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.

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.