Adaptive introgression across species boundaries in Heliconius butterflies

On the use of genome-wide data to model and date the time of anthropogenic hybridisation: An example from the Scottish wildcat

While hybridisation has long been recognised as an important natural phenomenon in evolution, the conservation of taxa subject to introgressive hybridisation from domesticated forms is a subject of intense debate. Hybridisation of Scottish wildcats and domestic cats is a good example in this regard. Here, we developed a modelling framework to determine the timescale of introgression using approximate Bayesian computation (ABC). Applying the model to ddRAD-seq data from 129 individuals, genotyped at 6546 loci, we show that a population of wildcats genetically distant from domestic cats is still present in Scotland. These individuals were found almost exclusively within the captive breeding programme. Most wild-living cats sampled were introgressed to some extent. The demographic model predicts high levels of gene-flow between domestic cats and Scottish wildcats (13% migrants per generation) over a short timeframe, the posterior mean for the onset of hybridisation (T₁ ) was 3.3 generations (~10 years) before present. Although the model had limited power to detect signals of ancient admixture, we found evidence that significant recent hybridisation may have occurred subsequent to the founding of the captive breeding population (T₂ ). The model consistently predicts T₁ after T₂ , estimated here to be 19.3 generations (~60 years) ago, highlighting the importance of this population as a resource for conservation management. Additionally, we evaluate the effectiveness of current methods to classify hybrids. We show that an optimised 35 SNP panel is a better predictor of the ddRAD-based hybrid score in comparison with a morphological method.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

The role of climate change in pollinator decline across the Northern Hemisphere is underestimated

Pollinator biodiversity loss occurs at unprecedented rates globally, with particularly sharp declines documented in the North Temperate Zone. There is currently no consensus on the main drivers of the decline. Although climate change is expected to drive biodiversity loss in the future, current warming is often suggested to have positive impacts on pollinator assemblages in higher latitudes. Consequently, pollinator conservation initiatives in Europe and the USA tend to lack climate adaptation initiatives, an omission of which may be risky if climate change has significant negative impacts on pollinators. To gain an understanding of the impacts of climate change on pollinator biodiversity in the Northern Hemisphere, we conducted a literature review on genetic, species and community level diversity. Our findings suggest that global heating most likely causes homogenization of pollinator assemblages at all levels of pollinator biodiversity, making them less resilient to future stochasticity. Aspects of biodiversity that are rarely measured (e.g. genetic diversity, β-diversity, species evenness) tend to be most affected, while some dimensions of climate change, such as fluctuations in winter weather conditions, changes in the length of the vegetational season and increased frequency of extreme weather events, that seldom receive attention in empirical studies, tend to be particularly detrimental to pollinators. Negative effects of global heating on pollinator biodiversity are most likely exacerbated by homogenous and fragmented landscapes, widespread across Europe and the US, which limit opportunities for range-shifts and reduce micro-climatic buffering. This suggests the need for conservation initiatives to focus on increasing landscape connectivity and heterogeneity at multiple spatial scales.

Exchanging fluids The sociocultural implications of microbial, cultural, and ethnic admixture in Latin America

Knowledge of evolutionary influences on patterns of human mating, social interactions, and differential health is increasing, yet these insights have rarely been applied to historical analyses of human population dynamics. The genetic and evolutionary forces behind biases in interethnic mating and in the health of individuals of different ethnic groups in Latin America and the Caribbean since the European colonization of America are still largely ignored. We discuss how historical and contemporary sociocultural interactions and practices are strongly influenced by population-level evolutionary forces. Specifically, we discuss the historical implications of functional (de facto) polygyny, sex-biased admixture, and assortative mating in Latin America. We propose that these three evolutionary mechanisms influenced mating patterns, shaping the genetic and cultural landscape across Latin America and the Caribbean. Further, we discuss how genetic differences between the original populations that migrated at different times into Latin America contributed to their accommodation to and survival in the different local ecologies and interethnic interactions. Relevant medical and social implications follow from the genetic and cultural changes reviewed.

Genomic analyses overturn two long-standing homoploid hybrid speciation hypotheses

The importance of hybridization in generating biological diversity has been historically controversial. Previously, inference about hybridization was limited by dependence on morphological data; with the advent of the next-generation sequencing tools for nonmodel organisms, the evolutionary significance of hybridization is more evident. Here, we test classic hypotheses of hybrid origins of two species in the Phlox pilosa complex. Morphological intermediacy motivated the hypotheses that Phlox amoena lighthipei and Phlox pilosa deamii were independent homoploid hybrid lineages derived from P. amoena amoena and P. pilosa pilosa. We use double-digest restriction site-associated DNA sequencing of individuals from throughout the range of these taxa to conduct the most thorough analysis of evolutionary history in this system to date. Surprisingly, we find no support for the hybrid origin of P. pilosa deamii or P. amoena lighthipei. Our data do identify a history of admixture in individuals collected at a contemporary hybrid zone between the putative parent lineages. We show that three very different evolutionary histories, only one of which involves hybrid origin, have produced intermediate or recombinant morphological traits between P. amoena amoena and P. pilosa pilosa. Although morphological data are still an efficient means of generating hypotheses about past gene flow, genomic data are now the standard of evidence for elucidating evolutionary history.

CARACTERIZACIÓN GENÉTICA DE LA POBLACIÓN DE Heliconius sara (Nymphalidae) EN LA ISLA GORGONA, COLOMBIA

La estructura genética de poblaciones de mariposas con distribución en islas y sus pares continentales ha sido poco documentada para el neotrópico. Este estudio presenta la caracterización de una población de Heliconius sara con distribución en la Isla Gorgona, ubicada en la región del Pacífico Oriental Colombiano. Para esto se examinaron secuencias parciales de un marcador mitocondrial incluyendo información obtenida del GenBank. Se comparó la diversidad y estructura genética con sus conespecíficos continentales y también con congéneres, con los que comparte un ancestro común cercano en el clado Sapho-Sara. Para el análisis de diversidad y estructura genética se realizó un análisis molecular de varianza. Este análisis muestra que la distancia entre la población de la isla y sus pares en el continente es consistente con la variación intraespecífica observada en otras especies del género Heliconius. Para la reconstrucción de la genealogía y datación reciente en el Pleistoceno superior del grupo monofilético de secuencias de H. sara, se realizó un análisis de inferencia bayesiana, así como una de máxima verosimilitud. Del análisis demográfico se seleccionó un modelo histórico de flujo asimétrico desde la isla hacia el continente que sugiere baja resistencia de la discontinuidad geográfica a la dispersión de esta mariposa diurna desde la isla. Este es el primer estudio en examinar un posible evento de aislamiento de una población insular de mariposas en Colombia.

Detecting adaptive introgression in human evolution using convolutional neural networks

Studies in a variety of species have shown evidence for positively selected variants introduced into a population via introgression from another, distantly related population—a process known as adaptive introgression. However, there are few explicit frameworks for jointly modelling introgression and positive selection, in order to detect these variants using genomic sequence data. Here, we develop an approach based on convolutional neural networks (CNNs). CNNs do not require the specification of an analytical model of allele frequency dynamics and have outperformed alternative methods for classification and parameter estimation tasks in various areas of population genetics. Thus, they are potentially well suited to the identification of adaptive introgression. Using simulations, we trained CNNs on genotype matrices derived from genomes sampled from the donor population, the recipient population and a related non-introgressed population, in order to distinguish regions of the genome evolving under adaptive introgression from those evolving neutrally or experiencing selective sweeps. Our CNN architecture exhibits 95% accuracy on simulated data, even when the genomes are unphased, and accuracy decreases only moderately in the presence of heterosis. As a proof of concept, we applied our trained CNNs to human genomic datasets—both phased and unphased—to detect candidates for adaptive introgression that shaped our evolutionary history.

Reticulate evolution in Conidae: Evidence of nuclear and mitochondrial introgression

Conidae is a hyperdiverse family of marine snails that has many hallmarks of adaptive radiation. Hybridization and introgression may contribute to such instances of rapid diversification by generating novel gene combinations that facilitate exploitation of distinct niches. Here we evaluated whether or not these mechanisms may have contributed to the evolutionary history of a subgenus of Conidae (Virroconus). Several observations hint at evidence of past introgression for members of this group, including incongruence between phylogenetic relationships inferred from mitochondrial gene sequences and morphology and widespread sympatry of many Virroconus species in the Indo-West Pacific. We generated and analyzed transcriptome data of Virroconus species to (i) infer a robust nuclear phylogeny, (ii) assess mitochondrial and nuclear gene tree discordance, and (iii) formally test for introgression of nuclear loci. We identified introgression of mitochondrial genomes and nuclear gene regions between ancestors of one pair of Virroconus species, and mitochondrial introgression between another pair. We also found evidence of adaptive introgression of conotoxin venom loci between a third pair of species. Together, our results demonstrate that hybridization and introgression impacted the evolutionary history of Virroconus and hence may have contributed to the adaptive radiation of Conidae.

Interspecific Gene Flow and the Evolution of Specialization in Black and White Rhinoceros

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros.

Positive selection plays a major role in shaping signatures of differentiation across the genomic landscape of two independent Ficedula flycatcher species pairs

A current debate within population genomics surrounds the relevance of patterns of genomic differentiation between closely related species for our understanding of adaptation and speciation. Mounting evidence across many taxa suggests that the same genomic regions repeatedly develop elevated differentiation in independent species pairs. These regions often coincide with high gene density and/or low recombination, leading to the hypothesis that the genomic differentiation landscape mostly reflects a history of background selection, and reveals little about adaptation or speciation. A comparative genomics approach with multiple independent species pairs at a timescale where gene flow and ILS are negligible permits investigating whether different evolutionary processes are responsible for generating lineage-specific versus shared patterns of species differentiation. We use whole-genome resequencing data of 195 individuals from four Ficedula flycatcher species comprising two independent species pairs: collared and pied flycatchers, and red-breasted and taiga flycatchers. We found that both shared and lineage-specific F_ST peaks could partially be explained by selective sweeps, with recurrent selection likely to underlie shared signatures of selection, whereas indirect evidence supports a role of recombination landscape evolution in driving lineage-specific signatures of selection. This work therefore provides evidence for an interplay of positive selection and recombination to genomic landscape evolution.

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.

Disentangling Population History and Character Evolution among Hybridizing Lineages

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Range and niche expansion through multiple interspecific hybridization: a genotyping by sequencing analysis of Cherleria (Caryophyllaceae)

BACKGROUND

Cherleria (Caryophyllaceae) is a circumboreal genus that also occurs in the high mountains of the northern hemisphere. In this study, we focus on a clade that diversified in the European High Mountains, which was identified using nuclear ribosomal (nrDNA) sequence data in a previous study. With the nrDNA data, all but one species was monophyletic, with little sequence variation within most species. Here, we use genotyping by sequencing (GBS) data to determine whether the nrDNA data showed the full picture of the evolution in the genomes of these species.

RESULTS

The overall relationships found with the GBS data were congruent with those from the nrDNA study. Most of the species were still monophyletic and many of the same subclades were recovered, including a clade of three narrow endemic species from Greece and a clade of largely calcifuge species. The GBS data provided additional resolution within the two species with the best sampling, C. langii and C. laricifolia, with structure that was congruent with geography. In addition, the GBS data showed significant hybridization between several species, including species whose ranges did not currently overlap.

CONCLUSIONS

The hybridization led us to hypothesize that lineages came in contact on the Balkan Peninsula after they diverged, even when those lineages are no longer present on the Balkan Peninsula. Hybridization may also have helped lineages expand their niches to colonize new substrates and different areas. Not only do genome-wide data provide increased phylogenetic resolution of difficult nodes, they also give evidence for a more complex evolutionary history than what can be depicted by a simple, branching phylogeny.

Unveiling cryptic diversity among Müllerian co-mimics: insights from the Western Palaearctic Syntomis moths (Lepidoptera: Erebidae: Arctiinae)

Accurate species delimitation is of primary importance in biodiversity assessments and in reconstructing patterns and processes in the diversification of life. However, the discovery of cryptic species in virtually all taxonomic groups unveiled significant gaps in our knowledge of biodiversity. Mimicry complexes are good candidates to source for cryptic species. Indeed, members of mimicry complexes undergo selective pressures on their habitus, which results in strong resemblance even between distantly related species. In this study, we used a multi-locus genetic approach to investigate the presence of cryptic diversity within a group of mimetic day-flying moths whose systematics has long been controversial, the Euro-Anatolian Syntomis. Results showed incongruence between species boundaries and the currently accepted taxonomy of this group. Both mitochondrial and nuclear markers indicate the presence of four, well-distinct genetic lineages. The genetic distance and time of divergence between the Balkan and Italian populations of S. marjana are the same as those found between S. phegea and S. ragazzii, the last two being well-distinct, broadly sympatrically occurring species. The divergence between the two lineages of S. marjana dates back to the Early Pleistocene, which coincided with substantial changes in climatic conditions and vegetation cover in Southern Europe that have likely induced geographic and ecological vicariance. Syntomis populations belonging to the taxa kruegeri (s. str.), albionica and quercii are now considered a separate species from marjana s. str. and are thus distinguished as Syntomis quercii Verity, 1914, bona sp., stat. nov. Our results show that the species richness of mimicry complexes inhabiting temperate regions might still be severely underestimated.

Population Transcriptomics Reveals Gene Flow and Introgression Between Two Non-sister Alpine Gentians

Distributional shifts driven by Quaternary climatic oscillations have been suggested to cause interspecific hybridization and introgression. In this study, we aimed to test this hypothesis by using population transcriptomes and coalescent modeling of two alpine none-sister gentians. Previous studies suggested that historical hybridizations occurred between Gentiana siphonantha and G. straminea in the high-altitude Qinghai-Tibet Plateau although both species are not sister to each other with the most recent divergence. In the present study, we sequenced transcriptomes of 33 individuals from multiple populations of G. siphonantha and G. straminea. The two species are well delimited by nuclear genomic SNPs while phylogenetic analyses of plastomes clustered one G. straminea individual into the G. siphonantha group. Further population structure analyses of the nuclear SNPs suggested that two populations of G. siphonantha were admixed with around 15% ancestry from G. straminea. These analyses suggested genetic introgressions from G. straminea to G. siphonantha. In addition, our coalescent-based modeling results revealed that gene flow occurred between the two species since Last Glacier Maximum after their initial divergence, which might have leaded to the observed introgressions. Our results underscore the significance of transcriptome population data in determining timescale of interspecific gene flow and direction of the resulting introgression.

The origin of the parrotfish species Scarus compressus in the Tropical Eastern Pacific: region-wide hybridization between ancient species pairs

Background

In the Tropical Eastern Pacific (TEP), four species of parrotfishes with complex phylogeographic histories co-occur in sympatry on rocky reefs from Baja California to Ecuador: Scarus compressus, S. ghobban, S. perrico, and S. rubroviolaceus. The most divergent, S. perrico, separated from a Central Indo-Pacific ancestor in the late Miocene (6.6 Ma). We tested the hypothesis that S. compressus was the result of ongoing hybridization among the other three species by sequencing four nuclear markers and a mitochondrial locus in samples spanning 2/3 of the latitudinal extent of the TEP.

Results

A Structure model indicated that K = 3 fit the nuclear data and that S. compressus individuals had admixed genomes. Our data could correctly detect and assign pure adults and F1 hybrids with > 0.90 probability, and correct assignment of F2s was also high in some cases. NewHybrids models revealed that 89.8% (n = 59) of the S. compressus samples were F1 hybrids between either S. perrico × S. ghobban or S. perrico × S. rubroviolaceus. Similarly, the most recently diverged S. ghobban and S. rubroviolaceus were hybridizing in small numbers, with half of the admixed individuals assigned to F1 hybrids and the remainder likely > F1 hybrids. We observed strong mito-nuclear discordance in all hybrid pairs. Migrate models favored gene flow between S. perrico and S. ghobban, but not other species pairs.

Conclusions

Mating between divergent species is giving rise to a region-wide, multispecies hybrid complex, characterized by a high frequency of parental and F1 genotypes but a low frequency of > F1 hybrids. Trimodal structure, and evidence for fertility of both male and female F1 hybrids, suggest that fitness declines sharply in later generation hybrids. In contrast, the hybrid population of the two more recently diverged species had similar frequencies of F1 and > F1 hybrids, suggesting accelerating post-mating incompatibility with time. Mitochondrial genotypes in hybrids suggest that indiscriminate mating by male S. perrico is driving pre-zygotic breakdown, which may reflect isolation of this endemic species for millions of years resulting in weak selection for conspecific mate recognition. Despite overlapping habitat use and high rates of hybridization, species boundaries are maintained by a combination of pre- and post-mating processes in this complex.

Gene Flow in Volant Vertebrates: Species Biology, Ecology and Climate Change

Gene flow, the exchange of genetic material between populations is an important biological process, which shapes and maintains biodiversity. The successful movement of individuals between populations depends on multiple factors determined by species biology and the environment. One of the most important factors regulating gene flow is the ability to move, and flight allows individuals to easily move across geographical barriers. Volant vertebrates are found on some of the remotest islands and contribute significantly to the biodiversity and ecosystem. The availability of next-generation sequencing data for non-model animals has substantially improved our understanding of gene flow and its consequences, allowing us to look at fine-scale patterns. However, most of our understanding regarding gene flow comes from the temperate regions and the Neotropics. The lack of studies from species-rich Asia is striking. In this review, we outline the importance of gene flow and the factors affecting gene flow, especially for volant vertebrates. We especially discuss research studies from tropical biomes of South and Southeast Asia, highlight the lacuna in literature and provide an outline for future studies in this species-rich region.

Tracing species replacement in Iberian marbled newts

Secondary contact between closely related species can lead to the formation of hybrid zones, allowing for interspecific gene flow. Hybrid zone movement can take place if one of the species possesses a competitive advantage over the other, ultimately resulting in species replacement. Such hybrid zone displacement is predicted to leave a genomic footprint across the landscape in the form of asymmetric gene flow (or introgression) of selectively neutral alleles from the displaced to the advancing species. Hybrid zone movement has been suggested for marbled newts in the Iberian Peninsula, supported by asymmetric gene flow and a distribution relict (i.e., an enclave) of Triturus marmoratus in the range of T. pygmaeus. We developed a panel of nuclear and mitochondrial SNP markers to test for the presence of a T. marmoratus genomic footprint in the Lisbon peninsula, south of the enclave. We found no additional populations of T. marmoratus. Analysis with the software Structure showed no genetic traces of T. marmoratus in T. pygmaeus. A principal component analysis showed some variation within the local T. pygmaeus, but it is unclear if this represents introgression from T. marmoratus. The results may be explained by (a) species replacement without introgressive hybridization and (b) displacement with hybridization followed by the near-complete erosion of the footprint by purifying selection. We predict that testing for a genomic footprint north of the reported enclave would confirm that species replacement in these marbled newts occurred with hybridization.

Convergence in sympatry: Evolution of blue-banded wing pattern in Morpho butterflies

Species interactions such as mimicry can promote trait convergence but disentangling this effect from those of shared ecology, evolutionary history, and niche conservatism is often challenging. Here by focusing on wing colour pattern variation within and between three butterfly species living in sympatry in a large proportion of their range, we tested the effect of species interactions on trait diversification. These butterflies display a conspicuous iridescent blue coloration on the dorsal side of their wings and a cryptic brownish colour on the ventral side. Combined with an erratic and fast flight, these colour patterns increase the difficulty of capture by predators and contribute to the high escape abilities of these butterflies. We hypothesize that, beyond their direct contribution to predator escape, these wing patterns can be used as signals of escape abilities by predators, resulting in positive frequency-dependent selection favouring convergence in wing pattern in sympatry. To test this hypothesis, we quantified dorsal wing pattern variations of 723 butterflies from the three species sampled throughout their distribution, including sympatric and allopatric situations and compared the phenotypic distances between species, sex and localities. We detected a significant effect of localities on colour pattern, and higher inter-specific resemblance in sympatry as compared to allopatry, consistent with the hypothesis of local convergence of wing patterns. Our results provide support to the existence of escape mimicry in the wild and stress the importance of estimating trait variation within species to understand trait variation between species, and to a larger extent, trait diversification at the macro-evolutionary scale.

Morning glory species co-occurrence is associated with asymmetrically decreased and cascading reproductive isolation

Hybridization between species can affect the strength of the reproductive barriers that separate those species. Two extensions of this effect are (1) the expectation that asymmetric hybridization or gene flow will have asymmetric effects on reproductive barrier strength and (2) the expectation that local hybridization will affect only local reproductive barrier strength and could therefore alter within‐species compatibility. We tested these hypotheses in a pair of morning glory species that exhibit asymmetric gene flow from highly selfing Ipomoea lacunosa into mixed‐mating Ipomoea cordatotriloba in regions where they co‐occur. Because of the direction of this gene flow, we predicted that reproductive barrier strength would be more strongly affected in I. cordatotriloba than I. lacunosa. We also predicted that changes to reproductive barriers in sympatric I. cordatotriloba populations would affect compatibility with allopatric populations of that species. We tested these predictions by measuring the strength of a reproductive barrier to seed set across the species’ ranges. Consistent with our first prediction, we found that sympatric and allopatric I. lacunosa produce the same number of seeds in crosses with I. cordatotriloba, whereas crosses between sympatric I. cordatotriloba and I. lacunosa are more successful than crosses between allopatric I. cordatotriloba and I. lacunosa. This difference in compatibility appears to reflect an asymmetric decrease in the strength of the barrier to seed set in sympatric I. cordatotriloba, which could be caused by I. lacunosa alleles that have introgressed into I. cordatotriloba. We further demonstrated that changes to sympatric I. cordatotriloba have decreased its ability to produce seeds with allopatric populations of the same species, in line with our second prediction. Thus, in a manner analogous to cascade reinforcement, we suggest that introgression associated with hybridization not only influences between‐species isolation but can also contribute to isolation within a species.

Maintenance of Species Differences in Closely Related Tetraploid Parasitic Euphrasia (Orobanchaceae) on an Isolated Island

Polyploidy is pervasive in angiosperm evolution and plays important roles in adaptation and speciation. However, polyploid groups are understudied due to complex sequence homology, challenging genome assembly, and taxonomic complexity. Here, we study adaptive divergence in taxonomically complex eyebrights (Euphrasia), where recent divergence, phenotypic plasticity, and hybridization blur species boundaries. We focus on three closely related tetraploid species with contrasting ecological preferences that are sympatric on Fair Isle, a small isolated island in the British Isles. Using a common garden experiment, we show a genetic component to the morphological differences present between these species. Using whole-genome sequencing and a novel k-mer approach we call "Tetmer", we demonstrate that the species are of allopolyploid origin, with a sub-genome divergence of approximately 5%. Using ∼2 million SNPs, we show sub-genome homology across species, with a very low sequence divergence characteristic of recent speciation. This genetic variation is broadly structured by species, with clear divergence of Fair Isle heathland Euphrasia micrantha, while grassland Euphrasia arctica and coastal Euphrasia foulaensis are more closely related. Overall, we show that tetraploid Euphrasia is a system of allopolyploids of postglacial species divergence, where adaptation to novel environments may be conferred by old variants rearranged into new genetic lineages.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Mimicry and mitonuclear discordance in nudibranchs: New insights from exon capture phylogenomics

Phylogenetic inference and species delimitation can be challenging in taxonomic groups that have recently radiated and where introgression produces conflicting gene trees, especially when species delimitation has traditionally relied on mitochondrial data and color pattern. Chromodoris, a genus of colorful and toxic nudibranch in the Indo-Pacific, has been shown to have extraordinary cryptic diversity and mimicry, and has recently radiated, ultimately complicating species delimitation. In these cases, additional genome-wide data can help improve phylogenetic resolution and provide important insights about evolutionary history. Here, we employ a transcriptome-based exon capture approach to resolve Chromodoris phylogeny with data from 2,925 exons and 1,630 genes, derived from 15 nudibranch transcriptomes. We show that some previously identified mimics instead show mitonuclear discordance, likely deriving from introgression or mitochondrial capture, but we confirm one "pure" mimic in Western Australia. Sister-species relationships and species-level entities were recovered with high support in both concatenated maximum likelihood (ML) and summary coalescent phylogenies, but the ML topologies were highly variable while the coalescent topologies were consistent across datasets. Our work also demonstrates the broad phylogenetic utility of 149 genes that were previously identified from eupulmonate gastropods. This study is one of the first to (a) demonstrate the efficacy of exon capture for recovering relationships among recently radiated invertebrate taxa, (b) employ genome-wide nuclear markers to test mimicry hypotheses in nudibranchs and (c) provide evidence for introgression and mitochondrial capture in nudibranchs.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation F_n individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

The conservation value of admixed phenotypes in a critically endangered species complex

In today's environmental crisis, conservationists are increasingly confronted with terminally endangered species whose last few surviving populations may be affected by allelic introgression from closely related species. Yet there is a worrying lack of evidence-based recommendations and solutions for this emerging problem. We analyzed genome-wide DNA markers and plumage variability in a critically endangered insular songbird, the Black-winged Myna (BWM, Acridotheres melanopterus). This species is highly threatened by the illegal wildlife trade, with its wild population numbering in the low hundreds, and its continued survival urgently depending on ex-situ breeding. Its three subspecies occur along a geographic gradient of melanism and are variably interpreted as three species. However, our integrative approach revealed that melanism poorly reflects the pattern of limited genomic differentiation across BWM subspecies. We also uncovered allelic introgression into the most melanistic subspecies, tertius, from the all-black congeneric Javan Myna (A. javanicus), which is native to the same islands. Based on our results, we recommend the establishment of three separate breeding programs to maintain subspecific traits that may confer local adaptation, but with the option of occasional cross-breeding between insurance populations in order to boost genetic diversity and increase overall viability prospects of each breeding program. Our results underscore the importance of evidence-based integrative approaches when determining appropriate conservation units. Given the rapid increase of terminally endangered organisms in need of ex-situ conservation, this study provides an important blueprint for similar programs dealing with phenotypically variable species.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

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

Identifying the traits causing reproductive isolation and the order in which they evolve is fundamental to understanding speciation. Here, we quantify prezygotic and intrinsic postzygotic isolation among allopatric, parapatric, and sympatric populations of the butterflies Heliconius elevatus and Heliconius pardalinus. Sympatric populations from the Amazon (H. elevatus and H. p. butleri) exhibit strong prezygotic isolation and rarely mate in captivity; however, hybrids are fertile. Allopatric populations from the Amazon (H. p. butleri) and Andes (H. p. sergestus) mate freely when brought together in captivity, but the female F1 hybrids are sterile. Parapatric populations (H. elevatus and H. p. sergestus) exhibit both assortative mating and sterility of female F1s. Assortative mating in sympatric populations is consistent with reinforcement in the face of gene flow, where the driving force, selection against hybrids, is due to disruption of mimicry and other ecological traits rather than hybrid sterility. In contrast, the lack of assortative mating and hybrid sterility observed in allopatric populations suggests that geographic isolation enables the evolution of intrinsic postzygotic reproductive isolation. Our results show how the types of reproductive barriers that evolve between species may depend on geography.

Thistledown velvet ants in the Desert Mimicry Ring and the evolution of white coloration: Müllerian mimicry, camouflage and thermal ecology

Adaptive coloration among animals is one of the most recognizable outcomes of natural selection. Here, we investigate evolutionary drivers of white coloration in velvet ants (Hymenoptera: Mutillidae), which has previously been considered camouflage with the fruit of creosote bush (Larrea tridentata). Our analyses indicate instead that velvet ants evolved white coloration millions of years before creosote bush was widespread in North America's hot deserts. Furthermore, velvet ants and the creosote fruit exhibit different spectral reflectance patterns, which appear distinct to potential insectivorous predators. While the white coloration in velvet ants likely did not evolve as camouflage, we find that white-coloured species remain cooler than their red/orange relatives, and therefore we infer the white coloration likely evolved in response to Neogene desertification. This study shows the importance of cross-disciplinary investigation and of testing multiple hypotheses when investigating evolutionary drivers of adaptive coloration.

The social supergene dates back to the speciation time of two Solenopsis fire ant species

Colony social organization of multiple Solenopsis fire ant species is determined by a supergene with two haplotypes SB and Sb, which are similar to X/Y sex chromosomes. The ancestral monogyne (single-queen) social form has been associated with homozygous SB/SB queens, while queens in colonies with the derived polygyne (multi-queen) social structure are heterozygous SB/Sb. By comparing 14 Solenopsis invicta genomes and the outgroup S. fugax, we dated the formation of the supergene to 1.1 (0.7-1.6) million years ago, much older than previous estimates, and close to the estimated time of speciation of the two socially polymorphic species S. invicta and S. richteri. We also used 12 S. invicta and S. richteri genomes to compare the evolutionary distances between these species and the distances between the social haplotypes, and found them to be similar. A phylogenetic analysis suggested that the monophyletic Sb clade is more closely related to S. richteri SB haplotypes than to S. invicta SB haplotypes. We conclude that the formation of the supergene occurred concomitantly with the process of speciation of the Solenopsis socially-polymorphic clade, and hypothesize that the Sb variant first arouse in one incipiently-speciating population and then introgressed into the other populations or species.

Speciation through chromosomal fusion and fission in Lepidoptera

Changes in chromosome numbers may strongly affect reproductive barriers, because individuals heterozygous for distinct karyotypes are typically expected to be at least partially sterile or to show reduced recombination. Therefore, several classic speciation models are based on chromosomal changes. One import mechanism generating variation in chromosome numbers is fusion and fission of existing chromosomes, which is particularly likely in species with holocentric chromosomes, i.e. chromosomes that lack a single centromere. Holocentric chromosomes evolved repeatedly across the tree of life, including in Lepidoptera. Although changes in chromosome numbers are hypothesized to be an important driver of the spectacular diversification of Lepidoptera, comparative studies across the order are lacking. We performed the first comprehensive literature survey of karyotypes for Lepidoptera species since the 1970s and tested if, and how, chromosomal variation might affect speciation. Even though a meta-analysis of karyological differences between closely related taxa did not reveal an effect on the degree of reproductive isolation, phylogenetic diversification rate analyses across the 16 best-covered genera indicated a strong, positive association of rates of chromosome number evolution and speciation. These findings suggest a macroevolutionary impact of varying chromosome numbers in Lepidoptera and likely apply to other taxonomic groups, especially to those with holocentric chromosomes. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Patterns of genomic divergence and introgression between Japanese stickleback species with overlapping breeding habitats

With only a few absolute geographic barriers in marine environments, the factors maintaining reproductive isolation among marine organisms remain elusive. However, spatial structuring in breeding habitat can contribute to reproductive isolation. This is particularly important for marine organisms that migrate to use fresh- or brackish water environments to breed. The Japanese Gasterosteus stickleback species, the Pacific Ocean three-spined stickleback (G. aculeatus) and the Japan Sea stickleback (G. nipponicus) overwinter in the sea, but migrate to rivers for spawning. Although they co-occur at several locations across the Japanese islands, they are reproductively isolated. Our previous studies in Bekanbeushi River showed that the Japan Sea stickleback spawns in the estuary, while the Pacific Ocean stickleback mainly spawns further upstream in freshwater. Overall genomic divergence was very high with many interspersed regions of introgression. Here, we investigated genomic divergence and introgression between the sympatric species in the much shorter Tokotan River, where they share spawning sites. The levels of genome-wide divergence were reduced and introgression was increased, suggesting that habitat isolation substantially contributes to a reduction in gene flow. We also found that genomic regions of introgression were largely shared between the two systems. Furthermore, some regions of introgression were located near loci with a heterozygote advantage for juvenile survival. Taken together, introgression may be partially driven by adaptation in this system. Although, the two species remain clearly genetically differentiated. Regions with low recombination rates showed especially low introgression. Speciation reversal is therefore likely prevented by barriers other than habitat isolation.

The roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena

BACKGROUND

Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri.

RESULTS

We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex.

CONCLUSIONS

Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

Defining Species When There is Gene Flow

Whatever one's definition of species, it is generally expected that individuals of the same species should be genetically more similar to each other than they are to individuals of another species. Here, we show that in the presence of cross-species gene flow, this expectation may be incorrect. We use the multispecies coalescent model with continuous-time migration or episodic introgression to study the impact of gene flow on genetic differences within and between species and highlight a surprising but plausible scenario in which different population sizes and asymmetrical migration rates cause a genetic sequence to be on average more closely related to a sequence from another species than to a sequence from the same species. Our results highlight the extraordinary impact that even a small amount of gene flow may have on the genetic history of the species. We suggest that contrasting long-term migration rate and short-term hybridization rate, both of which can be estimated using genetic data, may be a powerful approach to detecting the presence of reproductive barriers and to define species boundaries.[Gene flow; introgression; migration; multispecies coalescent; species concept; species delimitation.].

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first- and second-generation hybrids produced by controlled crosses between either two co-mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major-effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

In animals, bright colors often evolve to mimic other species when a resemblance is selectively favored. Understanding the proximate mechanisms underlying such color mimicry can give insights into how mimicry evolves-for example, whether color convergence evolves from a shared set of mechanisms or through the evolution of novel color production mechanisms. We studied color production mechanisms in poison frogs (Dendrobatidae), focusing on the mimicry complex of Ranitomeya imitator. Using reflectance spectrometry, skin pigment analysis, electron microscopy, and color modeling, we found that the bright colors of these frogs, both within and outside the mimicry complex, are largely structural and produced by iridophores but that color production depends crucially on interactions with pigments. Color variation and mimicry are regulated predominantly by iridophore platelet thickness and, to a lesser extent, concentration of the red pteridine pigment drosopterin. Compared with each of the four morphs of model species that it resembles, R. imitator displays greater variation in both structural and pigmentary mechanisms, which may have facilitated phenotypic divergence in this species. Analyses of nonmimetic dendrobatids in other genera demonstrate that these mechanisms are widespread within the family and that poison frogs share a complex physiological "color palette" that can produce diverse and highly reflective colors.

Gene Flow in the Müllerian Mimicry Ring of a Poisonous Papuan Songbird Clade (Pitohui; Aves)

Müllerian mimicry rings are remarkable symbiotic species assemblages in which multiple members share a similar phenotype. However, their evolutionary origin remains poorly understood. Although gene flow among species has been shown to generate mimetic patterns in some Heliconius butterflies, mimicry is believed to be due to true convergence without gene flow in many other cases. We investigated the evolutionary history of multiple members of a passerine mimicry ring in the poisonous Papuan pitohuis. Previous phylogenetic evidence indicates that the aposematic coloration shared by many, but not all, members of this genus is ancestral and has only been retained by members of the mimicry ring. Using a newly assembled genome and thousands of genomic DNA markers, we demonstrate gene flow from the hooded pitohui (Pitohui dichrous) into the southern variable pitohui (Pitohui uropygialis), consistent with shared patterns of aposematic coloration. The vicinity of putatively introgressed loci is significantly enriched for genes that are important in melanin pigment expression and toxin resistance, suggesting that gene flow may have been instrumental in the sharing of plumage patterns and toxicity. These results indicate that interspecies gene flow may be a more general mechanism in generating mimicry rings than hitherto appreciated.

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus)

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.

The evolution of sexual signaling is linked to odorant receptor tuning in perfume-collecting orchid bees

Sexual signaling is an important reproductive barrier known to evolve early during the formation of new species, but the genetic mechanisms that facilitate the divergence of sexual signals remain elusive. Here we isolate a gene linked to the rapid evolution of a signaling trait in a pair of nascent neotropical orchid bee lineages, Euglossa dilemma and E. viridissima. Male orchid bees acquire chemical compounds from their environment to concoct species-specific perfumes to later expose during courtship. We find that the two lineages acquire chemically distinct perfumes and are reproductively isolated despite low levels of genome-wide differentiation. Remarkably, variation in perfume chemistry coincides with rapid divergence in few odorant receptor (OR) genes. Using functional assays, we demonstrate that the derived variant of Or41 in E. dilemma is specific towards its species-specific major perfume compound, whereas the ancestral variant in E. viridissima is broadly tuned to multiple odorants. Our results show that OR evolution likely played a role in the divergence of sexual communication in natural populations.

Adaptive Changes in Hemoglobin Function in High-Altitude Tibetan Canids Were Derived via Gene Conversion and Introgression

A key question in evolutionary biology concerns the relative importance of different sources of adaptive genetic variation, such as de novo mutations, standing variation, and introgressive hybridization. A corollary question concerns how allelic variants derived from these different sources may influence the molecular basis of phenotypic adaptation. Here, we use a protein-engineering approach to examine the phenotypic effect of putatively adaptive hemoglobin (Hb) mutations in the high-altitude Tibetan wolf that were selectively introgressed into the Tibetan mastiff, a high-altitude dog breed that is renowned for its hypoxia tolerance. Experiments revealed that the introgressed coding variants confer an increased Hb–O₂ affinity in conjunction with an enhanced Bohr effect. We also document that affinity-enhancing mutations in the β-globin gene of Tibetan wolf were originally derived via interparalog gene conversion from a tandemly linked β-globin pseudogene. Thus, affinity-enhancing mutations were introduced into the β-globin gene of Tibetan wolf via one form of intragenomic lateral transfer (ectopic gene conversion) and were subsequently introduced into the Tibetan mastiff genome via a second form of lateral transfer (introgression). Site-directed mutagenesis experiments revealed that the increased Hb–O₂ affinity requires a specific two-site combination of amino acid replacements, suggesting that the molecular underpinnings of Hb adaptation in Tibetan mastiff (involving mutations that arose in a nonexpressed gene and which originally fixed in Tibetan wolf) may be qualitatively distinct from functionally similar changes in protein function that could have evolved via sequential fixation of de novo mutations during the breed’s relatively short duration of residency at high altitude.

What nature separated, and human joined together: About a spontaneous hybridization between two allopatric dogwood species (Cornus controversa and C. alternifolia)

In this study, possible hybridization between two allopatric species, Cornus controversa and Cornus alternifolia, was explored using molecular and morphological approaches. Scanning electron microscope analyses of the adaxial and the abaxial leaf surfaces yielded a few new not yet described characters typical for the particular species and intermediate for hybrids. With the use of 14 Random Amplified Polymorphic DNA and 5 Amplified Fragment Length Polymorphism primer combinations, 44 fragments species specific to C. controversa and 51 species specific to C. alternifolia were obtained. Most of these bands were also found in putative hybrids. All clustering analyses based on binary data combined from both methods confirmed a separate and intermediate status of the hybrids. Hybrid index estimates for hybrids C1-C5 indicated that all were the first generation of offspring (F₁). Chloroplast intergenic spacers (trnF-trnL and psbC-trnS) were used to infer the hybridization direction. Based on the assumption of maternal inheritance of chloroplast DNA, C. controversa seems to be the maternal parent of the hybrid. Internal transcribed spacer sequences of the five hybrids analyzed here indicated higher similarity with the sequences of C. controversa (all shared the majority of its single nucleotide polymorphisms). Sequence analysis of PI-like genes fully confirmed the hybrid origin of C1-C5 hybrids. Our results also showed that two specimens in the C. alternifolia group, A1 and A3, are not free of introgression. They are probably repeated backcrosses toward C. alternifolia. Furthermore, molecular data seem to point not only to unidirectional introgression toward C. controversa (the presence of hybrids) but to bidirectional introgression as well, since the presence of markers specific for C. controversa in the profiles of C. alternifolia specimen A3 was observed.

Melanism evolution in the cat family is influenced by intraspecific communication under low visibility

Melanism in the cat family has been associated with functions including camouflage, thermoregulation and parasite resistance. Here we investigate a new hypothesis proposing that the evolution of melanism in cats has additionally been influenced by communication functions of body markings. To evaluate this hypothesis, we assembled a species-level data set of morphological (body marks: white marks on the backs of ears) and ecological (circadian activity: arrhythmic/nocturnal, and environmental preference: open/closed) characteristics that could be associated with communication via body markings, and combined these data with a dated molecular phylogeny. Next, we tested the association between melanism and communication, first by relating species’ body marks with their ecological conditions, using a Bayesian implementation of the threshold model. Second, to explore the evolution of characteristics potentially influencing melanism in cat species, we modeled their evolution relative to melanism using models of coordinated vs. independent character changes. Our results suggest that white marks are associated with intraspecific communication between individuals that have non-melanistic phenotypes, as well as towards melanistic individuals (without white marks). The absence of white marks in a melanistic individual tends to be a limiting condition for intraspecific visual communication at night, resulting in an evolutionary dilemma for these species, i.e. to be almost invisible at night, but not to communicate visually. The comparative analysis of several evolutionary models indicated more support for the evolution of melanism being coordinated with the evolution of arrhythmic activity and white marks on the backs of ears.

The Genetic Basis of Scale-Loss Phenotype in the Rapid Radiation of Takifugu Fishes

Rapid radiation associated with phenotypic divergence and convergence provides an opportunity to study the genetic mechanisms of evolution. Here we investigate the genus Takifugu that has undergone explosive radiation relatively recently and contains a subset of closely-related species with a scale-loss phenotype. By using observations during development and genetic mapping approaches, we show that the scale-loss phenotype of two Takifugu species, T. pardalis Temminck & Schlegel and T. snyderi Abe, is largely controlled by an overlapping genomic segment (QTL). A search for candidate genes underlying the scale-loss phenotype revealed that the QTL region contains no known genes responsible for the evolution of scale-loss phenotype in other fishes. These results suggest that the genes used for the scale-loss phenotypes in the two Takifugu are likely the same, but the genes used for the similar phenotype in Takifugu and distantly related fishes are not the same. Meanwhile, Fgfrl1, a gene predicted to function in a pathway known to regulate bone/scale development was identified in the QTL region. Since Fgfr1a1, another memebr of the Fgf signaling pathway, has been implicated in scale loss/scale shape in fish distantly related to Takifugu, our results suggest that the convergence of the scale-loss phenotype may be constrained by signaling modules with conserved roles in scale development.

An amphibian species pushed out of Britain by a moving hybrid zone

Classical theory states that hybrid zones will be stable in troughs of low population density where dispersal is hampered. Yet, evidence for moving hybrid zones is mounting. One possible reason that moving zones have been underappreciated is that they may drive themselves into oblivion and with just the superseding species remaining, morphological and genetic signals of past species replacement may be difficult to appreciate. Using genetic data (32 diagnostic single nucleotide polymorphisms) from a clinal hybrid zone of the common toad (Bufo bufo) and the spined toad (Bufo spinosus) in France for comparison, alleles of the latter species were documented in common toads in the south of Great Britain, at frequencies in excess of 10%. Because long distance dispersal across the Channel is unlikely, the conclusion reached was that the continental toad hybrid zone which previously extended into Britain, moved southwards and extirpated B. spinosus. Species distribution models for the mid-Holocene and the present support that climate has locally changed in favour of B. bufo. The system bears resemblance with the demise of Homo neanderthalensis and the rise of Homo sapiens and provides an example that some paleoanthropologists demanded in support of a hominin "leaky replacement" scenario. The toad example is informative just because surviving pure B. spinosus and an extant slowly moving interspecific hybrid zone are available for comparison.

Eukaryote hybrid genomes

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve.