The genetic basis of an adaptive radiation: warning colour in two Heliconius species

Mimetic colour pattern races of Heliconius butterflies provide a striking example of adaptive radiation and numerous crossing experiments have investigated the genetics of these racial differences. However, colour pattern differentiation between closely related Heliconius species has not been previously studied. Here we present data from crosses between H. erato cyrbia and its sister species, H. himera. The genetic architecture underlying colour pattern divergence between these species is identical to that observed between races of H. erato. As in inter-racial crosses, colour pattern differences resulted from segregation at a few major loci. Evidence from 1321 offspring in 4 F₁, 17 backcross, 7 F₂ and 21 further crosses showed that two major loci controlled most of the colour pattern differences between H. erato and H. himera. There were strong interactions between these loci in their patterns of expression and evidence for other loci with relatively minor phenotypic effects. More importantly, based on patterns of expression within broods and linkage with Aconitase, we conclude that these major loci were homologous with those known to be responsible for colour pattern differences within H. erato. Our crosses also permit a re-evaluation of the relationships between colour pattern races of H. erato. This suggests that H. e. hydara, which occurs across a major mtDNA break, is the ancestral phenotype from which other races have evolved. Based on this assumption, we find no evidence to support the recent suggestion that apparently homologous colour pattern alleles have arisen multiple times.     

Characterization of microsatellite loci in neotropical Heliconius butterflies

The Heliconius butterflies offer exceptional opportunities for the study of the ecology and evolution of mimicry. Despite previous reports of difficulties in the development of microsatellite loci in Lepidoptera, we characterize 15 polymorphic loci in H. erato that show promise for genetic mapping and population studies in this and other species. Levels of variation were high, in both numbers and size ranges of alleles. The loci showed broad amplification success across the genus and in two other genera. All loci that amplified in a population of H. melpomene were polymorphic.     

Convergent evolution in the genetic basis of Müllerian mimicry in heliconius butterflies

The neotropical butterflies Heliconius melpomene and H. erato are Müllerian mimics that display the same warningly colored wing patterns in local populations, yet pattern diversity between geographic regions. Linkage mapping has previously shown convergent red wing phenotypes in these species are controlled by loci on homologous chromosomes. Here, AFLP bulk segregant analysis using H. melpomene crosses identified genetic markers tightly linked to two red wing-patterning loci. These markers were used to screen a H. melpomene BAC library and a tile path was assembled spanning one locus completely and part of the second. Concurrently, a similar strategy was used to identify a BAC clone tightly linked to the locus controlling the mimetic red wing phenotypes in H. erato. A methionine rich storage protein (MRSP) gene was identified within this BAC clone, and comparative genetic mapping shows red wing color loci are in homologous regions of the genome of H. erato and H. melpomene. Subtle differences in these convergent phenotypes imply they evolved independently using somewhat different developmental routes, but are nonetheless regulated by the same switch locus. Genetic mapping of MRSP in a third related species, the “tiger” patterned H. numata, has no association with wing patterning and shows no evidence for genomic translocation of wing-patterning loci.     

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.     

Wing patterning genes and coevolution of Müllerian mimicry in Heliconius butterflies: Support from phylogeography,cophylogeny, and divergence times

Examples of long‐term coevolution are rare among free‐living organisms. Müllerian mimicry in Heliconius butterflies had been suggested as a key example of coevolution by early genetic studies. However, research over the last two decades has been dominated by the idea that the best‐studied comimics, H. erato and H. melpomene, did not coevolve at all. Recently sequenced genes associated with wing color pattern phenotype offer a new opportunity to resolve this controversy. Here, we test the hypothesis of coevolution between H. erato and H. melpomene using Bayesian multilocus analysis of five color pattern genes and five neutral genetic markers. We first explore the extent of phylogenetic agreement versus conflict between the different genes. Coevolution is then tested against three aspects of the mimicry diversifications: phylogenetic branching patterns, divergence times, and, for the first time, phylogeographic histories. We show that all three lines of evidence are compatible with strict coevolution of the diverse mimicry wing patterns, contrary to some recent suggestions. Instead, these findings tally with a coevolutionary diversification driven primarily by the ecological force of Müllerian mimicry.     

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

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

Multiple Mechanisms of Photoreceptor Spectral Tuning in Heliconius Butterflies

The evolution of color vision is often studied through the lens of receptor gain relative to an ancestor with fewer spectral classes of photoreceptor. For instance, in Heliconius butterflies, a genus-specific UVRh opsin duplication led to the evolution of UV color discrimination in Heliconius erato females, a rare trait among butterflies. However, color vision evolution is not well understood in the context of loss. In Heliconius melpomene and Heliconius ismenius lineages, the UV2 receptor subtype has been lost, which limits female color vision in shorter wavelengths. Here, we compare the visual systems of butterflies that have either retained or lost the UV2 photoreceptor using intracellular recordings, ATAC-seq, and antibody staining. We identify several ways these butterflies modulate their color vision. In H. melpomene, chromatin reorganization has downregulated an otherwise intact UVRh2 gene, whereas in H. ismenius, pseudogenization has led to the truncation of UVRh2. In species that lack the UV2 receptor, the peak sensitivity of the remaining UV1 photoreceptor cell is shifted to longer wavelengths. Across Heliconius, we identify the widespread use of filtering pigments and co-expression of two opsins in the same photoreceptor cells. Multiple mechanisms of spectral tuning, including the molecular evolution of blue opsins, have led to the divergence of receptor sensitivities between species. The diversity of photoreceptor and ommatidial subtypes between species suggests that Heliconius visual systems are under varying selection pressures for color discrimination. Modulating the wavelengths of peak sensitivities of both the blue- and remaining UV-sensitive photoreceptor cells suggests that Heliconius species may have compensated for UV receptor loss.     

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome—that is, the genetic architecture of reproductive isolation—can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male‐transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.     

Cryptic genetic and wing pattern diversity in a mimetic Heliconius butterfly

Despite rampant colour pattern diversity in South America, Heliconius erato exhibits a ‘postman’ wing pattern throughout most of Central America. We examined genetic variation across the range of H. erato, including dense sampling in Central America, and discovered a deep genetic break, centred on the mountain range that runs through Costa Rica. This break is characterized by a novel mitochondrial lineage, which is nearly fixed in northern Central America, that branches basal to all previously described mitochondrial diversity in the species. Strong genetic differentiation also appears in Z‐linked and autosomal markers, and it is further associated with a distinct, but subtle, shift in wing pattern phenotype. Comparison of clines in wing phenotype, mtDNA and nuclear markers indicate they are all centred on the mountains dividing Costa Rica, but that cline width differs among data sets. Phylogeographical analyses, accounting for this new diversity, rewrite our understanding of mimicry evolution in this system. For instance, these results suggest that H. erato originated west of the Andes, perhaps in Central America, and as many as 1 million years before its co‐mimic, H. melpomene. Overall our data indicate that neutral genetic markers and colour pattern loci are congruent and converge on the same hypothesis—H. erato originated in northwest South America or Central America with a ‘postman’ phenotype and then radiated into the wealth of colour patterns present today.     

Signatures of selection in loci governing major colour patterns in <Emphasis Type="Italic">Heliconius</Emphasis> butterflies and related species

Background  

Protein-coding change is one possible genetic mechanism underlying the evolution of adaptive wing colour pattern variation in Heliconius butterflies. Here we determine whether 38 putative genes within two major Heliconius patterning loci, HmYb and HmB, show evidence of positive selection. Ratios of nonsynonymous to synonymous nucleotide changes (ω) were used to test for selection, as a means of identifying candidate genes within each locus that control wing pattern.     

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.     

Metamorphosis of a Butterfly-Associated Bacterial Community

Butterflies are charismatic insects that have long been a focus of biological research. They are also habitats for microorganisms, yet these microbial symbionts are little-studied, despite their likely importance to butterfly ecology and evolution. In particular, the diversity and composition of the microbial communities inhabiting adult butterflies remain uncharacterized, and it is unknown how the larval (caterpillar) and adult microbiota compare. To address these knowledge gaps, we used Illumina sequencing of 16S rRNA genes from internal bacterial communities associated with multiple life stages of the neotropical butterfly Heliconius erato. We found that the leaf-chewing larvae and nectar- and pollen-feeding adults of H. erato contain markedly distinct bacterial communities, a pattern presumably rooted in their distinct diets. Larvae and adult butterflies host relatively small and similar numbers of bacterial phylotypes, but few are common to both stages. The larval microbiota clearly simplifies and reorganizes during metamorphosis; thus, structural changes in a butterfly''s bacterial community parallel those in its own morphology. We furthermore identify specific bacterial taxa that may mediate larval and adult feeding biology in Heliconius and other butterflies. Although male and female Heliconius adults differ in reproductive physiology and degree of pollen feeding, bacterial communities associated with H. erato are not sexually dimorphic. Lastly, we show that captive and wild individuals host different microbiota, a finding that may have important implications for the relevance of experimental studies using captive butterflies.     

Genetic Evidence for Hybrid Trait Speciation in Heliconius Butterflies

Homoploid hybrid speciation is the formation of a new hybrid species without change in chromosome number. So far, there has been a lack of direct molecular evidence for hybridization generating novel traits directly involved in animal speciation. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has been proposed as a hybrid species, and its color pattern can be recreated by introgression of the H. m. melpomene red band into the genetic background of the yellow banded H. cydno cordula. This hybrid color pattern is also involved in mate choice and leads to reproductive isolation between H. heurippa and its close relatives. Here, we provide molecular evidence for adaptive introgression by sequencing genes across the Heliconius red band locus and comparing them to unlinked wing patterning genes in H. melpomene, H. cydno, and H. heurippa. 670 SNPs distributed among 29 unlinked coding genes (25,847bp) showed H. heurippa was related to H. c. cordula or the three species were intermixed. In contrast, among 344 SNPs distributed among 13 genes in the red band region (18,629bp), most showed H. heurippa related with H. c. cordula, but a block of around 6,5kb located in the 3′ of a putative kinesin gene grouped H. heurippa with H. m. melpomene, supporting the hybrid introgression hypothesis. Genealogical reconstruction showed that this introgression occurred after divergence of the parental species, perhaps around 0.43Mya. Expression of the kinesin gene is spatially restricted to the distal region of the forewing, suggesting a mechanism for pattern regulation. This gene therefore constitutes the first molecular evidence for adaptive introgression during hybrid speciation and is the first clear candidate for a Heliconius wing patterning locus.     

Estimating the relative fitness of local adaptive peaks: the aerodynamic costs of flight in mimetic passion-vine butterflies Heliconius

In a related paper, we demonstrated that mimetic Heliconius butterflies have converged in wing-beat frequency and degree of asymmetry in the wing motion, whereas sister species are dissimilar in these same traits. Warning signals of sympatric, distasteful species converge in evolutionary models in order to educate their predators more efficiently that the signal is associated with unprofitable prey. Barring other constraints, the behaviours of the different co-mimetic pairs should ultimately converge on that behaviour which minimizes the energetic cost of flight. We estimated the energetic cost of each mimic''s flight behaviour in order to predict the difference in height of each fitness peak and the direction of convergent selection qualitatively. Following adjustments for body mass, mimetic Heliconius melpomene and Heliconius erato required more aerodynamic power than Heliconius cydno and Heliconius sapho. This difference was attributed to the slower flight speeds and higher wing-beat frequencies of H. melpomene and H. erato. Consequently, H. melpomene and H. erato expended more energy per unit distance per unit body mass than H. cydno and H. sapho. However, differences in body mass may equalize energy budgets and stabilize the sympatric coexistence of the two pairs of co-mimics.     

Evaluating female remating rates in light of spermatophore degradation in Heliconius butterflies: pupal‐mating monandry versus adult‐mating polyandry

1. Butterflies are frequently used in comparative studies of sexual selection because of their diverse mating systems. In Heliconius, the two major clades in the genus are characterised by contrasting pupal‐mating and adult‐mating strategies. Adult‐mating females are considered to be promiscuous whereas pupal‐mating females are thought to be monandrous. 2. Counting spermatophores in female Lepidoptera is a common method for assessing patterns of female remating. However, in pupal‐mating Heliconius butterflies spermatophores can become completely degraded, potentially leading to underestimation of female remating rates. 3. We qualitatively characterised the different states of spermatophore degradation, and showed that complete degradation takes approximately 3 weeks in captive‐bred H. erato females. 4. We counted spermatophores and/or assayed spermatophore degradation in > 500 Heliconius females across 28 species sampled from natural populations. Among pupal‐maters these observations yielded a few rare observations of double mating by recently eclosed females, but generally indicated a lack of rematings. In contrast, approximately 25% of sampled adult‐mating females remated at least once. 5. Using a novel statistical analysis, we estimated remating rates from patterns of spermatophore degradation or from counts stratified by age, as indicated by wing‐wear. This analysis showed no statistically significant evidence for remating for the pupal‐mating H. erato whereas significant remating rates were detected for adult‐mating species. 6. The present results support the established view of Heliconius mating systems in which pupal‐maters are largely monandrous, whereas adult‐maters are polyandrous.     

STRONG NATURAL SELECTION IN A WARNING-COLOR HYBRID ZONE

Frequency-dependent selection on warning color can maintain narrow hybrid zones between unpalatable prey taxa. To measure such selection, we transferred marked Heliconius erato (Lepidoptera: Nymphalidae) in both directions across a 10-km-wide hybrid zone between Peruvian races differing in color pattern. These experimental H. erato were released at four sites, along with control H. erato of the phenotype native to each site. Survival of experimental butterflies was significantly lower than that of controls at two sites and overall. Most selection, measured as differences in survival, occurred soon after release. Selection against foreign morphs was 52% (confidence limits: 25–71 %) and was probably due to bird attacks on unusual warning-color morphs (more than 10% of the recaptures had beak marks). Since only three major loci determine the color-pattern differences, this suggests an average selection coefficient of 0.17 per locus, sufficient to maintain the narrow clines in H. erato.     

A major locus controls a biologically active pheromone component in Heliconius melpomene

Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal-rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species’ difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced F_ST. A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.     

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.     

Ultrastructure and morphogenesis of the wing scales in Heliconius erato phyllis (Lepidoptera: Nymphalidae): What silvery/brownish surfaces can tell us about the development of color patterning?

Usually the literature on Heliconius show three types of scales, classified based on the correlation between color and ultrastructure: type I – white and yellow, type II – black, and type III – orange and red. The ultrastructure of the scales located at the silvery/brownish surfaces of males/females is for the first time described in this paper. Besides, we describe the ontogeny of pigmentation, the scale morphogenesis and the maturation timing of scales fated to different colors in Heliconius erato phyllis. The silvery/brownish surfaces showed ultrastructurally similar scales to the type I, II and III. The ontogeny of pigmentation follows the sequence red, black, silvery/brownish and yellow. The maturation of yellow-fated scales, however, occurred simultaneously with the red-fated scales, before the pigmentation becomes visible. In spite of the scales at the silvery/brownish surfaces being ultrastructurally similar to the yellow, red and black scales, they mature after them; this suggests that the maturation timing does not show a relationship with the scale ultrastructure, with the deposition timing of the yellow pigment. The analysis of H. erato phyllis scale morphogenesis, as well as the scales ultrastructure and maturation timing, provided new findings into the developmental architecture of color pattern in Heliconius.