- Research article
- Open Access
Two sisters in the same dress: Heliconius cryptic species
© Giraldo et al; licensee BioMed Central Ltd. 2008
- Received: 06 August 2008
- Accepted: 28 November 2008
- Published: 28 November 2008
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.
Here we present phenotypic, behavioral and genetic evidence for the coexistence of two sympatric 'cryptic' species near Florencia in the eastern Andes of Colombia that share the same orange rayed colour pattern. These represent H. melpomene malleti and a novel taxon in the H. cydno group, here designated as novel race of Heliconius timareta, Heliconius timareta florencia. No-choice mating experiments show that these sympatric forms have strong assortative mating (≈96%) despite great similarity in colour pattern, implying enhanced divergence in pheromonal signals.
We hypothesize that these species might have resulted from recent convergence in colour pattern, perhaps facilitated by hybrid introgression of wing pattern genes.
- Reproductive Isolation
- Cryptic Species
- Colour Pattern
- Discal Cell
- Triose Phosphate Isomerase
Ecological selection is known to play an important role in speciation [1, 2]. Where ecological traits under divergent selection also affect the species recognition system, speciation can be rapid [3, 4]. Young species pairs provide evidence for this, such as stickleback morphs found in glacial lakes in which body size diverges through adaptation to different niches and is correlated with mate choice [5, 6]. Similarly, in Darwin's finches adaptive variation in beak size also promotes reproductive isolation through pleiotropic effects on song [7, 8]. In many butterflies colour pattern similarly plays a role in mate recognition [9, 10] and is also under ecological selection for signaling to predators, crypsis and thermoregulation .
In Heliconius butterflies shifts in colour pattern have been shown to play a major role in speciation [12, 13]. Closely related species and sub-species typically differ in colour pattern  and are adapted to local Müllerian mimicry rings in which distasteful species converge on a common pattern [15, 16]. The pattern differences between related forms lead to strong assortative mating [12, 13]. Frequency dependent selection contributes to maintaining the mimetic patterns and also causes post-mating isolation [17, 18]. Thus, between closely related species such as H. cydno and H. melpomene, rare hybrid individuals are likely to be strongly selected against, as their pattern will not be recognized by predators [12, 15]. These species overlap extensively across Central America and the Andes [19, 20]. Through their geographical range, H. melpomene mimics H. erato with a black background and red, yellow and orange marks, whereas H. cydno mimics species of the H. sapho group with a black-blue background and white and yellow marks [12, 21]. Although H. melpomene and H. cydno occasionally hybridize in nature (less than 0.1%), reproductive isolation is strong, with ecological isolation and colour pattern associated mate choice playing the major role [12, 22, 23].
The phylogeny of Heliconius supports a key role for pattern change in speciation, with almost all sister species differing in colour pattern . Exceptions such as H. sara and H. leucadia are far more genetically divergent than species in the H. melpomene clade, implying that speciation was relatively ancient . It was therefore a surprise when a putative 'cryptic' species was identified in the H. cydno/melpomene species complex, implying either mimetic convergence between closely related species, or speciation without pattern change . H. tristero was described as a new species from the southeastern Andes of Colombia with a cydno-like mtDNA haplotype, but a black, red and yellow "postman" pattern similar to sympatric H. melpomene mocoa . This species was initially met with scepticism, mainly due to the fact that only two individuals of H. tristero were collected, combined with the likelihood of mtDNA introgression between H. cydno and H. melpomene . It was therefore suggested that the H. tristero specimens most likely represent rare hybrids between H. cydno and H. melpomene.
Sampling and species assignment
A total of 196 adult butterflies were collected between 2002–2006 in Sucre (01°48'12¨N 75°39'19¨W) near Florencia, Colombia. For comparison of larval and adult morphology, an additional 118 butterflies representing H. cydno weymeri (N = 21), H. c. cordula (N = 19), H. c. zelinde (N = 12), H. melpomene vulcanus (N = 17), H. m. melpomene (N = 26) and H. m. mocoa (N = 23) were collected from various sites in Colombia. Some of the individuals from Florencia were maintained alive in separate insectaries at La Vega, Cundinamarca for mate and host choice experiments. Wings were removed from the remaining specimens and bodies preserved in DMSO or Ethanol (96%) for subsequent phenotypic (N = 131) and molecular analysis (N = 142) and stored in the permanent collection of the Instituto de Genética of Universidad de los Andes in Colombia (the collection is now part of the recently created Natural History Museum ANDES). DNA extractions were made from one-third of the thorax of each individual using the DNeasy tissue Kit (QIAGEN), following the manufacturers' protocol. Digital images were obtained by scanning the wings and quantitative measurements taken using TpsDig image software . One-way Analyses of Variance (ANOVA) on phenotypic measurements were performed with SPSS 11.0.4 software .
No-choice mating experiments
Once diagnostic characters for the two morphotypes had been established we carried out mate choice experiments to investigate reproductive isolation. No-choice mating trials were performed in 2 m × 2 m × 2 m insectaries in La Vega, Cundinamarca between June 2004 and October 2005. These experiments are a simulation of a natural situation where males encounter females singly, and estimate the reluctance of both sexes to mate inter-specifically. A virgin female (one to three days old) of each morphotype was presented to ten mature males (more than 10 days old) in a single insectary for two days. Successful and failed matings were recorded every 30 minutes between 6 am to 2 pm. Males were used only once. In order to detect matings that were unobserved, females were checked for the presence of a spermatophore in their reproductive tract. A binomial mating probability Pixj and 95% confidence intervals were obtained for each combination of i-type female and j-type male using Maximum Likelihood as previously described .
Host plant choice and larval morphology
Females of both Florencia populations (H. m. malleti and H. cydno cognate) and other populations of the two species were kept in individual insectaries with known host plants, Passiflora edulis, P. maliformis, P. ligularis, P. arborea, P. quadrangularis, P. oerstedii [30, 31], and as controls two species used by H. erato, P. suberosa, and P. rubra . The number of eggs laid per plant by each female was recorded twice a week. Offspring from the same females were used for analysis of larval morphology. H. melpomene and H. cydno larvae are distinguishable by cephalic colour  which is pale yellow in H. melpomene and orange in H. cydno. Pictures of larvae raised from wild females were taken under similar light conditions with a colour standard. Between four and seven larvae were analysed from each of 21 females collected in Florencia and compared to those from other localities. Pictures were processed using Scion Image (Scion Corporation, Frederick, MD, USA) and four RGB indexes calculated, R' = r/(r+g+b), G' = g/(r+g+b), B' = b/(r+g+b) and LM = R'-G' . Statistical significance of differences in indices was tested with one-way variance analysis using SPSS 11.0.4 software .
We sequenced a region of nuclear DNA spanning the 3' end of the Z-linked Triose phosphate isomerase (Tpi) gene, from thirteen female specimens of H. m. malleti and seven H. cydno cognate from Sucre Florencia, Caquetá and nine specimens of H. timareta from Ecuador (additional file 1). The Polymerase Chain Reaction (PCR) was performed using primers and conditions outlined elsewhere . The PCR products were electrophoretically separated on 1.5% low melting point agarose (Invitrogen), and the bands were cut from the gel and dissolved in gelase (Invitrogen). Clean PCR products were sequenced using the DNA sequence Kit (Big Dye 3.1, PE Applied Biosystems), in an ABI Prism 3100 Genetic Analyzer (PE Applied Biosystems).
We also sequenced fragments of two mitochondrial genes (COI-COII) from seventeen H.m. malleti, seven cydno-like and seven H. timareta individuals. In addition to sequences obtained here, we also included CO and Tpi sequences from GenBank (additional file 1). PAUP* v4.0b10  was used to search for a maximum parsimony tree, using a heuristic search with TBR branch swapping; bootstrap values were calculated with 5000 replicates using the same search conditions.
MrModeltest v2.2  was used to determine the most appropriate model of nucleotide substitution based on hierarchical likelihood ratio tests. For the mtDNA COI/COII data set MrModeltest identified the GTR+I+G model [37, 38], and for Tpi, the GTR+G model [37, 38]. Bayesian phylogenetic analyses were performed with MrBayes v3.1  following the analytical recommendations of the authors . The base frequency parameters determined by MrModeltest were used for analysis in MrBayes v3.1, with remaining parameters estimated using the GTR+I+G model for COI/COII and GTR+G for Tpi. Four differentially heated Markov chains were initiated from random trees, run for 106 generations and sampled every 100 cycles. Likelihood values were plotted against number of generations to determine the points at which stationary was reached. All trees sampled before these points were discarded and the remaining tree samples were used to generate a 50% majority rule consensus tree (n = 12101 for COI/COII and n = 9981 for Tpi). The posterior probability of each clade is provided by the percentage of trees identifying the clade [39, 40].
A total of 13 microsatellite loci were amplified and scored as described previously . Hardy-Weinberg equilibrium and linkage disequilibrium and their significance were tested for at each locus using Arlequin v2000 . Two Bayesian model-based clustering algorithms implemented in the programs STRUCTURE 2.2  and BAPS 4  were used to test the hypothesis that the Florencia population consisted of two clusters. We determined the number of ancestral clusters, K, using an ad hoc statistic ΔK based on the rate of change in the log probability of data for K between 1 and 5 in multiple runs . Each run consisted of 106 iterations, after a burning period of 104 iterations. To use these programs, Hardy-Weinberg and linkage equilibrium are assumed, and the software differentiates mixed populations on the basis of allele frequencies at each locus. Finally, population differentiation (FST) and genetic distances (DA) among the two Florencia types were calculated with Arlequin v2000 .
Microsatellite Data Analyses
Individuals from Florencia were screened for variation at thirteen microsatellite loci. The H. cydno cognate was less variable, with a mean observed heterozygosity of 0.39 as compared to 0.66 in H. m. malleti. All loci were in Hardy-Weinberg equilibrium within the two morphotypes, except Hm19, which showed significant heterozygote deficit for both forms, most likely due to null alleles which are known to be present at this locus . All subsequent analyses are presented with this locus removed. The microsatellite FST value between the two morphotypes in Florencia was large and significant (FST = 0.24226, P < 0.0001).
Five individuals had clearly intermediate assignment probabilities, suggesting possible hybrid genotypes. In support of this hypothesis, these individuals also had some evidence for colour pattern introgression, with either somewhat intermediate red line phenotypes (Figure 2) or melanic scales in the forewing band. These individuals represented around 3.5% of the total sample (5/142) with three of the five being consistent with an F1 genotype (Figure 5). Five additional individuals in the H. m. malleti cluster showed posterior probabilities less than 0.95 suggesting that these might represent backcross hybrids.
No-choice mate experiments
Results from the no-choice mating trials
H. cydno- type
H. cydno- type
Host plant preference
Nineteen females from Florencia (10 H. m. malleti and 9 H. cydno-type) were assessed for host preference. Females assigned as H. m. malleti used primarily two plants (471 eggs in total on P. oerstedii, 71%, and P. ligularis, 28%), while the H. cydno cognate females oviposited on many species (729 eggs in total on P. edulis 48%, P. ligularis 15%, P. oerstedii 22%, and P. quadrangularis, P. arborea and P. maliformis 14%).
The colour patterns of H. melpomene and H. cydno are traits under strong ecological selection that also contribute to speciation [12, 23, 26, 49, 50]. Divergence in mimetic pattern contributes to reproductive isolation by assortative mating, due to the use of colour as a mate recognition signal and by frequency dependent mimicry selection against rare colour pattern hybrids [12, 23, 26, 29]. The fact that such a major role for mimicry in the reproductive isolation of currently hybridizing species has been clearly demonstrated has led to an expectation that mimicry is unlikely between closely related species with incomplete reproductive isolation . The results of the present study clearly show that this is not always the case. The butterflies collected in Florencia represent two mimetic species in the H. melpomene and H. cydno clades respectively. Despite the extreme phenotypic similarity between the butterflies studied here, the data indicate clear concordance between nuclear, mtDNA and phenotypic markers in assigning over 90% of the individuals sampled to one or the other morphotype. These represent a population of the widespread race H. m. malleti, and a novel entity related to H. cydno. Nonetheless, there is clear evidence for ongoing interspecific hybridization with around 3.5% of individuals sampled representing clearly identifiable hybrids.
The strong concordance between markers and relative scarcity of hybrids implies strong reproductive isolation. The mating experiments described here demonstrate strong assortative mating, and host choice experiments further imply some degree of ecological isolation. The strength of assortative mating in our experiments is actually greater than that between the phenotypically very divergent H. melpomene melpomene and H. cydno cordula (≈ 82%, Mavárez et al. 2006). It therefore seems likely that there has been increased divergence in mating signals apart from colour – most likely pheromonal – to allow the Florencia species to coexist. As colour is used as a cue in mate finding, we would predict that males are likely to be attracted to the pattern of the "wrong" species, but that this is compensated for by divergence in other mating signals.
H. melpomene and H. cydno are known to differ in habitat preference and host use [30, 51]. Our host choice data imply that this ecological difference is maintained between the Florencia species, with the H. cydno cognate more of a host generalist similar to other populations of H. cydno. It seems likely that this corresponds to a preference for forest habitats, as is the case for other H. cydno populations, leading to ecological isolation. Although crosses between the Florencia forms have not been carried out, there is preliminary evidence that the H. cydno cognate morphotype is compatible with other populations of H. cydno, while the Florencia H. m. malleti shows hybrid sterility with H. cydno (Giraldo and Linares, unpub). It seems likely that the Florencia species show similar female hybrid sterility as compared to other sympatric H. melpomene and H. cydno populations.
In order to clarify the role of colour pattern in speciation we need to consider the order of divergence in different factors. Divergence might have occurred initially in factors other than colour pattern, such as habitat preference and pheromonal mating cues. Divergence in mimicry then occurred subsequent to speciation in most populations of the H. cydno group, except those in the eastern Andes such as that studied here. Alternatively, speciation might have been initially triggered by divergence in colour pattern with mimetic convergence a derived state acquired subsequent to speciation in the Florencia region. The first explanation is more parsimonious with respect to colour pattern, but may actually be less likely. The derived position of the H. cydno cognate form in the mtDNA phylogeny suggests that it has evolved from a more H. cydno like ancestor (Figure 3). Furthermore, the observation that intraspecific races of H. melpomene show both divergent colour patterns and strong assortative mating implies that mimicry and associated mate preferences are the first steps in divergence in this group .
Under the second scenario, an initial divergence in colour pattern associated with adaptation to different Müllerian mimicry rings became associated with further changes in ecology and hybrid inviability, eventually leading to speciation . Subsequently, strong mimetic selection in the sympatric Florencia population has led to convergence of H. cydno onto the H. melpomene pattern. This might have occurred through adaptive introgression of colour pattern genes between the species [21, 23, 26, 52–54]. It seems plausible that the lack of suitable mimicry models in the H. sapho group has led this east Andean population to secondarily converge on the H. melpomene pattern.
The new H. cydno cognate from Florencia resembles H. timareta, a known relative of H. cydno from eastern Ecuador. H. timareta is polymorphic with one orange rayed form that is very similar to the taxon from Florencia (Figure 6). Furthermore, a recently described subspecies, H. t. timoratus Lamas (1998) from Perú is a near perfect mimic of H. m. malleti, although the two have so far not been collected in sympatry (Figure 6). Due to the colour pattern similarity with H. t. timoratus of Perú and the east Andean distribution of both forms, here we propose the H. cydno cognate from Florencia as a new northern race of H. timareta, which forms a mimicry ring with H. m. malleti and H. e. lativitta in the southern foothills of Colombia. We propose to name this form after Florencia, the town where it was discovered. From a geographical point of view (Figure 6) the H. cydno clade in the eastern Andes is represented by, from north to south, H. cydno cordula, H. heurippa, H. timareta florencia and H. tristero in Colombia; H. timareta timareta in Ecuador and H. timareta timoratus and another new H. cydno-like taxon in Perú (Mavarez et al., in prep.).
H. timareta florencia is anatomically similar to H. timareta and H. cydno races, with a FW area of 5,26 to 6,64 cm2 (additional file 2). The FW has two principal pattern elements in a black background, 1. An irregular yellow postmedial band extending proximodistally from distal end of discal cell to R2–R3 fork, and laterally from subcostal to CU1a, and 2. A red "Dennis" element extending from the basal end of the discal cell to the CU1b-discal cell fork, which is generally shorter than the H. m. malleti FW element (the approximate basal fourth of the FW). The ventral FW is similar with a slight reduction in the yellow band area and the distinguishing red line with length 2,65 to 5,92 mm (Figure 1b and Figure 2). The HW also has a black background colour and an androconial distribution as in H. timareta. A reddish dennis-ray element  is similar in form to the dennis-ray in H. m. malleti, with a narrow red longitudinal bar D (Dennis) and six red radiate marks on the discal part of the upper side . The ventral HW underside is similar with narrower red areas. Males and females are phenotypically similar. H. t. florencia is registered in ZooBank with the unique digital identifier LSID: 49656F41-E817-4FD6-B28D-6AA362CC5268. In order to comply with the International Code of Zoological Nomenclature, paper copies of this electronic article have been deposited in the following libraries: The Entomology Library of the Natural History Museum, London, UK; The Balfour Library, Department of Zoology, University of Cambridge, UK; The Genetics Library, Department of Genetics, University of Cambridge, UK; Biblioteca Universidad Nacional de Colombia, Hemeroteca Nacional, Sede Bogotá; Biblioteca Universidad de los Andes, Bogotá, Colombia.
Type material: The holotype male specimen is deposited in the permanent collection of the Natural History Museum ANDES of Universidad de los Andes in Colombia. The data label reads: Andes-E-11517, Colombia, Caquetá, Florencia, Quebrada las Doraditas, 01/09/2008. The red holotype label reads: Holotipo, Andes-E-11517, Heliconius timareta florencia.
Obviously the study of cryptic species has a long history and, especially with the advent of modern molecular techniques many previously undescribed taxa have been detected . However, cryptic species often use sensory modalities that humans do not readily perceive, such as pheromones [57, 58], toxicity resistance  and imperceptible song differences . The surprising aspect of this study is therefore the discovery of cryptic species in well-studied taxa where speciation is commonly triggered by bright visual signals. It is increasingly becoming clear that tropical biodiversity is severely underestimated and that combining morphological and DNA sequence analysis is a powerful tool for species discovery. Finally, our mate and host choice experiments suggest that the sympatric species in Florencia combine the same ecological differences known to occur in other parts of the range of H. cydno and H. melpomene with enhanced pheromonal signals to compensate for the lack of colour pattern signals. Our results therefore highlight the fact that speciation is the combined result of divergence along multiple phenotypic axes.
We would like to thank Maribel González and Carlos Arias at the Smithsonian Tropical Research Institute for laboratory help with microsatellites and mtDNA amplification and Fernando Alda for help in microsatellite scoring and analysis. We also thank Mathieu Joron and Jesus Mavarez for ongoing discussion of cryptic species in the eastern Andes and for the original suggestion that these might represent sub-species of H. timareta. James Mallet, Gerardo Lamas and Jean François Le Crom for photographs of H. timareta from Ecuador and Perú and H. tristero. Ministerio del Ambiente of Ecuador for collecting permitions. This work was carried out with grants from the Smithsonian Tropical Research Institute (NG), Instituto Colombiano para el Desarrollo de la Ciencia y la Tecnología Francisco José de Caldas COLCIENCIAS 7155-CO (CS), Biotechnology and Biological Sciences Research Council and the Royal Society (CDJ), Banco de la República (ML and NG) and private donations (ML) from Continautos S.A., Proficol El Carmen S.A., Didacol S.A. and F. Arango, Colombia.
- Schluter D: The Ecology of Adaptive Radiation. 2000, New York: Oxford University PressGoogle Scholar
- Coyne JA, Orr HA: Speciation. 2004, Sunderland, Mass.: Sinauer AssociatesGoogle Scholar
- Rundle HD, Nagel L, Wenrick Boughman J, Schluter D: Natural selection and parallel speciation in sympatric sticklebacks. Science. 2000, 287: 306-308. 10.1126/science.287.5451.306.View ArticlePubMedGoogle Scholar
- Feder J, Xie X, Rull J, Velez S, Forbes A, Leung B, Dambroski H, Filchak KE, Aluja M: Mayr, Dobzhansky, and Bush: a golden braid of biogeography, inversions, differential introgression, clines and sympatric speciation?. Proc Natl Acad Sci USA. 2005, 102: 6573-6580. 10.1073/pnas.0502099102.PubMed CentralView ArticlePubMedGoogle Scholar
- McKinnon JS, Rundle HD: Speciation in nature: the threespine stickleback model systems. Trends Ecol Evol. 2002, 17: 480-488. 10.1016/S0169-5347(02)02579-X.View ArticleGoogle Scholar
- Nagel L, Schluter D: Body size, natural selection, and speciation in sticklebacks. Evolution. 1998, 52: 209-218. 10.2307/2410936.View ArticleGoogle Scholar
- Huber SK, Podos J: Beak morphology and song features covary in a population of Darwin's finches (Geospiza fortis). Biol J Linn Soc. 2006, 88: 489-498. 10.1111/j.1095-8312.2006.00638.x.View ArticleGoogle Scholar
- Huber SK, De León LF, Hendry AP, Bermingham E: Reproductive isolation of sympatric morphs in a population of Darwin's finches. Proc R Soc Lond B. 2007, 274: 1709-1714. 10.1098/rspb.2007.0224.View ArticleGoogle Scholar
- Fordyce JA, Nice CC, Forister ML, Shapiro AM: The significance of wing pattern diversity in the Lycaenidae: mate discrimination by two recently diverged species. J Evol Biol. 2002, 15: 871-879. 10.1046/j.1420-9101.2002.00432.x.View ArticleGoogle Scholar
- Kemp DJ, Rutowski RL: Condition dependence, quantitative genetics, and the potential signal content of iridescent ultraviolet butterfly coloration. Evolution. 2007, 61: 168-183. 10.1111/j.1558-5646.2007.00014.x.View ArticlePubMedGoogle Scholar
- Nijhout HF: The Development and Evolution of Butterfly Wing Patterns. 1991, Washington, DC: Smithsonian Institution PressGoogle Scholar
- Jiggins CD, Naisbit RE, Coe RL, Mallet J: Reproductive isolation caused by colour pattern mimicry. Nature. 2001, 411: 302-305. 10.1038/35077075.View ArticlePubMedGoogle Scholar
- Jiggins CD, Estrada C, Rodrigues A: Mimicry and the evolution of pre-mating isolation in Heliconius melpomene. J Evol Biol. 2004, 17: 680-691. 10.1111/j.1420-9101.2004.00675.x.View ArticlePubMedGoogle Scholar
- Brown KS: The biology of Heliconius and related genera. Ann Rev Entomol. 1981, 26: 427-456. 10.1146/annurev.en.26.010181.002235.View ArticleGoogle Scholar
- Mallet J, McMillan WO, Jiggins CD: Mimicry and warning colour at the boundary between races and species. Endless Forms: Species and Speciation. Edited by: Howard DJ, Berlocher SH. 1998, New York: Oxford University Press, 390-403.Google Scholar
- Mallet J, Gilbert LE: Why are there so many mimicry rings? Correlations between habitat, behaviour and mimicry in Heliconius butterflies. Biol J Linn Soc. 1995, 55: 159-180.Google Scholar
- Mallet J, Joron M: Evolution of diversity in warning colour and mimicry: Polymorphisms, shifting balance and speciation. Annu Rev Ecol Syst. 1999, 30: 201-233. 10.1146/annurev.ecolsys.30.1.201.View ArticleGoogle Scholar
- Kapan D: Three-butterfly system provides a field test of Müllerian mimicry. Nature. 2001, 409: 338-340. 10.1038/35053066.View ArticlePubMedGoogle Scholar
- Brower AVZ: Parallel race formation and the evolution of mimicry in Heliconius butterflies: a phylogenetic hypothesis from mitochondrial DNA sequences. Evolution. 1996, 50: 195-221. 10.2307/2410794.View ArticleGoogle Scholar
- Brown KS: Ecologia Geográfica e Evoluçâo nas Florestas Neotropicais. 1979, Campinas, Brazil: Universidade Estadual de CampinasGoogle Scholar
- Linares M: Origin of Neotropical mimetic biodiversity from a three-way hybrid zone of Heliconius butterflies. Proceedings of the International Symposium on Biodiversity and Systematics in Tropical Ecosystems, Bonn 1994. Edited by: Ulrich H. 1997, Bonn: Zoologishunginstitut und Museum Alexander KoenigGoogle Scholar
- Jiggins CD, Emelianov I, Mallet J: Assortative mating and speciation as pleiotropic effects of ecological adaptation: examples in moths and butterflies. Insect Evolutionary Ecology. Edited by: Fellowes M, Holloway G, Rolff J. 2005, London, UK.: Royal Entomological Society, 451-473.Google Scholar
- Jiggins CD: Ecological speciation in mimetic butterflies. Bioscience. 2008, 58: 541-548. 10.1641/B580610.View ArticleGoogle Scholar
- Beltrán M, Jiggins CD, Brower AVZ, Bermingham E, Mallet J: Do pollen feeding and pupal-mating have a single origin in Heliconius ? Inferences from multilocus sequence data. Biol J Linn Soc. 2007, 92: 221-239. 10.1111/j.1095-8312.2007.00830.x.View ArticleGoogle Scholar
- Brower AVZ: A new species of Heliconius (Lepidoptera: Nymphalidae), from Sourtheastern Colombia, revealed by analysis or mitochondrial DNA sequences. Zool J Linn Soc. 116: 317-332. 10.1111/j.1096-3642.1996.tb00126.x.Google Scholar
- Mallet J: Rapid speciation, hibridization and adaptive radiation in the Heliconius melpomene group. Speciation and Ecology. Edited by: Butlin RK, Bridle JR, Schluter D. 2008, Sheffield, UKGoogle Scholar
- Rohlf FJ: TPSdig: digitize landmarks from image files, scanner, or video. Stony Brook. 2001, New York: Departament of Evolutionary Biology, University of New YorkGoogle Scholar
- SPSS Inc: SPSS Base 10.0 for Windows User's guide. SPSS Inc. 1999, Chicago, IL: SPSS Inc, 11.0.4Google Scholar
- Naisbit RE, Jiggins CD, Mallet J: Disruptive sexual selection against hybrids contributes to speciation between Heliconius cydno and Heliconius melpomene. Proc R Soc B. 2001, 268: 1849-1854. 10.1098/rspb.2001.1753.PubMed CentralView ArticlePubMedGoogle Scholar
- Smiley JT: Plant chemistry and the evolution of host specificity: new evidence from Heliconius and Passiflora. Science. 1978, 201: 745-747. 10.1126/science.201.4357.745.View ArticlePubMedGoogle Scholar
- Benson WW, Brown KS, Gilbert LE: Coevolution of plants and herbivores: passion flower butterflies. Evolution. 1975, 29: 659-680. 10.2307/2407076.View ArticleGoogle Scholar
- Jiggins CD, McMillan WO, Mallet J: Host plant adaptation has not played a role in the recent speciation of Heliconius erato and Heliconius himera (Lepidoptera; Nymphalidae). Ecol Entomol. 1997, 22: 361-365. 10.1046/j.1365-2311.1997.00067.x.View ArticleGoogle Scholar
- Endler JA: On the measurement and classification of colour in studies of animal colour patterns. Biol J Linn Soc. 1990, 41: 315-352. 10.1111/j.1095-8312.1990.tb00839.x.View ArticleGoogle Scholar
- Beltrán M, Jiggins CD, Bull V, McMillan WO, Bermingham E, Mallet J: Phylogenetic discordance at the species boundary: gene genealogies in Heliconius butterflies. Mol Biol Evol. 2002, 19: 2176-2190.View ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). 2000, Sunderland, Massachusetts: Sinauer Associates, 4Google Scholar
- Nylander JAA: MrModeltest 2.0. Program distributed by the author. 2004, Evolutionary Biology Centre, Uppsala UniversityGoogle Scholar
- Rodríguez FJ, Oliver JL, Marín A, Medina JR: The general stochastic model of nucleotide substitution. J Theor Biol. 1990, 142: 485-501. 10.1016/S0022-5193(05)80104-3.View ArticlePubMedGoogle Scholar
- Yang Z: Estimating the pattern of the nucleotide sustitution. J mol Evol. 1994, 39: 105-111.PubMedGoogle Scholar
- Huelsenbeck JP, Ronquist FR: MrBayes: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP, Mark P: MrBayes 3.1 Manual. 2005, 1-66.Google Scholar
- Mavárez J, Gonzalez MA: A set of microsatellite loci for Heliconius melpomene and close relatives. Mol Ecol Notes. 2006, 6: 20-23. 10.1111/j.1471-8286.2005.01115.x.View ArticleGoogle Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin Ver.2000: A software for population genetics data analysis. Genetics and Biometry Laboratory. 2000, University of Geneva. SwitzerlandGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of Population Structure Using Multilocus Genotype Data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Corander J, Waldmann P, Sillanpää MJ: Bayesian analysis of genetic differentiation between populations. Genetics. 2003, 163: 367-374.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Linares M: The ghost of mimicry past: laboratory reconstruction of an extinct butterfly "race". Heredity. 1997, 78: 628-635. 10.1038/hdy.1997.102.View ArticleGoogle Scholar
- Linares M: The genetics of mimetic coloration in the butterfly Heliconius cydno weymeri. J Hered. 1996, 87: 142-149.View ArticleGoogle Scholar
- Jiggins CD, Mavárez J, Beltrán M, MacMillan WO, Johnston JS, Bermingham E: A genetic map of the mimetic butterfly, Heliconius melpomene. Genetics. 2005, 171: 557-570. 10.1534/genetics.104.034686.PubMed CentralView ArticlePubMedGoogle Scholar
- Kronforst MR, Kapan DD, Gilbert LE: Parallel genetic architecture of parallel adaptive radiations in mimetic Heliconius butterflies. Genetics. 2006, 174: 535-539. 10.1534/genetics.106.059527.PubMed CentralView ArticlePubMedGoogle Scholar
- Joron M, Papa R, Beltrán M, Chamberlain N, Mavárez J, Baxter S, Abanto M, Bermingham E, Humphray SJ, Rogers J, et al: A conserved supergene locus controls colour pattern diversity in Heliconius butterflies. PLoS BioL. 2006, 4: e303-10.1371/journal.pbio.0040303.PubMed CentralView ArticlePubMedGoogle Scholar
- Estrada C, Jiggins CD: Patterns of pollen feeding and habitat preference among Heliconius species. Ecol Entomol. 2002, 27: 448-456. 10.1046/j.1365-2311.2002.00434.x.View ArticleGoogle Scholar
- Linares M: Adaptive microevolution through hybridization and biotic destruction in the neotropics. PhD Dissertation. 1989, University of TexasGoogle Scholar
- Gilbert LE: Adaptive novelty through introgression in Heliconius wing patterns: evidence for shared genetic "tool box" from synthetic hybrid zones and a theory of diversification. Ecology and Evolution Taking Flight: Butterflies as Model Systems. Edited by: Boggs CL, Watt WB, Ehrlich PR. 2003, Chicago, University of Chicago Press, 281-318.Google Scholar
- Mallet J, Beltrán M, Neukirchen W, Linares M: Natural hybridization in Heliconiine butterflies: the species boundary is a continuum. BMC Evol Biol. 2007, 7: 28-10.1186/1471-2148-7-28.PubMed CentralView ArticlePubMedGoogle Scholar
- Emsley MG: The geographical distribution of the color-pattern components of Heliconius erato and Heliconius melpomene with genetical evidence for the systematic relationship between the two species. Zoologica, New York. 1964, 49: 245-286.Google Scholar
- Bickford D, Lohman D, Sodhi NS, Ng PKL, Meier R, Winker K, Ingram KK, Das I: Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2006, 22: 148-155. 10.1016/j.tree.2006.11.004.View ArticlePubMedGoogle Scholar
- Kozlov MV, Zhu J, Philipp P, Francke W, Zvereva EL, Hansson BS, Löfstedt C: Pheromone specificity in Eriocrania semipurpurella (Stephens) and E. sangii (Wood) (Lepidoptera: Eriocraniidae) based on chirality of semiochemicals. J Chem Ecol. 1996, 22: 431-454. 10.1007/BF02033647.View ArticlePubMedGoogle Scholar
- Byers JR, Struble DL: Identification of sex pheromones of two sibling species in dingy cutworm complex Feltia jaculifera (Lepidoptera: Noctuidae). J Chem Ecol. 1990, 16: 2981-2992. 10.1007/BF00979489.View ArticlePubMedGoogle Scholar
- Sturmbauer C, Opadiya GB, Niederstaetter H, Reidmann A, Dallinger R: Mitochondrial DNA reveals cryptic oligochaete species differing in cadmium resistance. Mol Biol Evol. 1999, 16: 967-974.View ArticlePubMedGoogle Scholar
- Henry CS: Singing and cryptic speciation in insects. Trends Ecol Evol. 1994, 9: 388-392. 10.1016/0169-5347(94)90061-2.View ArticlePubMedGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1660. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.