Localization of Müllerian mimicry genes on a dense linkage map of Heliconius erato

We report a dense genetic linkage map of Heliconius erato, a neotropical butterfly that has undergone a remarkable adaptive radiation in warningly colored mimetic wing patterns. Our study exploited natural variation segregating in a cross between H. erato etylus and H. himera to localize wing color pattern loci on a dense linkage map containing amplified fragment length polymorphisms (AFLP), microsatellites, and single-copy nuclear loci. We unambiguously identified all 20 autosomal linkage groups and the sex chromosome (Z). The map spanned a total of 1430 Haldane cM and linkage groups varied in size from 26.3 to 97.8 cM. The average distance between markers was 5.1 cM. Within this framework, we localized two major color pattern loci to narrow regions of the genome. The first gene, D, responsible for red/orange elements, had a most likely placement in a 6.7-cM region flanked by two AFLP markers on the end of a large 87.5-cM linkage group. The second locus, Sd, affects the melanic pattern on the forewing and was found within a 6.3-cM interval between flanking AFLP loci. This study complements recent linkage analysis of H. erato's comimic, H. melpomene, and forms the basis for marker-assisted physical mapping and for studies into the comparative genetic architecture of wing-pattern mimicry in Heliconius.     

Convergent, modular expression of ebony and tan in the mimetic wing patterns of Heliconius butterflies

The evolution of pigmentation in vertebrates and flies has involved repeated divergence at a small number of genes related to melanin synthesis. Here, we study insect melanin synthesis genes in Heliconius butterflies, a group characterised by its diversity of wing patterns consisting of black (melanin), and yellow and red (ommochrome) pigmented scales. Consistent with their respective biochemical roles in Drosophila melanogaster, ebony is upregulated in non-melanic wing regions destined to be pigmented red whilst tan is upregulated in melanic regions. Wing regions destined to be pigmented yellow, however, are downregulated for both genes. This pattern is conserved across multiple divergent and convergent phenotypes within the Heliconii, suggesting a conserved mechanism for the development of black, red and yellow pattern elements across the genus. Linkage mapping of five melanin biosynthesis genes showed that, in contrast to other organisms, these genes do not control pattern polymorphism. Thus, the pigmentation genes themselves are not the locus of evolutionary change but lie downstream of a wing pattern regulatory factor. The results suggest a modular system in which particular combinations of genes are switched on whenever red, yellow or black pattern elements are favoured by natural selection for diverse and mimetic wing patterns.     

Shared and divergent expression domains on mimetic Heliconius wings

SUMMARY Heliconius butterfly wing patterns show repeated convergence between species and have adaptive value in mimicry and mate choice, offering an opportunity to connect adaptive changes in phenotype with their underlying genotypes. Here we study forewing ommochrome pigmentation in Heliconius melpomene . We clone two new ommochrome pathway genes for the Lepidoptera, karmoisin and kynurenine formamidase ( kf  ), and analyze the expression patterns of all known ommochrome genes across pupal wing development. In combination with published work, this generates the first comparative gene expression data for the co-mimics Heliconius erato and H. melpomene . In both species cinnabar expression correlates with the forewing band, but the expression pattern of vermillion differs significantly between the mimics. This demonstrates that both shared and divergent expression patterns are associated with mimetic phenotypes between Heliconius species. Two genes not studied in H. erato, scarlet and possibly kf , also show enhanced expression in the forewing band of H. melpomene , implying co-ordinated upregulation of several members of this biosynthetic pathway during pattern formation.     

Comparative population genetics of a mimicry locus among hybridizing Heliconius butterfly species

The comimetic Heliconius butterfly species pair, H. erato and H. melpomene, appear to use a conserved Mendelian switch locus to generate their matching red wing patterns. Here we investigate whether H. cydno and H. pachinus, species closely related to H. melpomene, use this same switch locus to generate their highly divergent red and brown color pattern elements. Using an F2 intercross between H. cydno and H. pachinus, we first map the genomic positions of two novel red/brown wing pattern elements; the G locus, which controls the presence of red vs brown at the base of the ventral wings, and the Br locus, which controls the presence vs absence of a brown oval pattern on the ventral hind wing. The results reveal that the G locus is tightly linked to markers in the genomic interval that controls red wing pattern elements of H. erato and H. melpomene. Br is on the same linkage group but approximately 26 cM away. Next, we analyze fine-scale patterns of genetic differentiation and linkage disequilibrium throughout the G locus candidate interval in H. cydno, H. pachinus and H. melpomene, and find evidence for elevated differentiation between H. cydno and H. pachinus, but no localized signature of association. Overall, these results indicate that the G locus maps to the same interval as the locus controlling red patterning in H. melpomene and H. erato. This, in turn, suggests that the genes controlling red pattern elements may be homologous across Heliconius, supporting the hypothesis that Heliconius butterflies use a limited suite of conserved genetic switch loci to generate both convergent and divergent wing patterns.     

Sexual Size Dimorphism in the Color Pattern Elements of Two Mimetic Heliconius Butterflies

Sexual dichromatism and sexual dimorphism of body size are reasonably well studied in butterflies. Sexual size dimorphism of color pattern elements, however, is much less explored. The object of this study is Heliconius, a genus of butterflies well known for the coevolution between mate color preferences and mimicry. Given the sexual role of wing coloration, we investigated the existence of sexual size dimorphism in the wing color elements of a mimetic pair—Heliconius erato phyllis Fabricius and Heliconius besckei Ménétriés—and analyzed the allometric patterns of these traits. Correlation between size of elements in the dorsal and ventral wing surfaces were also estimated. In both species, three out of four elements were larger in males, but the non-dimorphic element was not the same. With regard to the allometric patterns, our most important finding was that smaller males of one species have proportionally larger yellow bars. This is the first study specifically concerning quantitative sexual dimorphism in the coloration of this well-known genus of butterflies and it opens new prospects to investigate sex-related natural selection and sexual selection of color pattern elements.     

Multi-Allelic Major Effect Genes Interact with Minor Effect QTLs to Control Adaptive Color Pattern Variation in Heliconius erato

Recent studies indicate that relatively few genomic regions are repeatedly involved in the evolution of Heliconius butterfly wing patterns. Although this work demonstrates a number of cases where homologous loci underlie both convergent and divergent wing pattern change among different Heliconius species, it is still unclear exactly how many loci underlie pattern variation across the genus. To address this question for Heliconius erato, we created fifteen independent crosses utilizing the four most distinct color pattern races and analyzed color pattern segregation across a total of 1271 F2 and backcross offspring. Additionally, we used the most variable brood, an F2 cross between H. himera and the east Ecuadorian H. erato notabilis, to perform a quantitative genetic analysis of color pattern variation and produce a detailed map of the loci likely involved in the H. erato color pattern radiation. Using AFLP and gene based markers, we show that fewer major genes than previously envisioned control the color pattern variation in H. erato. We describe for the first time the genetic architecture of H. erato wing color pattern by assessing quantitative variation in addition to traditional linkage mapping. In particular, our data suggest three genomic intervals modulate the bulk of the observed variation in color. Furthermore, we also identify several modifier loci of moderate effect size that contribute to the quantitative wing pattern variation. Our results are consistent with the two-step model for the evolution of mimetic wing patterns in Heliconius and support a growing body of empirical data demonstrating the importance of major effect loci in adaptive change.     

Genomic islands of divergence in hybridizing Heliconius butterflies identified by large-scale targeted sequencing

Heliconius butterflies represent a recent radiation of species, in which wing pattern divergence has been implicated in speciation. Several loci that control wing pattern phenotypes have been mapped and two were identified through sequencing. These same gene regions play a role in adaptation across the whole Heliconius radiation. Previous studies of population genetic patterns at these regions have sequenced small amplicons. Here, we use targeted next-generation sequence capture to survey patterns of divergence across these entire regions in divergent geographical races and species of Heliconius. This technique was successful both within and between species for obtaining high coverage of almost all coding regions and sufficient coverage of non-coding regions to perform population genetic analyses. We find major peaks of elevated population differentiation between races across hybrid zones, which indicate regions under strong divergent selection. These 'islands' of divergence appear to be more extensive between closely related species, but there is less clear evidence for such islands between more distantly related species at two further points along the 'speciation continuum'. We also sequence fosmid clones across these regions in different Heliconius melpomene races. We find no major structural rearrangements but many relatively large (greater than 1 kb) insertion/deletion events (including gain/loss of transposable elements) that are variable between races.     

Wnt signaling underlies evolution and development of the butterfly wing pattern symmetry systems

Most butterfly wing patterns are proposed to be derived from a set of conserved pattern elements known as symmetry systems. Symmetry systems are so-named because they are often associated with parallel color stripes mirrored around linear organizing centers that run between the anterior and posterior wing margins. Even though the symmetry systems are the most prominent and diverse wing pattern elements, their study has been confounded by a lack of knowledge regarding the molecular basis of their development, as well as the difficulty of drawing pattern homologies across species with highly derived wing patterns. Here we present the first molecular characterization of symmetry system development by showing that WntA expression is consistently associated with the major basal, discal, central, and external symmetry system patterns of nymphalid butterflies. Pharmacological manipulations of signaling gradients using heparin and dextran sulfate showed that pattern organizing centers correspond precisely with WntA, wingless, Wnt6, and Wnt10 expression patterns, thus suggesting a role for Wnt signaling in color pattern induction. Importantly, this model is supported by recent genetic and population genomic work identifying WntA as the causative locus underlying wing pattern variation within several butterfly species. By comparing the expression of WntA between nymphalid butterflies representing a range of prototypical symmetry systems, slightly deviated symmetry systems, and highly derived wing patterns, we were able to infer symmetry system homologies in several challenging cases. Our work illustrates how highly divergent morphologies can be derived from modifications to a common ground plan across both micro- and macro-evolutionary time scales.     

Evolution of the wave: aerodynamic and aposematic functions of butterfly wing motion

Many unpalatable butterfly species use coloration to signal their distastefulness to birds, but motion cues may also be crucial to ward off predatory attacks. In previous research, captive passion-vine butterflies Heliconius mimetic in colour pattern were also mimetic in motion. Here, I investigate whether wing motion changes with the flight demands of different behaviours. If birds select for wing motion as a warning signal, aposematic butterflies should maintain wing motion independently of behavioural context. Members of one mimicry group (Heliconius cydno and Heliconius sapho) beat their wings more slowly and their wing strokes were more asymmetric than their sister-species (Heliconius melpomene and Heliconius erato, respectively), which were members of another mimicry group having a quick and steady wing motion. Within mimicry groups, wing beat frequency declined as its role in generating lift also declined in different behavioural contexts. In contrast, asymmetry of the stroke was not associated with wing beat frequency or behavioural context-strong indication that birds process and store the Fourier motion energy of butterfly wings. Although direct evidence that birds respond to subtle differences in butterfly wing motion is lacking, birds appear to generalize a motion pattern as much as they encounter members of a mimicry group in different behavioural contexts.     

Interspecific sexual attraction because of convergence in warning colouration: is there a conflict between natural and sexual selection in mimetic species?

When species converge in their colour patterns because of mimicry, and those patterns are also used in mate recognition, there is a probability of conflicting selection pressures. Closely related species that mimic one another are particularly likely to face such confusion because of similarities in their courtship behaviour and ecology. We conducted experiments in greenhouse conditions to study interspecific attraction between two mimetic butterfly species, Heliconius erato and Heliconius melpomene. Both species spent considerable time approaching and courting females of the co-mimic species. Experiments using wing models demonstrated the importance of colour pattern in this interspecific attraction. Although males of H. melpomene were attracted to their co-mimics as much as to their own females, H. erato males were more efficient at distinguishing conspecifics, possibly using wing odours. Although preliminary, these results suggest that the use of additional cues may have evolved in H. erato to reduce the cost of convergence in visual signals with H. melpomene. Overall, our results showed that there might be a cost of mimetic convergence because of a reduction in the efficiency of species recognition. Such cost may contribute to explain the apparently stable diversity in Müllerian mimetic patterns in many tropical butterfly assemblages.     

Homologies in the colour patterns of the genus Heliconius (Lepidoptera: Nymphalidae)

The colour patterns of Heliconius butterflies are composed from a relatively simple set of pattern elements whose homologues are recognizable throughout the genus. Although Heliconius colour patterns look quite different from those of most nymphalids, these pattern elements are seen to derive from the generalized nymphalid groundplan. The differences arise primarily from the loss or positional shift of certain pattern elements, a high degree of fusion between individual pattern elements, and, in the forewing, asymmetries of the pattern elements relative to the wing-cell midline. The scheme of homologies we present is consistent with what is currently known about the comparative morphology and developmental physiology of colour pattern formation in Lepidoptera, and provides a framework for the interpretation of developmental, evolutionary and genetic studies in Heliconius.     

PHENOTYPIC CORRELATION STRUCTURE AMONG ELEMENTS OF THE COLOR PATTERN IN PRECIS COENIA (LEPIDOPTERA: NYMPHALIDAE)

The component parts of butterfly wing patterns are arranged in sets of serially homologous pattern elements, repeated from wing cell to wing cell. Measurements were made on the sizes and positions of these elements on two successive, independent, sets of specimens in order to elucidate the phenotypic correlation structure among pattern elements. That portion of the correlation between measures due to overall size variation was accounted for through two alternate methods: multiple regression on two vein length measures, which represent wing size, and a Wright-style factor analysis. The sizes of pattern elements belonging to a homologous series were found to be significantly correlated whereas those of non-homologous elements varied independently. The degree of correlation among homologs varied, and, in the case of eyespot sizes, appeared to be inversely related to the degree of their morphological divergence. Although not correlated in size, the positions of non-homologous elements that lie within the same wing cell are moderately correlated. The results support current developmental models for the ontogeny of butterfly color pattern.     

Evolution of a mimicry supergene from a multilocus architecture

The origin and evolution of supergenes have long fascinated evolutionary biologists. In the polymorphic butterfly Heliconius numata, a supergene controls the switch between multiple different forms, and results in near-perfect mimicry of model species. Here, we use a morphometric analysis to quantify the variation in wing pattern observed in two broods of H. numata with different alleles at the supergene locus, 'P'. Further, we genotype the broods to associate the variation we capture with genetic differences. This allows us to begin mapping the quantitative trait loci that have minor effects on wing pattern. In addition to finding loci on novel chromosomes, our data, to our knowledge, suggest for the first time that ancestral colour-pattern loci, known to have major effects in closely related species, may contribute to the wing patterns displayed by H. numata, despite the large transfer of effects to the supergene.     

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.     

UV photoreceptors and UV-yellow wing pigments in Heliconius butterflies allow a color signal to serve both mimicry and intraspecific communication

Mimetic wing coloration evolves in butterflies in the context of predator confusion. Unless butterfly eyes have adaptations for discriminating mimetic color variation, mimicry also carries a risk of confusion for the butterflies themselves. Heliconius butterfly eyes, which express recently duplicated ultraviolet (UV) opsins, have such an adaptation. To examine bird and butterfly color vision as sources of selection on butterfly coloration, we studied yellow wing pigmentation in the tribe Heliconiini. We confirmed, using reflectance and mass spectrometry, that only Heliconius use 3-hydroxy-DL-kynurenine (3-OHK), which looks yellow to humans but reflects both UV- and long-wavelength light, whereas butterflies in related genera have chemically unknown yellow pigments mostly lacking UV reflectance. Modeling of these color signals reveals that the two UV photoreceptors of Heliconius are better suited to separating 3-OHK from non-3-OHK spectra compared with the photoreceptors of related genera or birds. The co-occurrence of potentially enhanced UV vision and a UV-reflecting yellow wing pigment could allow unpalatable Heliconius private intraspecific communication in the presence of mimics. Our results are the best available evidence for the correlated evolution of a color signal and color vision. They also suggest that predator visual systems are error prone in the context of mimicry.     

Mutants highlight the modular control of butterfly eyespot patterns

SUMMARY The eyespots on butterfly wings are thought to be serially homologous pattern elements. Yet eyespots differ greatly in number, shape, color, and size, within and among species. To what extent do these serially homologues have separate developmental identities, upon which selection acts to create diversity? We examined x‐ray–induced mutations for the eyespots of the nymphalid butterfly Bicyclus anynana that highlight the modular control of these serially homologous wing pattern elements. These mutations reduce or eliminate individual eyespots, or groups of eyespots, with no further effect on the wing color pattern. The collection of mutants highlights a greater potential developmental repertoire than that observed across the genus Bicyclus. We studied in detail one such mutation, of codominant effect, that causes the elimination of two adjacent eyespots on the ventral hindwing. By analyzing the expression of genes known to be involved in eyespot formation, we found an alteration in the differentiation of the “organizing” cells at the eyespot's center. No such cells differentiate in the wing subdivisions lacking the two eyespots in the mutants. We propose several developmental models, based on wing compartmentalization in Drosophila, that provide the first framework for thinking about the molecular evolution of butterfly wing pattern modularity.     

Elements of butterfly wing patterns

The color patterns on the wings of butterflies are unique among animal color patterns in that the elements that make up the overall pattern are individuated. Unlike the spots and stripes of vertebrate color patterns, the elements of butterfly wing patterns have identities that can be traced from species to species, and typically across genera and families. Because of this identity it is possible to recognize homologies among pattern elements and to study their evolution and diversification. Individuated pattern elements evolved from non-individuated precursors by compartmentalization of the wing into areas that became developmentally autonomous with respect to color pattern formation. Developmental compartmentalization led to the evolution of serially repeated elements and the emergence of serial homology. In these compartments, serial homologues were able to acquire site-specific developmental regulation and this, in turn, allowed them to diverge morphologically. Compartmentalization of the wing also reduced the developmental correlation among pattern elements. The release from this developmental constraint, we believe, enabled the great evolutionary radiation of butterfly wing patterns. During pattern evolution, the same set of individual pattern elements is arranged in novel ways to produce species-specific patterns, including such adaptations as mimicry and camouflage.     

A simulation study of mutations in the genetic regulatory hierarchy for butterfly eyespot focus determination

The color patterns on the wings of butterflies have been an important model system in evolutionary developmental biology. A recent computational model tested genetic regulatory hierarchies hypothesized to underlie the formation of butterfly eyespot foci [Evans, T.M., Marcus, J.M., 2006. A simulation study of the genetic regulatory hierarchy for butterfly eyespot focus determination. Evol. Dev. 8, 273-283]. The computational model demonstrated that one proposed hierarchy was incapable of reproducing the known patterns of gene expression associated with eyespot focus determination in wild-type butterflies, but that two slightly modified alternative hierarchies were capable of reproducing all of the known gene expressions patterns. Here we extend the computational models previously implemented in Delphi 2.0 to two mutants derived from the squinting bush brown butterfly (Bicyclus anynana). These two mutants, comet and Cyclops, have aberrantly shaped eyespot foci that are produced by different mechanisms. The comet mutation appears to produce a modified interaction between the wing margin and the eyespot focus that results in a series of comet-shaped eyespot foci. The Cyclops mutation causes the failure of wing vein formation between two adjacent wing-cells and the fusion of two adjacent eyespot foci to form a single large elongated focus in their place. The computational approach to modeling pattern formation in these mutants allows us to make predictions about patterns of gene expression, which are largely unstudied in butterfly mutants. It also suggests a critical experiment that will allow us to distinguish between two hypothesized genetic regulatory hierarchies that may underlie all butterfly eyespot foci.