Locus-specific introgression in young hybrid swarms: Drift may dominate selection

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Introgressive hybridization in the west Pacific pen shells (genus Atrina): Restricted interspecies gene flow within the genome

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

Rapid introgression of non-native alleles following hybridization between a native Anolis lizard species and a cryptic invader across an urban landscape

Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced-representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non-native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non-native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non-native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non-native genetic material even in the absence of ongoing immigration, indicating that selection favouring non-native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non-native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long-term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual-based simulations under a stepping-stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

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.