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.

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.

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.

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.

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.

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.

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.

Complex modular architecture around a simple toolkit of wing pattern genes.

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

An introgressed wing pattern acts as a mating cue.

Heliconius butterflies provide good examples of both homoploid hybrid speciation and ecological speciation. In particular, examples of adaptive introgression have been detected among the subspecies of Heliconius timareta, which acquired red color pattern elements from H. melpomene. We tested whether the introgression of red wing pattern elements into H. timareta florencia might also be associated with incipient reproductive isolation (RI) from its close relative, H. timareta subsp. nov., found in the eastern Andes. No choice experiments show a 50% reduction in mating between females of H. t. subsp. nov. and males of H .t. florencia, but not in the reciprocal direction. In choice experiments using wing models, males of H. timareta subsp. nov. approach and court red phenotypes less than their own, whereas males of H. t. florencia prefer models with a red phenotype. Intrinsic postzygotic isolation was not detected in crosses between these H. timareta races. These results suggest that a color pattern trait gained by introgression is triggering RI between H. timareta subsp. nov. and H. t. florencia.

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.

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.

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.

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.

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.

Assortative mating preferences among hybrids offers a route to hybrid speciation.

Homoploid speciation generates species without a change in chromosome number via introgressive hybridization and has been considered rare in animals. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has a color pattern that can be recreated by introgression of the H. melpomene red band into an H. cydno genetic background. Wild H. heurippa males show assortative mating based on color pattern and we here investigate the origin of this preference by studying first-generation backcross hybrids between H. melpomene and H. cydno that resemble H. heurippa. These hybrids show assortative mating preferences, showing a strong preference for their own color pattern over that of either parental species. This is consistent with a genetic basis to wing pattern preference and implies, first, that assortative mating preferences would facilitate the initial establishment of a homozygous hybrid color pattern by increasing the likelihood that early generation hybrids mate among themselves. Second, once established such a lineage would inherit assortative mating preferences that would lead to partial reproductive isolation from parental lineages.

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.

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.