Bidirectional adaptive introgression between two ecologically divergent sparrow species

On the origin of an insular hybrid butterfly lineage

A new species can form through hybridization between species. Hybrid speciation in animals has been intensely debated, partly because hard evidence for the process has been difficult to obtain. Here, we report the discovery of a European hybrid butterfly lineage, a finding that can be considered surprising given the intense and long-term study of European butterflies. The lineage we describe is mainly inhabiting an island in the Baltic Sea and was previously designated as a subspecies (horkei) of one of the parental species (Aricia artaxerxes). By analyzing whole-genome resequencing data and developing a novel cluster analysis based on historical recombination events (Fisher junctions), we determine that horkei originated by hybridization between the nonsister species A. artaxerxes and A. agestis. This hybridization event occurred approximately 54,000 years ago, predating the last glaciation of the current distribution range. Horkei must therefore have persisted long enough to be able to colonize its current range, despite that this area lies between the current distributions of the parental species. The hybrid origin, the maintenance of genomic integrity across times of dramatic climate change, and the expression of a combination of parental traits suggest that horkei could be in the process of hybrid speciation.

Ongoing introgression of a secondary sexual plumage trait in a stable avian hybrid zone

Hybrid zones are dynamic systems where natural selection, sexual selection, and other evolutionary forces can act on reshuffled combinations of distinct genomes. The movement of hybrid zones, individual traits, or both are of particular interest for understanding the interplay between selective processes. In a hybrid zone involving two lek-breeding birds, secondary sexual plumage traits of Manacus vitellinus, including bright yellow collar and olive belly color, have introgressed ~50 km asymmetrically across the genomic center of the zone into populations more genetically similar to Manacus candei. Males with yellow collars are preferred by females and are more aggressive than parental M. candei, suggesting that sexual selection was responsible for the introgression of male traits. We assessed the spatial and temporal dynamics of this hybrid zone using historical (1989-1994) and contemporary (2017-2020) transect samples to survey both morphological and genetic variation. Genome-wide single nucleotide polymorphism data and several male phenotypic traits show that the genomic center of the zone has remained spatially stable, whereas the olive belly color of male M. vitellinus has continued to introgress over this time period. Our data suggest that sexual selection can continue to shape phenotypes dynamically, independent of a stable genomic transition between species.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

Revisiting the Baltimore-Bullock's Oriole hybrid zone reveals changing plumage colour in Bullock's Orioles

Hybrid zones are dynamic areas where populations of two or more interbreeding species may change through an influx of novel genetic material resulting from hybridization or selection on standing genetic variation. Documenting changes in populations through time, however, is challenging because repeated samples are often missing or because long-term storage can affect trait morphologies, especially colour traits that may fade through time. We document a change in carotenoid-based orange breast feathers of Bullock's Orioles (Icterus bullockii) from the Great Plains hybrid zone, USA. Contemporary Bullock's Orioles are more orange than historic individuals from the same location sampled approximately 60 years ago. Spectrophotometry revealed that contemporary Bullock's Orioles resemble orange colour profiles of Baltimore Orioles (I. galbula), the species with which they hybridize. Fading or changes in diet hypotheses do not appear to explain the shift in colour we report for Bullock's Orioles. We propose that these changes in colour are facilitated through introgression with Baltimore Orioles, and favoured by females that choose brighter, more orange males. Our study highlights the long memory of natural history collections and how they offer new insights to the dynamic roll of hybrid zones in trait evolution between interacting species.

Hybrid enrichment of adaptive variation revealed by genotype-environment associations in montane sedges

The role of hybridization in diversification is complex and may result in many possible outcomes. Not only can hybridization produce new lineages, but those lineages may contain unique combinations of adaptive genetic variation derived from parental taxa that allow hybrid‐origin lineages to occupy unique environmental space relative to one (or both) parent(s). We document such a case of hybridization between two sedge species, Carex nova and Carex nelsonii (Cyperaceae), that occupy partially overlapping environmental space in the southern Rocky Mountains, USA. In the region hypothesized to be the origin of the hybrid lineage, one parental taxon (C. nelsonii) is at the edge of its environmental tolerance. Hybrid‐origin individuals display mixed ancestry between the parental taxa—of nearly 7000 unlinked loci sampled, almost 30% showed evidence of excess ancestry from one parental lineage—approximately half displayed a genomic background skewed towards one parent, and half skewed towards the other. To test whether excess ancestry loci may have conferred an adaptive advantage to the hybrid‐origin lineage, we conducted genotype–environment association analyses on different combinations of loci—with and without excess ancestry—and with multiple contrasts between the hybrids and parental taxa. Loci with skewed ancestry showed significant environmental associations distinguishing the hybrid lineage from one parent (C. nelsonii), whereas loci with relatively equal representation of parental ancestries showed no such environmental associations. Moreover, the overwhelming majority of candidate adaptive loci with respect to environmental gradients also had excess ancestry from a parental lineage, implying these loci have facilitated the persistence of the hybrid lineage in an environment unsuitable to at least one parent.

Bidirectional introgression between Mus musculus domesticus and Mus spretus

Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles—including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.

Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus)

Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.

Prevalence and Adaptive Impact of Introgression

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

Sperm length divergence as a potential prezygotic barrier in a passerine hybrid zone

The saltmarsh sparrow Ammospiza caudacuta and Nelson's sparrow A. nelsoni differ in ecological niche, mating behavior, and plumage, but they hybridize where their breeding distributions overlap. In this advanced hybrid zone, past interbreeding and current backcrossing result in substantial genomic introgression in both directions, although few hybrids are currently produced in most locations. However, because both species are nonterritorial and have only brief male-female interactions, it is difficult to determine to what extent assortative mating explains the low frequency of hybrid offspring. Since females often copulate with multiple males, a role of sperm as a postcopulatory prezygotic barrier appears plausible. Here, we show that sperm length differs between the two species in the hybrid zone, with low among-male variation consistent with strong postcopulatory sexual selection on sperm cells. We hypothesize that divergence in sperm length may constitute a reproductive barrier between species, as sperm length co-evolves with the size of specialized female sperm storage tubules. Sperm does not appear to act as a postzygotic barrier, as sperm from hybrids was unexceptional.

Signatures of Introgression across the Allele Frequency Spectrum

The detection of introgression from genomic data is transforming our view of species and the origins of adaptive variation. Among the most widely used approaches to detect introgression is the so-called ABBA-BABA test or D-statistic, which identifies excess allele sharing between nonsister taxa. Part of the appeal of D is its simplicity, but this also limits its informativeness, particularly about the timing and direction of introgression. Here we present a simple extension, D frequency spectrum or DFS, in which D is partitioned according to the frequencies of derived alleles. We use simulations over a large parameter space to show how DFS carries information about various factors. In particular, recent introgression reliably leads to a peak in DFS among low-frequency derived alleles, whereas violation of model assumptions can lead to a lack of signal at low frequencies. We also reanalyze published empirical data from six different animal and plant taxa, and interpret the results in the light of our simulations, showing how DFS provides novel insights. We currently see DFS as a descriptive tool that will augment both simple and sophisticated tests for introgression, but in the future it may be usefully incorporated into probabilistic inference frameworks.

Genic introgression from an invasive exotic fungal forest pathogen increases the establishment potential of a sibling native pathogen

Significant hybridization between the invasive North American fungal plant pathogen Heterobasidion irregulare and its Eurasian sister species H. annosum is ongoing in Italy. Whole genomes of nine natural hybrids were sequenced, assembled and compared with those of three genotypes each of the two parental species. Genetic relationships among hybrids and their level of admixture were determined. A multi-approach pipeline was used to assign introgressed genomic blocks to each of the two species. Alleles that introgressed from H. irregulare to H. annosum were associated with pathways putatively related to saprobic processes, while alleles that introgressed from the native to the invasive species were mainly linked to gene regulation. There was no overlap of allele categories introgressed in the two directions. Phenotypic experiments documented a fitness increase in H. annosum genotypes characterized by introgression of alleles from the invasive species, supporting the hypothesis that hybridization results in putatively adaptive introgression. Conversely, introgression from the native into the exotic species appeared to be driven by selection on genes favoring genome stability. Since the introgression of specific alleles from the exotic H. irregulare into the native H. annosum increased the invasiveness of the latter species, we propose that two invasions may be co-occurring: the first one by genotypes of the exotic species, and the second one by alleles belonging to the exotic species. Given that H. irregulare represents a threat to European forests, monitoring programs need to track not only exotic genotypes in native forest stands, but also exotic alleles introgressed in native genotypes.

Patterns of introgression vary within an avian hybrid zone

Background

Exploring hybrid zone dynamics at different spatial scales allows for better understanding of local factors that influence hybrid zone structure. In this study, we tested hypotheses about drivers of introgression at two spatial scales within the Saltmarsh Sparrow (Ammospiza caudacuta) and Nelson’s Sparrow (A. nelsoni) hybrid zone. Specifically, we evaluated the influence of neutral demographic processes (relative species abundance), natural selection (exogenous environmental factors and genetic incompatibilities), and sexual selection (assortative mating) in this mosaic hybrid zone. By intensively sampling adults (n = 218) and chicks (n = 326) at two geographically proximate locations in the center of the hybrid zone, we determined patterns of introgression on a fine scale across sites of differing habitat. We made broadscale comparisons of patterns from the center with those of prior studies in the southern edge of the hybrid zone.

Results

A panel of fixed SNPs (135) identified from ddRAD sequencing was used to calculate a hybrid index and determine genotypic composition/admixture level of the populations. Another panel of polymorphic SNPs (589) was used to assign paternity and reconstruct mating pairs to test for sexual selection. On a broad-scale, patterns of introgression were not explained by random mating within marshes. We found high rates of back-crossing and similarly low rates of recent-generation (F1/F2) hybrids in the center and south of the zone. Offspring genotypic proportions did not meet those expected from random mating within the parental genotypic distribution. Additionally, we observed half as many F1/F2 hybrid female adults than nestlings, while respective male groups showed no difference, in support of Haldane’s Rule. The observed proportion of interspecific mating was lower than expected when accounting for mate availability, indicating assortative mating was limiting widespread hybridization. On a fine spatial scale, we found variation in the relative influence of neutral and selective forces between inland and coastal habitats, with the smaller, inland marsh influenced primarily by neutral demographic processes, and the expansive, coastal marsh experiencing higher selective pressures in the form of natural (exogenous and endogenous) and sexual selection.

Conclusions

Multiple drivers of introgression, including neutral and selective pressures (exogenous, endogenous, and sexual selection), are structuring this hybrid zone, and their relative influence is site and context-dependent.

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.

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.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.

Ancestral Hybridization Facilitated Species Diversification in the Lake Malawi Cichlid Fish Adaptive Radiation

The adaptive radiation of cichlid fishes in East African Lake Malawi encompasses over 500 species that are believed to have evolved within the last 800,000 years from a common founder population. It has been proposed that hybridization between ancestral lineages can provide the genetic raw material to fuel such exceptionally high diversification rates, and evidence for this has recently been presented for the Lake Victoria region cichlid superflock. Here, we report that Lake Malawi cichlid genomes also show evidence of hybridization between two lineages that split 3-4 Ma, today represented by Lake Victoria cichlids and the riverine Astatotilapia sp. "ruaha blue." The two ancestries in Malawi cichlid genomes are present in large blocks of several kilobases, but there is little variation in this pattern between Malawi cichlid species, suggesting that the large-scale mosaic structure of the genomes was largely established prior to the radiation. Nevertheless, tens of thousands of polymorphic variants apparently derived from the hybridization are interspersed in the genomes. These loci show a striking excess of differentiation across ecological subgroups in the Lake Malawi cichlid assemblage, and parental alleles sort differentially into benthic and pelagic Malawi cichlid lineages, consistent with strong differential selection on these loci during species divergence. Furthermore, these loci are enriched for genes involved in immune response and vision, including opsin genes previously identified as important for speciation. Our results reinforce the role of ancestral hybridization in explosive diversification by demonstrating its significance in one of the largest recent vertebrate adaptive radiations.

Divergent selection and drift shape the genomes of two avian sister species spanning a saline-freshwater ecotone

The role of species divergence due to ecologically based divergent selection-or ecological speciation-in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline-freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes-the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)-to evaluate the influence of selective and demographic processes in shaping genome-wide patterns of divergence. Genome-wide divergence between these two recently diverged sister species was notably high (genome-wide F ST = 0.32). Against a background of high genome-wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism-all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Differential introgression of a female competitive trait in a hybrid zone between sex-role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male-male competition in maintaining or eroding species barriers. Although female-female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex-role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female-female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex-role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.