Abstract Hybridisation between wild and domestic taxa can favour the spread of domestic alleles into wild populations through backcrossing. The complex interplay of random genetic drift, recombination, and selection can shape the fate of introgressed alleles. Maladaptive domestic variants are likely to be purged by natural selection, but others may persist across generations. It has long been known that the Apennine Italian wolf population, exposed to large numbers of free-ranging dogs, has experienced extensive introgression. The unusually high frequency of black wolves observed in Italy, compared to other European populations, may parallel patterns documented in North American wolves, where the melanistic K B allele at the CBD103 gene, of domestic origin, has spread over thousands of years of introgression. We tested whether the K B mutation entered the peninsular Italian wolf population via hybridisation and spread through adaptive introgression. Genome-wide analyses of black and wild-type (grey-coated) Apennine wolves showed no clear signatures of recent dog ancestry in most melanistic animals. Our ancestry reconstruction approaches identified two distinct K B haplogroups of domestic origin, suggesting multiple introgression events. Notably, we found molecular evidence consistent with balancing selection on the K B haplotypes, whose functional role, nonetheless, warrants further research. Therefore, the microevolutionary genomic and ecological consequences of wolf-dog hybridisation in Italy should be carefully investigated to inform appropriate science-based conservation management strategies.

Anthropogenic hybridisation between wild and domestic taxa poses a significant threat to species integrity, including the endangered European wildcat. To enable reliable molecular assessment of admixture with domestic cats and to increase the accuracy of hybrid class assignment we optimised an existing reduced microfluidic 96 Single Nucleotide Polymorphism (SNP) panel. We selected SNPs from a genome-wide dataset for maximum FST between both taxa and replaced 60 SNPs from the previous 96 SNP panel. Comparison of both panels based on simulated hybrid genotypes and real-world genotypes proof the higher discriminatory power of the optimised panel, which allows for reliable assignment of F1 and F2 hybrids, as well as 1st and 2nd generation backcrosses. Additionally, we successfully tested the panel for both tissue and non-invasively collected hair samples, demonstrating the suitability of the new panel for implementation in wildcat monitoring programmes.

RK Wayne has arguably been the most influential geneticist of canids, famously promoting the conservation of wolves in his homeland, the United States. His influence has been felt in other countries and regions outside the contiguous United States, where he inspired others, also including former graduate students and research fellows of his, to use modern molecular techniques to examine the evolutionary biology of canids to inform the conservation and management of wolves. In this review, we focus on the implications of Wayne's work on wolves outside the United States. He envisioned a clear future for wolf conservation research, involving the study of wolves' ecological and genetic diversity, and the description of ecotypes requiring conservation. He also documented widespread hybridization among canids and introgression of DNA from domestic dogs to wolves, a process that started dozens of thousands of years ago. His work therefore calls for innovative studies, such as examining the potential fitness benefits of introgression. Inspired by his results, for example, on the purging of deleterious alleles in small populations, wolf researchers should use novel molecular tools to challenge other conservation genetics paradigms. Overall, RK Wayne's work constitutes a call for answers, which as scientists or citizens concerned with conservation matters, we are obliged to address, as we contribute to monitoring and maintaining biodiversity during our period of dramatic transformations of the biosphere.

Wolves (Canis lupus) are generally monitored by visual observations, camera traps, and DNA traces. In this study, we evaluated acoustic monitoring of wolf howls as a method for monitoring wolves, which may permit detection of wolves across longer distances than that permitted by camera traps. We analyzed acoustic data of wolves’ howls collected from both wild and captive ones. The analysis focused on individual and subspecies recognition. Furthermore, we aimed to determine the usefulness of acoustic monitoring in the field given the limited data for Eurasian wolves. We analyzed 170 howls from 16 individual wolves from 3 subspecies: Arctic (Canis lupus arctos), Eurasian (C. l. lupus), and Northwestern wolves (C. l. occidentalis). Variables from the fundamental frequency (f0) (lowest frequency band of a sound signal) were extracted and used in discriminant analysis, classification matrix, and pairwise post-hoc Hotelling test. The results indicated that Arctic and Eurasian wolves had subspecies identifiable calls, while Northwestern wolves did not, though this sample size was small. Identification on an individual level was successful for all subspecies. Individuals were correctly classified with 80%–100% accuracy, using discriminant function analysis. Our findings suggest acoustic monitoring could be a valuable and cost-effective tool that complements camera traps, by improving long-distance detection of wolves.

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 ≤ Qi ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Qi ≤ 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.

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.

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

The Sardinian dhole ( Cynotherium sardous )[1][1] was an iconic and unique canid species of canid that was endemic of Sardinia and Corsica until it became extinct at the end of the Late Pleistocene[2][2]–[5][3]. Given its peculiar dental morphology, small body size and high level of endemism, several canids have been proposed as possible ancestors of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor [3][4],[6][5]–[9][6]. Morphometric analyses[3][4],[6][5],[8][7]–[12][8] have failed to clarify the evolutionary relationship with other canids. We sequenced the genome of a ca 21,100 year old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found it represents a separate taxon from all other living canids from Eurasia, Africa and North America, and that the Sardinian and Asian dhole lineages diverged ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, that ended approximately 500-300 ka, when the landbridge between Sardinia and mainland Italy was broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids - probably a result of the long-term isolation - that could have contributed to the subsequent extinction of the Sardinian dhole. ### Competing Interest Statement The authors have declared no competing interest. [1]: #ref-1 [2]: #ref-2 [3]: #ref-5 [4]: #ref-3 [5]: #ref-6 [6]: #ref-9 [7]: #ref-8 [8]: #ref-12

Local adaptation of animals to the environment can abruptly become a burden when faced with rapid climatic changes such as those foreseen for the Italian peninsula over the next 70 years. Our study investigates the genetic structure of the Italian goat populations and links it with the environment and how genetics might evolve over the next 50 years. We used one of the largest national datasets including > 1000 goats from 33 populations across the Italian peninsula collected by the Italian Goat Consortium and genotyped with over 50 k markers. Our results showed that Italian goats can be discriminated in three groups reflective of the Italian geography and its geo-political situation preceding the country unification around two centuries ago. We leveraged the remarkable genetic and geographical diversity of the Italian goat populations and performed landscape genomics analysis to disentangle the relationship between genotype and environment, finding 64 SNPs intercepting genomic regions linked to growth, circadian rhythm, fertility, and inflammatory response. Lastly, we calculated the hypothetical future genotypic frequencies of the most relevant SNPs identified through landscape genomics to evaluate their long-term effect on the genetic structure of the Italian goat populations. Our results provide an insight into the past and the future of the Italian local goat populations, helping the institutions in defining new conservation strategy plans that could preserve their diversity and their link to local realities challenged by climate change.

The Sardinian dhole (Cynotherium sardous)1 was an iconic and unique canid species that was endemic to Sardinia and Corsica until it became extinct at the end of the Late Pleistocene.2-5 Given its peculiar dental morphology, small body size, and high level of endemism, several extant canids have been proposed as possible relatives of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor.3,6-9 Morphometric analyses3,6,8-12 have failed to clarify the evolutionary relationship with other canids.We sequenced the genome of a ca-21,100-year-old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found that it represents a separate taxon from all other living canids from Eurasia, Africa, and North America, and that the Sardinian dhole lineage diverged from the Asian dhole ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, which ended approximately 500-300 ka, when the land bridge between Sardinia and mainland Italy was already broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids-probably a result of the long-term isolation-that could have contributed to the subsequent extinction of the Sardinian dhole.

Background The low cost and rapidity of microsatellite analysis have led to the development of several markers for many species. Because in non-invasive genetics it is recommended to genotype individuals using few loci, generally a subset of markers is selected. The choice of different marker panels by different research groups studying the same population can cause problems and bias in data analysis. A priority issue in conservation genetics is the comparability of data produced by different labs with different methods. Here, we compared data from previous and ongoing studies to identify a panel of microsatellite loci efficient for the long-term monitoring of Apennine brown bears ( Ursus arctos marsicanus ), aiming at reducing genotyping uncertainty and allowing reliable individual identifications overtimes. Results We examined all microsatellite markers used up to now and identified 19 candidate loci. We evaluated the efficacy of 13 of the most commonly used loci analyzing 194 DNA samples belonging to 113 distinct bears selected from the Italian national biobank. We compared data from 4 different marker subsets on the basis of genotyping errors, allelic patterns, observed and expected heterozygosity, discriminatory powers, number of mismatching pairs, and probability of identity. The optimal marker set was selected evaluating the low molecular weight, the high discriminatory power, and the low occurrence of genotyping errors of each primer. We calibrated allele calls and verified matches among genotypes obtained in previous studies using the complete set of 13 STRs (Short Tandem Repeats), analyzing six invasive DNA samples from distinct individuals. Differences in allele-sizing between labs were consistent, showing a substantial overlap of the individual genotyping. Conclusions The proposed marker set comprises 11 Ursus specific markers with the addition of cxx20, the canid-locus less prone to genotyping errors, in order to prevent underestimation (maximizing the discriminatory power) and overestimation (minimizing the genotyping errors) of the number of Apennine brown bears. The selected markers allow saving time and costs with the amplification in multiplex of all loci thanks to the same annealing temperature. Our work optimizes the available resources by identifying a shared panel and a uniform methodology capable of improving comparisons between past and future studies.

The grey wolf (Canis lupus) is one of the most widely distributed mammals in which a variety of distinct populations have been described. However, given their currently fragmented distribution and recent history of human-induced population decline, little is known about the events that led to their differentiation. Based on the analysis of whole canid genomes, we examined the divergence times between Southern European wolf populations and their ancient demographic history. We found that all present-day Eurasian wolves share a common ancestor ca 36 000 years ago, supporting the hypothesis that all extant wolves derive from a single population that subsequently expanded after the Last Glacial Maximum. We also estimated that the currently isolated European populations of the Iberian Peninsula, Italy and the Dinarics-Balkans diverged very closely in time, ca 10 500 years ago, and maintained negligible gene flow ever since. This indicates that the current genetic and morphological distinctiveness of Iberian and Italian wolves can be attributed to their isolation dating back to the end of the Pleistocene, predating the recent human-induced extinction of wolves in Central Europe by several millennia.

We revise the taxonomy of the Sylvia cantillans complex, a group of phenotypically distinct warblers with mainly parapatric distributions around a large part of the Mediterranean basin. We redefine the species limits using a combination of mitochondrial and nuclear markers and we objectively link available names to the genetically defined lineages by genotyping the surviving type specimens. In addition, the study of archival documents clarifies the exact composition of type series and provides further evidence for the identification of lost types. These results support the recognition of three species-level taxa: Moltoni’s warbler, Sylvia subalpina (north-central Italy, Corsica, Sardinia and the Balearics); the western subalpine warbler, S. iberiae (North Africa, Iberia, southern France and extreme north-west Italy); and the eastern subalpine warbler, S. cantillans, with subspecies S. cantillans cantillans (southern Italy, Sicily) and S. cantillans albistriata (Balkans, Greece, western Turkey).

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.

Next Generation Sequencing (NGS) and related technologies have revolutionized the field of conservation and population genetics, providing novel tools and the capacity to discover thousands of new Single Nucleotide Polymorphisms (SNPs) for the analysis of population parameters. However, gathering NGS data for organisms with very large genomes, such as amphibians, remains challenging because it is still unclear how the current methods perform. Here, we use the Genotyping-by-Sequencing (GBS) approach to generate SNP data for the genotyping of two amphibian species that are of conservation concern, the Sardinian brook salamander (Euproctus platycephalus) and the Italian stream frog (Rana italica). Both E. platycephalus and R. italica have very large genomes (5.53 Gb and >20 Gb, respectively) so genomic data are not available for either of them. We used 95 individual samples and one Illumina lane for each species, with an additional lane for E. platycephalus. After filtering, we obtained 961 and 854 high-coverage SNPs for E. platycephalus and R. italica, respectively. Our results suggest that GBS can serve as a reliable and cost-effective method for genotyping large amphibian genomes, including non-model species.

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.