Does sexual conflict contribute to the evolution of novel warning patterns?

Why warning patterns are so diverse is an enduring evolutionary puzzle. Because predators associate particular patterns with unpleasant experiences, an individual's predation risk should decrease as the local density of its warning pattern increases, promoting pattern monomorphism. Distasteful Heliconius butterflies are known for their diversity of warning patterns. Here, we explore whether interlocus sexual conflict can contribute to their diversification. Male Heliconius use warning patterns as mating cues, but mated females may suffer costs if this leads to disturbance, favouring novel patterns. Using simulations, we show that under our model conditions drift alone is unlikely to cause pattern diversification, but that sexual conflict can assist such a process. We also find that genetic architecture influences the evolution of male preferences, which track changes in warning pattern due to sexual selection. When male attraction imposes costs on females, this affects the speed at which novel pattern alleles increase. In two experiments, females laid fewer eggs with males present. However, although males in one experiment showed less interest in females with manipulated patterns, we found no evidence that female colouration mitigates sex-specific costs. Overall, male attraction to conspecific warning patterns may impose an unrecognized cost on Heliconius females, but further work is required to determine this experimentally.

Light environment influences mating behaviours during the early stages of divergence in tropical butterflies.

Speciation is facilitated when traits under divergent selection also act as mating cues. Fluctuations in sensory conditions can alter signal perception independently of adaptation to the broader sensory environment, but how this fine-scale variation may constrain or promote behavioural isolation has received little attention. The warning patterns of Heliconius butterflies are under selection for aposematism and act as mating cues. Using computer vision, we extracted behavioural data from 1481 h of video footage, for 387 individuals. We show that the putative hybrid species H. heurippa and its close relative H. timareta linaresi differ in their response to divergent warning patterns, but that these differences are strengthened with increased local illuminance. Trials with live individuals reveal low-level assortative mating that is sufficiently explained by differences in visual attraction. Finally, results from hybrid butterflies are consistent with linkage between a major warning pattern gene and the corresponding behaviour, though the differences in behaviour we observe are unlikely to cause rapid reproductive isolation as predicted under a model of hybrid trait speciation. Overall, our results reveal that the contribution of ecological mating cues to reproductive isolation may depend on the immediate sensory conditions during which they are displayed to conspecifics.

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection.

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping-by-sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies.

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole-genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.

The Scent Chemistry of Heliconius Wing Androconia.

Neotropical Heliconius butterflies are members of various mimicry rings characterized by diverse colour patterns. In the present study we investigated whether a similar diversity is observed in the chemistry of volatile compounds present in male wing androconia. Recent research has shown that these androconia are used during courting of females. Three to five wild-caught male Heliconius individuals of 17 species and subspecies were analyzed by GC/MS. Most of the identified compounds originate from common fatty acids precursors, including aldehydes, alcohols, acetates or esters preferentially with a C18 and C20 chain, together with some alkanes. The compounds occurred in species-specific mixtures or signatures. For example, octadecanal is characteristic for H. melpomene, but variation in composition between the individuals was observed. Cluster analysis of compound occurrence in individual bouquets and analyses based on biosynthetic motifs such as functional group, chain length, or basic carbon-backbone modification were used to reveal structural patterns. Mimetic pairs contain different scent bouquets, but also some compounds in common, whereas sympatric species, both mimetic and non-mimetic, have more distinct compound compositions. The compounds identified here may play a role in mate choice thus helping maintain species integrity in a butterfly genus characterized by pervasive interspecific gene flow.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification.

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Wing patterning gene redefines the mimetic history of Heliconius butterflies.

The mimetic butterflies Heliconius erato and Heliconius melpomene have undergone parallel radiations to form a near-identical patchwork of over 20 different wing-pattern races across the Neotropics. Previous molecular phylogenetic work on these radiations has suggested that similar but geographically disjunct color patterns arose multiple times independently in each species. The neutral markers used in these studies, however, can move freely across color pattern boundaries, and therefore might not represent the history of the adaptive traits as accurately as markers linked to color pattern genes. To assess the evolutionary histories across different loci, we compared relationships among races within H. erato and within H. melpomene using a series of unlinked genes, genes linked to color pattern loci, and optix, a gene recently shown to control red color-pattern variation. We found that although unlinked genes partition populations by geographic region, optix had a different history, structuring lineages by red color patterns and supporting a single origin of red-rayed patterns within each species. Genes closely linked (80-250 kb) to optix exhibited only weak associations with color pattern. This study empirically demonstrates the necessity of examining phenotype-determining genomic regions to understand the history of adaptive change in rapidly radiating lineages. With these refined relationships, we resolve a long-standing debate about the origins of the races within each species, supporting the hypothesis that the red-rayed Amazonian pattern evolved recently and expanded, causing disjunctions of more ancestral patterns.

Adaptive introgression of a visual preference gene.

Visual preferences are important drivers of mate choice and sexual selection, but little is known of how they evolve at the genetic level. In this study, we took advantage of the diversity of bright warning patterns displayed by Heliconius butterflies, which are also used during mate choice. Combining behavioral, population genomic, and expression analyses, we show that two Heliconius species have evolved the same preferences for red patterns by exchanging genetic material through hybridization. Neural expression of regucalcin1 correlates with visual preference across populations, and disruption of regucalcin1 with CRISPR-Cas9 impairs courtship toward conspecific females, providing a direct link between gene and behavior. Our results support a role for hybridization during behavioral evolution and show how visually guided behaviors contributing to adaptation and speciation are encoded within the genome.

Genetic evidence for hybrid trait speciation in heliconius butterflies.

Homoploid hybrid speciation is the formation of a new hybrid species without change in chromosome number. So far, there has been a lack of direct molecular evidence for hybridization generating novel traits directly involved in animal speciation. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has been proposed as a hybrid species, and its color pattern can be recreated by introgression of the H. m. melpomene red band into the genetic background of the yellow banded H. cydno cordula. This hybrid color pattern is also involved in mate choice and leads to reproductive isolation between H. heurippa and its close relatives. Here, we provide molecular evidence for adaptive introgression by sequencing genes across the Heliconius red band locus and comparing them to unlinked wing patterning genes in H. melpomene, H. cydno, and H. heurippa. 670 SNPs distributed among 29 unlinked coding genes (25,847bp) showed H. heurippa was related to H. c. cordula or the three species were intermixed. In contrast, among 344 SNPs distributed among 13 genes in the red band region (18,629bp), most showed H. heurippa related with H. c. cordula, but a block of around 6,5kb located in the 3' of a putative kinesin gene grouped H. heurippa with H. m. melpomene, supporting the hybrid introgression hypothesis. Genealogical reconstruction showed that this introgression occurred after divergence of the parental species, perhaps around 0.43Mya. Expression of the kinesin gene is spatially restricted to the distal region of the forewing, suggesting a mechanism for pattern regulation. This gene therefore constitutes the first molecular evidence for adaptive introgression during hybrid speciation and is the first clear candidate for a Heliconius wing patterning locus.

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone.

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.

Speciation by hybridization in Heliconius butterflies.

Speciation is generally regarded to result from the splitting of a single lineage. An alternative is hybrid speciation, considered to be extremely rare, in which two distinct lineages contribute genes to a daughter species. Here we show that a hybrid trait in an animal species can directly cause reproductive isolation. The butterfly species Heliconius heurippa is known to have an intermediate morphology and a hybrid genome, and we have recreated its intermediate wing colour and pattern through laboratory crosses between H. melpomene, H. cydno and their F1 hybrids. We then used mate preference experiments to show that the phenotype of H. heurippa reproductively isolates it from both parental species. There is strong assortative mating between all three species, and in H. heurippa the wing pattern and colour elements derived from H. melpomene and H. cydno are both critical for mate recognition by males.

Deep mitochondrial divergence within a Heliconius butterfly species is not explained by cryptic speciation or endosymbiotic bacteria.

BACKGROUND: Cryptic population structure can be an indicator of incipient speciation or historical processes. We investigated a previously documented deep break in the mitochondrial haplotypes of Heliconius erato chestertonii to explore the possibility of cryptic speciation, and also the possible presence of endosymbiont bacteria that might drive mitochondrial population structure. RESULTS: Among a sample of 315 individuals from 16 populations of western Colombia, two principal mtDNA clades were detected with 2.15% divergence and we confirmed this structure was weakly associated with geography. The first mtDNA clade included 87% of individuals from northern populations and was the sister group of H. erato members of Andes western, while the second clade contained most individuals from southern populations (78%), which shared haplotypes with an Ecuadorian race of H. erato. In contrast, analysis using AFLP markers showed H. e. chestertonii to be a genetically homogeneous species with no association between mitochondrial divergence and AFLP structure. The lack of congruence between molecular markers suggests that cryptic speciation is not a plausible explanation for the deep mitochondrial divergence in H. e chestertonii. We also carried out the first tests for the presence of endosymbiontic bacteria in Heliconius, and identified two distinct lineages of Wolbachia within H. e. chestertonii. However, neither of the principal mitochondrial clades of H. e. chestertonii was directly associated with the patterns of infection. CONCLUSIONS: We conclude that historical demographic processes are the most likely explanation for the high mitochondrial differentiation in H. e. chestertonii, perhaps due to gene flow between Cauca valley H. e. chestertonii and west Pacific slope populations of H. erato.

Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range.

In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter- and intraspecific levels in mimetic Heliconius butterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified "indicator" compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.

Interplay between Developmental Flexibility and Determinism in the Evolution of Mimetic Heliconius Wing Patterns.

To what extent can we predict how evolution occurs? Do genetic architectures and developmental processes canalize the evolution of similar outcomes in a predictable manner? Or do historical contingencies impose alternative pathways to answer the same challenge? Examples of Müllerian mimicry between distantly related butterfly species provide natural replicates of evolution, allowing us to test whether identical wing patterns followed parallel or novel trajectories. Here, we explore the role that the signaling ligand WntA plays in generating mimetic wing patterns in Heliconius butterflies, a group with extraordinary mimicry-related wing pattern diversity. The radiation is relatively young, and numerous cases of wing pattern mimicry have evolved within the last 2.5-4.5 Ma. WntA is an important target of natural selection and is one of four major effect loci that underlie much of the pattern variation in the group. We used CRISPR/Cas9 targeted mutagenesis to generate WntA-deficient wings in 12 species and a further 10 intraspecific variants, including three co-mimetic pairs. In all tested butterflies, WntA knockouts affect pattern broadly and cause a shift among every possible scale cell type. Interestingly, the co-mimics lacking WntA were very different, suggesting that the gene networks that pattern a wing have diverged considerably among different lineages. Thus, although natural selection channeled phenotypic convergence, divergent developmental contexts between the two major Heliconius lineages opened different developmental routes to evolve resemblance. Consequently, even under very deterministic evolutionary scenarios, our results underscore a surprising unpredictability in the developmental paths underlying convergence in a recent radiation.

Gene flow and the genealogical history of Heliconius heurippa.

BACKGROUND: The neotropical butterfly Heliconius heurippa has a hybrid colour pattern, which also contributes to reproductive isolation, making it a likely example of hybrid speciation. Here we used phylogenetic and coalescent-based analyses of multilocus sequence data to investigate the origin of H. heurippa. RESULTS: We sequenced a mitochondrial region (CoI and CoII), a sex-linked locus (Tpi) and two autosomal loci (w and sd) from H. heurippa and the putative parental species, H. cydno and H. melpomene. These were analysed in combination with data from two previously sequenced autosomal loci, Dll and Inv. H. heurippa was monophyletic at mtDNA and Tpi, but showed a shared distribution of alleles derived from both parental lineages at all four autosomal loci. Estimates of genetic differentiation showed that H. heurippa is closer to H. cydno at mtDNA and three autosomal loci, intermediate at Tpi, and closer to H. melpomene at Dll. Using coalescent simulations with the Isolation-Migration model (IM), we attempted to establish the incidence of gene flow in the origin of H. heurippa. This analysis suggested that ongoing introgression is frequent between all three species and variable in extent between loci. CONCLUSION: Introgression, which is a necessary precursor of hybrid speciation, seems to have also blurred the coalescent history of these species. The origin of Heliconius heurippa may have been restricted to introgression of few colour pattern genes from H. melpomene into the H. cydno genome, with little evidence of genomic mosaicism.

Two sisters in the same dress: Heliconius cryptic species.

BACKGROUND: Sister species divergence and reproductive isolation commonly results from ecological adaptation. In mimetic Heliconius butterflies, shifts in colour pattern contribute to pre- and post-mating reproductive isolation and are commonly correlated with speciation. Closely related mimetic species are therefore not expected, as they should lack several important sources of reproductive isolation. RESULTS: Here we present phenotypic, behavioral and genetic evidence for the coexistence of two sympatric 'cryptic' species near Florencia in the eastern Andes of Colombia that share the same orange rayed colour pattern. These represent H. melpomene malleti and a novel taxon in the H. cydno group, here designated as novel race of Heliconius timareta, Heliconius timareta florencia. No-choice mating experiments show that these sympatric forms have strong assortative mating ( approximately 96%) despite great similarity in colour pattern, implying enhanced divergence in pheromonal signals. CONCLUSION: We hypothesize that these species might have resulted from recent convergence in colour pattern, perhaps facilitated by hybrid introgression of wing pattern genes.

A hybrid zone provides evidence for incipient ecological speciation in Heliconius butterflies.

In Heliconius butterflies, it has been proposed that speciation occurs through a combination of divergence in ecological habitat preferences and mimetic colour patterns. Here we test this hypothesis by investigating a parapatric form of the widespread species Heliconius erato. Mendelian (colour patterns) and molecular genetic data permit us to address hypotheses about introgression and genetic differentiation between different populations. Combined analysis of colour pattern, microsatellite loci and mitochondrial DNA showed that Heliconius erato venus and Heliconius erato chestertonii form a bimodal hybrid zone implying partial reproductive isolation. In a sample of 121 individuals collected in sympatry, 25% were hybrids representing a significant deficit of heterozygotes compared to the Hardy-Weinberg expectation. Seven microsatellite loci, analysed for a subset of these individuals, showed marked differentiation between the parental taxa, and unambiguously identified two genotypic clusters concordant with our phenotypic classification of individuals. Mitochondrial DNA analysis showed H. erato venus as a monophyletic group well differentiated from H. erato chestertonii, implying a lack of historical introgression between the populations. Heliconius erato chestertonii is therefore an incipient species that maintains its integrity despite high levels of hybridization. Moreover, H. erato chestertonii is found at higher altitudes than other races of H. erato and has a distinct colour pattern and mimetic relationship. Hence, there are now two examples of parapatric incipient species related to H. erato, H. himera and H. erato chestertonii, both of which are associated with higher altitudes, more arid habitats and distinct mimetic relationships. This implies that parapatric habitat adaptation is a likely cause of speciation in this group.

Regulation of testosterone degradation in Comamonas testosteroni.

Recently, we have identified a gene encoding a LuxR-type factor, TeiR (Testosterone-inducible Regulator), which positively regulates steroid degradation in Comamonas testosteroni. Herein, we demonstrate that TeiR interacts in vivo with steroid catabolic gene promoters. The presence of testosterone induces a significant TeiR protein increase at the early logarithmic phase of growth. Interestingly, it is not until the early stationary phase where the activation of a steroid-inducible gene promoter is observed, indicating that testosterone might not be the true inductor of the steroid degradation pathway. In addition, beta-galactosidase expression driven by a testosterone-inducible promoter is prematurely activated in cells cultured in medium supplemented with ethyl acetate extracts obtained from the early stationary phase cell-free supernatants of C. testosteroni grown in presence of testosterone. Complementation experiments of C. testosteroni wild type performed with teiR deletion constructs indicate that extra-copies of deleted-TeiR exert a dominant negative effect on the wild-type TeiR protein. While, when C. testosteroni teiR mutants were used to carry out complementation assays only the full length gene can overcome the teiR mutant phenotype. Altogether these findings indicate that TeiR regulates steroid catabolic genes interacting with their promoters and suggest that this interaction requires the presence of a testosterone-derived metabolite to induce the system.

Natural hybridization in heliconiine butterflies: the species boundary as a continuum.

BACKGROUND: To understand speciation and the maintenance of taxa as separate entities, we need information about natural hybridization and gene flow among species. RESULTS: Interspecific hybrids occur regularly in Heliconius and Eueides (Lepidoptera: Nymphalidae) in the wild: 26-29% of the species of Heliconiina are involved, depending on species concept employed. Hybridization is, however, rare on a per-individual basis. For one well-studied case of species hybridizing in parapatric contact (Heliconius erato and H. himera), phenotypically detectable hybrids form around 10% of the population, but for species in sympatry hybrids usually form less than 0.05% of individuals. There is a roughly exponential decline with genetic distance in the numbers of natural hybrids in collections, both between and within species, suggesting a simple "exponential failure law" of compatibility as found in some prokaryotes. CONCLUSION: Hybridization between species of Heliconius appears to be a natural phenomenon; there is no evidence that it has been enhanced by recent human habitat disturbance. In some well-studied cases, backcrossing occurs in the field and fertile backcrosses have been verified in insectaries, which indicates that introgression is likely, and recent molecular work shows that alleles at some but not all loci are exchanged between pairs of sympatric, hybridizing species. Molecular clock dating suggests that gene exchange may continue for more than 3 million years after speciation. In addition, one species, H. heurippa, appears to have formed as a result of hybrid speciation. Introgression may often contribute to adaptive evolution as well as sometimes to speciation itself, via hybrid speciation. Geographic races and species that coexist in sympatry therefore form part of a continuum in terms of hybridization rates or probability of gene flow. This finding concurs with the view that processes leading to speciation are continuous, rather than sudden, and that they are the same as those operating within species, rather than requiring special punctuated effects or complete allopatry. Although not qualitatively distinct from geographic races, nor "real" in terms of phylogenetic species concepts or the biological species concept, hybridizing species of Heliconius are stably distinct in sympatry, and remain useful groups for predicting morphological, ecological, behavioural and genetic characteristics.