Background: Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri. Results: We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex. Conclusions: Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

Selection shapes genetic diversity around target mutations, yet little is known about how selection on specific loci affects the genetic trajectories of populations, including their genome-wide patterns of diversity and demographic responses. Adaptive introgression provides a way to assess how adaptive evolution at one locus impacts whole-genome biology. Here we study the patterns of genetic variation and geographic structure in a neotropical butterfly, Heliconius numata, and its closely related allies in the so-called melpomene-silvaniform subclade. H. numata is known to have evolved a supergene via the introgression of an adaptive inversion about 2.2 My ago, triggering a polymorphism maintained by balancing selection. This locus controls a polymorphism of wing patterns involved in mimicry associations with distinct groups of co-mimics, and butterflies show disassortative mate preferences and heterozygote advantage at this locus. We contrasted patterns of genetic diversity and structure 1) among extant polymorphic and monomorphic populations of H. numata, 2) between H. numata and its close relatives, and 3) between ancestral lineages in a phylogenetic framework. We show that H. numata populations which carry the introgressed inversions in a balanced polymorphism show markedly distinct patterns of diversity compared to all other taxa. They show the highest diversity and demographic estimates in the entire clade, as well as a remarkably low level of geographic structure and isolation by distance across the entire Amazon basin. By contrast, monomorphic populations of H. numata as well as its sister species and their ancestral lineages all show the lowest effective population sizes and genetic diversity in the clade, and higher levels of geographical structure across the continent. This suggests that the large effective population size of polymorphic populations could be a property associated with harbouring the supergene. Our results are consistent with the hypothesis that the adaptive introgression of the inversion triggered a shift from directional to balancing selection and a change in gene flow due to disassortative mating, causing a general increase in genetic diversity and the homogenisation of genomes at the continental scale.

A population of Parides anchises nephalion was studied during seven months in SE Brazil. The population size was about 10-20 individuals (with theoretic maximum near 60 individuals), with small variations in some months. Sex ratio was male biased, with males dominating in all months. The age structure was not stable, with an increase in new individuals before the population peak in December. The residence time was 14.1 ± 8.2 days for males and 9.0 ± 3.6 to females, with the maximum registered of 30 days. Males can travel distances of up to 400 m, but most individuals were always recaptured in the same site. The mean forewing length was greater in females. The population features agree with those found in other species of Parides in other neotropical sites.

No abstract

This article documents the addition of 299 microsatellite marker loci and nine pairs of single-nucleotide polymorphism (SNP) EPIC primers to the Molecular Ecology Resources (MER) Database. Loci were developed for the following species: Alosa pseudoharengus, Alosa aestivalis, Aphis spiraecola, Argopecten purpuratus, Coreoleuciscus splendidus, Garra gotyla, Hippodamia convergens, Linnaea borealis, Menippe mercenaria, Menippe adina, Parus major, Pinus densiflora, Portunus trituberculatus, Procontarinia mangiferae, Rhynchophorus ferrugineus, Schizothorax richardsonii, Scophthalmus rhombus, Tetraponera aethiops, Thaumetopoea pityocampa, Tuta absoluta and Ugni molinae. These loci were cross-tested on the following species: Barilius bendelisis, Chiromantes haematocheir, Eriocheir sinensis, Eucalyptus camaldulensis, Eucalyptus cladocalix, Eucalyptus globulus, Garra litaninsis vishwanath, Garra para lissorhynchus, Guindilla trinervis, Hemigrapsus sanguineus, Luma chequen. Guayaba, Myrceugenia colchagüensis, Myrceugenia correifolia, Myrceugenia exsucca, Parasesarma plicatum, Parus major, Portunus pelagicus, Psidium guayaba, Schizothorax richardsonii, Scophthalmus maximus, Tetraponera latifrons, Thaumetopoea bonjeani, Thaumetopoea ispartensis, Thaumetopoea libanotica, Thaumetopoea pinivora, Thaumetopoea pityocampa ena clade, Thaumetopoea solitaria, Thaumetopoea wilkinsoni and Tor putitora. This article also documents the addition of nine EPIC primer pairs for Euphaea decorata, Euphaea formosa, Euphaea ornata and Euphaea yayeyamana.