Dataset Information

Natural hybridization in heliconiine butterflies: the species boundary as a continuum.

ABSTRACT:

Background

To understand speciation and the maintenance of taxa as separate entities, we need information about natural hybridization and gene flow among species.

Results

Interspecific hybrids occur regularly in Heliconius and Eueides (Lepidoptera: Nymphalidae) in the wild: 26-29% of the species of Heliconiina are involved, depending on species concept employed. Hybridization is, however, rare on a per-individual basis. For one well-studied case of species hybridizing in parapatric contact (Heliconius erato and H. himera), phenotypically detectable hybrids form around 10% of the population, but for species in sympatry hybrids usually form less than 0.05% of individuals. There is a roughly exponential decline with genetic distance in the numbers of natural hybrids in collections, both between and within species, suggesting a simple "exponential failure law" of compatibility as found in some prokaryotes.

Conclusion

Hybridization between species of Heliconius appears to be a natural phenomenon; there is no evidence that it has been enhanced by recent human habitat disturbance. In some well-studied cases, backcrossing occurs in the field and fertile backcrosses have been verified in insectaries, which indicates that introgression is likely, and recent molecular work shows that alleles at some but not all loci are exchanged between pairs of sympatric, hybridizing species. Molecular clock dating suggests that gene exchange may continue for more than 3 million years after speciation. In addition, one species, H. heurippa, appears to have formed as a result of hybrid speciation. Introgression may often contribute to adaptive evolution as well as sometimes to speciation itself, via hybrid speciation. Geographic races and species that coexist in sympatry therefore form part of a continuum in terms of hybridization rates or probability of gene flow. This finding concurs with the view that processes leading to speciation are continuous, rather than sudden, and that they are the same as those operating within species, rather than requiring special punctuated effects or complete allopatry. Although not qualitatively distinct from geographic races, nor "real" in terms of phylogenetic species concepts or the biological species concept, hybridizing species of Heliconius are stably distinct in sympatry, and remain useful groups for predicting morphological, ecological, behavioural and genetic characteristics.

SUBMITTER: Mallet J

PROVIDER: S-EPMC1821009 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundA key to understanding the origins of species is determining the evolutionary processes that drive the patterns of genomic divergence during speciation. New genomic technologies enable the study of high-resolution genomic patterns of divergence across natural speciation continua, where taxa pairs with different levels of reproductive isolation can be used as proxies for different stages of speciation. Empirical studies of these speciation continua can provide valuable insights into how genomes diverge during speciation.MethodsWe examine variation across a handful of genomic regions in parapatric and allopatric populations of Heliconius butterflies with varying levels of reproductive isolation. Genome sequences were mapped to 2.2-Mb of the H. erato genome, including 1-Mb across the red color pattern locus and multiple regions unlinked to color pattern variation.ResultsPhylogenetic analyses reveal a speciation continuum of pairs of hybridizing races and incipient species in the Heliconius erato clade. Comparisons of hybridizing pairs of divergently colored races and incipient species reveal that genomic divergence increases with ecological and reproductive isolation, not only across the locus responsible for adaptive variation in red wing coloration, but also at genomic regions unlinked to color pattern.DiscussionWe observe high levels of divergence between the incipient species H. erato and H. himera, suggesting that divergence may accumulate early in the speciation process. Comparisons of genomic divergence between the incipient species and allopatric races suggest that limited gene flow cannot account for the observed high levels of divergence between the incipient species.ConclusionsOur results provide a reconstruction of the speciation continuum across the H. erato clade and provide insights into the processes that drive genomic divergence during speciation, establishing the H. erato clade as a powerful framework for the study of speciation.

Project description:We investigate the influence of three common definitions of the solute/solvent dielectric boundary (DB) on the accuracy of the electrostatic solvation energy ?Gel computed within the Poisson Boltzmann and the generalized Born models of implicit solvation. The test structures include small molecules, peptides and small proteins; explicit solvent ?Gel are used as accuracy reference. For common atomic radii sets BONDI, PARSE (and ZAP9 for small molecules) the use of van der Waals (vdW) DB results, on average, in considerably larger errors in ?Gel than the molecular surface (MS) DB. The optimal probe radius ?w for which the MS DB yields the most accurate ?Gel varies considerably between structure types. The solvent accessible surface (SAS) DB becomes optimal at ?w ~ 0.2 Å (exact value is sensitive to the structure and atomic radii), at which point the average accuracy of ?Gel is comparable to that of the MS-based boundary. The geometric equivalence of SAS to vdW surface based on the same atomic radii uniformly increased by ?w gives the corresponding optimal vdW DB. For small molecules, the optimal vdW DB based on BONDI + 0.2 Å radii can yield ?Gel estimates at least as accurate as those based on the optimal MS DB. Also, in small molecules, pairwise charge-charge interactions computed with the optimal vdW DB are virtually equal to those computed with the MS DB, suggesting that in this case the two boundaries are practically equivalent by the electrostatic energy criteria. In structures other than small molecules, the optimal vdW and MS dielectric boundaries are not equivalent: the respective pairwise electrostatic interactions in the presence of solvent can differ by up to 5 kcal/mol for individual atomic pairs in small proteins, even when the total ?Gel are equal. For small proteins, the average decrease in pairwise electrostatic interactions resulting from the switch from optimal MS to optimal vdW DB definition can be mimicked within the MS DB definition by doubling of the solute dielectric constant. However, the use of the higher interior dielectric does not eliminate the large individual deviations between pairwise interactions computed within the two DB definitions. It is argued that while the MS based definition of the dielectric boundary is more physically correct in some types of practical calculations, the choice is not so clear in some other common scenarios.