Toward a genome-wide approach for detecting hybrids: informative SNPs to detect introgression between domestic cats and European wildcats (Felis silvestris)

Genome-wide variant analyses reveal new patterns of admixture and population structure in Australian dingoes

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

Detecting introgressive hybridization to maintain genetic integrity in endangered large waterbird: a case study in milky stork

Hybridization between milky stork (Mycteria cinerea) and painted stork (M. leucocephala) occurs frequently in captivity. Dusit Zoo is a captive breeding facility where storks with phenotypically ambiguous patterns have recently been observed, and their status remaining inconclusive. Here, we used a combination of phenotypic characters and genetic markers (cytochrome b and 14 microsatellite markers) to distinguish and identify hybrids from the two parental species (n = 114). Haplotype analysis revealed asymmetric mtDNA introgression from M. cinerea to M. leucocephala, with twelve morphologically classified M. leucocephala individuals carrying heterospecific mtDNA. Comprehensive biparental genetic assessments identified 33% of all three genetic clusters as admixed individuals, of which most were either F2 hybrids, backcrosses with M. leucocephala, or hybrids of unknown generation, implying weak premating isolation with the absence of intrinsic postzygotic isolation between parentals. Morphological analysis demonstrated that the absence or indistinctness of a black bar across the breast is the most noticeable trait to identify these hybrids. The endangered M. cinerea was found to have genomic contamination from M. leucocephala and vice versa, with at least 41 hybrid individuals being identified. These findings provide critical information for detecting hybrids and identifying suitable breeding stocks with genetic purity for future reintroduction and conservation management.

Genomic and machine learning-based screening of aquaculture-associated introgression into at-risk wild North American Atlantic salmon (Salmo salar) populations

The negative genetic impacts of gene flow from domestic to wild populations can be dependent on the degree of domestication and exacerbated by the magnitude of pre-existing genetic differences between wild populations and the domestication source. Recent evidence of European ancestry within North American aquaculture Atlantic salmon (Salmo salar) has elevated the potential impact of escaped farmed salmon on often at-risk wild North American salmon populations. Here, we compare the ability of single nucleotide polymorphism (SNP) and microsatellite (SSR) marker panels of different sizes (7-SSR, 100-SSR and 220K-SNP) to detect introgression of European genetic information into North American wild and aquaculture populations. Linear regression comparing admixture predictions for a set of individuals common to the three datasets showed that the 100-SSR panel and 7-SSR panels replicated the full 220K-SNP-based admixture estimates with low accuracy (r² of .64 and .49, respectively). Additional tests explored the effects of individual sample size and marker number, which revealed that ~300 randomly selected SNPs could replicate the 220K-SNP admixture predictions with greater than 95% fidelity. We designed a custom SNP panel (301-SNP) for European admixture detection in future monitoring work and then developed and tested a python package, salmoneuadmix (https://github.com/CNuge/SalmonEuAdmix), which uses a deep neural network to make de novo estimates of individuals' European admixture proportion without the need to conduct complete admixture analysis utilizing baseline samples. The results demonstrate the mobilization of targeted SNP panels and machine learning in support of at-risk species conservation and management.

Historic samples reveal loss of wild genotype through domestic chicken introgression during the Anthropocene

Human activities have precipitated a rise in the levels of introgressive gene flow among animals. The investigation of conspecific populations at different time points may shed light on the magnitude of human-mediated introgression. We used the red junglefowl Gallus gallus, the wild ancestral form of the chicken, as our study system. As wild junglefowl and domestic chickens readily admix, conservationists fear that domestic introgression into junglefowl may compromise their wild genotype. By contrasting the whole genomes of 51 chickens with 63 junglefowl from across their natural range, we found evidence of a loss of the wild genotype across the Anthropocene. When comparing against the genomes of junglefowl from approximately a century ago using rigorous ancient-DNA protocols, we discovered that levels of domestic introgression are not equal among and within modern wild populations, with the percentage of domestic ancestry around 20-50%. We identified a number of domestication markers in which chickens are deeply differentiated from historic junglefowl regardless of breed and/or geographic provenance, with eight genes under selection. The latter are involved in pathways dealing with development, reproduction and vision. The wild genotype is an allelic reservoir that holds most of the genetic diversity of G. gallus, a species which is immensely important to human society. Our study provides fundamental genomic infrastructure to assist in efforts to prevent a further loss of the wild genotype through introgression of domestic alleles.

Development of InDel markers for interspecific hybridization between hill pigeons and feral pigeons based on whole-genome re-sequencing

Interspecific hybridization occurs among birds, and closely related sister taxa tend to hybridize at a high rate. Genomic hybridization markers are useful for understanding the patterns and processes of hybridization and for conserving endangered species in captivity and the wild. In this study, we developed genomic hybridization markers for the F1 progeny of the sister taxa feral pigeons (Columba livia var. domestica) and endangered hill pigeons (Columba rupestris) (family Columbidae). Using whole-genome re-sequencing data, we performed genome-wide analysis for insertion/deletion (InDel) polymorphisms and validated using primers. We conducted polymerase chain reaction (PCR) and agarose gel electrophoresis to identify species-specific InDels. We produced eight F1 hybrids of hill and feral pigeons, and their samples were tested by re-performing analyses and sequencing using 11 species-specific InDel polymorphisms. Eight InDel markers simultaneously amplified two DNA fragments from all F1 hybrids, and there was no abnormality in the sequencing results. The application of genomic tools to detect hybrids can play a crucial role in the assessment of hybridization frequency in the wild. Moreover, systematic captive propagation efforts with hybrids can help control the population decline of hill pigeons.

Genetics of randomly bred cats support the cradle of cat domestication being in the Near East

Cat domestication likely initiated as a symbiotic relationship between wildcats (Felis silvestris subspecies) and the peoples of developing agrarian societies in the Fertile Crescent. As humans transitioned from hunter-gatherers to farmers ~12,000 years ago, bold wildcats likely capitalized on increased prey density (i.e., rodents). Humans benefited from the cats' predation on these vermin. To refine the site(s) of cat domestication, over 1000 random-bred cats of primarily Eurasian descent were genotyped for single-nucleotide variants and short tandem repeats. The overall cat population structure suggested a single worldwide population with significant isolation by the distance of peripheral subpopulations. The cat population heterozygosity decreased as genetic distance from the proposed cat progenitor's (F.s. lybica) natural habitat increased. Domestic cat origins are focused in the eastern Mediterranean Basin, spreading to nearby islands, and southernly via the Levantine coast into the Nile Valley. Cat population diversity supports the migration patterns of humans and other symbiotic species.

Species-informative SNP markers for characterising freshwater prawns of genus Macrobrachium in Cameroon

Single Nucleotide Polymorphisms (SNPs) are now popular for a myriad of applications in animal and plant species including, ancestry assignment, conservation genetics, breeding, and traceability of animal products. The objective of this study was to develop a customized cost-effective SNP panel for genetic characterisation of Macrobrachium species in Cameroon. The SNPs identified in a previous characterization study were screened as viable candidates for the reduced panel. Starting from a full set of 1,814 SNPs, a total of 72 core SNPs were chosen using conventional approaches: allele frequency differentials, minor allele frequency profiles, and Wright’s Fst statistics. The discriminatory power of reduced set of informative SNPs were then tested using the admixture analysis, principal component analysis, and discriminant analysis of principal components. The panel of prioritised SNP markers (i.e., N = 72 SNPs) distinguished Macrobrachium species with 100% accuracy. However, large sample size is needed to identify more informative SNPs for discriminating genetically closely related species, including M. macrobrachion versus M. vollenhovenii and M. sollaudii versus M. dux. Overall, the findings in this study show that we can accurately characterise Macrobrachium using a small set of core SNPs which could be useful for this economically important species in Cameroon. Given the results obtained in this study, a larger independent validation sample set will be needed to confirm the discriminative capacity of this SNP panel for wider commercial and research applications.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as “pure” polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.

Whole genome resequencing data enables a targeted SNP panel for conservation and aquaculture of Oreochromis cichlid fishes

Cichlid fish of the genus Oreochromis form the basis of the global tilapia aquaculture and fisheries industries. Broodstocks for aquaculture are often collected from wild populations, which in Africa may be from locations containing multiple Oreochromis species. However, many species are difficult to distinguish morphologically, hampering efforts to maintain good quality farmed strains. Additionally, non-native farmed tilapia populations are known to be widely distributed across Africa and to hybridize with native Oreochromis species, which themselves are important for capture fisheries. The morphological identification of these hybrids is particularly unreliable. Here, we describe the development of a single nucleotide polymorphism (SNP) genotyping panel from whole-genome resequencing data that enables targeted species identification in Tanzania. We demonstrate that an optimized panel of 96 genome-wide SNPs based on F_ST outliers performs comparably to whole genome resequencing in distinguishing species and identifying hybrids. We also show this panel outperforms microsatellite-based and phenotype-based classification methods. Case studies indicate several locations where introduced aquaculture species have become established in the wild, threatening native Oreochromis species. The novel SNP markers identified here represent an important resource for assessing broodstock purity in hatcheries and helping to conserve unique endemic biodiversity.

Consequences of Hybridization in Mammals: A Systematic Review

Hybridization, defined as breeding between two distinct taxonomic units, can have an important effect on the evolutionary patterns in cross-breeding taxa. Although interspecific hybridization has frequently been considered as a maladaptive process, which threatens species genetic integrity and survival via genetic swamping and outbreeding depression, in some cases hybridization can introduce novel adaptive variation and increase fitness. Most studies to date focused on documenting hybridization events and analyzing their causes, while relatively little is known about the consequences of hybridization and its impact on the parental species. To address this knowledge gap, we conducted a systematic review of studies on hybridization in mammals published in 2010-2021, and identified 115 relevant studies. Of 13 categories of hybridization consequences described in these studies, the most common negative consequence (21% of studies) was genetic swamping and the most common positive consequence (8%) was the gain of novel adaptive variation. The total frequency of negative consequences (49%) was higher than positive (13%) and neutral (38%) consequences. These frequencies are biased by the detection possibilities of microsatellite loci, the most common genetic markers used in the papers assessed. As negative outcomes are typically easier to demonstrate than positive ones (e.g., extinction vs hybrid speciation), they may be over-represented in publications. Transition towards genomic studies involving both neutral and adaptive variation will provide a better insight into the real impacts of hybridization.

Informing conservation strategies with museum genomics: Long-term effects of past anthropogenic persecution on the elusive European wildcat

Like many carnivore species, European wildcats (Felis silvestris) have suffered severe anthropogenic population declines in the past, resulting in a strong population bottleneck at the beginning of the 20th century. In Germany, the species has managed to survive its near extinction in small isolated areas and is currently recolonizing former habitats owing to legal protection and concerted conservation efforts. Here, we SNP-genotyped and mtDNA-sequenced 56 historical and 650 contemporary samples to assess the impact of massive persecution on genetic diversity, population structure, and hybridization dynamics of wildcats. Spatiotemporal analyses suggest that the presumed postglacial differentiation between two genetically distinct metapopulations in Germany is in fact the result of the anthropogenic bottleneck followed by re-expansion from few secluded refugia. We found that, despite the bottleneck, populations experienced no severe genetic erosion, nor suffered from elevated inbreeding or showed signs of increased hybridization with domestic cats. Our findings have significant implications for current wildcat conservation strategies, as the data analyses show that the two presently recognized wildcat population clusters should be treated as a single conservation unit. Although current populations appear under no imminent threat from genetic factors, fostering connectivity through the implementation of forest corridors will facilitate the preservation of genetic diversity and promote long-term viability. The present study documents how museum collections can be used as essential resource for assessing long-term anthropogenic effects on natural populations, for example, regarding population structure and the delineation of appropriate conservation units, potentially informing todays' species conservation.

A comparison of microsatellites and genome-wide SNPs for the detection of admixture brings the first molecular evidence for hybridization between Mustela eversmanii and M. putorius (Mustelidae, Carnivora)

Introgressive hybridization can pose a serious threat to endangered species which have an overlapping distribution such as in the case of two polecat species, Mustela eversmanii and M. putorius, in Europe. The population size of steppe polecat is known to continuously shrink, whereas its sister species, the European polecat, is still somehow widespread. In this study, we perform an analysis using microsatellite (SSR) and genomic (SNP) data sets to identify natural hybrids between polecats. Four populations were genotyped for eight polymorphic SSR loci, and thousands of unlinked SNPs were generated using a reduced-representation sequencing approach, RADseq, to characterize the genetic make-up of allopatric populations and to identify hybrids in the sympatric area. We applied standard population genetic analyses to characterize the populations based on their SSR allelic frequency. Only a single sample out of 48 sympatric samples showed exact intermediacy that we identified as an F1 hybrid. Additionally, one specimen was indicated in the genomic data sets as backcrossed. Other backcrosses, indicated by SSRs, were not validated by SNPs, which highlights the higher efficacy of the genomic method to identify backcrossed individuals. The low frequency of hybridization suggests that the difference in habitat preference of the two species may act as a barrier to admixture. Therefore, it is apparently unlikely that polecat populations are threatened by significant introgression. The two species showed a clear genetic differentiation using both techniques. We found higher genetic diversity values in the sympatric steppe polecat population than in the other studies on polecat populations. Although M. putorius is a hunted species in most countries, genetic diversity values indicate worse conditions in Europe than in the protected sibling species M. eversmanii. Suspending hunting and providing protected status of the former seems to be reasonable and timely.

Genome-Wide SNPs Detect Hybridisation of Marsupial Gliders (Petaurus breviceps breviceps × Petaurus norfolcensis) in the Wild

Petaurus breviceps and Petaurus norfolcensis have produced hybrids in captivity, however there are no reported cases of Petaurus hybridisation in the wild. This study uses morphological data, mitochondrial DNA, and nuclear genome-wide SNP markers to confirm P. breviceps breviceps × P. norfolcensis hybridisation within their natural range on the central coast of New South Wales, Australia. Morphological data identified a potential hybrid that was confirmed with next-generation sequencing technology and 10,111 genome-wide SNPs. Both STRUCTURE and NewHybrid analyses identified the hybrid as a P. norfolcensis backcross, which suggests an initial F1 hybrid was fertile. The mitochondrial DNA matched that of a P. b. breviceps, indicating that a P. b. breviceps female initially mated with a P. norfolcensis male to produce a fertile female offspring. Our study is an important example of how genome-wide SNPs can be used to identify hybrids where the distribution of congeners overlaps. Hybridisation between congeners is likely to become more frequent as climate changes and habitats fragment, resulting in increased interactions and competition for resources and mates.

Updating the Bibliography of Interbreeding among Canis in North America

This bibliography provides a collection of references that documents the evolution of studies evidencing interbreeding among Canis species in North America. Over the past several decades, advances in biology and genomic technology greatly improved our ability to detect and characterize species interbreeding, which has significance for understanding species in a changing landscape as well as for endangered species management. This bibliography includes a discussion within each category of interbreeding, the timeline of developing evidence, and includes a review of past research conducted on experimental crosses. Research conducted in the early 20th century is rich with detailed records and photographs of hybrid offspring development and behavior. With the progression of molecular methods, studies can estimate historical demographic parameters and detect chromosomal patterns of ancestry. As these methods continue to increase in accessibility, the field will gain a deeper and richer understanding of the evolutionary history of North American Canis.

Comparative Performance of Popular Methods for Hybrid Detection using Genomic Data

Interspecific hybridization is an important evolutionary phenomenon that generates genetic variability in a population and fosters species diversity in nature. The availability of large genome scale datasets has revolutionized hybridization studies to shift from the observation of the presence or absence of hybrids to the investigation of the genomic constitution of hybrids and their genome-specific evolutionary dynamics. Although a handful of methods have been proposed in an attempt to identify hybrids, accurate detection of hybridization from genomic data remains a challenging task. In addition to methods that infer phylogenetic networks or that utilize pairwise divergence, site pattern frequency based and population genetic clustering approaches are popularly used in practice, though the performance of these methods under different hybridization scenarios has not been extensively examined. Here, we use simulated data to comparatively evaluate the performance of four tools that are commonly used to infer hybridization events: the site pattern frequency based methods HyDe and the D-statistic (i.e., the ABBA-BABA test) and the population clustering approaches structure and ADMIXTURE. We consider single hybridization scenarios that vary in the time of hybridization and the amount of incomplete lineage sorting (ILS) for different proportions of parental contributions (γ); introgressive hybridization; multiple hybridization scenarios; and a mixture of ancestral and recent hybridization scenarios. We focus on the statistical power to detect hybridization and the false discovery rate (FDR) for comparisons of the D-statistic and HyDe, and the accuracy of the estimates of γ as measured by the mean squared error for HyDe, structure, and ADMIXTURE. Both HyDe and the D-statistic are powerful for detecting hybridization in all scenarios except those with high ILS, although the D-statistic often has an unacceptably high FDR. The estimates of γ in HyDe are impressively robust and accurate whereas structure and ADMIXTURE sometimes fail to identify hybrids, particularly when the proportional parental contributions are asymmetric (i.e., when γ is close to 0). Moreover, the posterior distribution estimated using structure exhibits multimodality in many scenarios, making interpretation difficult. Our results provide guidance in selecting appropriate methods for identifying hybrid populations from genomic data.

Biased assessment of ongoing admixture using STRUCTURE in the absence of reference samples

Detection of hybridization and introgression is important in ecological research as in conservation and evolutionary biology. STRUCTURE is one of the most popular software to study introgression and allows estimating what proportion of the genome of each individual belongs to each ancestral population, even in cases where no reference sample from the ancestral nonadmixed populations is previously identified. In spite of its frequent use, some studies have indicated that ancestry estimates may not always be reliable. We simulated population data under different conditions with regard to the genetic differentiation between ancestral populations, number of loci considered, number of alleles per marker and hybridization rate, and analysed data with STRUCTURE. When reference samples were not included, the comparison of the known degree of admixture for each simulated individual and the value estimated with STRUCTURE revealed a strong underestimation of the level of introgression, classifying many admixed individuals as nonadmixed. This derives from an inaccurate estimation of the ancestral allele frequencies. When samples from the nonadmixed ancestral population were included as reference in the analyses, the bias in the estimations was reduced. The most accurate estimates were obtained when potentially admixed samples were few in relation to reference samples. Thus, whenever possible, a very large proportion of nonadmixed reference samples should be included in admixture assessments and different approaches should be combined. The misestimate of the amount of introgression can impair our understanding of the evolutionary history of species and misguide conservation efforts.

Protein metabolism and physical fitness are physiological determinants of body condition in Southern European carnivores

The physiological significance of biometric body condition indices (bBCI) is poorly understood. We hypothesized that bBCI are composite metrics of nutritional physiology, physical fitness and health. To test this hypothesis, we first compared the performance of eight bBCI, using 434 Southern European carnivores from six species as a model system; and then identified, by non-destructive methods, the hematology and serum biochemistry correlates of three selected bBCI. Fulton’s K Index, Major Axis Regression Residuals and Scaled Mass Index were the only bBCI insensitive to the effect of sex and age. The most informative physiological parameters in explaining the variation of these bBCI were the albumin (Effect Size (ES) = − 1.66 to − 1.76), urea (ES = 1.61 to 1.85) and total bilirubin (ES = − 1.62 to − 1.79). Hemoglobin and globulins (positive) and cholesterol (negative) were moderately informative (0.9 <|ES|< 1.5). This study shows that most bBCI do not control for the effect of age and sex in Southern European carnivores. Our results support that bBCI are composite measures of physiologic processes, reflecting a positive gradient from protein-poor to protein-rich diets, accompanied by increased physical fitness. Biometric body condition indices allow the integration of ecologically relevant physiological aspects in an easily obtained metric.

Harmonizing hybridization dissonance in conservation

A dramatic increase in the hybridization between historically allopatric species has been induced by human activities. However, the notion of hybridization seems to lack consistency in two respects. On the one hand, it is inconsistent with the biological species concept, which does not allow for interbreeding between species, and on the other hand, it is considered either as an evolutionary process leading to the emergence of new biodiversity or as a cause of biodiversity loss, with conservation implications. In the first case, we argue that conservation biology should avoid the discussion around the species concept and delimit priorities of conservation units based on the impact on biodiversity if taxa are lost. In the second case, we show that this is not a paradox but an intrinsic property of hybridization, which should be considered in conservation programmes. We propose a novel view of conservation guidelines, in which human-induced hybridization may also be a tool to enhance the likelihood of adaptation to changing environmental conditions or to increase the genetic diversity of taxa affected by inbreeding depression. The conservation guidelines presented here represent a guide for the development of programmes aimed at protecting biodiversity as a dynamic evolutionary system.

Species assignment in forensics and the challenge of hybrids

Forensic identification of species is in growing demand, particularly from law enforcement authorities in the areas of wildlife, fisheries and hunting as well as food authentication. Within the non-human forensic genetics expanding applications' field, the major current difficulties result from the lack of standards and genetic databases as well as the poor or absent taxonomic definition of several groups. Here we focus on a forensically important and overlooked problem in species identification: the exclusive use of uniparental markers, a common practice in current genetic barcoding methodologies, may lead to incorrect or impossible assignment whenever hybrids can occur (frequently, not only in domesticates, but also in the wild). For example, if one of these cases involves a mammal, and mitochondrial DNA alone is used (which in instances may be the only type of DNA sequence available in databases), the sample will be wrongfully assigned to the female parental species, completely missing the detection of a possible hybrid animal. The importance of this issue in the forensic contributions to food authentication, wildlife and conservation genetics is analyzed. We present a cautionary guidance on the forensic reporting of results avoiding this error.

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.

Identification of Ancestry Informative Marker (AIM) Panels to Assess Hybridisation between Feral and Domestic Sheep

Simple Summary

Once present in the entirety of Europe, mouflon (wild sheep) became extinct due to intense hunting, but remnant populations survived and became feral on the Mediterranean islands of Corsica and Sardinia. Although now protected by regional laws, Sardinian mouflon is threatened by crossbreeding with domestic sheep causing genetic hybridisation. The spread of domestic genes can be detrimental for wild populations as it dilutes the genetic features that characterise them. This work aimed to identify diagnostic tools that could be applied to monitor the level of hybridisation between mouflon and domestic sheep. Tens of thousands of genetic markers known as single nucleotide polymorphisms (SNPs) were screened and we identified the smallest number of SNPs necessary to discriminate between pure mouflon and sheep. We produced four SNP panels of different sizes which were able to assess the hybridisation level of a mouflon and we verified that the SNP panels efficacy is independent of the domestic sheep breed involved in the hybrid. The implementation of these results into actual diagnostic tools will help the conservation of this unique and irreplaceable mouflon population, and the methodology applied can easily be transferred to other case studies of interest.

Abstract

Hybridisation of wild populations with their domestic counterparts can lead to the loss of wildtype genetic integrity, outbreeding depression, and loss of adaptive features. The Mediterranean island of Sardinia hosts one of the last extant autochthonous European mouflon (Ovis aries musimon) populations. Although conservation policies, including reintroduction plans, have been enforced to preserve Sardinian mouflon, crossbreeding with domestic sheep has been documented. We identified panels of single nucleotide polymorphisms (SNPs) that could act as ancestry informative markers able to assess admixture in feral x domestic sheep hybrids. The medium-density SNP array genotyping data of Sardinian mouflon and domestic sheep (O. aries aries) showing pure ancestry were used as references. We applied a two-step selection algorithm to this data consisting of preselection via Principal Component Analysis followed by a supervised machine learning classification method based on random forest to develop SNP panels of various sizes. We generated ancestry informative marker (AIM) panels and tested their ability to assess admixture in mouflon x domestic sheep hybrids both in simulated and real populations of known ancestry proportions. All the AIM panels recorded high correlations with the ancestry proportion computed using the full medium-density SNP array. The AIM panels proposed here may be used by conservation practitioners as diagnostic tools to exclude hybrids from reintroduction plans and improve conservation strategies for mouflon 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.

Range-wide patterns of human-mediated hybridisation in European wildcats

Hybridisation between wild taxa and their domestic congeners is a significant conservation issue. Domestic species frequently outnumber their wild relatives in population size and distribution and may therefore genetically swamp the native species. The European wildcat (Felis silvestris) has been shown to hybridise with domestic cats (Felis catus). Previously suggested spatially divergent introgression levels have not been confirmed on a European scale due to significant differences in the applied methods to assess hybridisation of the European wildcat. We analysed 926 Felis spp. samples from 13 European countries, using a set of 86 selected ancestry-informative SNPs, 14 microsatellites, and ten mitochondrial and Y-chromosome markers to study regional hybridisation and introgression patterns and population differentiation. We detected 51 hybrids (four F1 and 47 F2 or backcrosses) and 521 pure wildcats throughout Europe. The abundance of hybrids varied considerably among studied populations. All samples from Scotland were identified as F2 hybrids or backcrosses, supporting previous findings that the genetic integrity of that wildcat population has been seriously compromised. In other European populations, low to moderate levels of hybridisation were found, with the lowest levels being in Central and Southeast Europe. The occurrence of distinct maternal and paternal markers between wildcat and domestic cat suggests that there were no severe hybridisation episodes in the past. The overall low (< 1%) prevalence of F1 hybrids suggests a low risk of hybridisation for the long-term genetic integrity of the wildcat in most of Europe. However, regionally elevated introgression rates confirm that hybridisation poses a potential threat. We propose regional in-depth monitoring of hybridisation rates to identify factors driving hybridisation so as to develop effective strategies for conservation.

Assessing Risks to Wildlife from Free-Roaming Hybrid Cats: The Proposed Introduction of Pet Savannah Cats to Australia as a Case Study

Simple Summary

The domestic cat, Felis catus, is often cross-bred with other species in the cat family to produce hybrid or ‘designer’ cats that are sought by people as pets. However, hybrid cats are often surrendered to wildlife shelters, or released, which leads to concern that they may establish free-roaming populations and damage native wildlife. In 2008, the Australian government rejected an application, on precautionary grounds, to import savannah cats (hybrids of the domestic cat and serval Leptailurus serval) into the country. We review the limited information informing this decision and then present a framework that identifies the native mammal species likely to have been most at risk of predation from savannah cats if importation and establishment had occurred. Assuming that savannah cats hunt similar prey to those that are hunted by both parent species, we estimate that 91% of Australia’s extant terrestrial mammal fauna would likely face some risk of predation from savannah cats, including 93% of non-volant mammal species that have threatened conservation status. The framework results strongly validate the decision to ban savannah cats from Australia. We suggest that our framework approach could be adapted to assess the likely risks that are posed by the arrival of other hybrid cats or hybrids of other animals.

Abstract

Hybrid cats—created by crossing different species within the family Felidae—are popular pets, but they could potentially threaten native species if they escape and establish free-roaming populations. To forestall this possibility, the Australian government imposed a specific ban on importation of the savannah cat, a hybrid created by crossing the domestic cat Felis catus and serval Leptailurus serval, in 2008. We develop a decision–framework that identifies those species of non-volant native mammals in Australia that would likely have been susceptible to predation by savannah cats if importation and establishment had occurred. We assumed that savannah cats would hunt ecologically similar prey to those that are depredated by both the domestic cat and the serval, and categorised native mammals as having different levels of susceptibility to predation by savannah cats based on their size, habitat range, and behaviour. Using this framework, we assessed savannah cats as likely to add at least 28 extant native mammal species to the 168 that are known already to be susceptible to predation by the domestic cat, posing a risk to 91% of Australia’s extant non-volant terrestrial mammal species (n = 216) and to 93% of threatened mammal species. The framework could be generalised to assess risks from any other hybrid taxa.

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.

Distinguishing the victim from the threat: SNP-based methods reveal the extent of introgressive hybridization between wildcats and domestic cats in Scotland and inform future in situ and ex situ management options for species restoration

The degree of introgressive hybridization between the Scottish wildcat and domestic cat has long been suspected to be advanced. Here, we use a 35-SNP-marker test, designed to assess hybridization between wildcat and domestic cat populations in Scotland, to assess a database of 295 wild-living and captive cat samples, and test the assumptions of the test using 3,097 SNP markers generated independently in a subset of the data using ddRAD. We discovered that despite increased genetic resolution provided by these methods, wild-living cats in Scotland show a complete genetic continuum or hybrid swarm structure when judged against reference data. The historical population of wildcats, although hybridized, clearly groups at one end of this continuum, as does the captive population of wildcats. The interpretation of pelage scores against nuclear genetic data continues to be problematic. This is probably because of a breakdown in linkage equilibrium between wildcat pelage genes as the two populations have become increasingly mixed, meaning that pelage score or SNP score alone is poor diagnostic predictors of hybrid status. Until better tools become available, both should be used jointly, where possible, when making management decisions about individual cats. We recommend that the conservation community in Scotland must now define clearly what measures are to be used to diagnose a wildcat in the wild in Scotland, if future conservation action is to be effective.

Host specificity driving genetic structure and diversity in ectoparasite populations: Coevolutionary patterns in Apodemus mice and their lice

A degree of host specificity, manifested by the processes of host-parasite cospeciations and host switches, is assumed to be a major determinant of parasites' evolution. To understand these patterns and formulate appropriate ecological hypotheses, we need better insight into the coevolutionary processes at the intraspecific level, including the maintenance of genetic diversity and population structure of parasites and their hosts. Here, we address these questions by analyzing large-scale molecular data on the louse Polyplax serrata and its hosts, mice of the genus Apodemus, across a broad range of European localities. Using mitochondrial DNA sequences and microsatellite data, we demonstrate the general genetic correspondence of the Apodemus/Polyplax system to the scenario of the postglacial recolonization of Europe, but we also show several striking discrepancies. Among the most interesting are the evolution of different degrees of host specificity in closely related louse lineages in sympatry, or decoupled population structures of the host and parasites in central Europe. We also find strong support for the prediction that parasites with narrower host specificity possess a lower level of genetic diversity and a deeper pattern of interpopulation structure as a result of limited dispersal and smaller effective population size.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Hindering the illegal trade in dog and cat furs through a DNA-based protocol for species identification

In Western countries dogs and cats are the most popular pets, and people are increasingly opposed to their rearing for the fur industry. In 2007, a Regulation of the European Union (EU) banned the use and trade of dog and cat furs, but an official analytical protocol to identify them as source species was not provided, and violations of law are still frequent in all Member States. In this paper we report on the development and validation of a simple and affordable DNA method for species detection in furs to use as an effective tool to combat illegal trade in fur products. A set of mitochondrial primers was designed for amplification of partial cytochrome b, control region and ND1 gene in highly degraded samples, like furs and pelts. Our amplification workflow involved the use of a non-specific primer pair to perform a first test to identify the species through sequencing, then the application of species-specific primer pairs to use in singleplex end-point PCRs as confirmation tests. The advantage of this two-step procedure is twofold: on the one hand it minimises the possibility of negative test results from degraded samples, since failure of amplification with a first set of primers can be offset by successful amplification of the second, and on the other it adds confidence and reliability to final authentication of species. All designed primers were validated on a reference collection of tissue samples, obtaining solid results in terms of specificity, sensitivity, repeatability and reproducibility. Application of the protocol on real caseworks from seized furs yielded successful results also from old and dyed furs, suggesting that age and chemical staining do not necessarily affect positive amplifications. Major pros of this approach are: (1) sensitive and informative primer sets for detection of species; (2) short PCR amplicons for the analysis of poor quality DNA; (3) binding primers that avoid contamination from human DNA; (4) user-friendly protocol for any laboratory equipped for analysis of low-copy-number DNA. Our molecular procedure proved to be a good starting point for enforcing the EU Regulation against dog and cat fur trade in forensic contexts where source attribution is essential to the assignment of responsibilities.

Applications and efficiencies of the first cat 63K DNA array

The development of high throughput SNP genotyping technologies has improved the genetic dissection of simple and complex traits in many species including cats. The properties of feline 62,897 SNPs Illumina Infinium iSelect DNA array are described using a dataset of over 2,000 feline samples, the most extensive to date, representing 41 cat breeds, a random bred population, and four wild felid species. Accuracy and efficiency of the array’s genotypes and its utility in performing population-based analyses were evaluated. Average marker distance across the array was 37,741 Kb, and across the dataset, only 1% (625) of the markers exhibited poor genotyping and only 0.35% (221) showed Mendelian errors. Marker polymorphism varied across cat breeds and the average minor allele frequency (MAF) of all markers across domestic cats was 0.21. Population structure analysis confirmed a Western to Eastern structural continuum of cat breeds. Genome-wide linkage disequilibrium ranged from 50–1,500 Kb for domestic cats and 750 Kb for European wildcats (Felis silvestris silvestris). Array use in trait association mapping was investigated under different modes of inheritance, selection and population sizes. The efficient array design and cat genotype dataset continues to advance the understanding of cat breeds and will support monogenic health studies across feline breeds and populations.

hybriddetective: A workflow and package to facilitate the detection of hybridization using genomic data in r

The ability to detect and characterize hybridization in nature has long been of interest to many fields of biology and often has direct implications for wildlife management and conservation. The capacity to identify the presence of hybridization, and quantify the numbers of individuals belonging to different hybrid classes, permits inference on the magnitude of, and timescale over which, hybridization has been or is occurring. Here, we present an r package and associated workflow developed for the detection, with estimates of efficiency and accuracy, of multigenerational hybrid individuals using genetic or genomic data in conjunction with the program newhybrids. This package includes functions for the identification and testing of diagnostic panels of markers, the simulation of multigenerational hybrids, and the quantification and visualization of the efficiency and accuracy with which hybrids can be detected. Overall, this package delivers a streamlined hybrid analysis platform, providing improvements in speed, ease of use and repeatability over current ad hoc approaches. The latest version of the package and associated documentation are available on GitHub (https://github.com/bwringe/hybriddetective).

Assessing patterns of genetic admixture between sheep breeds: Case study in Algeria

In developing countries, cross-breeding between local breeds and indigene or exotic breeds represents one of the main threats to the livestock diversity, leading to genetic dilution and loss of unique allelic combination underlying essential local adaptive traits. In this study, two Algerian sheep breeds, known to be highly admixed, were considered as a case study, to demonstrate how combination of different methodologies coupled with the use of specific softwares can be efficient to assess the spatial structuration of a hybrid zone, even in a case of extreme admixture. A fine sampling covering distribution areas of both breeds was implemented in order to study the admixture area and adjacent zones from a phenotypic (i.e., 19 quantitative traits were considered) and a genetic point of view (i.e., 21 microsatellites markers were used). Both approaches gave concordant patterns, highlighting areas with sheep most differentiated (or less admixed) for each breed. In detail, the region of Biskra appeared as the most preserved for the Ouled-Djellal breed and the northwest of Laghouat was identified as the most preserved area for the Rembi breed. The approach proposed in the study offers a low-cost solution to identify the most representative flocks of a breed, allowing the implementation of efficient conservation plans.

European wildcat populations are subdivided into five main biogeographic groups: consequences of Pleistocene climate changes or recent anthropogenic fragmentation?

Extant populations of the European wildcat are fragmented across the continent, the likely consequence of recent extirpations due to habitat loss and over‐hunting. However, their underlying phylogeographic history has never been reconstructed. For testing the hypothesis that the European wildcat survived the Ice Age fragmented in Mediterranean refuges, we assayed the genetic variation at 31 microsatellites in 668 presumptive European wildcats sampled in 15 European countries. Moreover, to evaluate the extent of subspecies/population divergence and identify eventual wild × domestic cat hybrids, we genotyped 26 African wildcats from Sardinia and North Africa and 294 random‐bred domestic cats. Results of multivariate analyses and Bayesian clustering confirmed that the European wild and the domestic cats (plus the African wildcats) belong to two well‐differentiated clusters (average Ф_ST = 0.159, r_st = 0.392, P > 0.001; Analysis of molecular variance [AMOVA]). We identified from c. 5% to 10% cryptic hybrids in southern and central European populations. In contrast, wild‐living cats in Hungary and Scotland showed deep signatures of genetic admixture and introgression with domestic cats. The European wildcats are subdivided into five main genetic clusters (average Ф_ST = 0.103, r_st = 0.143, P > 0.001; AMOVA) corresponding to five biogeographic groups, respectively, distributed in the Iberian Peninsula, central Europe, central Germany, Italian Peninsula and the island of Sicily, and in north‐eastern Italy and northern Balkan regions (Dinaric Alps). Approximate Bayesian Computation simulations supported late Pleistocene–early Holocene population splittings (from c. 60 k to 10 k years ago), contemporary to the last Ice Age climatic changes. These results provide evidences for wildcat Mediterranean refuges in southwestern Europe, but the evolution history of eastern wildcat populations remains to be clarified. Historical genetic subdivisions suggest conservation strategies aimed at enhancing gene flow through the restoration of ecological corridors within each biogeographic units. Concomitantly, the risk of hybridization with free‐ranging domestic cats along corridor edges should be carefully monitored.