First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in 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.     

A set of microsatellite markers for Heliconius melpomene and closely related species

The butterflies in the genus Heliconius offer an exceptional opportunity for the study of the ecology and genetics of an adaptive radiation due to their extensive intra‐ and interspecific variation in wing colour patterns and mimetic associations. Here, we characterize 22 polymorphic microsatellite loci in Heliconius melpomene that have been shown to be useful for linkage mapping and population studies in this and other species. Levels of variation were high, although heterozygosity deficiencies were found in most loci, probably due to null alleles. The loci showed broad amplification success on six other species across the genus.     

Cryptic differences in colour among Müllerian mimics: how can the visual capacities of predators and prey shape the evolution of wing colours?

Antagonistic interactions between predators and prey often lead to co‐evolution. In the case of toxic prey, aposematic colours act as warning signals for predators and play a protective role. Evolutionary convergence in colour patterns among toxic prey evolves due to positive density‐dependent selection and the benefits of mutual resemblance in spreading the mortality cost of educating predators over a larger prey assemblage. Comimetic species evolve highly similar colour patterns, but such convergence may interfere with intraspecific signalling and recognition in the prey community, especially for species involved in polymorphic mimicry. Using spectrophotometry measures, we investigated the variation in wing coloration among comimetic butterflies from distantly related lineages. We focused on seven morphs of the polymorphic species Heliconius numata and the seven corresponding comimetic species from the genus Melinaea. Significant differences in the yellow, orange and black patches of the wing were detected between genera. Perceptions of these cryptic differences by bird and butterfly observers were then estimated using models of animal vision based on physiological data. Our results showed that the most strikingly perceived differences were obtained for the contrast of yellow against a black background. The capacity to discriminate between comimetic genera based on this colour contrast was also evaluated to be higher for butterflies than for birds, suggesting that this variation in colour, likely undetectable to birds, might be used by butterflies for distinguishing mating partners without losing the benefits of mimicry. The evolution of wing colour in mimetic butterflies might thus be shaped by the opposite selective pressures exerted by predation and species recognition.     

Pervasive genetic associations between traits causing reproductive isolation in Heliconius butterflies

Ecological speciation proceeds through the accumulation of divergent traits that contribute to reproductive isolation, but in the face of gene flow traits that characterize incipient species may become disassociated through recombination. Heliconius butterflies are well known for bright mimetic warning patterns that are also used in mate recognition and cause both pre- and post-mating isolation between divergent taxa. Sympatric sister taxa representing the final stages of speciation, such as Heliconius cydno and Heliconius melpomene, also differ in ecology and hybrid fertility. We examine mate preference and sterility among offspring of crosses between these species and demonstrate the clustering of Mendelian colour pattern loci and behavioural loci that contribute to reproductive isolation. In particular, male preference for red patterns is associated with the locus responsible for the red forewing band. Two further colour pattern loci are associated, respectively, with female mating outcome and hybrid sterility. This genetic architecture in which ‘speciation genes’ are clustered in the genome can facilitate two controversial models of speciation, namely divergence in the face of gene flow and hybrid speciation.     

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.     

A hybrid zone provides evidence for incipient ecological speciation in Heliconius butterflies

In Heliconius butterflies, it has been proposed that speciation occurs through a combination of divergence in ecological habitat preferences and mimetic colour patterns. Here we test this hypothesis by investigating a parapatric form of the widespread species Heliconius erato. Mendelian (colour patterns) and molecular genetic data permit us to address hypotheses about introgression and genetic differentiation between different populations. Combined analysis of colour pattern, microsatellite loci and mitochondrial DNA showed that Heliconius erato venus and Heliconius erato chestertonii form a bimodal hybrid zone implying partial reproductive isolation. In a sample of 121 individuals collected in sympatry, 25% were hybrids representing a significant deficit of heterozygotes compared to the Hardy-Weinberg expectation. Seven microsatellite loci, analysed for a subset of these individuals, showed marked differentiation between the parental taxa, and unambiguously identified two genotypic clusters concordant with our phenotypic classification of individuals. Mitochondrial DNA analysis showed H. erato venus as a monophyletic group well differentiated from H. erato chestertonii, implying a lack of historical introgression between the populations. Heliconius erato chestertonii is therefore an incipient species that maintains its integrity despite high levels of hybridization. Moreover, H. erato chestertonii is found at higher altitudes than other races of H. erato and has a distinct colour pattern and mimetic relationship. Hence, there are now two examples of parapatric incipient species related to H. erato, H. himera and H. erato chestertonii, both of which are associated with higher altitudes, more arid habitats and distinct mimetic relationships. This implies that parapatric habitat adaptation is a likely cause of speciation in this group.     

Two sisters in the same dress: <Emphasis Type="Italic">Heliconius</Emphasis> cryptic species

Background  

Sister species divergence and reproductive isolation commonly results from ecological adaptation. In mimetic Heliconius butterflies, shifts in colour pattern contribute to pre- and post-mating reproductive isolation and are commonly correlated with speciation. Closely related mimetic species are therefore not expected, as they should lack several important sources of reproductive isolation.     

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.     

Genetic constraints on wing pattern variation in Lycaeides butterflies: A case study on mapping complex,multifaceted traits in structured populations

Patterns of phenotypic variation within and among species can be shaped and constrained by trait genetic architecture. This is particularly true for complex traits, such as butterfly wing patterns, that consist of multiple elements. Understanding the genetics of complex trait variation across species boundaries is difficult, as it necessitates mapping in structured populations and can involve many loci with small or variable phenotypic effects. Here, we investigate the genetic architecture of complex wing pattern variation in Lycaeides butterflies as a case study of mapping multivariate traits in wild populations that include multiple nominal species or groups. We identify conserved modules of integrated wing pattern elements within populations and species. We show that trait covariances within modules have a genetic basis and thus represent genetic constraints that can channel evolution. Consistent with this, we find evidence that evolutionary changes in wing patterns among populations and species occur in the directions of genetic covariances within these groups. Thus, we show that genetic constraints affect patterns of biological diversity (wing pattern) in Lycaeides, and we provide an analytical template for similar work in other systems.     

Hybrid zones and the speciation continuum in Heliconius butterflies

Tropical butterflies in the genus Heliconius have long been models in the study of the stages of speciation. Heliconius are unpalatable to predators, and many species are notable for multiple geographic populations with striking warning colour pattern differences associated with Müllerian mimicry. A speciation continuum is evident in Heliconius hybrid zones. Examples range from hybrid zones across which (a) there is little genetic differentiation other than at mimicry loci, but where hybrids are common, (b) to ‘bimodal‘ hybrid zones with strong genetic divergence and few hybrids, (c) through to ‘good’ sympatric species, with hybridization extremely rare or absent. Now, in this issue of Molecular Ecology, Arias et al. ( 2012 ) have found an intermediate case in Colombian Heliconius cydno showing evidence for assortative mating and molecular differences, but where hybrids are abundant.     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

Ecological and genetic associations across a Heliconius hybrid zone

Differences in habitat use can bridge early and late stages of speciation by initiating assortative mating. Heliconius colour pattern races might select habitats over which each pattern confers a relative fitness advantage because signal efficacy of wing patterns can vary by environment. Thus habitat preferences could serve to promote the evolution of mimetic colour patterns for mate choice. Here I compare colour pattern genotype and phenotype frequencies to environmental variation across the H. erato hydara x H. erato erato hybrid zone in French Guiana to determine whether races exhibit habitat preferences. I found that genotype and phenotype frequencies correspond to differences in land cover moreso than to other environmental factors. Temporal shifts in colour pattern genotypes, phenotypes and land cover also were associated at individual sample sites, which further suggests that H. erato races differ in habitat use and that habitat preferences may promote speciation among Heliconius butterflies.     

Mimicry and the evolution of premating isolation in Heliconius melpomene Linnaeus

Ecological divergence can cause speciation if adaptive traits have pleiotropic effects on mate choice. In Heliconius butterflies, mimetic patterns play a role in mate detection between sister species, as well as signalling to predators. Here we show that male butterflies from four recently diverged parapatric populations of Heliconius melpomene are more likely to approach and court their own colour patterns as compared with those of other races. A few exceptions, where males were more attracted to patterns other than their own, suggest that some mimetic patterns are sub-optimal in mate choice. Genotype frequencies in hybrid zones between races of H. melpomene suggest that mating is random, so reinforcement is unlikely to have played a role in intra-specific divergence. In summary, co-evolved divergence of colour pattern and mate preference occurs rapidly and is likely the first step in Heliconius speciation.     

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.     

Mimicry: developmental genes that contribute to speciation

Despite renewed interest in the role of natural selection as a catalyst for the origin of species, the developmental and genetic basis of speciation remains poorly understood. Here we describe the genetics of Müllerian mimicry in Heliconius cydno and H. melpomene (Lepidoptera: Nymphalidae), sister species that recently diverged to mimic other Heliconius. This mimetic shift was a key step in their speciation, leading to pre- and postmating isolation. We identify 10 autosomal loci, half of which have major effects. At least eight appear to be homologous with genes known to control pattern differences within each species. Dominance has evolved under the influence of identifiable "modifier" loci rather than being a fixed characteristic of each locus. Epistasis is found at many levels: phenotypic interaction between specific pairs of genes, developmental canalization due to polygenic modifiers so that patterns are less sharply defined in hybrids, and overall fitness through ecological selection against nonmimetic hybrid genotypes. Most of the loci are clustered into two genomic regions or "supergenes," suggesting color pattern evolution is constrained by preexisting linked elements that may have arisen via tandem duplication rather than having been assembled by natural selection. Linkage, modifiers, and epistasis affect the strength of mimicry as a barrier to gene flow between these naturally hybridizing species and may permit introgression in genomic regions unlinked to those under disruptive selection. Müllerian mimics in Heliconius use different genetic architectures to achieve the same mimetic patterns, implying few developmental constraints. Therefore, although developmental and genomic constraints undoubtedly influence the evolutionary process, their effects are probably not strong in comparison with natural selection.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

A study of wing morphology and fluctuating asymmetry in interspecific hybrids between Drosophila buzzatii and D. koepferae

In this work we investigate the effect of interspecific hybridization on wing morphology using geometric morphometrics in the cactophilic sibling species D. buzzatii and D. koepferae. Wing morphology in F₁ hybrids exhibited an important degree of phenotypic plasticity and differs significantly from both parental species. However, the pattern of morphological variation between hybrids and the parental strains varied between wing size and wing shape, across rearing media, sexes, and crosses, suggesting a complex genetic architecture underlying divergence in wing morphology. Even though there was significant fluctuating asymmetry for both, wing size and shape in F₁ hybrids and both parental species, there was no evidence of an increased degree of fluctuating asymmetry in hybrids as compared to parental species. These results are interpreted in terms of developmental stability as a function of a balance between levels of heterozygosity and the disruption of coadaptation as an indirect consequence of genomic divergence.     

Ecologically relevant cryptic species in the highly polymorphic Amazonian butterfly Mechanitis mazaeus s.l. (Lepidoptera: Nymphalidae; Ithomiini)

The understanding of mimicry has relied on a strong biosystematic framework ever since early naturalists first recognized this textbook example of natural selection. We follow in this tradition, applying new biosystematics information to resolve problems in an especially difficult genus of tropical butterflies. Mechanitis species are important components of Neotropical mimetic communities. However, their colour pattern variability has presented challenges for systematists, and has made it difficult to study the very mimicry they so nicely illustrate. The South American Mechanitis mazaeus and relatives have remained particularly intractable. Recent systematists have recognized one highly polytypic species, whereas earlier work recognized the melanic Andean foothill races as a distinct species: Mechanitis messenoides. Recent molecular evidence suggests M. mazaeus and M. messenoides are genetically well differentiated, but evidence of morphological and ecological differences indicative of separate species was still lacking. Thus, it remains to be conclusively demonstrated whether this is an extreme case of a polymorphic mimetic species, or whether distinct co‐mimetic lineages are involved. Here we provide evidence that M. mazaeus and M. messenoides are ecologically distinct and identify consistent morphological differences in both adult and immature stages. These ecological and morphological differences are correlated with mitochondrial sequence data. In spite of some overlap in almost all traits, wing shape, adult colour pattern, and larval colour pattern differ between the two species, in addition to clutch size and larval host use in local sympatry. Although three well‐differentiated mitochondrial DNA (mtDNA) haplogroups were identified within these two species, one for M. mazaeus and two within M. messenoides, no morphological or ecological differences were found between two mtDNA haplogroups, both of which appear to belong to M. messenoides. We conclude that M. mazaeus and M. messenoides are distinct although highly polymorphic species, each with multiple sympatric co‐mimetic forms, and suggest that further work is needed to clarify the identity of other phenotypes and subspecies of Mechanitis. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 540–560.