What, if anything, are hybrids: enduring truths and challenges associated with population structure and gene flow

Simulating the efficacy of wolf-dog hybridization management with individual-based modeling

Introgressive hybridization between wolves and dogs is a conservation concern due to its potentially deleterious long-term evolutionary consequences. European legislation requires that wolf-dog hybridization be mitigated through effective management. We developed an individual-based model (IBM) to simulate the life cycle of gray wolves that incorporates aspects of wolf sociality that affect hybridization rates (e.g., the dissolution of packs after the death of one/both breeders) with the goal of informing decision-making on management of wolf-dog hybridization. We applied our model by projecting hybridization dynamics in a local wolf population under different mate choice and immigration scenarios and contrasted results of removal of admixed individuals with their sterilization and release. In several scenarios, lack of management led to complete admixture, whereas reactive management interventions effectively reduced admixture in wolf populations. Management effectiveness, however, strongly depended on mate choice and number and admixture level of individuals immigrating into the wolf population. The inclusion of anthropogenic mortality affecting parental and admixed individuals (e.g., poaching) increased the probability of pack dissolution and thus increased the probability of interbreeding with dogs or admixed individuals and boosted hybridization and introgression rates in all simulation scenarios. Recognizing the necessity of additional model refinements (appropriate parameterization, thorough sensitivity analyses, and robust model validation) to generate management recommendations applicable in real-world scenarios, we maintain confidence in our model's potential as a valuable conservation tool that can be applied to diverse situations and species facing similar threats.

A Geometric Morphometrics Approach to the Study of Natural Variations and Hybrid Detection in Populations of Alnus incana (L.) Moench and Alnus rohlenae Vít, Douda and Mandák

Landmark-based geometric morphometrics (GM) was used to examine, for the first time, spontaneous hybridization between Alnus incana (L.) Moench and Alnus rohlenae Vít, Douda and Mandák, and to assess inter- and intrapopulation variability in leaf shape, leaf size and venation in natural populations in Serbia (Western Balkans). Two geographically distant (30 km) and two close (1.2 km) populations were selected to examine hybridization. The variability in leaf shapes was assessed by canonical variate analysis and linear discriminant analysis performed on the symmetric component of variation. Covariation between the symmetric component of shape variation and the number of pairs of secondary leaf veins was investigated with partial least squares analysis. Static allometry was examined for the first time in the genus Alnus Mill. A higher proportion of A. incana leaves was classified as A. rohlenae in geographically close populations, which is in accordance with the hypothesis about spontaneous hybridization. No single leaf of A. rohlenae was classified as A. incana, indicating that putative hybrids can only be found in grey alder populations. This study demonstrates that GM is a powerful tool for species delimitation and hybrid detection in the genus Alnus and it can be used for preliminary screening in hybrid zones.

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

The puzzle of plant hybridisation: a high propensity to hybridise but few hybrid zones reported

An interesting conundrum was recently revealed by R. Abbott when he found that the number of hybrid zones reported in the literature for plants is very low, given the propensity of plants to hybridise. In another literature survey on hybrid zones performed over the period 1970-2022, we found that the number of hybrid zones reported for vertebrates was 2.3 times greater than that reported for vascular plants, even though there are about six times more vascular plant species than vertebrates. Looking at the number of papers reporting hybrid zones, there are 4.9 times more on vertebrates than on vascular plants. These figures support the relevance of this conundrum. In this paper we aim to shed light on this question by providing a structured discussion of the causes that may underlie this conundrum. We propose six non-mutually exclusive factors, namely lack or deficit of spatial structure, lack or deficit of genetic structure, effects of hybridisation between non-closely related species, lability of plant hybrid zones over time, botanists' perception of hybridisation, and deficit of population genetic data. There does not appear to be a single factor that explains our puzzle, which applies to all cases of plants where hybridisation is detected but no hybrid zone is reported. It is argued that some plant features suggest that the puzzle is not, at least entirely, due to insufficient knowledge of the specific cases, a hypothesis that should be addressed with a wider range of empirical data across different taxonomic groups.

Disentangling the causes of high polymorphism sharing in sympatric Petunia species from subtropical highland grasslands: insights from nuclear diversity

Genetic polymorphism sharing between closely related and sympatric plant species could result from common ancestry, ancient or recent hybridization. Here we analyzed four Petunia species from the subtropical highland grasslands in southern South America based on nuclear diversity to disentangle the causes of high polymorphism sharing between them. We genotyped microsatellite loci, employed population genetic methods to estimate variability, species limits, and ancient and recent gene flow, and assigned individuals to genetic and taxonomic groups. Finally, we modeled evolutionary processes to determine the impact of Quaternary climate changes on species phylogenetic relationships. Our results indicated that genetic diversity was strongly influenced by expansion and habitat fragmentation during the Quaternary cycles. The extensive polymorphism sharing is mainly due to species' common ancestry, and we did not discard ancient hybridization. Nowadays, niche differentiation is the primary driver for maintaining genetic and morphological limits between the four analysed Petunia species and there is no recent gene flow between them.

Origins and evolution of biological novelty

Understanding the origins and impacts of novel traits has been a perennial interest in many realms of ecology and evolutionary biology. Here, we build on previous evolutionary and philosophical treatments of this subject to encompass novelties across biological scales and eco-evolutionary perspectives. By defining novelties as new features at one biological scale that have emergent effects at other biological scales, we incorporate many forms of novelty that have previously been treated in isolation (such as novelty from genetic mutations, new developmental pathways, new morphological features, and new species). Our perspective is based on the fundamental idea that the emergence of a novelty, at any biological scale, depends on its environmental and genetic context. Through this lens, we outline a broad array of generative mechanisms underlying novelty and highlight how genomic tools are transforming our understanding of the origins of novelty. Lastly, we present several case studies to illustrate how novelties across biological scales and systems can be understood based on common mechanisms of change and their environmental and genetic contexts. Specifically, we highlight how gene duplication contributes to the evolution of new complex structures in visual systems; how genetic exchange in symbiosis alters functions of both host and symbiont, resulting in a novel organism; and how hybridisation between species can generate new species with new niches.

Hybridization boosters diversification in a Neotropical Bulbophyllum (Orchidaceae) group

Genetic data shows that cryptic hybrids are more common than previously thought and that hybridization and introgression are widespread processes. Regardless, studies on hybridization are scarce for the highly speciose Bulbophyllum. The genus presents more than 2,200 species and many examples of recent radiations, in which hybridization is expected to be frequent. Currently, only four natural Bulbophyllum hybrids are recognized, all of them recently described based on morphological evidence. Here we test whether genomic evidence supports the hybrid status of two Neotropical Bulbophyllum species, while also evaluating the impact of this phenomenon on the genomes of the putative parental species. We also assess if there is evidence of hybridization among B. involutum and B. exaltatum, sister species that diverged recently. We leverage the power of next-generation sequence data, associated with model-based analysis for three systems putatively constituted by two parental species and one hybrid. All taxa belong to the Neotropical B. sect. Didactyle clade. We found evidence of hybridization in all studied systems. Despite the occurrence of hybridization, there are no signs of backcrossing. Because of the high propensity of hybridization across many taxa, the common occurrence of hybridization during the evolutionary history of B. sect. Didactyle means it is time to account for and examine its evolutionary role in these orchids.

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

Hybridization between closely related songbirds is related to human habitat disturbance

Human habitat disturbances can promote hybridization between closely related, but typically reproductively isolated, species. We explored whether human habitat disturbances are related to hybridization between two closely related songbirds, black-capped and mountain chickadees, using both genomic and citizen science data sets. First, we genotyped 409 individuals from across both species' ranges using reduced-representation genome sequencing and compared measures of genetic admixture to a composite measure of human landscape disturbance. Then, using eBird observations, we compared human landscape disturbance values for sites where phenotypically diagnosed hybrids were observed to locations where either parental species was observed to determine whether hybrid chickadees are reported in more disturbed areas. We found that hybridization between black-capped and mountain chickadees positively correlates with human habitat disturbances. From genomic data, we found that (1) hybrid index (HI) significantly increased with habitat disturbance, (2) more hybrids were sampled in disturbed habitats, (3) mean HIs were higher in disturbed habitats versus wild habitats, and (4) hybrids were detected in habitats with significantly higher disturbance values than parentals. Using eBird data, we found that both hybrid and black-capped chickadees were significantly more disturbance-associated than mountain chickadees. Surprisingly, we found that nearly every black-capped chickadee we sampled contained some proportion of hybrid ancestry, while we detected very few mountain chickadee backcrosses. Our results highlight that hybridization between black-capped and mountain chickadees is widespread, but initial hybridization is rare (few F1s were detected). We conclude that human habitat disturbances can erode pre-zygotic reproductive barriers between chickadees and that post-zygotic isolation is incomplete. Understanding what becomes of recently hybridizing species following large-scale habitat disturbances is a new, but pressing, consideration for successfully preserving genetic biodiversity in a rapidly changing world.

Complex population history affects admixture analyses in nine-spined sticklebacks

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree Plathymenia reticulata (n = 233 individuals). Genetic structure analysis of P. reticulata was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F CT = 0.222 and F ST = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F2) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in P. reticulata. The low genetic diversity of the P. reticulata forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

Identification of Wolf-Dog Hybrids in Europe – An Overview of Genetic Studies

Significant development of genetic tools during the last decades provided opportunities for more detailed analyses and deeper understanding of species hybridization. New genetic markers allowed for reliable identification of admixed individuals deriving from recent hybridization events (a few generations) and those originating from crossings up to 19 generations back. Implementation of microsatellites (STRs) together with Bayesian clustering provided abundant knowledge regarding presence of admixed individuals in numerous populations and helped understand the problematic nature of studying hybridization (i.a., defining a reliable thresholds for recognizing individuals as admixed or obtaining well-grounded results representing actual proportion of hybrids in a population). Nevertheless, their utilization is limited to recent crossbreeding events. Single Nucleotide Polymorphisms (SNPs) proved to be more sensible tools for admixture analyses furnishing more reliable knowledge, especially for older generation backcrosses. Small sets of Ancestry Informative Markers (AIMs) of both types of markers were effective enough to implement in monitoring programs, however, SNPs seem to be more appropriate because of their ability to identify admixed individuals up to 3rd generations. The main aim of this review is to summarize abundant knowledge regarding identification of wolf-dog hybrids in Europe and discuss the most relevant problems relating to the issue, together with advantages and disadvantages of implemented markers and approaches.

Genetic structure in Louisiana Iris species reveals patterns of recent and historical admixture

PREMISE

When divergent lineages come into secondary contact, reproductive isolation may be incomplete, thus providing an opportunity to investigate how speciation is manifested in the genome. The Louisiana Irises (Iris, series Hexagonae) comprise a group of three or more ecologically and reproductively divergent lineages that can produce hybrids where they come into contact. We estimated standing genetic variation to understand the current distribution of population structure in the Louisiana Irises.

METHODS

We used genotyping-by-sequencing techniques to sample the genomes of Louisiana Iris species across their ranges. We sampled 20 populations (n = 632 individuals) across 11,249 loci and used Entropy and PCA models to assess population genetic data.

RESULTS

We discovered evidence for interspecific gene flow in parts of the range. Our analysis revealed patterns of population structure at odds with widely accepted nominal taxonomy. We discovered undescribed hybrid populations, designated as belonging to the I. brevicaulis lineage. Iris nelsonii shared significant ancestry with only one of the purported parent species, I. fulva, evidence inconsistent with a hybrid origin.

CONCLUSIONS

This study provides several key findings important to the investigation of standing genetic variation in the Louisiana Iris species complex. Compared to the other nominal species, I. brevicaulis contains a large amount of genetic diversity. In addition, we discovered a previously unknown hybrid zone between I. brevicaulis and I. hexagona along the Texas coast. Finally, our results do not support the long-standing hypothesis that I. nelsonii has mixed ancestry from three parental taxa.

Displaced clines in an avian hybrid zone (Thamnophilidae: Rhegmatorhina) within an Amazonian interfluve

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi.

A holocenic and dynamic hybrid zone between two cactophilic Drosophila species in a coastal lowland plain of the Brazilian Atlantic Forest

Hybridization and introgression are processes that contribute to shaping biological diversity. The factors promoting the formation of these processes are multiples but poorly explored in a biogeographical and ecological context. In the southeast coastal plain of the Brazilian Atlantic Forest, a hybrid zone was described between two closely related cactophilic species, Drosophila antonietae and D. serido. Here, we revisited and analysed specimens from this hybrid zone to evaluate its temporal and spatial dynamic. We examined allopatric and sympatric populations of the flies using independent sources of data such as mitochondrial and nuclear sequences, microsatellite loci, morphometrics of wings and male genitalia, and climatic niche models. We also verified the emergence of the flies from necrotic tissues of collected cacti to verify the role of host association for the population dynamics. Our results support the existence of a hybrid zone due to secondary contact and limited to the localities where the two species are currently in contact. Furthermore, we detected asymmetric bidirectional introgression and the maintenance of the species integrity, ecological association and morphological characters, suggesting selection and limited introgression. Considering our paleomodels, probably this hybrid zone is recent and the contact occurred during the Holocene to the present day, favoured by range expansion of their populations due to expansion of open and dry areas in eastern South America during palaeoclimatic and geomorphological events.

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.

Extent and direction of introgressive hybridization of mule and white-tailed deer in western Canada

Hybridization of mule deer (Odocoileus hemionus) and white-tailed deer (O. virginianus) appears to be a semi-regular occurrence in western North America. Previous studies confirmed the presence of hybrids in a variety of sympatric habitats, but their developing molecular resources limited identification to the earliest, most admixed generations. For this reason, estimates of hybrid production in wild populations often rely on anecdotal reports. As well, white-tailed deer populations' continued encroachment into historically mule deer-occupied habitats due to changes in land use, habitat homogenization, and a warming climate may increase opportunities for interspecific encounters. We sought to quantify the prevalence and extent of hybrid deer in the prairies of western Canada using a SNP assay with enhanced discriminating power. By updating the available molecular resources, we sought to identify and characterize previously cryptic introgression. We also investigated the influence of various parameters on hybridity by way of logistic regression. We observed overall hybridization rates of ~1.0%, slightly lower than that reported by previous studies, and found white-tailed-like hybrids to be more common than their mule deer-like counterparts. Here, we build upon past studies of hybridization in North American deer by increasing hybrid detection power, expanding sample sizes, demonstrating a new molecular resource applicable to future research and observing asymmetrical directionality of introgression.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Asymmetric introgression reveals the genetic architecture of a plumage trait

Genome-wide variation in introgression rates across hybrid zones offers a powerful opportunity for studying population differentiation. One poorly understood pattern of introgression is the geographic displacement of a trait implicated in lineage divergence from genome-wide population boundaries. While difficult to interpret, this pattern can facilitate the dissection of trait genetic architecture because traits become uncoupled from their ancestral genomic background. We studied an example of trait displacement generated by the introgression of head plumage coloration from personata to alba subspecies of the white wagtail. A previous study of their hybrid zone in Siberia revealed that the geographic transition in this sexual signal that mediates assortative mating was offset from other traits and genetic markers. Here we show that head plumage is associated with two small genetic regions. Despite having a simple genetic architecture, head plumage inheritance is consistent with partial dominance and epistasis, which could contribute to its asymmetric introgression.

Hybrid zones are windows into the evolutionary process. Semenov et al. find that the head plumage differences between white wagtail subspecies have a simple genetic basis involving two small genetic regions, in which partially dominant and epistatic interactions help to explain how this sexual signal has become decoupled from other plumage traits.

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Genetic Admixture in the Population of Wild Apple (Malus Sieversii) from the Tien Shan Mountains, Kazakhstan

There is growing attention given to gene flow between crops and the wild relatives as global landscapes have been rapidly converted into agricultural farm fields over the past century. Crop-to-wild introgression may advance the extinction risks of rare plants through demographic swamping and/or genetic swamping. Malus sieversii, the progenitor of the apple, is exclusively distributed along the Tien Shan mountains. Habitat fragmentation and hybridization between M. sieversii and the cultivated apples have been proposed to be the causal mechanism of the accelerated extinction risk. We examined the genetic diversity pattern of eleven wild and domesticated apple populations and assessed the gene flow between M. sieversii and the cultivated apples in Kazakhstan using thirteen nuclear microsatellite loci. On average, apple populations harbored fairly high within-population diversity, whereas population divergences were very low suggesting likely influence of human-mediated dispersal. Assignment results showed a split pattern between the cultivated and wild apples and frequent admixture among the apple populations. Coupled with the inflated contemporary migration rates, the admixture pattern might be the signature of increased human intervention within the recent past. Our study highlighted the prevalent crop to wild gene flow of apples occurring in Kazakhstan, proposing an accelerated risk of genetic swamping.

Genome-Wide DArTSeq Genotyping and Phenotypic Based Assessment of Within and Among Accessions Diversity and Effective Sample Size in the Diverse Sorghum, Pearl Millet, and Pigeonpea Landraces

Germplasm should be conserved in such a way that the genetic integrity of a given accession is maintained. In most genebanks, landraces constitute a major portion of collections, wherein the extent of genetic diversity within and among landraces of crops vary depending on the extent of outcrossing and selection intensity infused by farmers. In this study, we assessed the level of diversity within and among 108 diverse landraces and wild accessions using both phenotypic and genotypic characterization. This included 36 accessions in each of sorghum, pearl millet, and pigeonpea, conserved at ICRISAT genebank. We genotyped about 15 to 25 individuals within each accession, totaling 1,980 individuals using the DArTSeq approach. This resulted in 45,249, 19,052, and 8,211 high-quality single nucleotide polymorphisms (SNPs) in pearl millet, sorghum, and pigeonpea, respectively. Sorghum had the lowest average phenotypic (0.090) and genotypic (0.135) within accession distances, while pearl millet had the highest average phenotypic (0.227) and genotypic (0.245) distances. Pigeonpea had an average of 0.203 phenotypic and 0.168 genotypic within accession distances. Analysis of molecular variance also confirms the lowest variability within accessions of sorghum (26.3%) and the highest of 80.2% in pearl millet, while an intermediate in pigeonpea (57.0%). The effective sample size required to capture maximum variability and to retain rare alleles while regeneration ranged from 47 to 101 for sorghum, 155 to 203 for pearl millet, and 77 to 89 for pigeonpea accessions. This study will support genebank curators, in understanding the dynamics of population within and among accessions, in devising appropriate germplasm conservation strategies, and aid in their utilization for crop improvement.

Riverscape properties contribute to the origin and structure of a hybrid zone in a Neotropical freshwater fish

Understanding the structure of hybrid zones provides valuable insights about species boundaries and speciation, such as the evolution of barriers to gene flow and the strength of selection. In river networks, studying evolutionary processes in hybrid zones can be especially challenging, given the influence of past and current river properties along with biological species-specific traits. Here, we suggest that a natural hybrid zone between two divergent lineages of the sexually dimorphic Neotropical fish Nematocharax venustus was probably established by secondary contact as a result of a river capture event between the Contas and Pardo river basins. This putative river capture is supported by hydrogeological evidence of elbows of capture, wind gaps and geological faults. The morphological (colour pattern) and genetic (mtDNA and RADseq) variation reveal a clinal transition between parental lineages along the main river, with predominance of F2 hybrids at the centre of the hybrid zone, absence of early generation backcrosses and different levels of hybridization in the tributaries. We highlight that different sources of information are crucial for understanding how the riverscape spatial history influences the connectivity between and within river systems and, consequently, the dynamics of gene flow between freshwater lineages/species.

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental-like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.

The protected tree Dimorphandra wilsonii (Fabaceae) is a population of inter-specific hybrids: recommendations for conservation in the Brazilian Cerrado/Atlantic Forest ecotone

BACKGROUNDS AND AIMS

Dimorphandra wilsonii Rizzini, a critically endangered and protected tree, has a restricted distribution in the ecotone between the Cerrado and the Atlantic Forest in south-eastern Brazil. In this area, it co-occurs with D. mollis Benth., a common tree from the Cerrado, and D. exaltata Schott., a rare tree from the Atlantic Forest. Previous studies of D. wilsonii indicated heterozygosity excess at the individual level. Field observation of some intermediate phenotypes between D. wilsonii and both congeners suggests hybridization of D. wilsonii with D. mollis and/or D. exaltata. Here, we tested the hypothesis that D. wilsonii may have originated from hybridization between D. exaltata and D. mollis. We also performed cytogenetic analysis to examine if the heterozygosity excess could be explained by polyploidy in D. wilsonii.

METHODS

We evaluated the genetic diversity and population structure of D. wilsonii using 11 nuclear simple sequence repeats (SSRs) genotyped in 152 individuals sampled across the taxon's range. We performed comparative genetic analyses using overlapping SSR markers between D. wilsonii and previously published SSR data in D. mollis and D. exaltata to subsequently perform a series of allelic comparisons, multivariate and Bayesian analysis.

KEY RESULTS

Our results suggest that D. wilsonii individuals are most likely to correspond to F1 hybrids between D. exaltata and D. mollis. Cytogenetic analysis indicated that D. wilsonii is diploid with the same chromosome number as D. mollis (2n = 2x = 28).

CONCLUSIONS

Our study raises questions about the taxonomic status and the evolutionary future of D. wilsonii. We suggest that the conservation and management strategy for D. wilsonii should be revised and that it should take into account both parental Dimorphandra species in the ecotone, with special emphasis on the threatened D. exaltata. Finally, this study highlights the value of genetic information for the design of conservation strategies.

What could be the fate of secondary contact zones between closely related plant species?

Interspecific hybridization has been fundamental in plant evolution. Nevertheless, the fate of hybrid zones throughout the generations remains poorly addressed. We analyzed a pair of recently diverged, interfertile, and sympatric Petunia species to ask what fate the interspecific hybrid population has met over time. We analyzed the genetic diversity in two generations from two contact sites and evaluated the effect of introgression. To do this, we collected all adult plants from the contact zones, including canonicals and intermediary colored individuals, and compared them with purebred representatives of both species based on seven highly informative microsatellite loci. We compared the genetic diversity observed in the contact zones with what is seen in isolated populations of each species, considering two generations of these annual species. Our results have confirmed the genetic differentiation between the species and the hybrid origin of the majority of the intermediary colored individuals. We also observed a differentiation related to genetic variability and inbreeding levels among the populations. Over time, there were no significant differences per site related to genetic diversity or phenotype composition. We found two stable populations kept by high inbreeding and backcross rates that influence the genetic diversity of their parental species through introgression.

Environment-dependent introgression from Quercus dentata to a coastal ecotype of Quercus mongolica var. crispula in northern Japan

Introgression from one species in a specific environment to another may facilitate colonization of the environment by the recipient species. However, such environment-dependent introgression has been clarified in limited plant taxa. In northern Japan, there are two interfertile oak species: Quercus dentata (Qd) in coastal areas and Q. mongolica var. crispula (Qc) in inland areas. However, at higher latitudes where Qd is rare, a coastal Qc ecotype with Qd-like traits is distributed in the coastal areas. We distinguished inland Qc, coastal Qc, and coastal Qd populations based on genome-wide genotypes and multitrait phenotypes and verified introgression from coastal Qd to coastal Qc using reduced library sequencing. Genotypes and phenotypes differed among the populations, and coastal Qc was intermediate between inland Qc and coastal Qd. The ABBA-BABA test showed introgression from coastal Qd to coastal Qc. In coastal Qc, we found various stages of introgression after the first generation of backcross but detected no genomic regions where introgression was enhanced. Overall, we show evidence for introgression from a coastal species to an ecotype of an inland species, which has colonized the coastal environment. It remains unclear whether introgressed alleles are selected in the coastal environment.

Hybridizing salamanders experience accelerated diversification

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

Variable outcomes of hybridization between declining Alosa alosa and Alosa fallax

Hybridization dynamics between co-occurring species in environments where human-mediated changes take place are important to quantify for furthering our understanding of human impacts on species evolution and for informing management. The allis shad Alosa alosa (Linnaeus, 1758) and twaite shad Alosa fallax (Lacépède, 1803), two clupeids sister species, have been severely impacted by human activities across Europe. The shrinkage of A. alosa distribution range along with the decline of the remaining populations' abundance threatens its persistence. The main objective was to evaluate the extent of hybridization and introgression between those interacting species. We developed a set of 77 species-specific SNP loci that allowed a better resolution than morphological traits as they enabled the detection of hybrids up to the third generation. Variable rates of contemporary hybridization and introgression patterns were detected in 12 studied sites across the French Atlantic coast. Mitochondrial markers revealed a cyto-nuclear discordance almost invariably involving A. alosa individuals with an A. fallax mitochondrial DNA and provided evidence of historical asymmetric introgression. Overall, contemporary and historical introgression revealed by nuclear and mitochondrial markers strongly suggests that a transfer of genes occurs from A. fallax toward A. alosa genome since at least four generations. Moreover, the outcomes of introgression greatly depend on the catchments where local processes are thought to occur. Undoubtedly, interspecific interaction and gene flow should not be overlooked when considering the management of those species.

Tall Pinus luzmariae trees with genes from P. herrerae

CONTEXT

Pinus herrerae and P. luzmariae are endemic to western Mexico, where they cover an area of more than 1 million hectares. Pinus herrerae is also cultivated in field trials in South Africa and South America, because of its considerable economic importance as a source of timber and resin. Seed quality, afforestation success and desirable traits may all be influenced by the presence of hybrid trees in seed stands.

AIMS

We aimed to determine the degree of hybridization between P. herrerae and P. luzmariae in seed stands of each species located in the Sierra Madre Occidental, Durango, Mexico.

METHODS

AFLP molecular markers from samples of 171 trees across five populations were analyzed with STRUCTURE and NewHybrids software to determine the degree of introgressive hybridization. The accuracy of STRUCTURE and NewHybrids in detecting hybrids was quantified using the software Hybridlab 1.0. Morphological analysis of 131 samples from two populations of P. herrerae and two populations of P. luzmariae was also conducted by Random Forest classification. The data were compared by Principal Coordinate Analysis (PCoA) in GenAlex 6.501.

RESULTS

Hybridization between Pinus herrerae and P. luzmariae was observed in all seed stands under study and resulted in enhancement of desirable silvicultural traits in the latter species. In P. luzmariae, only about 16% molecularly detected hybrids correspond to those identified on a morphological basis. However, the morphology of P. herrerae is not consistent with the molecularly identified hybrids from one population and is only consistent with 3.3 of those from the other population.

CONCLUSIONS

This is the first report of hybrid vigour (heterosis) in Mexican pines. Information about hybridization and introgression is essential for developing effective future breeding programs, successful establishment of plantations and management of natural forest stands. Understanding how natural hybridization may influence the evolution and adaptation of pines to climate change is a cornerstone to sustainable forest management including adaptive silviculture.

A standardized approach to empirically define reliable assignment thresholds and appropriate management categories in deeply introgressed populations

Anthropogenic hybridization is recognized as a major threat to the long-term survival of natural populations. While identifying F1 hybrids might be simple, the detection of older admixed individuals is far from trivial and it is still debated whether they should be targets of management. Examples of anthropogenic hybridization have been described between wolves and domestic dogs, with numerous cases detected in the Italian wolf population. After selecting appropriate wild and domestic reference populations, we used empirical and simulated 39-autosomal microsatellite genotypes, Bayesian assignment and performance analyses to develop a workflow to detect different levels of wolf x dog admixture. Membership proportions to the wild cluster (q_iw) and performance indexes identified two q-thresholds which allowed to efficiently classify the analysed genotypes into three assignment classes: pure (with no or negligible domestic ancestry), older admixed (with a marginal domestic ancestry) and recent admixed (with a clearly detectable domestic ancestry) animals. Based on their potential to spread domestic variants, such classes were used to define three corresponding management categories: operational pure, introgressed and operational hybrid individuals. Our multiple-criteria approach can help wildlife managers and decision makers in more efficiently targeting the available resources for the long-term conservation of species threatened by anthropogenic hybridization.

Recombining Your Way Out of Trouble: The Genetic Architecture of Hybrid Fitness under Environmental Stress

Hybridization between species can either promote or impede adaptation. But we know very little about the genetic basis of hybrid fitness, especially in nondomesticated organisms, and when populations are facing environmental stress. We made genetically variable F2 hybrid populations from two divergent Saccharomyces yeast species. We exposed populations to ten toxins and sequenced the most resilient hybrids on low coverage using ddRADseq to investigate four aspects of their genomes: 1) hybridity, 2) interspecific heterozygosity, 3) epistasis (positive or negative associations between nonhomologous chromosomes), and 4) ploidy. We used linear mixed-effect models and simulations to measure to which extent hybrid genome composition was contingent on the environment. Genomes grown in different environments varied in every aspect of hybridness measured, revealing strong genotype–environment interactions. We also found selection against heterozygosity or directional selection for one of the parental alleles, with larger fitness of genomes carrying more homozygous allelic combinations in an otherwise hybrid genomic background. In addition, individual chromosomes and chromosomal interactions showed significant species biases and pervasive aneuploidies. Against our expectations, we observed multiple beneficial, opposite-species chromosome associations, confirmed by epistasis- and selection-free computer simulations, which is surprising given the large divergence of parental genomes (∼15%). Together, these results suggest that successful, stress-resilient hybrid genomes can be assembled from the best features of both parents without paying high costs of negative epistasis. This illustrates the importance of measuring genetic trait architecture in an environmental context when determining the evolutionary potential of genetically diverse hybrid populations.

Increased genetic marker density reveals high levels of admixture between red deer and introduced Japanese sika in Kintyre, Scotland

Hybridization is a natural process at species range boundaries, but increasing numbers of species are hybridizing due to direct or indirect human activities. In such cases of anthropogenic hybridization, subsequent introgression can threaten the survival of native species. To date, many such systems have been studied with too few genetic markers to assess the level of threat resulting from advanced backcrossing. Here, we use 44,999 single nucleotide polymorphisms (SNPs) and the ADMIXTURE program to study two areas of Scotland where a panel of 22 diagnostic microsatellites previously identified introgression between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon). In Kintyre, we reclassify 26% of deer from the pure species categories to the hybrid category whereas in the NW Highlands we only reclassify 2%. As expected, the reclassified individuals are mostly advanced backcrosses. We also investigate the ability of marker panels selected on different posterior allele frequency criteria to find hybrids assigned by the full marker set and show that in our data, ancestry informative markers (i.e. those that are highly differentiated between the species, but not fixed) are better than diagnostic markers (those markers that are fixed between the species) because they are more evenly distributed in the genome. Diagnostic loci are concentrated on the X chromosome to the detriment of autosomal coverage.

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south-western and south-eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south-western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south-eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150-km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2-27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population-level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.

Agricultural insect hybridization and implications for pest management

Biological invasions, the expansion of agricultural frontiers, and climate change favor encounters of divergent lineages of animals and plants, increasing the likelihood of hybridization. However, hybridization of insect species and its consequences for agroecosystems have not received sufficient attention. Gene exchange between distinct and distant genetic pools can improve the survival and reproduction of insect pests, and threaten beneficial insects in disturbed agricultural environments. Hybridization may be the underlying explanation for the recurrent pest outbreaks and control failures in putative hybrid zones, as suspected for bollworm, corn borer, whiteflies, and stink bugs. Reliable predictions of the types of changes that can be expected in pest insect genomes and fitness, and of their impacts on the fate of species and populations remain elusive. Typical steps in pest management, such as insect identification, pest monitoring, and control are likely affected by gene flow and adaptive introgression mediated by hybridization, and we do not have ways to respond to or mitigate the problem. To address the adverse effects of farming intensification and global trade, we must ensure that current integrated pest management programs incorporate up-to-date monitoring and diagnostic tools. The rapid identification of hybrids, quantification of levels of introgression, and in-depth knowledge of what genes have been transferred may help to explain and predict insect population outbreaks and control failures in the future. © 2019 Society of Chemical Industry.

Genomic signatures of introgression between commercial and native bumblebees, Bombus terrestris, in western Iberian Peninsula-Implications for conservation and trade regulation

Human-mediated introductions of species may have profound impacts on native ecosystems. One potential impact with largely unforeseen consequences is the potential admixture of introduced with autochthonous species through hybridization. Throughout the world, bumblebees have been deliberately introduced for crop pollination with known negative impacts on native pollinators. Given the likely allochthonous origin of commercial bumblebees used in Portugal (subspecies Bombus terrestris terrestris and B. t. dalmatinus), our aim was to assess their putative introgression with the native Iberian subspecies B. terrestris lusitanicus. We analysed one mitochondrial gene, cytochrome c oxidase subunit I (COX1) and genomic data involving thousands of genome-wide restriction-site-associated DNA markers (RAD-seq). In the mitochondrial COX1 analyses, we detected one relatively common haplotype in commercial bumblebees, also present in wild samples collected nearby the greenhouses where the commercial hives are used. In the RAD-seq analysis, we found a clear genetic differentiation between native and commercial lineages. Furthermore, we detected candidate hybrids in the wild, as well as putatively escaped commercial bumblebees, some of which being potentially fertile males. Although we cannot assess directly the fitness effects of introgressed alleles, there is a risk of maladaptive allele introgression to the local bumblebee subspecies, which can negatively impact autochthon populations. One immediate recommendation to farmers is for the proper disposal of hive boxes, after their use in greenhouses, so as to minimize the risk of escapees contaminating native populations. On the other hand, the feasibility of using local subspecies B. t. lusitanicus, preferably with local production, should be evaluated.

Asymmetric introgression between fishes in the Red River basin of Texas is associated with variation in water quality

When ecologically divergent taxa encounter one another, hybrid zones can form when reproductive isolation is incomplete. The location of such hybrid zones can be influenced by environmental variables, and an ecological context can provide unique insights into the mechanisms by which species diverge and are maintained. Two ecologically differentiated species of small benthic fishes, the endemic and imperiled prairie chub, Macrhybopsis australis, and the shoal chub, Macrhybopsis hyostoma, are locally sympatric within the upper Red River Basin of Texas. We integrated population genomic data and environmental data to investigate species divergence and the maintenance of species boundaries in these two species. We found evidence of advanced-generation asymmetric hybridization and introgression, with shoal chub alleles introgressing more frequently into prairie chubs than the reciprocal. Using a Bayesian Genomic Cline framework, patterns of genomic introgression were revealed to be quite heterogeneous, yet shoal chub alleles were found to have likely selectively introgressed across species boundaries significantly more often than prairie chub alleles, potentially explaining some of the observed asymmetry in hybridization. These patterns were remarkably consistent across two sampled geographic regions of hybridization. Several environmental variables were found to significantly predict individual admixture, suggesting ecological isolation might maintain species boundaries.

Genomic consequences of a recent three-way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human-mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make-up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double-digested restriction site-associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple-way admixture with population genomic data.