No abstract

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

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

Among animals, evidence for homoploid hybrid speciation (HHS, i.e. the creation of a hybrid lineage without a change in chromosome number) was limited until recently to the virgin chub, Gila seminuda, and some controversial data in support of hybrid status for the red wolf, Canis rufus. This scarcity of evidence, together with pessimistic attitudes among zoologists about the evolutionary importance of hybridisation, prompted the view that HHS is extremely rare among animals, especially as compared with plants. However, in recent years, the literature on animal HHS has expanded to include several new putative examples in butterflies, ants, flies and fishes. We argue that this evidence suggests that HHS is far more common than previously thought and use it to provide insights into some of the genetic and ecological aspects associated with this type of speciation among animals.

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.

The mimetic butterflies Heliconius erato and Heliconius melpomene have undergone parallel radiations to form a near-identical patchwork of over 20 different wing-pattern races across the Neotropics. Previous molecular phylogenetic work on these radiations has suggested that similar but geographically disjunct color patterns arose multiple times independently in each species. The neutral markers used in these studies, however, can move freely across color pattern boundaries, and therefore might not represent the history of the adaptive traits as accurately as markers linked to color pattern genes. To assess the evolutionary histories across different loci, we compared relationships among races within H. erato and within H. melpomene using a series of unlinked genes, genes linked to color pattern loci, and optix, a gene recently shown to control red colorpattern variation. We found that although unlinked genes partition populations by geographic region, optix had a different history, structuring lineages by red color patterns and supporting a single origin of red-rayed patterns within each species. Genes closely linked (80-250 kb) to optix exhibited only weak associations with color pattern. This study empirically demonstrates the necessity of examining phenotype-determining genomic regions to understand the history of adaptive change in rapidly radiating lineages. With these refined relationships, we resolve a long-standing debate about the origins of the races within each species, supporting the hypothesis that the red-rayed Amazonian pattern evolved recently and expanded, causing disjunctions of more ancestral patterns.Müllerian mimicry | population genetics | phylogeography R esearchers typically rely on neutrally evolving loci to generate a phylogenetic and population genetic history of adaptive divergence. The rationale is that these markers provide an unbiased view of the relationships among divergent phenotypes and a better understanding of the evolutionary processes generating variation. However, the genome is a complicated mosaic shaped by an interplay of mutation, drift, selection, and recombination. Recombination allows different regions of the genome to experience alternative restrictions to gene flow, and thus develop different evolutionary trajectories. The closer a genetic marker is to the alleles responsible for adaptive differences, the more likely that it will trace the history of phenotypic change.Understanding how phenotypic variation is generated in nature is greatly enhanced by studying groups that are actively undergoing diversification. By deciphering the history of such diverse phenotypes we gain a clearer understanding of the evolutionary process, including the tempo and mode of phenotypic change. Heliconius butterflies present one of the most striking examples of a recent phenotypic radiation. The 40 species in the genus exhibit hundreds of wing patterns that are involved in Müllerian mimicry complexes, where distasteful species converge on a shared warning signal to avoid predation. Th...

We compare the historical demographies of two Mü llerian comimetic butterfly species: Heliconius erato and Heliconius melpomene. These species show an extensive parallel geographic divergence in their aposematic wing phenotypes. Recent studies suggest that this coincident mosaic results from simultaneous demographic processes shaped by extrinsic forces over Pleistocene climate fluctuations. However, DNA sequence variation at two rapidly evolving unlinked nuclear loci, Mannose phosphate isomerase (Mpi) and Triose phosphate isomerase (Tpi), show that the comimetic species have quite different quaternary demographies. In H. erato, despite ongoing lineage sorting across the Andes, nuclear genealogical estimates showed little geographical structure, suggesting high historical gene flow. Coalescent-based demographic analysis revealed population growth since the Pliocene period. Although these patterns suggest vicariant population subdivision associated with the Andean orogeny, they are not consistent with hypotheses of Pleistocene population fragmentation facilitating allopatric wing phenotype radiation in H. erato. In contrast, nuclear genetic diversity, , in H. melpomene was reduced relative to its comimic and revealed three phylogeographical clades. The pattern of coalescent events within regional clades was most consistent with population growth in relatively isolated populations after a recent period of restricted population size. These different demographic histories suggest that the wingpattern radiations were not coincident in the two species. Instead, larger effective population size (Ne) in H. erato, together with profound population change in H. melpomene, supports an earlier hypothesis that H. erato diversified first as the model species of this remarkable mimetic association.

Homoploid hybrid speciation (HHS) is the establishment of a novel species through introgressive hybridization without a change in chromosome number. We discuss different routes by which this might occur and propose a novel term, 'hybrid trait speciation', which combines the idea that hybridization can generate adaptive novelty with the 'magic trait' model of ecological speciation. Heliconius butterflies contain many putative examples of hybrid colour patterns, but only recently has the HHS hypothesis been tested explicitly in this group. Molecular data has shown evidence for gene flow between many distinct species. Furthermore, the colour pattern of Heliconius heurippa can be recreated in laboratory crosses between Heliconius melpomene and Heliconius cydno and, crucially, plays a role in assortative mating between the three species. Nonetheless, although the genome of H. heurippa shows evidence for hybridization, it is not a mosaic of the two parental species. Instead, ongoing hybridization has likely blurred any signal of the original speciation event. We argue that where hybridization leads to novel adaptive traits that also cause reproductive isolation, it is likely to trigger speciation.