Geographic contrasts between pre- and postzygotic barriers are consistent with reinforcement in Heliconius butterflies

Complexity of Chemical Emissions Increases Concurrently with Sexual Maturity in Heliconius Butterflies

Pheromone communication is widespread among animals. Since it is often involved in mate choice, pheromone production is often tightly controlled. Although male sex pheromones (MSPs) and anti-aphrodisiacs have been studied in some Heliconius butterfly species, little is known about the factors affecting their production and release in these long-lived butterflies. Here, we investigate the effect of post-eclosion age on chemical blends from pheromone-emitting tissues in Heliconius atthis and Heliconius charithonia, exhibiting respectively free-mating and pupal-mating strategies that are hypothesised to differently affect the timing of their pheromone emissions. We focus on two different tissues: the wing androconia, responsible for MSPs used in courtship, and the genital tip, the production site for anti-aphrodisiac pheromones that affect post-mating behaviour. Gas chromatography-mass spectrometric analysis of tissue extracts from virgin males and females of both species from day 0 to 8 post-eclosion demonstrates the following. Some ubiquitous fatty acid precursors are already detectable at day 0. The complexity of the chemical blends increases with age regardless of tissue or sex. No obvious difference in the time course of blend production was evident between the two species, but female tissues in H. charithonia were more affected by age than in H. atthis. We suggest that compounds unique to male androconia and genitals and whose amount increases with age are potential candidates for future investigation into their roles as pheromones. While this analysis revealed some of the complexity in Heliconius chemical ecology, the effects of other factors, such as the time of day, remain unknown.

Hybrid speciation driven by multilocus introgression of ecological traits

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

Pervasive mimicry in flight behavior among aposematic butterflies

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Major patterns in the introgression history of Heliconius butterflies

Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.

On the impermanence of species: The collapse of genetic incompatibilities in hybridizing populations

Species pairs often become genetically incompatible during divergence, which is an important source of reproductive isolation. An idealized picture is often painted where incompatibility alleles accumulate and fix between diverging species. However, recent studies have shown both that incompatibilities can collapse with ongoing hybridization, and that incompatibility loci can be polymorphic within species. This paper suggests some general rules for the behavior of incompatibilities under hybridization. In particular, we argue that redundancy of genetic pathways can strongly affect the dynamics of intrinsic incompatibilities. Since fitness in genetically redundant systems is unaffected by introducing a few foreign alleles, higher redundancy decreases the stability of incompatibilities during hybridization, but also increases tolerance of incompatibility polymorphism within species. We use simulations and theories to show that this principle leads to two types of collapse: in redundant systems, exemplified by classical Dobzhansky-Muller incompatibilities, collapse is continuous and approaches a quasi-neutral polymorphism between broadly sympatric species, often as a result of isolation-by-distance. In nonredundant systems, exemplified by co-evolution among genetic elements, incompatibilities are often stable, but can collapse abruptly with spatial traveling waves. As both types are common, the proposed principle may be useful in understanding the abundance of genetic incompatibilities in natural populations.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Exploitation of an ancestral pheromone biosynthetic pathway contributes to diversification in Heliconius butterflies

During courtship, male butterflies of many species produce androconial secretions containing male sex pheromones (MSPs) that communicate species identity and affect female choice. MSPs are thus likely candidates as reproductive barriers, yet their role in speciation remains poorly studied. Although Heliconius butterflies are a model system in speciation, their MSPs have not been investigated from a macroevolutionary perspective. We use GC/MS to characterize male androconial secretions in 33 of the 69 species in the Heliconiini tribe. We found these blends to be species-specific, consistent with a role in reproductive isolation. We detected a burst in blend diversification rate at the most speciose genus, Heliconius; a consequence of Heliconius and Eueides species using a fatty acid (FA) metabolic pathway to unlock more complex blends than basal Heliconiini species, whose secretions are dominated by plant-like metabolites. A comparison of 10 sister species pairs demonstrates a striking positive correlation between blend dissimilarity and range overlap, consistent with character displacement or reinforcement in sympatry. These results demonstrate for the first time that MSP diversification can promote reproductive isolation across this group of butterflies, showcasing how implementation of an ancestral trait, the co-option of the FA metabolic pathway for pheromone production, can facilitate rapid speciation.

Head and Tail Oxidized Terpenoid Esters from Androconia of Heliconius erato Butterflies

Heliconius erato is a neotropical butterfly species that is part of a complex mimicry ring, with colorful wing patterns. For intraspecific communication, males use pheromones that are released from two different scent-emitting structures. Scent glands located near the abdominal claspers of males, containing antiaphrodisiac pheromones, release a highly complex mixture of compounds that is transferred to females during mating, rendering them unattractive to other males. On the other hand, androconia, scent-emitting scale areas on the wings of male butterflies, release a structurally more restricted set of compounds that likely serves an aphrodisiac role. We report here on two structurally related compounds that are the major androconial constituents, produced in high amounts and are not volatile due to their high molecular mass. Their structures were established by extensive analysis of mass, infrared, and NMR spectra, as well as microderivatization reactions of the natural extract. After establishing synthetic access, the compounds were unequivocally identified as two unusual head and tail oxidized terpenoids, (4E,8E,12E)-4,8,12-trimethyl-16-oxoheptadeca-4,8,12-trien-1-yl oleate (1) and stearate (2). Although behavioral assays are necessary to fully comprehend their role in the chemical communication of the species, hypotheses for their use by the butterflies are also discussed.

Prevalence and Adaptive Impact of Introgression

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

Sexual selection and introgression in avian hybrid zones: Spotlight on Manacus

Hybrid zones offer a window into the processes and outcomes of evolution, from species formation or fusion to genomic underpinnings of specific traits and isolating mechanisms. Sexual selection is believed to be an important factor in speciation processes, and hybrid zones present special opportunities to probe its impact. The manakins (Aves, Pipridae) are a promising group in which to study the interplay of sexual selection and natural hybridization: they show substantial variation across the family in the strength of sexual selection they experience, they readily hybridize within and between genera, and they appear to have formed hybrid species, a rare event in birds. A hybrid zone between two manakins in the genus Manacus is unusual in that plumage and behavioral traits of one species have introgressed asymmetrically into populations of the second species through positive sexual selection, then apparently stalled at a river barrier. This is one of a handful of documented examples of asymmetric sexual trait introgression with a known selective mechanism. It offers opportunities to examine reproductive isolation, introgression, plumage color evolution, and natural factors enhancing or constraining the effects of sexual selection in real time. Here, we review previous work in this system, propose new hypotheses for observed patterns, and recommend approaches to test them.

Synteny-Based Genome Assembly for 16 Species of Heliconius Butterflies, and an Assessment of Structural Variation across the Genus

Heliconius butterflies (Lepidoptera: Nymphalidae) are a group of 48 neotropical species widely studied in evolutionary research. Despite the wealth of genomic data generated in past years, chromosomal level genome assemblies currently exist for only two species, Heliconius melpomene and Heliconius erato, each a representative of one of the two major clades of the genus. Here, we use these reference genomes to improve the contiguity of previously published draft genome assemblies of 16 Heliconius species. Using a reference-assisted scaffolding approach, we place and order the scaffolds of these genomes onto chromosomes, resulting in 95.7-99.9% of their genomes anchored to chromosomes. Genome sizes are somewhat variable among species (270-422 Mb) and in one small group of species (Heliconius hecale, Heliconius elevatus, and Heliconius pardalinus) expansions in genome size are driven mainly by repetitive sequences that map to four small regions in the H. melpomene reference genome. Genes from these repeat regions show an increase in exon copy number, an absence of internal stop codons, evidence of constraint on nonsynonymous changes, and increased expression, all of which suggest that at least some of the extra copies are functional. Finally, we conducted a systematic search for inversions and identified five moderately large inversions fixed between the two major Heliconius clades. We infer that one of these inversions was transferred by introgression between the lineages leading to the erato/sara and burneyi/doris clades. These reference-guided assemblies represent a major improvement in Heliconius genomic resources that enable further genetic and evolutionary discoveries in this genus.

The limits to parapatric speciation 3: evolution of strong reproductive isolation in presence of gene flow despite limited ecological differentiation

Gene flow tends to impede the accumulation of genetic divergence. Here, we determine the limits for the evolution of postzygotic reproductive isolation in a model of two populations that are connected by gene flow. We consider two selective mechanisms for the creation and maintenance of a genetic barrier: local adaptation leads to divergence among incipient species due to selection against migrants, and Dobzhansky–Muller incompatibilities (DMIs) reinforce the genetic barrier through selection against hybrids. In particular, we are interested in the maximum strength of the barrier under a limited amount of local adaptation, a challenge that many incipient species may initially face. We first confirm that with classical two-locus DMIs, the maximum amount of local adaptation is indeed a limit to the strength of a genetic barrier. However, with three or more loci and cryptic epistasis, this limit holds no longer. In particular, we identify a minimal configuration of three epistatically interacting mutations that is sufficient to confer strong reproductive isolation.

This article is part of the theme issue ‘Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers’.