Genetic differentiation without mimicry shift in a pair of hybridizing Heliconius species (Lepidoptera: Nymphalidae)

Butterflies in the genus Heliconius have undergone rapid adaptive radiation for warning patterns and mimicry, and are excellent models to study the mechanisms underlying diversification. In Heliconius, mimicry rings typically involve distantly related species, whereas closely related species often join different mimicry rings. Genetic and behavioural studies have n how reproductive isolation in many pairs of Heliconius taxa is largely mediated by natural and sexual selection on wing colour patterns. However, recent studies have uncovered new cases in which pairs of closely related species are near‐perfect mimics of each other. Here, we provide morphometric and genetic evidence for the coexistence of two closely related, hybridizing co‐mimetic species on the eastern slopes of the Andes, H. melpomene amaryllis and H. timareta ssp. nov. , which is described here as H. timareta thelxinoe . A joint analysis of multilocus genotyping and geometric morphometrics of wing shape shows a high level of differentiation between the two species, with only limited gene flow and mixing. Some degree of genetic mixing can be detected, but putative hybrids were rare, only one of 175 specimens being a clear hybrid. In contrast, we found phenotypic differentiation between populations of H. timareta thelxinoe , possibly indicative of strong selection for local mimicry in different communities. In this pair of species, the absence of breakdown of genetic isolation despite near‐identical wing patterns implies that factors other than wing patterns keep the two taxa apart, such as chemical or behavioural signals, or ecological adaptation along a strong altitudinal gradient. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 830–847.     

Pupal Melanization in Heliconius erato phyllis (Lepidoptera; Nymphalidae): Genetic and Environmental Effects

This paper deals with estimates of heritability and the effect of two different kinds of environment on the melanization of pupae in the butterfly Heliconius erato phyllis. The results for heritability showed greater differences when the offspring were regressed on male or female values for pupal colour. Pupal colour followed a system of discrete scores, from 2 for the light pupae to 5 for the darkest one. Estimate of the heritability (h²) was 0.44 when the average score of the siblings were regressed on the mid-parental value. When regressed on the male parent, h² was larger (0.54), being very low when regressed on the female parent (0.09). Estimate by the analysis of variance was also 0.44 when both male and female sibs were included in the analysis (for male sibs, h²=0.43 and for female sibs, 0.47). To test for the effect of the environment, each brood was divided in three, one being the control, the second subjected to a black environment when reaching the fifth instar and the third to a white environment. Those in the black environment originated dark pupae that scored 5 or 4; for those in the white environment, there was no difference with the controls. Caterpillars when entering the prepupal stage were also subjected to the black treatment (early and late prepupal stage, respectively, for treatments called black A and black B). A significant effect was observed only when early prepupal stage was subjected to the black environment.     

Behavioral and Physiological Differences between Two Parapatric Heliconius Species1

The behavior and physiology of two parapatric sibling species, i Heliconius erato cyrbia Godt. and H. himera Hew., were investigated to assess if environmental adaptation enabled stable morphological, genetic, and ecological differences to exist in the face of hybridization. Morning and evening activity, egg production, and larval development time of H. himera and H. erato in insectaries were recorded; individuals were collected in allopatry and in sympatry from a hybrid zone in which the species overlapped. Studies were performed at ambient conditions within the natural range of H. himera. H. himera was considerably more active than H, erato flying earlier in the morning and later in the evening, even when both species were collected in sympatry. Similarly, H. himera laid more eggs, and the hatched larvae developed more rapidly. The results suggest that physiological constraints are an important selective force that may have been important in speciation and counteracts hybridization in the maintenance of the H. himeral H. erato contact zone. Ecological selection, arising from adaptation to low temperatures, may help explain the competitive exclusion of H. erato by H. himera in the drier, cooler montane habitat favored by the latter species.     

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

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

Molecular evolution of a long wavelength-sensitive opsin in mimetic Heliconius butterflies (Lepidoptera: Nymphalidae)

This study examines the pattern of opsin nucleotide and amino acid substitution among mimetic species 'rings' of Heliconius butterflies that are characterized by divergent wing colour patterns. A long wavelength opsin gene, OPS1 , was sequenced from each of seven species of Heliconius and one species of Dryas (Lepidoptera: Nymphalidae). A parsimony analysis of OPS1 nucleotide and amino acid sequences resulted in a phylogeny that was consistent with that presented by Brower & Egan in 1997, which was based on mitochondrial cytochrome oxidase I and II as well as nuclear wingless genes. Nodes in the OPS1 phylogeny were well supported by bootstrap analysis and decay indices. An analysis of specific sites within the gene indicates that the accumulation of amino acid substitutions has occurred independently of the morphological diversification of Heliconius wing colour patterns. Amino acid substitutions were examined with respect to their location within the opsin protein and their possible interactions with the chromophore and the G-protein. Of the 15 amino acid substitutions identified among the eight species, one nonconservative replacement (A226Q) was identified in a position that may be involved in binding with the G-protein.     

Comparative effects of inbreeding and outbreeding on the ontogeny of Heliconius erato phyllis (Lepidoptera: Nymphalidae)

Inbreeding depression has been reported in various groups of organisms, including insects. Estimates of inbreeding consequences were obtained by comparing 12 life‐history and morphological traits among nine inbred families (F = 0.25) and 16 outbred families (F = 0) of the Neotropical butterfly Heliconius erato phyllis. A Student's t‐test showed statistically significant differences for pupal weight and right forewing area, both in males and in females, between inbred and outbred families. Survival during development, from egg hatching to adulthood, also differed significantly between inbred and outbred families. The average number of haploid lethal equivalents was 0.17 for pupal weight, 0.15 for forewing area and 0.71 for survival from hatching to adulthood. The results of this study confirm that the consequences of inbreeding are more deleterious to life history traits than to morphological ones.     

Accumulation of Dobzhansky-Muller incompatibilities within a spatially structured population

Abstract.— A simple, deterministic analysis predicts that accumulation of Dobzhansky-Muller incompatibilities by a spatially structured population strongly depends on the number of negative interactions of an allele. If an allele can be incompatible with alleles at only one locus, incompatibilities accumulate linearly with time. In contrast, if an allele can participate in multiple pairwise incompatibilities with alleles at different loci, the expected number of incompatibilities eventually increases quadratically.     

Patterns of divergence in the olive sunbird Cyanomitra olivacea (Aves: Nectariniidae) across the African rainforest–savanna ecotone

In the debate over modes of vertebrate diversification in tropical rainforests, two competing hypotheses of speciation predominate: those that emphasize the role of geographical isolation during glacial periods and those that stress the role of ecology and diversifying selection across ecotones or environmental gradients. To investigate the relative roles of selection versus isolation in refugia, we contrasted genetic and morphologic divergence of the olive sunbird (Cyanomitra olivacea) at 18 sites (approximately 200 individuals) across the forest–savanna ecotone of Central Africa in a region considered to have harboured three hypothesized refugia during glacial periods. Habitats were characterized using bioclimatic and satellite remote‐sensing data. We found relatively high levels of gene flow between ecotone and forest populations and between refugia. Consistent with a pattern of divergence‐with‐gene‐flow, we found morphological characters to be significantly divergent across the gradient [forest versus ecotone (mean ± SD): wing length 60.47 ± 1.81 mm versus 62.18 ± 1.35 mm; tarsus length 15.51 ± 0.82 mm versus 16.00 ± 0.57 mm; upper mandible length 21.77 ± 1.09 mm versus 23.19 ± 0.98 mm, respectively]. Within‐habitat comparisons across forest and ecotone sites showed no significant differences in morphology. The results show that divergence in morphological traits is tied to environmental variables across the gradient and is occurring despite gene flow. The pattern of divergence‐with‐gene‐flow found is similar to that described for other rainforest species across the gradient. These results suggest that neither refugia, nor isolation‐by‐distance have played a major role in divergence in the olive sunbird, although ecological differences along the forest and savanna ecotone may impose significant selection pressures on the phenotype and potentially be important in diversification. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 821–835.     

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.     

Herbivore handling of a Plants trichome: the case of Heliconius charithonia (L.) (Lepidoptera: Nymphalidae) and Passiflora lobata (Killip) Hutch. (Passifloraceae)

Trichomes reduce herbivore attack on plants by physically and/or chemically inhibiting movement or other activities. Despite evidence that herbivores are negatively affected by trichomes there also reports of insect counter-adaptations that circumvent the plants defense. This paper reports on a study that investigated the likely mechanisms employed by larvae of the nymphalid butterfly, Heliconius charithonia (L.), that allow it to feed on a host that is presumably protected by hooked trichomes (Passiflora lobata (Killip) Hutch). Evidence were gathered using data from direct observations of larval movement and behavior, faeces analysis, scanning electron microscopy of plant surface and experimental analysis of larval movement on plants with and without trichomes (manually removed). The latter involved a comparison with a non specialist congener, Heliconius pachinus Salvin. Observations showed that H. charithonia larvae are capable of freeing themselves from entrapment on trichome tips by physical force. Moreover, wandering larvae lay silk mats on the trichomes and remove their tips by biting. In fact, trichome tips were found in the faeces. Experimental removal of trichomes aided in the movement of the non specialist but had no noticeable effect on the specialist larvae. These results support the suggestion that trichomes are capable of deterring a non specialist herbivore (H. pachinus). The precise mechanisms that allow the success of H. charithonia are not known, but I suggest that a blend of behavioral as well as physical resistance mechanisms is involved. Future studies should ascertain whether larval integument provides physical resistance to trichomes.     

Erebia serotina Descimon & de Lesse 1953 (Lepidoptera: Nymphalidae), a recurrent hybrid between two distantly related species

Erebia serotina was described in 1953 as a scarce, low-elevation endemic Pyrenean species flying late in the season. At least 34 individuals are known from various locations. However, the absence of females suggests a hybrid origin, and E. epiphron and either E. pronoe or E. manto have been proposed as possible parents. Electrophoretic analysis of five allozyme loci and sequencing of three mitochondrial DNA segments and one nuclear gene now demonstrate that E. serotina results from the cross between E. epiphron females and E. pronoe males. We have used our and previously published sequence data to generate a molecular phylogenetic tree of the genus Erebia which shows that these two species are only distantly related. The question of why they happen to hybridize on a seemingly routine basis is thus raised.     

Biometrical evidence for adaptations of the salivary glands to pollen feeding in Heliconius butterflies (Lepidoptera: Nymphalidae)

The Neotropical genus Heliconius (Nymphalidae) is unique among butterflies for its pollen-feeding behaviour. With the application of saliva, they extract amino acids from pollen grains on the outside of the proboscis. We predicted that the salivary glands of pollen-feeding Heliconiinae would show adaptations to this derived feeding behaviour. A biometrical analysis of the salivary glands revealed that pollen-feeding butterflies of the genus Heliconius have disproportionately longer and more voluminous salivary glands than nonpollen-feeding Nymphalidae. The first two components in the principal component analysis explained approximately 95% of the total variance. The size-dependent factor score coefficients of body length and salivary gland parameters were predominately represented on axis 1. They significantly discriminated pollen-feeding from nonpollen-feeding heliconiines on that axis. Factor score coefficients for the volume of the secretory region of the salivary glands separated heliconiines from the outgroup species. The detailed biometrical analysis of salivary glands features thus provides strong evidence that the secretory regions of the salivary glands are larger in pollen-feeding butterflies. We concluded that pollen feeding is associated with a high production of salivary fluid.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 604–612.     

Embryogenesis of Heliconius erato (Lepidoptera,Nymphalidae): A contribution to the anatomical development of an evo‐devo model organism

This study reports on the embryogenesis of Heliconius erato phyllis between blastoderm formation and the prehatching larval stage. Syncytial blastoderm formation occurred approximately 2 h after egg laying (AEL) and at about 4 h, the cellular blastoderm was formed. The germ band arose from the entire length of the blastoderm, and rapidly became compacted occupying approximately two‐thirds of the egg length. At about 7 h AEL, protocephalon and protocorm differentiation occurred. Continued proliferation of the germ band was followed by penetration into the yolk mass, forming a C‐shaped embryo at about 10 h. Approximately 12 h AEL, the gnathal, thoracic and abdominal segments became visible. The primordium of the mouthparts and thoracic legs formed as paired evaginations, while the prolegs formed as paired lobes. At about 30 h, the embryo reversed dorsoventrally. Approximately 32 h AEL, the protocephalon and gnathal segments fused, shifting the relative position of the rudimentary appendages in this region. At about 52 h, the embryo was U‐shaped in lateral view and at approximately 56 h, the bristles began evagination from the larval cuticle. Larvae hatched at about 72 h. We found that H. erato phyllis followed an embryonic pattern consistent with long‐germ embryogenesis. Thus, we believe that H. erato phyllis should be classified as a long‐germ lepidopteran. The study of H. erato phyllis embryogenesis provided a structural glimpse into the morphogenetic events that occur in the Heliconius egg period. This study could help future molecular approaches to understanding the evolution of Heliconius development.     

Highly polymorphic microsatellite loci for Speyeria idalia (Lepidoptera: Nymphalidae)

Four polymorphic microsatellite loci were identified in the butterfly Speyeria idalia. We constructed a phagemid library and screened approximately 120 000 inserts. Probing with GT15, we identified 36 positives and designed polymerase chain reaction (PCR) primers for 12 potential loci. Of those loci, only four consistently produced polymorphic, diallelic PCR products in the expected size range. These results are consistent with previous studies concerning the low frequency of microsatellite loci in the lepidoptera, although these four loci are highly polymorphic and therefore likely to provide information on the fine scale genetic structure among populations in this species.     

Hybridization among sympatric species of Rhododendron (Ericaceae) in Turkey: morphological and molecular evidence

Rhododendron (Ericaceae) is a large genus in which barriers to hybridization are especially weak, but many species are maintained in sympatry. Hybridization among four species of Rhododendron subsect. Pontica, which occur in sympatry in Turkey, was investigated. Material of R. ponticum, R. smirnovii, R. ungernii, and R. caucasicum and their putative hybrids was collected from the wild. Based on morphology, chloroplast DNA and nuclear ribosomal DNA restriction fragment length polymorphism (RFLP) profiles, each accession was identified as a species or hybrid combination. Five of the six possible hybrid combinations among the four species were detected. Rhododendron ponticum × R. smirnovii was represented by a single individual and R. caucasicum × R. smirnovii by one small group of hybrid plants. The combinations R. ponticum × R. ungernii and R. ungernii × R. smirnovii showed evidence of frequent backcrossing, while R. ponticum × R. caucasicum appeared unusual in that an intermediate hybrid type was abundant, whereas hybrids with phenotypes approaching either parent were rare. Possible explanations of this latter situation are discussed. The results suggest that natural hybridization among Rhododendron species is common and that ecological factors are important in maintaining integrity when species occur in sympatry.     

Genetic panmixia within a narrow contact zone between chromosomally and ecologically distinct black fly sibling species (Diptera: Simuliidae)

Hybrid zones are windows into the speciation process, and their study can give clues into the maintenance and breakdown of species boundaries. Using both genetic and ecological tools, we investigate lineage diversification across a contact zone characterized by chromosome rearrangements. We show that black fly sibling species, Simulium arcticum sensu stricto (s.s.) and Simulium saxosum, lack genetic differentiation at both microsatellite and mtDNA loci in allopatry and sympatry, as well as exhibit high levels of gene flow and continuous chromosome variation in sympatry. Furthermore, hybrid frequencies at the contact zone are similar to those seen between races, rather than species. In contrast, S. arcticum s.s. and S. saxosum maintain ecological differences and distinct habitat associations ‐ the contact zone situated at the margin of suitable habitat for each sibling species. Moreover, gene flow occurs only in a narrow band along an ecological transition. Except for the contact zone, S. arcticum s.s. and S. saxosum hybrids do not occur elsewhere within the sibling species' ranges. Although S. arcticum s.s. and S. saxosum maintain the potential to interbreed freely, we conclude that habitat associations and, perhaps, chromosome systems prevent expansion of ranges and assimilation of lineages.     

The genomics of incompatibility factors and sex determination in hybridizing species of Cottus (Pisces)

Cottus rhenanus and Cottus perifretum have formed hybrid lineages and narrow hybrid zones that can be best explained through the action of natural selection. However, the underlying selective forces as well as their genomic targets are not well understood. This study identifies genomic regions in the parental species that cause hybrid incompatibilities and tests whether these manifest in a sex-specific manner to learn about processes that affect natural hybridization in Cottus. Interspecific F₂ crosses were analyzed for 255 markers for genetic mapping and to detect transmission distortion as a sign for genetic incompatibilities. The Cottus map consists of 24 linkage groups with a total length of 1575.4 cM. A male heterogametic (XY) sex determination region was found on different linkage groups in the two parental species. Genetic incompatibilities were incomplete, varied among individuals and populations and were not associated with the heterogametic sex. The variance between populations and individuals makes it unlikely that there are species-specific incompatibility loci that could affect the gene pool of natural hybrids in a simple and predictable way. Conserved synteny with sequenced fish genomes permits to genetically study the Cottus genome through the transfer of genomic information from the model fish species. Homology relationships of candidate genomic regions in Cottus indicate that sex determination is not based on the same genomic regions found in other fish species. This suggests a fast evolutionary turnover of the genetic basis of sex determination that, together with the small size of the heterogametic regions, may contribute to the absence of fitness effects related to the Haldane''s rule.     

Hybrid zone origins, species boundaries, and the evolution of wing-pattern diversity in a polytypic species complex of North American admiral butterflies (Nymphalidae: Limenitis)

Hybrid zones present opportunities to study the effects of gene flow, selection, and recombination in natural populations and, thus, provide insights into the genetic and phenotypic changes that occur early in speciation. Here we investigate a hybrid zone between mimetic (Limenitis arthemis astyanax) and nonmimetic (Limenitis arthemis arthemis) populations of admiral butterflies using DNA sequence variation from mtDNA and seven nuclear gene loci. We find three distinct mitochondrial clades within this complex, and observe a strong overall concordance between wing-pattern phenotypes and mitochondrial variation. Nuclear gene genealogies, in contrast, revealed no evidence of exclusivity for either wing-pattern phenotype, suggesting incomplete barriers to gene exchange and/or insufficient time for lineage sorting. Coalescent simulations indicate that gene flow between these two subspecies is highly asymmetric, with the majority of migration occurring from mimetic into nonmimetic populations. Selective sweeps of alleles responsible for mimetic phenotypes may have occurred more than once when mimetic and nonmimetic Limenitis occurred together in the presence of the model (Battus philenor).