Disentangling Timing of Admixture, Patterns of Introgression, and Phenotypic Indicators in a Hybridizing Wolf Population

A first exploratory comparison of the behaviour of wolves (Canis lupus) and wolf-dog hybrids in captivity

Extensive introgression of genes from domesticated taxa may be a serious threat for the genomic integrity and adaptability of wild populations. Grey wolves (Canis lupus) are especially vulnerable to this phenomenon, but there are no studies yet assessing the potential behavioural effects of dog-introgression in wolves. In this study, we conducted a first systematic comparison of admixed (N = 11) and non-admixed (N = 14) wolves in captivity, focusing on their reaction to unfamiliar humans and novel objects, and the cohesiveness of their social groups. When exposed to unfamiliar humans in the experimental task, wolves were more vigilant, fearful and aggressive than admixed wolves, and less likely to approach humans, but also more likely to spend time in human proximity. When exposed to novel objects, wolves were more aggressive than admixed wolves, less likely to spend time in object proximity, and more likely to interact with objects, but also less vigilant and as fearful as admixed wolves. Finally, social networks were more cohesive in wolves than in admixed wolves. Although caution is needed when comparing groups of captive individuals with different life experiences, our study suggests that dog admixture may lead to important behavioural changes in wolves, with possible implications for conservation strategies.

Genetic swamping of the critically endangered Scottish wildcat was recent and accelerated by disease

The European wildcat population in Scotland is considered critically endangered as a result of hybridization with introduced domestic cats,1,2 though the time frame over which this gene flow has taken place is unknown. Here, using genome data from modern, museum, and ancient samples, we reconstructed the trajectory and dated the decline of the local wildcat population from viable to severely hybridized. We demonstrate that although domestic cats have been present in Britain for over 2,000 years,3 the onset of hybridization was only within the last 70 years. Our analyses reveal that the domestic ancestry present in modern wildcats is markedly over-represented in many parts of the genome, including the major histocompatibility complex (MHC). We hypothesize that introgression provides wildcats with protection against diseases harbored and introduced by domestic cats, and that this selection contributes to maladaptive genetic swamping through linkage drag. Using the case of the Scottish wildcat, we demonstrate the importance of local ancestry estimates to both understand the impacts of hybridization in wild populations and support conservation efforts to mitigate the consequences of anthropogenic and environmental change.

The population bottleneck of the Iberian wolf impacted genetic diversity but not admixture with domestic dogs: A temporal genomic approach

After decades of intense persecution, the Iberian wolf subspecies faced a severe bottleneck in the 1970s that considerably reduced its range and population size, nearly leading to its extinction in central and southern Iberian Peninsula. Such population decline could have impacted the genetic diversity of Iberian wolves through different processes, namely genetic drift and dynamics of hybridization with domestic dogs. By contrasting the genomes of 68 contemporary with 54 historical samples spanning the periods before and immediately after the 1970s bottleneck, we found evidence of its impact on genetic diversity and dynamics of wolf-dog hybridization. Our genome-wide assessment revealed that wolves and dogs form two well-differentiated genetic groups in Iberia and that hybridization rates did not increase during the bottleneck. However, an increased number of hybrid individuals was found over time during the population re-expansion, particularly at the edge of the wolf range. We estimated a low percentage of dog ancestry (~1.4%) in historical samples, suggesting that dog introgression was not a key driver for wolf extinction in central and southern Iberia. Our findings also unveil a significant decline in genetic diversity in contemporary samples, with the highest proportion of homozygous segments in the genome being recently inherited. Overall, our study provides unprecedented insight into the impact of a sharp decline on the Iberian wolf genome and refines our understanding of the ecological and evolutionary drivers of wolf-dog hybridization in the wild.

Possible origins and implications of atypical morphologies and domestication-like traits in wild golden jackals (Canis aureus)

Deciphering the origins of phenotypic variations in natural animal populations is a challenging topic for evolutionary and conservation biologists. Atypical morphologies in mammals are usually attributed to interspecific hybridisation or de-novo mutations. Here we report the case of four golden jackals (Canis aureus), that were observed during a camera-trapping wildlife survey in Northern Israel, displaying anomalous morphological traits, such as white patches, an upturned tail, and long thick fur which resemble features of domesticated mammals. Another individual was culled under permit and was genetically and morphologically examined. Paternal and nuclear genetic profiles, as well as geometric morphometric data, identified this individual as a golden jackal rather than a recent dog/wolf-jackal hybrid. Its maternal haplotype suggested past introgression of African wolf (Canis lupaster) mitochondrial DNA, as previously documented in other jackals from Israel. When viewed in the context of the jackal as an overabundant species in Israel, the rural nature of the surveyed area, the abundance of anthropogenic waste, and molecular and morphological findings, the possibility of an individual presenting incipient stages of domestication should also be considered.

Citizen science facilitates first ever genetic detection of wolf-dog hybridization in Indian savannahs

Human demographic expansion has confined wildlife to fragmented habitats, often in proximity to human-modified landscapes. Such interfaces facilitate increased interactions between feral or domesticated animals and wildlife, posing a high risk to wild species. This is especially relevant for free-ranging dogs (Canis lupus familiaris) and wild canids like gray wolves (Canis lupus) and golden jackals (Canis aureus). Wolf-dog hybridization may lead to a significant reduction of specific adaptations in wolves that could result in the decline of wolf populations. Detection and genetic discrimination of hybrids between dogs and wolves are challenging because of their complex demographic history and close ancestry. Citizen scientists identified two phenotypically different-looking individuals and subsequently collected non-invasive samples that were used by geneticists to test wolf-dog hybridization. Genomic data from shed hair samples of suspected hybrid individuals using double-digest restriction-site-associated DNA (ddRAD) sequencing resulted in 698 single nucleotide polymorphism (SNP) markers. We investigated the genetic origin of these two individuals analyzed with genetically known dogs, wolves, and other canid species including jackals and dholes (Cuon alpinus). Our results provide the first genetic evidence of one F2 hybrid and the other individual could be a complex hybrid between dogs and wolves. Our results re-iterate the power of next-generation sequencing (NGS) for non-invasive samples as an efficient tool for detecting hybrids. Our results suggest the need for more robust monitoring of wolf populations and highlight the tremendous potential for collaborative approaches between citizens and conservation scientists to detect and monitor threats to biodiversity.

How to Choose? Comparing Different Methods to Count Wolf Packs in a Protected Area of the Northern Apennines

Despite a natural rewilding process that caused wolf populations in Europe to increase and expand in the last years, human-wolf conflicts still persist, threatening the long-term wolf presence in both anthropic and natural areas. Conservation management strategies should be carefully designed on updated population data and planned on a wide scale. Unfortunately, reliable ecological data are difficult and expensive to obtain and often hardly comparable through time or among different areas, especially because of different sampling designs. In order to assess the performance of different methods to estimate wolf (Canis lupus L.) abundance and distribution in southern Europe, we simultaneously applied three techniques: wolf howling, camera trapping and non-invasive genetic sampling in a protected area of the northern Apennines. We aimed at counting the minimum number of packs during a single wolf biological year and evaluating the pros and cons for each technique, comparing results obtained from different combinations of these three methods and testing how sampling effort may affect results. We found that packs' identifications could be hardly comparable if methods were separately used with a low sampling effort: wolf howling identified nine, camera trapping 12 and non-invasive genetic sampling eight packs. However, increased sampling efforts produced more consistent and comparable results across all used methods, although results from different sampling designs should be carefully compared. The integration of the three techniques yielded the highest number of detected packs, 13, although with the highest effort and cost. A common standardised sampling strategy should be a priority approach to studying elusive large carnivores, such as the wolf, allowing for the comparison of key population parameters and developing shared and effective conservation management plans.

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

Limited movement of an avian hybrid zone in relation to regional variation in magnitude of climate change

Studies of natural hybrid zones can provide documentation of range shifts in response to climate change and identify loci important to reproductive isolation. Using a temporal (36-38 years) comparison of the black-capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee hybrid zone, we investigated movement of the western portion of the zone (western Missouri) and assessed whether loci and pathways underpinning reproductive isolation were similar to those in the eastern portion of the hybrid zone. Using 92 birds sampled along the hybrid zone transect in 2016 and 68 birds sampled between 1978 and 1980, we generated 11,669 SNPs via ddRADseq. These SNPs were used to assess movement of the hybrid zone through time and to evaluate variation in introgression among loci. We demonstrate that the interface has moved ~5 km to the northwest over the last 36-38 years, that is, at only one-fifth the rate at which the eastern portion (e.g., Pennsylvania, Ohio) of the hybrid zone has moved. Temperature trends over the last 38 years reveal that eastern areas have warmed 50% more than western areas in terms of annual mean temperature, possibly providing an explanation for the slower movement of the hybrid zone in Missouri. Our results suggest hybrid zone movement in broadly distributed species, such as chickadees, will vary between areas in response to local differences in the impacts of climate change.

“Guess Who’s Coming to Dinner”: Molecular Tools to Reconstruct multilocus Genetic Profiles from Wild Canid Consumption Remains

Non-invasive genetic sampling is a practical tool to monitor pivotal ecological parameters and population dynamic patterns of endangered species. It can be particularly suitable when applied to elusive carnivores such as the Apennine wolf (Canis lupus italicus) and the European wildcat (Felis silvestris silvestris), which can live in overlapping ecological contexts and sometimes share their habitats with their domestic free-ranging relatives, increasing the risk of anthropogenic hybridisation. In this case study, we exploited all the ecological and genetic information contained in a single biological canid faecal sample, collected in a forested area of central Italy, to detect any sign of trophic interactions between wolves and European wildcats or their domestic counterparts. Firstly, the faecal finding was morphologically examined, showing the presence of felid hair and claw fragment remains. Subsequently, total genomic DNA contained in the hair and claw samples was extracted and genotyped, through a multiple-tube approach, at canid and felid diagnostic panels of microsatellite loci. Finally, the obtained individual multilocus genotypes were analysed with reference wild and domestic canid and felid populations to assess their correct taxonomic status using Bayesian clustering procedures. Assignment analyses classified the genotype obtained from the endothelial cells present on the hair sample as a wolf with slight signals of dog ancestry, showing a qi = 0.954 (C.I. 0.780–1.000) to the wolf cluster, and the genotype obtained from the claw as a domestic cat, showing a qi = 0.996 (95% C.I. = 0.982–1.000) to the domestic cat cluster. Our results clearly show how a non-invasive multidisciplinary approach allows the cost-effective identification of both prey and predator genetic profiles and their taxonomic status, contributing to the improvement of our knowledge about feeding habits, predatory dynamics, and anthropogenic hybridisation risk in threatened species.

Genomic and morphological data shed light on the complexities of shared ancestry between closely related duck species

Causes for genomic and morphological similarities among recently radiated species are often multifaceted and are further convoluted among species that readily interbreed. Here, we couple genomic and morphological trait comparisons to test the extent that ancestry and gene flow explain the retention of mallard-like traits within a sister species, the Mexican duck. First, we confirm that these taxa remain genetically structured, and that Mexican ducks exhibit an isolation-by-distance pattern. Despite the assumption of wide-spread hybridization, we found only a few late-stage hybrids, all from the southwestern USA. Next, assessing 23 morphological traits, we developed a genetically-vetted morphological key that is > 97% accurate in distinguishing across sex-age cohorts of Mexican ducks, mallards, and hybrids. During key development, we determined that 25% of genetically pure, immature male Mexican ducks of the northern population naturally displayed mallard-like traits in their formative plumage. In fact, applying this key to 55 museum specimens, we identified that only four of the 14 specimens originally classified as phenotypic hybrids were truly hybrids. We discuss how genomic and morphological comparisons shed light into the mechanism(s) underlying the evolution of complex phenotypic traits in recent radiations, and how misunderstanding the true morphological diversity within Mexican ducks resulted in taxonomic revisions that hindered conservation efforts.

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.

A reduced SNP panel to trace gene flow across southern European wolf populations and detect hybridization with other Canis taxa

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ Q_i ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Q_i ≤ 0.827 to parental groups). An assay design and test with invasive and non-invasive canid samples performed successfully for 178 SNPs. By separating natural population admixture from inter-specific hybridization, our reduced panel can help advance evolutionary research, monitoring, and timely conservation management.

Settlement of a stable wolf pack in a highly anthropic area of Pisan hills: Relationship with animal husbandry and hunting in a human-wolf coexistence perspective

Anthropic areas play a pivot role for main wolf conservation challenges. Wolf presence in the higher Pisan hills has been well documented while wolf settlement in the lower Pisan hills is still uncertain. In this study, long-term information on wolf presence in a highly anthropic area of the lower Pisan hills was collected by using non-invasive monitoring techniques. Furthermore, both the relationship of this predator with human activity and the impact of hunting on wolf presence have been investigated. The results obtained indicate the presence of a stable and reproductive wolf pack composed by both Italian wolf and hybrids individuals in the municipalities of Crespina Lorenzana and Casciana Terme Lari. A high impact of wolf on livestock was recorded in this area since no prevention systems were adopted by farmers. Wolf appears not to have a negative impact on wild boar population. Similarly, wild boar drive hunting does not appear to affect the wolf pack presence in the area. Thereby wolf may play a key role as controller of wild population. Prevention strategies improvement becomes instrumental to promote wolf-human coexistence. Further investigation to monitor pack hybridization level and turnover and to assess the impact of packs on wild population and livestock in anthropic areas is desirable.

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.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

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.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples

Background

Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples.

Results

Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations.

Conclusions

We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07761-5.

Genomic evidence for the Chinese mountain cat as a wildcat conspecific (Felis silvestris bieti) and its introgression to domestic cats

The Qinghai-Tibet Plateau endemic Chinese mountain cat has a controversial taxonomic status, whether it is a true species or a wildcat (Felis silvestris) subspecies and whether it has contributed to cat (F. s. catus) domestication in East Asia. Here, we sampled F. silvestris lineages across China and sequenced 51 nuclear genomes, 55 mitogenomes, and multilocus regions from 270 modern or museum specimens. Genome-wide analyses classified the Chinese mountain cat as a wildcat conspecific F. s. bieti, which was not involved in cat domestication of China, thus supporting a single domestication origin arising from the African wildcat (F. s. lybica). A complex hybridization scenario including ancient introgression from the Asiatic wildcat (F. s. ornata) to F. s. bieti, and contemporary gene flow between F. s. bieti and sympatric domestic cats that are likely recent Plateau arrivals, raises the prospect of disrupted wildcat genetic integrity, an issue with profound conservation implications.

Systematic Review on Local Ancestor Inference From a Mathematical and Algorithmic Perspective

Genotypic data provide deep insights into the population history and medical genetics. The local ancestry inference (LAI) (also termed local ancestry deconvolution) method uses the hidden Markov model (HMM) to solve the mathematical problem of ancestry reconstruction based on genomic data. HMM is combined with other statistical models and machine learning techniques for particular genetic tasks in a series of computer tools. In this article, we surveyed the mathematical structure, application characteristics, historical development, and benchmark analysis of the LAI method in detail, which will help researchers better understand and further develop LAI methods. Firstly, we extensively explore the mathematical structure of each model and its characteristic applications. Next, we use bibliometrics to show detailed model application fields and list articles to elaborate on the historical development. LAI publications had experienced a peak period during 2006-2016 and had kept on moving in the following years. The efficiency, accuracy, and stability of the existing models were evaluated by the benchmark. We find that phased data had higher accuracy in comparison with unphased data. We summarize these models with their distinct advantages and disadvantages. The Loter model uses dynamic programming to obtain a globally optimal solution with its parameter-free advantage. Aligned bases can be used directly in the Seqmix model if the genotype is hard to call. This research may help model developers to realize current challenges, develop more advanced models, and enable scholars to select appropriate models according to given populations and datasets.

Whole-genome analyses provide no evidence for dog introgression in Fennoscandian wolf populations

Hybridization and admixture can threaten the genetic integrity of populations and be of particular concern to endangered species. Hybridization between grey wolves and dogs has been documented in many wolf populations worldwide and is a prominent example of human-mediated hybridization between a domesticated species and its wild relative. We analysed whole-genome sequences from >200 wolves and >100 dogs to study admixture in Fennoscandian wolf populations. A principal component analysis of genetic variation and admixture showed that wolves and dogs were well-separated, without evidence for introgression. Analyses of local ancestry revealed that wolves had <1% mixed ancestry, levels comparable to the degree of mixed ancestry in many dogs, and likely not resulting from recent wolf-dog hybridization. We also show that the founders of the Scandinavian wolf population were genetically inseparable from Finnish and Russian Karelian wolves, pointing at the geographical origin of contemporary Scandinavian wolves. Moreover, we found Scandinavian-born animals among wolves sampled in Finland, demonstrating bidirectional gene flow between the Scandinavian Peninsula and eastern countries. The low incidence of admixture between wolves and dogs in Fennoscandia may be explained by the fact that feral dogs are rare in this part of Europe and that careful monitoring and management act to remove hybrids before they backcross into wolf populations.

Improving Selection Detection with Population Branch Statistic on Admixed Populations

Detecting natural selection signals in admixed populations can be problematic since the source of the signal typically dates back prior to the admixture event. On one hand, it is now possible to study various source populations before a particular admixture thanks to the developments in ancient DNA (aDNA) in the last decade. However, aDNA availability is limited to certain geographical regions and the sample sizes and quality of the data might not be sufficient for selection analysis in many cases. In this study, we explore possible ways to improve detection of pre-admixture signals in admixed populations using a local ancestry inference approach. We used masked haplotypes for population branch statistic (PBS) and full haplotypes constructed following our approach from Yelmen et al. (2019) for cross-population extended haplotype homozygosity (XP-EHH), utilizing forward simulations to test the power of our analysis. The PBS results on simulated data showed that using masked haplotypes obtained from ancestry deconvolution instead of the admixed population might improve detection quality. On the other hand, XP-EHH results using the admixed population were better compared with the local ancestry method. We additionally report correlation for XP-EHH scores between source and admixed populations, suggesting that haplotype-based approaches must be used cautiously for recently admixed populations. Additionally, we performed PBS on real South Asian populations masked with local ancestry deconvolution and report here the first possible selection signals on the autochthonous South Asian component of contemporary South Asian populations.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

The relevance of gene flow with wild relatives in understanding the domestication process

The widespread use of genomic tools has allowed for a deeper understanding of the genetics and the evolutionary dynamics of domestication. Recent studies have suggested that multiple domestications and introgression are more common than previously thought. However, the ability to correctly infer the many aspects of domestication process depends on having an adequate representation of wild relatives. Cultivated maize (Zea mays ssp. mays) is one of the most important crops in the world, with a long and a relatively well-documented history of domestication. The current consensus points towards a single domestication event from teosinte Zea mays ssp. parviglumis from the Balsas Basin in Southwestern Mexico. However, the underlying diversity of teosintes from Z. mays ssp. parviglumis and Zea mays ssp. mexicana was not taken into account in early studies. We used 32 739 single nucleotide polymorphisms (SNPs) obtained from 29 teosinte populations and 43 maize landraces to explore the relationship between wild and cultivated members of Zea. We then inferred the levels of gene flow among teosinte populations and maize, the degree of population structure of Zea mays subspecies, and the potential domestication location of maize. We confirmed a strong geographic structure within Z. mays ssp. parviglumis and documented multiple gene flow events with other members of the genus, including an event between Z. mays ssp. mexicana and maize. Our results suggest that the likely ancestor of maize may have been domesticated in Jalisco or in the southern Pacific Coast and not in the Balsas Basin as previously thought. In this context, different populations of both teosinte subspecies have contributed to modern maize's gene pool. Our results point towards a long period of domestication marked by gene flow with wild relatives, confirming domestication as long and ongoing process.

Living a dog’s life: a putative gray wolf in a feral dog group

The gray wolf and the domestic dog are closely related species that can interbreed and produce fertile offspring. In settings where unrestrained dogs are present in the wild, hybridization can happen naturally. However, the behavior of the resulting hybrids and their ecological impact is largely understudied. In September–November 2018, a putative gray wolf was repeatedly camera-trapped in a group of 10 presumably feral dogs in a remote mountainous area (the Osogovo Mountain) along the border between Bulgaria and North Macedonia. The most feasible explanation for this individual’s atypical behavior is that it is of hybrid origin (assumption based on phenotype). To the best of our knowledge, this is the first documented observation of such a kind. A discussion of its recruitment and position in the group is presented, setting the basis for further investigation of the complex interaction between wolves, dogs and hybrids in the wild.

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.

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.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

Global Phylogeographic and Admixture Patterns in Grey Wolves and Genetic Legacy of An Ancient Siberian Lineage

The evolutionary relationships between extinct and extant lineages provide important insight into species' response to environmental change. The grey wolf is among the few Holarctic large carnivores that survived the Late Pleistocene megafaunal extinctions, responding to that period's profound environmental changes with loss of distinct lineages and phylogeographic shifts, and undergoing domestication. We reconstructed global genome-wide phylogeographic patterns in modern wolves, including previously underrepresented Siberian wolves, and assessed their evolutionary relationships with a previously genotyped wolf from Taimyr, Siberia, dated at 35 Kya. The inferred phylogeographic structure was affected by admixture with dogs, coyotes and golden jackals, stressing the importance of accounting for this process in phylogeographic studies. The Taimyr lineage was distinct from modern Siberian wolves and constituted a sister lineage of modern Eurasian wolves and domestic dogs, with an ambiguous position relative to North American wolves. We detected gene flow from the Taimyr lineage to Arctic dog breeds, but population clustering methods indicated closer similarity of the Taimyr wolf to modern wolves than dogs, implying complex post-divergence relationships among these lineages. Our study shows that introgression from ecologically diverse con-specific and con-generic populations was common in wolves' evolutionary history, and could have facilitated their adaptation to environmental change.

Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations

The survival of indigenous European wildcat (Felis silvestris silvestris) populations can be locally threatened by introgressive hybridization with free-ranging domestic cats. Identifying pure wildcats and investigating the ancestry of admixed individuals becomes thus a conservation priority. We analyzed 63k cat Single Nucleotide Polymorphisms (SNPs) with multivariate, Bayesian and gene-search tools to better evaluate admixture levels between domestic and wild cats collected in Europe, timing and ancestry proportions of their hybrids and backcrosses, and track the origin (wild or domestic) of the genomic blocks carried by admixed cats, also looking for possible deviations from neutrality in their inheritance patterns. Small domestic ancestry blocks were detected in the genomes of most admixed cats, which likely originated from hybridization events occurring from 6 to 22 generations in the past. We identified about 1,900 outlier coding genes with excess of wild or domestic ancestry compared to random expectations in the admixed individuals. More than 600 outlier genes were significantly enriched for Gene Ontology (GO) categories mainly related to social behavior, functional and metabolic adaptive processes (wild-like genes), involved in cognition and neural crest development (domestic-like genes), or associated with immune system functions and lipid metabolism (parental-like genes). These kinds of genomic ancestry analyses could be reliably applied to unravel the admixture dynamics in European wildcats, as well as in other hybridizing populations, in order to design more efficient conservation plans.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

Two decades of non-invasive genetic monitoring of the grey wolves recolonizing the Alps support very limited dog introgression

Potential hybridization between wolves and dogs has fueled the sensitive conservation and political debate underlying the recovery of the grey wolf throughout Europe. Here we provide the first genetic analysis of wolf-dog admixture in an area entirely recolonized, the northwestern Alps. As part of a long-term monitoring program, we performed genetic screening of thousands of non-invasive samples collected in Switzerland and adjacent territories since the return of the wolf in the mid-1990s. We identified a total of 115 individuals, only 2 of them showing significant signs of admixture stemming from past interbreeding with dogs, followed by backcrossing. This low rate of introgression (<2% accounting for all wolves ever detected over 1998–2017) parallels those from other European populations, especially in Western Europe (<7%). Despite potential hybridization with stray dogs, few founders and strong anthropogenic pressures, the genetic integrity of the Alpine population has remained intact throughout the entire recolonization process. In a context of widespread misinformation, this finding should reduce conflicts among the different actors involved and facilitate wolf conservation. Real-time genetic monitoring will be necessary to identify potential hybrids and support an effective management of this emblematic population.

Combining Bayesian genetic clustering and ecological niche modeling: Insights into wolf intraspecific genetic structure

The distribution of intraspecific genetic variation and how it relates to environmental factors is of increasing interest to researchers in macroecology and biogeography. Recent studies investigated the relationships between the environment and patterns of intraspecific genetic variation across species ranges but only few rigorously tested the relation between genetic groups and their ecological niches. We quantified the relationship of genetic differentiation (F ST) and the overlap of ecological niches (as measured by n-dimensional hypervolumes) among genetic groups resulting from spatial Bayesian genetic clustering in the wolf (Canis lupus) in the Italian peninsula. Within the Italian wolf population, four genetic clusters were detected, and these clusters showed different ecological niches. Moreover, different wolf clusters were significantly related to differences in land cover and human disturbance features. Such differences in the ecological niches of genetic clusters should be interpreted in light of neutral processes that hinder movement, dispersal, and gene flow among the genetic clusters, in order to not prematurely assume any selective or adaptive processes. In the present study, we found that both the plasticity of wolves-a habitat generalist-to cope with different environmental conditions and the occurrence of barriers that limit gene flow lead to the formation of genetic intraspecific genetic clusters and their distinct ecological niches.

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.

Wolf outside, dog inside? The genomic make-up of the Czechoslovakian Wolfdog

BACKGROUND

Genomic methods can provide extraordinary tools to explore the genetic background of wild species and domestic breeds, optimize breeding practices, monitor and limit the spread of recessive diseases, and discourage illegal crossings. In this study we analysed a panel of 170k Single Nucleotide Polymorphisms with a combination of multivariate, Bayesian and outlier gene approaches to examine the genome-wide diversity and inbreeding levels in a recent wolf x dog cross-breed, the Czechoslovakian Wolfdog, which is becoming increasingly popular across Europe.

RESULTS

Pairwise FST values, multivariate and assignment procedures indicated that the Czechoslovakian Wolfdog was significantly differentiated from all the other analysed breeds and also well-distinguished from both parental populations (Carpathian wolves and German Shepherds). Coherently with the low number of founders involved in the breed selection, the individual inbreeding levels calculated from homozygosity regions were relatively high and comparable with those derived from the pedigree data. In contrast, the coefficient of relatedness between individuals estimated from the pedigrees often underestimated the identity-by-descent scores determined using genetic profiles. The timing of the admixture and the effective population size trends estimated from the LD patterns reflected the documented history of the breed. Ancestry reconstruction methods identified more than 300 genes with excess of wolf ancestry compared to random expectations, mainly related to key morphological features, and more than 2000 genes with excess of dog ancestry, playing important roles in lipid metabolism, in the regulation of circadian rhythms, in learning and memory processes, and in sociability, such as the COMT gene, which has been described as a candidate gene for the latter trait in dogs.

CONCLUSIONS

In this study we successfully applied genome-wide procedures to reconstruct the history of the Czechoslovakian Wolfdog, assess individual wolf ancestry proportions and, thanks to the availability of a well-annotated reference genome, identify possible candidate genes for wolf-like and dog-like phenotypic traits typical of this breed, including commonly inherited disorders. Moreover, through the identification of ancestry-informative markers, these genomic approaches could provide tools for forensic applications to unmask illegal crossings with wolves and uncontrolled trades of recent and undeclared wolfdog hybrids.

Widespread, long-term admixture between grey wolves and domestic dogs across Eurasia and its implications for the conservation status of hybrids

Hybridisation between a domesticated species and its wild ancestor is an important conservation problem, especially if it results in the introgression of domestic gene variants into wild species. Nevertheless, the legal status of hybrids remains unregulated, partially because of the limited understanding of the hybridisation process and its consequences. The occurrence of hybridisation between grey wolves and domestic dogs is well documented from different parts of the wolf geographic range, but little is known about the frequency of hybridisation events, their causes and the genetic impact on wolf populations. We analysed 61K SNPs spanning the canid genome in wolves from across Eurasia and North America and compared that data to similar data from dogs to identify signatures of admixture. The haplotype block analysis, which included 38 autosomes and the X chromosome, indicated the presence of individuals of mixed wolf-dog ancestry in most Eurasian wolf populations, but less admixture was present in North American populations. We found evidence for male-biased introgression of dog alleles into wolf populations, but also identified a first-generation hybrid resulting from mating between a female dog and a male wolf. We found small blocks of dog ancestry in the genomes of 62% Eurasian wolves studied and melanistic individuals with no signs of recent admixed ancestry, but with a dog-derived allele at a locus linked to melanism. Consequently, these results suggest that hybridisation has been occurring in different parts of Eurasia on multiple timescales and is not solely a recent phenomenon. Nevertheless, wolf populations have maintained genetic differentiation from dogs, suggesting that hybridisation at a low frequency does not diminish distinctiveness of the wolf gene pool. However, increased hybridisation frequency may be detrimental for wolf populations, stressing the need for genetic monitoring to assess the frequency and distribution of individuals resulting from recent admixture.