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.     

A genetic linkage map of the mimetic butterfly Heliconius melpomene

Heliconius melpomene is a mimetic butterfly that exhibits great geographic variation in color pattern. We present here a genetic linkage map based on analysis of genetic markers in 73 individuals from a single F(2) family, offspring of a cross between H. m. cythera from western Ecuador and H. m. melpomene from French Guiana. A novel "three-step method" is described for the analysis of dominant markers in an F(2) cross, using outbred parental strains and taking advantage of the lack of crossing over in female Lepidoptera. This method is likely to prove useful for future mapping studies in outbred species with crossing over restricted to one sex, such as the Lepidoptera and Drosophila. The resulting linkage map has 21 linkage groups corresponding to the 21 chromosomes of H. melpomene and includes 219 AFLP markers, 23 microsatellites, 19 single-copy nuclear genes, and the color pattern switch genes Yb and Sb. The marker density is high, averaging >1/7 cM. The total map length is 1616 cM and the average chromosome length is 77 cM. The genome size of H. melpomene was estimated to be 292 Mb, giving a relationship of physical-to-map distance of 180 kb/cM. This map forms the basis for future comparative linkage analysis of color pattern evolution in Heliconius.     

Opsin cDNA sequences of a UV and green rhodopsin of the satyrine butterfly Bicyclus anynana

The cDNAs of an ultraviolet (UV) and long-wavelength (LW) (green) absorbing rhodopsin of the bush brown Bicyclus anynana were partially identified. The UV sequence, encoding 377 amino acids, is 76-79% identical to the UV sequences of the papilionids Papilio glaucus and Papilio xuthus and the moth Manduca sexta. A dendrogram derived from aligning the amino acid sequences reveals an equidistant position of Bicyclus between Papilio and Manduca. The sequence of the green opsin cDNA fragment, which encodes 242 amino acids, represents six of the seven transmembrane regions. At the amino acid level, this fragment is more than 80% identical to the corresponding LW opsin sequences of Dryas, Heliconius, Papilio (rhodopsin 2) and Manduca. Whereas three LW absorbing rhodopsins were identified in the papilionid butterflies, only one green opsin was found in B. anynana.     

Butterfly wing colours: scale beads make white pierid wings brighter

The wing-scale morphologies of the pierid butterflies Pieris rapae (small white) and Delias nigrina (common jezabel), and the heliconine Heliconius melpomene are compared and related to the wing-reflectance spectra. Light scattering at the wing scales determines the wing reflectance, but when the scales contain an absorbing pigment, reflectance is suppressed in the absorption wavelength range of the pigment. The reflectance of the white wing areas of P. rapae, where the scales are studded with beads, is considerably higher than that of the white wing areas of H. melpomene, which has scales lacking beads. The beads presumably cause the distinct matt-white colour of the wings of pierids and function to increase the reflectance amplitude. This will improve the visual discrimination between conspecific males and females.     

Butterfly wing pattern mutants: developmental heterochrony and co-ordinately regulated phenotypes

Butterfly wings are colored late in development, when pigments are synthesized in specialized wing scale cells in a fixed developmental succession. In this succession, colored pigments are deposited first and the remaining areas are later melanized black or brown. Here we studied the developmental changes underlying two wing pattern mutants, firstly melanic mutants of the swallowtail Papilio glaucus, in which the yellow background is turned black, and secondly a Spotty mutant of the satyrid Bicyclus anynana, which carries two additional eyespots. Despite the very different pattern changes in these two mutants, they are both associated with changes in rates of scale development and correspondingly, the final color pattern. In the melanic swallowtail, background scales originally destined to become yellow (normally developing early and synthesizing papiliochrome) show delayed development, fail to make papiliochrome, and subsequently melanize at the same time as scales in the wild-type black pattern. In the B. anynana eyespot, scale maturation begins with the central white focus, then progresses to the surrounding gold ring and later finishes with melanization of the black center. Mutants showing additional eyespots display accelerated rates of scale development (corresponding to new eyespots) in wing cells not normally occupied by eyespots. Thus by either delaying or accelerating rates of scale development, the final color, or position, of a wing pattern element can be changed. We propose that this heterochrony of scale development is a basic mechanism of color pattern formation on which developmental mutants act to change lepidopteran color patterns. Received: 20 April 2000 / Accepted: 19 July 2000     

Butterfly genomics eclosing

Technological and conceptual advances of the last decade have led to an explosion of genomic data and the emergence of new research avenues. Evolutionary and ecological functional genomics, with its focus on the genes that affect ecological success and adaptation in natural populations, benefits immensely from a phylogenetically widespread sampling of biological patterns and processes. Among those organisms outside established model systems, butterflies offer exceptional opportunities for multidisciplinary research on the processes generating and maintaining variation in ecologically relevant traits. Here we highlight research on wing color pattern variation in two groups of Nymphalid butterflies, the African species Bicyclus anynana (subfamily Satyrinae) and species of the South American genus Heliconius (subfamily Heliconiinae), which are emerging as important systems for studying the nature and origins of functional diversity. Growing genomic resources including genomic and cDNA libraries, dense genetic maps, high-density gene arrays, and genetic transformation techniques are extending current gene mapping and expression profiling analysis and enabling the next generation of research questions linking genes, development, form, and fitness. Efforts to develop such resources in Bicyclus and Heliconius underscore the general challenges facing the larger research community and highlight the need for a community-wide effort to extend ongoing functional genomic research on butterflies.     

ASSORTATIVE MATING PREFERENCES AMONG HYBRIDS OFFERS A ROUTE TO HYBRID SPECIATION

Homoploid speciation generates species without a change in chromosome number via introgressive hybridization and has been considered rare in animals. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has a color pattern that can be recreated by introgression of the H. melpomene red band into an H. cydno genetic background. Wild H. heurippa males show assortative mating based on color pattern and we here investigate the origin of this preference by studying first-generation backcross hybrids between H. melpomene and H. cydno that resemble H. heurippa . These hybrids show assortative mating preferences, showing a strong preference for their own color pattern over that of either parental species. This is consistent with a genetic basis to wing pattern preference and implies, first, that assortative mating preferences would facilitate the initial establishment of a homozygous hybrid color pattern by increasing the likelihood that early generation hybrids mate among themselves. Second, once established such a lineage would inherit assortative mating preferences that would lead to partial reproductive isolation from parental lineages.     

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.     

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.     

Adaptive introgression across species boundaries in Heliconius butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies.     

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.     

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.     

Concerted evolution and developmental integration in modular butterfly wing patterns

Developing organisms are thought to be modular in organization so that traits in different modules evolve independently whereas traits within a module change in a concerted manner. The eyespot pattern in Bicyclus anynana butterflies provides an ideal system where morphological modularity can be dissected and different levels of genetic integration analyzed. Several lines of evidence show that all eyespots in an individual butterfly are genetically integrated, suggesting that the whole pattern, rather than the separate eyespots, should be considered as a single character. However, despite the strong genetic correlations between the two eyespots on the dorsal forewing of B. anynana, there is great potential for independent changes. Here we use laboratory lines selected in different directions for the size of those eyespots to study correlated responses in the whole eyespot pattern. We show clear changes in eyespot size across all wing surfaces, which depend on eyespot position along the anterior-posterior axis. There are also changes in the number of extra eyespots and in eyespot color composition but no changes in eyespot position relative to wing margin. Our analysis of eyespot pattern modularity is discussed in the light of what is known about the cellular and genetic mechanisms of eyespot formation and the great potential for evolutionary diversification in butterfly wing patterns.     

A new mimetic species of Heliconius (Lepidoptera: Nymphalidae), from southeastern Colombia, revealed by cladistic analysis of mitochondrial DNA sequences

A new species of Heliconius and a new geographical race of Heliconius melpomene are described from the vicinity of Mocoa, Dpto. Putumayo, Colombia, based on molecular and morphological characters. The new species, H. tristero , is a close relative of H. cydno , a geographically differentiated species which lacks red coloration and engages in Müllerian mimicry with other blue and yellow Heliconius species in Central and northwestern South America. H. tristero has switched mimetic associations, instead mimicking the local, sympatric forms of two widespread mimetic species, H. erato and H. melpomene. This discovery provides evidence that the splinter species H. heurippa, H. tristero and H. timareta represent phenotypically divergent members of the H. cydno group that are endemic to successive river valleys on the eastern slope of the northern Andean Cordillera. The nominal taxon Heliconius amaryllis bellula Stichel, currently misapplied to both H. tristero and H. melpomene populations from the Mocoa region of Colombia, is considered here to represent a hybrid between H. heurippa and H. tristero. The Mocoa melpomene race is formally named Heliconius melpomene mocoa , new subspecies.     

Phylogenetic codivergence supports coevolution of mimetic Heliconius butterflies

The unpalatable and warning-patterned butterflies Heliconius erato and Heliconius melpomene provide the best studied example of mutualistic Müllerian mimicry, thought-but rarely demonstrated-to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of H. erato and H. melpomene, and this was initially hailed as one of the most striking known cases of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of H. erato and H. melpomene, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of H. erato and H. melpomene. Cophylogenetic historical reconstructions support repeated codivergence of mimetic populations, from the base of the sampled radiations. Pairwise distance correlation tests, based on our coalescent analyses plus recently published AFLP and wing colour pattern gene data, also suggest that the phylogenies of H. erato and H. melpomene show significant topological congruence. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of H. erato and H. melpomene occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. Our results suggest that differences in within-species genetic divergence are the result of a greater overall effective population size for H. erato relative to H. melpomene and do not imply incongruence in the timing of their phylogenetic radiations. Repeated codivergence between Müllerian co-mimics, predicted to exert mutual selection pressures, strongly suggests coevolution. Our results therefore support a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change.     

Patterns of pollen feeding and habitat preference among Heliconius species

Abstract 1. The ecological circumstances that precipitate speciation remain poorly understood. Here, a community of Heliconius butterflies in lowland Panama was studied to investigate patterns of pollen use, and more specifically the ecological changes associated with the recent divergence of Heliconius melpomene (Linnaeus) and H. cydno (Doubleday).
2. Considering the seven commonest Heliconius species in the community, 32 types of pollen or spore were encountered in pollen loads but only five pollen species were common. Systematic exploitation of pollen was therefore confined to a small proportion of the flowers visited.
3. Most of the variation in pollen load composition between individuals was explained by differences in collecting locality. The exception was Psiguria , which was used in all habitats by the melpomene / hecale clade far more than by the erato / sapho clade. This may suggest an ancestral switch within Heliconius towards increased reliance on Psiguria pollen.
4. Heliconius cydno and H. melpomene differed significantly in pollen load composition for three of the five most commonly collected pollen species. This is most probably explained by differences in habitat preference; H. melpomene and its co-mimic H. erato are found in open habitat while H. cydno and its co-mimic H. sapho are found in closed-canopy forest.
5. As melpomene and cydno are known to hybridise occasionally, such differences in adult microhabitat contribute to pre-mating isolation. Habitat divergence between H. cydno and H. melpomene , which is associated with changes in mimicry, must have played a role in their recent speciation.     

Phylogeny of Bicyclus (Lepidoptera: Nymphalidae) inferred from COI, COII, and EF-1alpha gene sequences

Despite the fact that Bicyclus anynana has become an important model species for wing-pattern developmental biology and studies of phenotypic plasticity, little is known of the evolutionary history of the genus Bicyclus and the position of B. anynana. Understanding the evolution of development as well as the evolution of plasticity can be attempted in this species-rich genus that displays a large range of wing patterns with variable degrees of phenotypic responses to the environment. A context to guide extrapolations from population genetic studies within B. anynana to those between closely related species has been long overdue. A phylogeny of 54 of the 80 known Bicyclus species is presented based on the combined 3000-bp sequences of two mitochondrial genes, cytochrome oxidase I and II, and the nuclear gene, elongation factor 1alpha. A series of tree topologies, constructed either from the individual genes or from the combined data, using heuristic searches under a variety of weighting schemes were compared under the best maximum-likelihood models fitted for each gene separately. The most likely tree topology to have generated the three data sets was found to be a tree resulting from a combined MP analysis with equal weights. Most phylogenetic signal for the analysis comes from silent substitutions at the third position, and despite the faster rate of evolution and higher levels of homoplasy of the mitochondrial genes relative to the nuclear gene, the latter does not show substantially stronger support for basal clades. Finally, moving branches from the chosen tree topology to other positions on the tree so as to comply better with a previous morphological study did not significantly affect tree length.     

Color pattern formation on the wing of a butterfly Pieris rapae. 2. Color determination and scale development

It has been shown that microcautery on the prospective apical black region of the early pupal forewing of a butterfly, Pieris rapae , causes alteration of the scale color on the adult wing and a delay in histogenesis of the pupal wing. From these results, it has been assumed that the developmental delay of scale cells in the pupal wing alters their developmental fate and the hypothesis that different color fates of scales are determined by differences in the developmental timetables between scale cells is proposed. In this study, we attempted to find the developmental timetables of individual scales expressing specific color to test this hypothesis. It was found that the holes on the upper surface of a scale become larger as they develop and the hole sizes of scales in the white region are always larger than in the black region on the same wings either during pupal period or after eclosion. This suggests that the scale hole size is a good index that reflects developmental rate of the scale and a difference in the hole size between adult scales is attributed to a difference in the developmental timetables when their ancestral scale precursor cells were in the pupal period. A comparison of the hole sizes between adult scales in different color regions suggested that normal white scales were in a more advanced state than were the black ones but white scales induced by microcautery were in a less advanced state than black ones on the same wing. This supports our hypothesis.