Hybrid zones can be valuable tools for studying evolution and identifying genomic regions responsible for adaptive divergence and underlying phenotypic variation. Hybrid zones between subspecies of Heliconius butterflies can be very narrow and are maintained by strong selection acting on colour pattern. The co-mimetic species H. erato and H. melpomene have parallel hybrid zones where both species undergo a change from one colour pattern form to another. We use restriction associated DNA sequencing to obtain several thousand genome wide sequence markers and use these to analyse patterns of population divergence across two pairs of parallel hybrid zones in Peru and Ecuador. We compare two approaches for analysis of this type of data; alignment to a reference genome and de novo assembly, and find that alignment gives the best results for species both closely (H. melpomene) and distantly (H. erato, ~15% divergent) related to the reference sequence. Our results confirm that the colour pattern controlling loci account for the majority of divergent regions across the genome, but we also detect other divergent regions apparently unlinked to colour pattern differences. We also use association mapping to identify previously unmapped colour pattern loci, in particular the Ro locus. Finally, we identify within our sample a new cryptic population of H. timareta in Ecuador, which occurs at relatively low altitude and is mimetic with H. melpomene malleti.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present “haplotagging,” a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Hybrid zones can be valuable tools for studying evolution and identifying genomic regions responsible for adaptive divergence and underlying phenotypic variation. Hybrid zones between subspecies of Heliconius butterflies can be very narrow and are maintained by strong selection acting on color pattern. The comimetic species, H. erato and H. melpomene, have parallel hybrid zones in which both species undergo a change from one color pattern form to another. We use restriction-associated DNA sequencing to obtain several thousand genome-wide sequence markers and use these to analyze patterns of population divergence across two pairs of parallel hybrid zones in Peru and Ecuador. We compare two approaches for analysis of this type of data-alignment to a reference genome and de novo assembly-and find that alignment gives the best results for species both closely (H. melpomene) and distantly (H. erato,~15% divergent) related to the reference sequence. Our results confirm that the color pattern controlling loci account for the majority of divergent regions across the genome, but we also detect other divergent regions apparently unlinked to color pattern differences. We also use association mapping to identify previously unmapped color pattern loci, in particular the Ro locus. Finally, we identify a new cryptic population of H. timareta in Ecuador, which occurs at relatively low altitude and is mimetic with H. melpomene malleti.

One-sentence summary:Haplotagging, a novel linked-read sequencing technique that enables whole genome haplotyping in large populations, reveals the formation of a novel hybrid race in parallel hybrid zones of two co-mimicking Heliconius butterfly species through strikingly parallel divergences in their genomes. Short title:Haplotagging reveals parallel formation of hybrid races Abstract Genetic variation segregates as linked sets of variants, or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. And yet, genomic data often lack haplotype information, due to constraints in sequencing technologies. Here we present "haplotagging", a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabasesize haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high-and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species the geographic clines for the major wing pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the centre of the hybrid zone.We propose that shared warning signalling (Müllerian mimicry) may couple the cline shifts seen in both species, and facilitate the parallel co-emergence of a novel hybrid morph in both co-mimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

To the extent that geography correlates with particular environmental parameters, the geographical distribution of phylogenetically related social and nonsocial organisms should shed light on the conditions that lead to sociality versus nonsociality. Social spiders are notorious for being concentrated in tropical regions of the world, occupying a set of habitats more restricted than those available to the phylogenetic lineages in which they occur. Here we document a parallel pattern involving elevation in the spider genus Anelosimus in America and describe the biology of a newly discovered social species found at what appears to be the altitudinal edge of sociality in the genus. We show that this is a cooperative permanent-social species with highly female-biased sex ratios but colonies that are one to two orders of magnitude smaller than those of a low-elevation congener of similar body size. We suggest that the absence of subsocial Anelosimus species in the lowland rain forest may be due to an increased probability of maternal death in this habitat due to greater predation and/or precipitation, while absence of a sufficient supply of large insects at high elevations or latitudes may restrict social species to low- to midelevation tropical moist forests. We refer to these as the "maternal survival" and "prey size" hypotheses, respectively, and suggest that both in combination may explain the geographical distribution of sociality in the genus.

The origin and evolution of supergenes have long fascinated evolutionary biologists. In the polymorphic butterfly Heliconius numata, a supergene controls the switch between multiple different forms, and results in near-perfect mimicry of model species. Here, we use a morphometric analysis to quantify the variation in wing pattern observed in two broods of H. numata with different alleles at the supergene locus, 'P '. Further, we genotype the broods to associate the variation we capture with genetic differences. This allows us to begin mapping the quantitative trait loci that have minor effects on wing pattern. In addition to finding loci on novel chromosomes, our data, to our knowledge, suggest for the first time that ancestral colour-pattern loci, known to have major effects in closely related species, may contribute to the wing patterns displayed by H. numata, despite the large transfer of effects to the supergene.

Despite more than a century of biological research on the evolution and maintenance of mimetic signals, the relative frequencies of models and mimics necessary to establish and maintain Batesian mimicry in natural populations remain understudied. Here we investigate the frequency-dependent dynamics of imperfect Batesian mimicry, using predation experiments involving artificial butterfly models. We use two geographically distinct populations of butterflies that vary in their relative frequencies of a putatively defended model () and Batesian mimic (). We found that in Costa Rica, where both species share similar abundances, Batesian mimicry breaks down, and predators more readily attack artificial butterfly models of the presumed mimic, By contrast, in Ecuador, where (model) is significantly more abundant than (mimic), both species are equally protected from predation. Our results provide compelling experimental evidence that imperfect Batesian mimicry is frequency-dependent on the relative abundance of models and mimics in natural populations, and contribute to the growing body of evidence that complex dynamics, such as seasonality or the availability of alternative prey, influence the evolution of mimetic traits.

Bright, highly reflective iridescent colours can be seen across nature and are produced by the scattering of light from nanostructures. Heliconius butterflies have been widely studied for their diversity and mimicry of wing colour patterns. Despite iridescence evolving multiple times in this genus, little is known about the genetic basis of the colour and the development of the structures which produce it. Heliconius erato can be found across Central and South America, but only races found in western Ecuador and Colombia have developed blue iridescent colour. Here, we use crosses between iridescent and non-iridescent races of H. erato to study phenotypic variation in the resulting F 2 generation. Using measurements of blue colour from photographs, we find that iridescent structural colour is a quantitative trait controlled by multiple genes, with strong evidence for loci on the Z sex chromosome. Iridescence is not linked to the Mendelian colour pattern locus that also segregates in these crosses (controlled by the gene cortex ). Small-angle X-ray scattering data show that spacing between longitudinal ridges on the scales, which affects the intensity of the blue reflectance, also varies quantitatively in F 2 crosses.