First-generation linkage map of the warningly colored butterfly Heliconius erato.

We report the first genetic linkage map of Heliconius erato, a species that shows remarkable variation in its warningly colored wing patterns. We use crosses between H. erato and its sister species, H. himera, to place two major color pattern genes, D and Cr, on a linkage map containing AFLP, allozyme, microsatellite and single-copy nuclear loci. We identified all 21 linkage groups in an initial genetic screen of 22 progeny from an F1 female x male H. himera family. Of the 229 markers, 87 used to identify linkage groups were also informative in 35 progeny from a sibling backcross (H. himera female x F1 male). With these, and an additional 33 markers informative in the second family, we constructed recombinational maps for 19 of the 21 linkage groups. These maps varied in length from 18.1 to 431.1 centimorgans (cM) and yielded an estimated total length of 2400 cM. The average distance between markers was 23 cM, and eight of the 19 linkage groups, including the sex chromosome (Z) and the chromosome containing the Cr locus, contained two or more codominant anchor loci. Of the three potential candidate genes mapped here, Cubitus interruptus (Ci), Decapentaplegic (Dpp) and Wingless (Wg), only Ci was linked, although loosely, to a known Heliconius color pattern locus. This work is an important first step for constructing a denser genetic map of the H. erato color pattern radiation and for a comparative genomic study of the architecture of mimicry in Heliconius butterflies.

Historical demography of Mullerian mimicry in the neotropical Heliconius butterflies.

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, theta, 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 wing-pattern radiations were not coincident in the two species. Instead, larger effective population size (N(e)) 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.