Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies

VOC Characterization of Byasa hedistus (Lepidoptera: Papilionidae) and Its Visual and Olfactory Responses during Foraging and Courtship

Color and odor are crucial cues for butterflies during foraging and courtship. While most sexual dimorphic butterflies rely more on vision, our understanding of how butterflies with similar coloration use different signals remains limited. This study investigated the visual and olfactory behavioral responses of the similarly colored butterfly Byasa hedistus during foraging and courtship. While visiting artificial flowers of different colors, we found that B. hedistus exhibits an innate color preference, with a sequence of preferences for red, purple, and blue. The frequency of flower visits by B. hedistus significantly increased when honey water was sprayed on the artificial flowers, but it hardly visited apetalous branches with honey water. This proves that locating nectar sources by odor alone is difficult in the absence of floral color guides. During courtship, males are active while females hardly chase; only two models were observed: males chasing males and males chasing females. The courtship process includes four behaviors: slowing approach, straight chasing, hovering, and spinning. B. hedistus cannot distinguish between sexes based on color, as there is no significant difference in color and shape between them. Twenty-three VOCs (>1%) were identified in B. hedistus, with 21 shared by both sexes, while ketones are specific to males. These VOCs are principally represented by cineole, β-pinene, and linalool. When cineole was added to butterfly mimics, many butterflies were attracted to them, but the butterflies did not seem to distinguish between males and females. This suggests that cineole may be the feature VOC for identifying conspecific groups. Adding β-pinene and linalool to mimics induced numerous butterflies to chase, hover, spin around, and attempt to mate with them. This suggests that β-pinene and linalool are crucial cues indicating the presence of females.

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.

Chemical species recognition in an adaptive radiation of Hawaiian Tetragnatha spiders (Araneae: Tetragnathidae)

Studies of adaptive radiations have played a central role in our understanding of reproductive isolation. Yet the focus has been on human-biased visual and auditory signals, leaving gaps in our knowledge of other modalities. To date, studies on chemical signals in adaptive radiations have focused on systems with multimodal signalling, making it difficult to isolate the role chemicals play in reproductive isolation. In this study we examine the use of chemical signals in the species recognition and adaptive radiation of Hawaiian Tetragnatha spiders by focusing on entire communities of co-occurring species, and conducting behavioural assays in conjunction with chemical analysis of their silks using gas chromatography-mass spectrometry. Male spiders significantly preferred the silk extracts of conspecific mates over those of sympatric heterospecifics. The compounds found in the silk extracts, long chain alkyl methyl ethers, were remarkably species-specific in the combination and quantity. The differences in the profile were greatest between co-occurring species and between closely related sibling species. Lastly, there were significant differences in the chemical profile between two populations of a particular species. These findings provide key insights into the role chemical signals play in the attainment and maintenance of reproductive barriers between closely related co-occurring species.

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.

Variation of Clasper Scent Gland Composition of Heliconius Butterflies from a Biodiversity Hotspot

Male Heliconius butterflies possess two pheromone emitting structures, wing androconia and abdominal clasper scent glands. The composition of the clasper scent gland of males of 17 Heliconius and Eueides species from an overlapping area in Ecuador, comprising three mimicry groups, was investigated by GC/MS. The chemical signal serves as an anti-aphrodisiac signal that is transferred from males to females during mating, indicating the mating status of the female to prevent them from harassment by other males. In addition, the odour may also serve in predator defence. There is potential for convergence driven by mimicry, although, such convergence might be detrimental for species recognition of the butterflies within the mimicry ring, making mating more difficult. More than 500 compounds were detected, consisting of volatile, semi-volatile or non-volatile compounds, including terpenes, fatty acid esters or aromatic compounds. Several novel esters were identified by GC/MS and GC/IR data, microderivatisation and synthesis, including butyl (Z)-3-dodecenoate and other (Z)-3-alkenoates, 3-oxohexyl citronellate and 5-methylhexa-3,5-dienyl (E)-2,3-dihydrofarnesoate. The secretions were found to be species specific, potentially allowing for species differentiation. Statistical analysis of the compounds showed differentiation by phylogenetic clade and species, but not by mimicry group.

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Evolutionary origins of sexual dimorphism: Lessons from female-limited mimicry in butterflies

The striking female-limited mimicry observed in some butterfly species is a text-book example of sexually dimorphic trait submitted to intense natural selection. Two main evolutionary hypotheses, based on natural and sexual selection respectively, have been proposed. Predation pressure favoring mimicry toward defended species could be higher in females because of their slower flight, and thus overcome developmental constraints favoring the ancestral trait that limits the evolution of mimicry in males but not in females. Alternatively, the evolution of mimicry in males could be limited by female preference for non-mimetic males. However, the evolutionary origin of female preference for non-mimetic males remains unclear. Here, we hypothesize that costly sexual interactions between individuals from distinct sympatric species might intensify because of mimicry, therefore promoting female preference for non-mimetic trait. Using a mathematical model, we compare the evolution of female-limited mimicry when assuming either alternative selective hypotheses. We show that the patterns of divergence of male and female trait from the ancestral traits can differ between these selection regimes. We specifically highlight that divergence in female trait is not a signature of the effect of natural selection. Our results also evidence why female-limited mimicry is more frequently observed in Batesian mimics.

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.

Chemical identity of cuticular lipid components in the mimetic swallowtail butterfly Papilio polytes

The swallowtail Papilio polytes shows Batesian and female-limited polymorphic mimicry. In Japan, P . polytes females have two different forms: the cyrus form is non-mimetic and resembles males, whereas the polytes form mimics Pachliopta aristolochiae and Byasa ( Atrophaneura ) alcinous as unpalatable models. During mating, P . polytes males use cuticular lipids to distinguish non-mimetic females from conspecific males and sympatric sister species. In this study, we investigated whether compositional differences in cuticular lipids exist between mimetic and non-mimetic females of P . polytes and between mimetic females and their model species. The mimetic and non-mimetic females had nearly identical cuticular lipid profiles, which differed from those of males. The two model species exhibited sexually dimorphic and species-specific cuticular lipid compositions, which were distinctly different from those of mimetic P . polytes females. These results strongly suggest that P . polytes females maintain the identity of cuticular lipid profiles regardless of the mimicry type, and this feature helps males recognize mimetic females as the correct mating partners.

Convergent morphology and divergent phenology promote the coexistence of Morpho butterfly species

The coexistence of closely-related species in sympatry is puzzling because ecological niche proximity imposes strong competition and reproductive interference. A striking example is the widespread wing pattern convergence of several blue-banded Morpho butterfly species with overlapping ranges of distribution. Here we perform a series of field experiments using flying Morpho dummies placed in a natural habitat. We show that similarity in wing colour pattern indeed leads to interspecific territoriality and courtship among sympatric species. In spite of such behavioural interference, demographic inference from genomic data shows that sympatric closely-related Morpho species are genetically isolated. Mark-recapture experiments in the two most closely-related species unravel a strong temporal segregation in patrolling activity of males. Such divergence in phenology reduces the costs of reproductive interference while simultaneously preserving the benefits of convergence in non-reproductive traits in response to common ecological pressures. Henceforth, the evolution of multiple traits may favour species diversification in sympatry by partitioning niche in different dimensions.

Identification and Composition of Clasper Scent Gland Components of the Butterfly Heliconius erato and Its Relation to Mimicry

The butterfly Heliconius erato occurs in various mimetic morphs. The male clasper scent gland releases an anti‐aphrodisiac pheromone and additionally contains a complex mixture of up to 350 components, varying between individuals. In 114 samples of five different mimicry groups and their hybrids 750 different compounds were detected by gas chromatography/mass spectrometry (GC/MS). Many unknown components occurred, which were identified using their mass spectra, gas chromatography/infrared spectroscopy (GC/IR)‐analyses, derivatization, and synthesis. Key compounds proved to be various esters of 3‐oxohexan‐1‐ol and (Z)‐3‐hexen‐1‐ol with (S)‐2,3‐dihydrofarnesoic acid, accompanied by a large variety of other esters with longer terpene acids, fatty acids, and various alcohols. In addition, linear terpenes with up to seven uniformly connected isoprene units occur, e. g. farnesylfarnesol. A large number of the compounds have not been reported before from nature. Discriminant analyses of principal components of the gland contents showed that the iridescent mimicry group differs strongly from the other, mostly also separated, mimicry groups. Comparison with data from other species indicated that Heliconius recruits different biosynthetic pathways in a species‐specific manner for semiochemical formation.

Neural basis of acoustic species recognition in a cryptic species complex

Sexual traits that promote species recognition are important drivers of reproductive isolation, especially among closely related species. Identifying neural processes that shape species differences in recognition is crucial for understanding the causal mechanisms of reproductive isolation. Temporal patterns are salient features of sexual signals that are widely used in species recognition by several taxa, including anurans. Recent advances in our understanding of temporal processing by the anuran auditory system provide an opportunity to investigate the neural basis of species-specific recognition. The anuran inferior colliculus consists of neurons that are selective for temporal features of calls. Of potential relevance are auditory neurons known as interval-counting neurons (ICNs) that are often selective for the pulse rate of conspecific advertisement calls. Here, we tested the hypothesis that ICNs mediate acoustic species recognition by exploiting the known differences in temporal selectivity in two cryptic species of gray treefrog (Hyla chrysoscelis and Hyla versicolor). We examined the extent to which the threshold number of pulses required to elicit behavioral responses from females and neural responses from ICNs was similar within each species but potentially different between the two species. In support of our hypothesis, we found that a species difference in behavioral pulse number thresholds closely matched the species difference in neural pulse number thresholds. However, this relationship held only for ICNs that exhibited band-pass tuning for conspecific pulse rates. Together, these findings suggest that differences in temporal processing of a subset of ICNs provide a mechanistic explanation for reproductive isolation between two cryptic treefrog species.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Evolutionary trade-offs between male secondary sexual traits revealed by a phylogeny of the hyperdiverse tribe Eumaeini (Lepidoptera: Lycaenidae)

Male butterflies in the hyperdiverse tribe Eumaeini possess an unusually complex and diverse repertoire of secondary sexual characteristics involved in pheromone production and dissemination. Maintaining multiple sexually selected traits is likely to be metabolically costly, potentially resulting in trade-offs in the evolution of male signals. However, a phylogenetic framework to test hypotheses regarding the evolution and maintenance of male sexual traits in Eumaeini has been lacking. Here, we infer a comprehensive, time-calibrated phylogeny from 379 loci for 187 species representing 91% of the 87 described genera. Eumaeini is a monophyletic group that originated in the late Oligocene and underwent rapid radiation in the Neotropics. We examined specimens of 818 of the 1096 described species (75%) and found that secondary sexual traits are present in males of 91% of the surveyed species. Scent pads and scent patches on the wings and brush organs associated with the genitalia were probably present in the common ancestor of Eumaeini and are widespread throughout the tribe. Brush organs and scent pads are negatively correlated across the phylogeny, exhibiting a trade-off in which lineages with brush organs are unlikely to regain scent pads and vice versa. In contrast, scent patches seem to facilitate the evolution of scent pads, although they are readily lost once scent pads have evolved. Our results illustrate the complex interplay between natural and sexual selection in the origin and maintenance of multiple male secondary sexual characteristics and highlight the potential role of sexual selection spurring diversification in this lineage.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

3-Acetoxy-fatty acid isoprenyl esters from androconia of the ithomiine butterfly Ithomia salapia

Male ithomiine butterflies (Nymphalidae: Danainae) have hairpencils on the forewings (i.e., androconia) that disseminate semiochemicals during courtship. While most ithomiines are known to contain derivatives of pyrrolizidine alkaloids, dihydropyrrolizines, or γ-lactones in these androconia, here we report on a new class of fatty acid esters identified in two subspecies, Ithomia salapia aquinia and I. s. derasa. The major components were identified as isoprenyl (3-methyl-3-butenyl) (Z)-3-acetoxy-11-octadecenoate, isoprenyl (Z)-3-acetoxy-13-octadecenoate (12) and isoprenyl 3-acetoxyoctadecanoate (11) by GC/MS and GC/IR analyses, microderivatizations, and synthesis of representative compounds. The absolute configuration of 12 was determined to be R. The two subspecies differed not only in the composition of the ester bouquet, but also in the composition of more volatile androconial constituents. While some individuals of I. s. aquinia contained ithomiolide A (3), a pyrrolizidine alkaloid derived γ-lactone, I. s. derasa carried the sesquiterpene α-elemol (8) in the androconia. These differences might be important for the reproductive isolation of the two subspecies, in line with previously reported low gene exchange between the two species in regions where they co-occur. Furthermore, the occurrence of positional isomers of unsaturated fatty acid derivatives indicates activity of two different desaturases within these butterflies, Δ9 and Δ11, which has not been reported before in male Lepidoptera.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Divergence of chemosensing during the early stages of speciation

Chemosensory communication is essential to insect biology, playing indispensable roles during mate-finding, foraging, and oviposition behaviors. These traits are particularly important during speciation, where chemical perception may serve to establish species barriers. However, identifying genes associated with such complex behavioral traits remains a significant challenge. Through a combination of transcriptomic and genomic approaches, we characterize the genetic architecture of chemoperception and the role of chemosensing during speciation for a young species pair of Heliconius butterflies, Heliconius melpomene and Heliconius cydno We provide a detailed description of chemosensory gene-expression profiles as they relate to sensory tissue (antennae, legs, and mouthparts), sex (male and female), and life stage (unmated and mated female butterflies). Our results untangle the potential role of chemical communication in establishing barriers during speciation and identify strong candidate genes for mate and host plant choice behaviors. Of the 252 chemosensory genes, HmOBP20 (involved in volatile detection) and HmGr56 (a putative synephrine-related receptor) emerge as strong candidates for divergence in pheromone detection and host plant discrimination, respectively. These two genes are not physically linked to wing-color pattern loci or other genomic regions associated with visual mate preference. Altogether, our results provide evidence for chemosensory divergence between H. melpomene and H. cydno, two rarely hybridizing butterflies with distinct mate and host plant preferences, a finding that supports a polygenic architecture of species boundaries.