The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system.

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

Dominance mechanisms in supergene alleles controlling butterfly wing pattern variation: insights from gene expression in Heliconius numata.

Loci under balancing selection, where multiple alleles are maintained, offer a relevant opportunity to investigate the role of natural selection in shaping genetic dominance: the high frequency of heterozygotes at these loci has been shown to enable the evolution of dominance among alleles. In the butterfly Heliconius numata, mimetic wing color variations are controlled by an inversion polymorphism of a circa 2 Mb genomic region (supergene P), with strong dominance between sympatric alleles. To test how differences in dominance observed on wing patterns correlate with variations in expression levels throughout the supergene region, we sequenced the complete transcriptome of heterozygotes at the prepupal stage and compared it to corresponding homozygotes. By defining dominance based on non-overlapping ranges of transcript expression between genotypes, we found contrasting patterns of dominance between the supergene and the rest of the genome; the patterns of transcript expression in the heterozygotes were more similar to the expression observed in the dominant homozygotes in the supergene region. Dominance also differed among the three subinversions of the supergene, suggesting possible epistatic interactions among their gene contents underlying dominance evolution. We found the expression pattern of the melanization gene cortex located in the P-region to predict wing pattern phenotype in the heterozygote. We also identify new candidate genes that are potentially involved in mimetic color pattern variations highlighting the relevance of transcriptomic analyses in heterozygotes to pinpoint candidate genes in non-recombining regions.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system.

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which multiple species share the same conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which vivid colour and pattern are produced in a Müllerian mimicry complex of poison frogs. We did this by first assembling a high-quality de novo genome assembly for the mimic poison frog Ranitomeya imitator. This assembled genome is 6.8 Gbp in size, with a contig N50 of 300 Kbp R. imitator and two colour morphs from both Ranitomeya fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes that are differentially expressed throughout development, many of them related to melanocyte development, melanin synthesis, iridophore development and guanine synthesis. Polytypic differences within species may be the result of differences in expression and/or timing of expression, whereas convergence for colour pattern between species does not appear to be due to the same changes in gene expression. In addition, we identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species. Finally, we hypothesize that genes in the keratin family are important for producing different structural colours within these frogs.

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.

Warning signals are under positive frequency-dependent selection in nature.

Positive frequency-dependent selection (FDS) is a selection regime where the fitness of a phenotype increases with its frequency, and it is thought to underlie important adaptive strategies resting on signaling and communication. However, whether and how positive FDS truly operates in nature remains unknown, which hampers our understanding of signal diversity. Here, we test for positive FDS operating on the warning color patterns of chemically defended butterflies forming multiple coexisting mimicry assemblages in the Amazon. Using malleable prey models placed in localities showing differences in the relative frequencies of warningly colored prey, we demonstrate that the efficiency of a warning signal increases steadily with its local frequency in the natural community, up to a threshold where protection stabilizes. The shape of this relationship is consistent with the direct effect of the local abundance of each warning signal on the corresponding avoidance knowledge of the local predator community. This relationship, which differs from purifying selection acting on each mimetic pattern, indicates that predator knowledge, integrated over the entire community, is saturated only for the most common warning signals. In contrast, among the well-established warning signals present in local prey assemblages, most are incompletely known to local predators and enjoy incomplete protection. This incomplete predator knowledge should generate strong benefits to life history traits that enhance warning efficiency by increasing the effective frequency of prey visible to predators. Strategies such as gregariousness or niche convergence between comimics may therefore readily evolve through their effects on predator knowledge and warning efficiency.

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.

Selection of Valid Reference Genes for Reverse Transcription Quantitative PCR Analysis in Heliconius numata (Lepidoptera: Nymphalidae).

Identifying the genetic basis of adaptive variation is challenging in non-model organisms and quantitative real time PCR. is a useful tool for validating predictions regarding the expression of candidate genes. However, comparing expression levels in different conditions requires rigorous experimental design and statistical analyses. Here, we focused on the neotropical passion-vine butterflies Heliconius, non-model species studied in evolutionary biology for their adaptive variation in wing color patterns involved in mimicry and in the signaling of their toxicity to predators. We aimed at selecting stable reference genes to be used for normalization of gene expression data in RT-qPCR analyses from developing wing discs according to the minimal guidelines described in Minimum Information for publication of Quantitative Real-Time PCR Experiments (MIQE). To design internal RT-qPCR controls, we studied the stability of expression of nine candidate reference genes (actin, annexin, eF1α, FK506BP, PolyABP, PolyUBQ, RpL3, RPS3A, and tubulin) at two developmental stages (prepupal and pupal) using three widely used programs (GeNorm, NormFinder and BestKeeper). Results showed that, despite differences in statistical methods, genes RpL3, eF1α, polyABP, and annexin were stably expressed in wing discs in late larval and pupal stages of Heliconius numata This combination of genes may be used as a reference for a reliable study of differential expression in wings for instance for genes involved in important phenotypic variation, such as wing color pattern variation. Through this example, we provide general useful technical recommendations as well as relevant statistical strategies for evolutionary biologists aiming to identify candidate-genes involved adaptive variation in non-model organisms.

Unexpected colour pattern variation in mimetic frogs: implication for the diversification of warning signals in the genus Ranitomeya.

Predation is expected to promote uniformity in the warning coloration of defended prey, but also mimicry convergence between aposematic species. Despite selection constraining both colour-pattern and population divergence, many aposematic animals display numerous geographically structured populations with distinct warning signal. Here, we explore the extent of phenotypic variation of sympatric species of Ranitomeya poison frogs and test for theoretical expectations on variation and convergence in mimetic signals. We demonstrate that both warning signal and mimetic convergence are highly variable and are negatively correlated: some localities display high variability and no mimicry while in others the phenotype is fixed and mimicry is perfect. Moreover, variation in warning signals is always present within localities, and in many cases this variation overlaps between populations, such that variation is continuous. Finally, we show that coloration is consistently the least variable element and is likely of greater importance for predator avoidance compared to patterning. We discuss the implications of our results in the context of warning signal diversification and suggest that, like many other locally adapted traits, a combination of standing genetic variation and founding effect might be sufficient to enable divergence in colour pattern.

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

The role of predators in maintaining the geographic organization of aposematic signals.

Selective predation of aposematic signals is expected to promote phenotypic uniformity. But while these signals may be uniform within a population, numerous species display impressive variations in warning signals among adjacent populations. Predators from different localities who learn to avoid distinct signals while performing intense selection on others are thus expected to maintain such a geographic organization. We tested this assumption by placing clay frog models, representing distinct color morphs of the Peruvian poison dart frog Ranitomeya imitator and a nonconspicuous frog, reciprocally between adjacent localities. In each locality, avian predators were able to discriminate between warning signals; the adjacent exotic morph experienced up to four times more attacks than the local one and two times more than the nonconspicuous phenotype. Moreover, predation attempts on the exotic morph quickly decreased to almost nil, suggesting rapid learning. This experiment offers direct evidence for the existence of different predator communities performing localized homogenizing selection on distinct aposematic signals.

Advergence in Müllerian mimicry: the case of the poison dart frogs of Northern Peru revisited.

Whether the evolution of similar aposematic signals in different unpalatable species (i.e. Müllerian mimicry) is because of phenotypic convergence or advergence continues to puzzle scientists. The poison dart frog Ranitomeya imitator provides a rare example in support of the hypothesis of advergence: this species was believed to mimic numerous distinct model species because of high phenotypic variability and low genetic divergence among populations. In this study, we test the evidence in support of advergence using a population genetic framework in two localities where R. imitator is sympatric with different model species, Ranitomeya ventrimaculata and Ranitomeya variabilis. Genetic analyses revealed incomplete sorting of mitochondrial haplotypes between the two model species. These two species are also less genetically differentiated than R. imitator populations on the basis of both mitochondrial and nuclear DNA comparisons. The genetic similarity between the model species suggests that they have either diverged more recently than R. imitator populations or that they are still connected by gene flow and were misidentified as different species. An analysis of phenotypic variability indicates that the model species are as variable as R. imitator. These results do not support the hypothesis of advergence by R. imitator. Although we cannot rule out phenotypic advergence in the evolution of Müllerian mimicry, this study reopens the discussion regarding the direction of the evolution of mimicry in the R. imitator system.

Evolution and genetic architecture of disassortative mating at a locus under heterozygote advantage.

The evolution of mate choice is a major topic in evolutionary biology because it is thought to be a key factor in trait and species diversification. Here, we aim at uncovering the ecological conditions and genetic architecture enabling the puzzling evolution of disassortative mating based on adaptive traits. This rare form of mate choice is observed for some polymorphic traits but theoretical predictions on the emergence and persistence of this behavior are largely lacking. Thus, we developed a mathematical model to specifically understand the evolution of disassortative mating based on mimetic color pattern in the polymorphic butterfly Heliconius numata. We confirm that heterozygote advantage favors the evolution of disassortative mating and show that disassortative mating is more likely to emerge if at least one allele at the trait locus is free from any recessive deleterious mutations. We modeled different possible genetic architectures underlying mate choice behavior, such as self-referencing alleles, or specific preference or rejection alleles. Our results showed that self-referencing or rejection alleles linked to the color pattern locus enable the emergence of disassortative mating. However, rejection alleles allow the emergence of disassortative mating only when the color pattern and preference loci are tightly linked.

Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms.

Chromosomal inversions are ubiquitous in genomes and often coordinate complex phenotypes, such as the covariation of behavior and morphology in many birds, fishes, insects or mammals1-11. However, why and how inversions become associated with polymorphic traits remains obscure. Here we show that despite a strong selective advantage when they form, inversions accumulate recessive deleterious mutations that generate frequency-dependent selection and promote their maintenance at intermediate frequency. Combining genomics and in vivo fitness analyses in a model butterfly for wing-pattern polymorphism, Heliconius numata, we reveal that three ecologically advantageous inversions have built up a heavy mutational load from the sequential accumulation of deleterious mutations and transposable elements. Inversions associate with sharply reduced viability when homozygous, which prevents them from replacing ancestral chromosome arrangements. Our results suggest that other complex polymorphisms, rather than representing adaptations to competing ecological optima, could evolve because chromosomal rearrangements are intrinsically prone to carrying recessive harmful mutations.

Balanced polymorphisms and their divergence in a Heliconius butterfly.

The evolution of mimicry in similarly defended prey is well described by the Müllerian mimicry theory, which predicts the convergence of warning patterns in order to gain the most protection from predators. However, despite this prediction, we can find great diversity of color patterns among Müllerian mimics such as Heliconius butterflies in the neotropics. Furthermore, some species have evolved the ability to maintain multiple distinct warning patterns in single populations, a phenomenon known as polymorphic mimicry. The adaptive benefit of these polymorphisms is questionable since variation from the most common warning patterns is expected to be disadvantageous as novel signals are punished by predators naive to them. In this study, we use artificial butterfly models throughout Central and South America to characterize the selective pressures maintaining polymorphic mimicry in Heliconius doris. Our results highlight the complexity of positive frequency-dependent selection, the principal selective pressure driving convergence among Müllerian mimics, and its impacts on interspecific variation of mimetic warning coloration. We further show how this selection regime can both limit and facilitate the diversification of mimetic traits.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata.

BACKGROUND: Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata. RESULTS: We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns. CONCLUSIONS: This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Genomic architecture and introgression shape a butterfly radiation.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Territorial aggressiveness on the arboreal ant Azteca alfari by Camponotus blandus in French Guiana due to behavioural constraints.

This study reports new information on interactions between two sympatric ant species, the plant-ant Azteca alfari (Dolichoderinae) living in association with the myrmecophyte Cecropia obtusa (Cecropiaceae) and Camponotus blandus (Formicinae), a ground-nesting, arboreal-foraging species. Workers of A. alfari forage only on the foliage and the upper parts of the trunk of their host Cecropia, while C. blandus nests in the ground but frequently forages and patrols pioneer tree foliage, including Cecropia. The activity pattern of A. alfari and the number of C. blandus on Cecropia obtusa was monitored hourly during a two-day period in a disturbed area in French Guiana. The maximum activity of C. blandus occurred between 8:30 and 12:30, at which time A. alfari had retreated within the domatia and were least present on the trunks. Even though aggressive confrontations were observed, C. blandus workers often initiate confrontations but do not prey on A. alfari nor exploit food bodies produced by Cecropia, the principal food source of A. alfari. Hence hostility appears to be the result of territoriality. Differences in their foraging rhythms are proposed as promoting resource and territory partitioning in this ant assemblage.

Major Improvements to the Heliconius melpomene Genome Assembly Used to Confirm 10 Chromosome Fusion Events in 6 Million Years of Butterfly Evolution.

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, we incorporated ∼20x coverage of Pacific Biosciences sequencing, and scaffolded the haploid genome using an assembly of this long-read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an N50 length of 2.1 Mb, an N50 number of 34, and with 99% of the genome placed, and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million yr.

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.