Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation.

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

Enhancer evolutionary co-option through shared chromatin accessibility input.

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

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.

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone.

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

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.

Revisiting the developmental and cellular role of the pigmentation gene yellow in Drosophila using a tagged allele.

Pigmentation is a diverse and ecologically relevant trait in insects. Pigment formation has been studied extensively at the genetic and biochemical levels. The temporality of pigment formation during animal development, however, is more elusive. Here, we examine this temporality, focusing on yellow, a gene involved in the formation of black melanin. We generated a protein-tagged yellow allele in the fruit fly Drosophila melanogaster, which allowed us to precisely describe Yellow expression pattern at the tissue and cellular levels throughout development. We found Yellow expressed in the pupal epidermis in patterns prefiguring black pigmentation. We also found Yellow expressed in a few central neurons from the second larval instar to adult stages, including a subset of neurons adjacent to the clock neurons marked by the gene Pdf. We then specifically examined the dynamics of Yellow expression domain and subcellular localization in relationship to pigment formation. In particular, we showed how a late step of re-internalization is regulated by the large low-density lipoprotein receptor-related protein Megalin. Finally we suggest a new function for Yellow in the establishment of sharp pigmentation pattern boundaries, whereby this protein may assume a structural role, anchoring pigment deposits or pigmentation enzymes in the cuticle.

Does divergent selection predict the evolution of mate preference and reproductive isolation in the tropical butterfly genus Melinaea (Nymphalidae: Ithomiini)?

Many studies have shown that speciation can be facilitated when a trait under divergent selection also causes assortative mating. In Müllerian mimetic butterflies, a change in wing colour pattern can cause reproductive isolation. However, colour pattern divergence does not always lead to reproductive isolation. Understanding how divergent selection affects speciation requires identifying the mechanisms that promote mate preference and/or choosiness. This study addresses whether shifts in wing colour pattern drives mate preference and reproductive isolation in the tropical butterfly genus Melinaea (Nymphalidae: Ithomiini), and focuses on five taxa that form a speciation continuum, from subspecies to fully recognized species. Using genetic markers, wing colour pattern quantification, male pheromone characterization and behavioural assays of mating preference, we characterize the extent of genetic and phenotypic differentiation between taxa and compare it to the level of reproductive isolation. We show strong premating isolation between the closely related species M. satevis and M. marsaeus, in addition to genetic and phenotypic (colour pattern and pheromones) differentiation. By contrast, M. menophilus and M. marsaeus consist of pairs of subspecies that differ for colour pattern but that cannot be differentiated genetically. Pheromonal differentiation of subspecies was significant only for M. marsaeus, although most individuals were indistinguishable. Melinaea menophilus and M. marsaeus also differ in the strength of assortative mating, suggesting that mate preference has evolved only in M. marsaeus, consistent with selection against maladaptive offspring, as subspecific 'hybrids' of M. marsaeus have intermediate, non-mimetic colour patterns, unlike those of M. menophilus which display either parental phenotypes. We conclude that a shift in colour pattern per se is not sufficient for reproductive isolation, but rather, the evolution of assortative mating may be caused by selection against maladaptive intermediate phenotypes. This study suggests that mate preference and assortative mating evolve when adaptive, and that even in the early stages of divergence, reproductive isolation can be nearly complete due to mating preferences.

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.

Refining mimicry: phenotypic variation tracks the local optimum.

Müllerian mimicry between chemically defended preys is a textbook example of natural selection favouring phenotypic convergence onto a shared warning signal. Studies of mimicry have concentrated on deciphering the ecological and genetic underpinnings of dramatic switches in mimicry association, producing a well-known mosaic distribution of mimicry patterns across geography. However, little is known about the accuracy of resemblance between natural comimics when the local phenotypic optimum varies. In this study, using analyses of wing shape, pattern and hue, we quantify multimodal phenotypic similarity between butterfly comimics sharing the so-called postman pattern in different localities with varying species composition. We show that subtle but consistent variation between populations of the localized species, Heliconius timareta thelxinoe, enhance resemblance to the abundant comimics which drive the mimicry in each locality. Those results suggest that rarer comimics track the changes in the phenotypic optimum caused by gradual changes in the composition of the mimicry community, providing insights into the process by which intraspecific diversity of mimetic pattern may arise. Furthermore, our results suggest a multimodal evolution of similarity, with coordinated convergence in different features of the phenotype such as wing outline, pattern and hue. Finally, multilocus genotyping allows estimating local hybridization rates between H. timareta and comimic H. melpomene in different populations, raising the hypothesis that mimicry refinement between closely related comimics may be enhanced by adaptive introgression at loci modifying the accuracy of resemblance.

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'.

Regulatory encoding of quantitative variation in spatial activity of a Drosophila enhancer.

Developmental enhancers control the expression of genes prefiguring morphological patterns. The activity of an enhancer varies among cells of a tissue, but collectively, expression levels in individual cells constitute a spatial pattern of gene expression. How the spatial and quantitative regulatory information is encoded in an enhancer sequence is elusive. To link spatial pattern and activity levels of an enhancer, we used systematic mutations of the yellow spot enhancer, active in developing Drosophila wings, and tested their effect in a reporter assay. Moreover, we developed an analytic framework based on the comprehensive quantification of spatial reporter activity. We show that the quantitative enhancer activity results from densely packed regulatory information along the sequence, and that a complex interplay between activators and multiple tiers of repressors carves the spatial pattern. Our results shed light on how an enhancer reads and integrates trans-regulatory landscape information to encode a spatial quantitative pattern.

Morpho morphometrics: Shared ancestry and selection drive the evolution of wing size and shape in Morpho butterflies.

Butterfly wings harbor highly diverse phenotypes and are involved in many functions. Wing size and shape result from interactions between adaptive processes, phylogenetic history, and developmental constraints, which are complex to disentangle. Here, we focus on the genus Morpho (Nymphalidae: Satyrinae, 30 species), which presents a high diversity of sizes, shapes, and color patterns. First, we generate a comprehensive molecular phylogeny of these 30 species. Next, using 911 collection specimens, we quantify the variation of wing size and shape across species, to assess the importance of shared ancestry, microhabitat use, and sexual selection in the evolution of the wings. While accounting for phylogenetic and allometric effects, we detect a significant difference in wing shape but not size among microhabitats. Fore and hindwings covary at the individual and species levels, and the covariation differs among microhabitats. However, the microhabitat structure in covariation disappears when phylogenetic relationships are taken into account. Our results demonstrate that microhabitat has driven wing shape evolution, although it has not strongly affected forewing and hindwing integration. We also found that sexual dimorphism of forewing shape and color pattern are coupled, suggesting a common selective force.

Evolution of dominance mechanisms at a butterfly mimicry supergene.

Genetic dominance in polymorphic loci may respond to selection; however, the evolution of dominance in complex traits remains a puzzle. We analyse dominance at a wing-patterning supergene controlling local mimicry polymorphism in the butterfly Heliconius numata. Supergene alleles are associated with chromosomal inversion polymorphism, defining ancestral versus derived alleles. Using controlled crosses and the new procedure, Colour Pattern Modelling, allowing whole-wing pattern comparisons, we estimate dominance coefficients between alleles. Here we show strict dominance in sympatry favouring mimicry and inconsistent dominance throughout the wing between alleles from distant populations. Furthermore, dominance among derived alleles is uncoordinated across wing-pattern elements, producing mosaic heterozygous patterns determined by a hierarchy in colour expression. By contrast, heterozygotes with an ancestral allele show complete, coordinated dominance of the derived allele, independently of colours. Therefore, distinct dominance mechanisms have evolved in association with supergene inversions, in response to strong selection on mimicry polymorphism.

Wing shape variation associated with mimicry in butterflies.

Mimetic resemblance in unpalatable butterflies has been studied by evolutionary biologists for over a century, but has largely focused on the convergence in wing color patterns. In Heliconius numata, discrete color-pattern morphs closely resemble comimics in the distantly related genus Melinaea. We examine the possibility that the shape of the butterfly wing also shows adaptive convergence. First, simple measures of forewing dimensions were taken of individuals in a cross between H. numata morphs, and showed quantitative differences between two of the segregating morphs, f. elegans and f. silvana. Second, landmark-based geometric morphometric and elliptical Fourier outline analyses were used to more fully characterize these shape differences. Extension of these techniques to specimens from natural populations suggested that, although many of the coexisting morphs could not be discriminated by shape, the differences we identified between f. elegans and f. silvana hold in the wild. Interestingly, despite extensive overlap, the shape variation between these two morphs is paralleled in their respective Melinaea comimics. Our study therefore suggests that wing-shape variation is associated with mimetic resemblance, and raises the intriguing possibility that the supergene responsible for controlling the major switch in color pattern between morphs also contributes to wing shape differences in H. numata.

Roles of dental development and adaptation in rodent evolution.

In paleontology, many changes affecting morphology, such as tooth shape in mammals, are interpreted as ecological adaptations that reflect important selective events. Despite continuing studies, the identification of the genetic bases and key ecological drivers of specific mammalian dental morphologies remains elusive. Here we focus on the genetic and functional bases of stephanodonty, a pattern characterized by longitudinal crests on molars that arose in parallel during the diversification of murine rodents. We find that overexpression of Eda or Edar is sufficient to produce the longitudinal crests defining stephanodonty in transgenic laboratory mice. Whereas our dental microwear analyses show that stephanodonty likely represents an adaptation to highly fibrous diet, the initial and parallel appearance of stephanodonty may have been facilitated by developmental processes, without being necessarily under positive selection. This study demonstrates how combining development and function can help to evaluate adaptive scenarios in the evolution of new morphologies.