Estimating Wolf Population Size and Dynamics by Field Monitoring and Demographic Models: Implications for Management and Conservation.

We estimated the current size and dynamics of the wolf population in Tuscany and investigated the trends and demographic drivers of population changes. Estimates were obtained by two different approaches: (i) mixed-technique field monitoring (from 2014 to 2016) that found the minimum observed pack number and estimated population size, and (ii) an individual-based model (run by Vortex software v. 10.3.8.0) with demographic inputs derived from a local intensive study area and historic data on population size. Field monitoring showed a minimum population size of 558 wolves (SE = 12.005) in 2016, with a density of 2.74 individuals/100 km2. The population model described an increasing trend with an average annual rate of increase λ = 1.075 (SE = 0.014), an estimated population size of about 882 individuals (SE = 9.397) in 2016, and a density of 4.29 wolves/100 km2. Previously published estimates of wolf population were as low as 56.2% compared to our field monitoring estimation and 34.6% in comparison to our model estimation. We conducted sensitivity tests to analyze the key parameters driving population changes based on juvenile and adult mortality rates, female breeding success, and litter size. Mortality rates played a major role in determining intrinsic growth rate changes, with adult mortality accounting for 62.5% of the total variance explained by the four parameters. Juvenile mortality was responsible for 35.8% of the variance, while female breeding success and litter size had weak or negligible effects. We concluded that reliable estimates of population abundance and a deeper understanding of the role of different demographic parameters in determining population dynamics are crucial to define and carry out appropriate conservation and management strategies to address human-wildlife conflicts.

Anthropogenic hybridization and its influence on the adaptive potential of the Sardinian wild boar (Sus scrofa meridionalis).

The wild boar (Sus scrofa meridionalis) arrived in Sardinia with the first human settlers in the early Neolithic with the potential to hybridize with the domestic pig (S. s. domesticus) throughout its evolution on the island. In this paper, we investigated the possible microevolutionary effects of such introgressive hybridization on the present wild boar population, comparing Sardinian wild specimens with several commercial pig breeds and Sardinian local pigs, along with a putatively unadmixed wild boar population from Central Italy, all genotyped with a medium density SNP chip. We first aimed at identifying hybrids in the population using different approaches, then examined genomic regions enriched for domestic alleles in the hybrid group, and finally we applied two methods to find regions under positive selection to possibly highlight instances of domestic adaptive introgression into a wild population. We found three hybrids within the Sardinian sample (3.1% out of the whole dataset). We reported 11 significant windows under positive selection with a method that looks for overly differentiated loci in the target population, compared with other two populations. We also identified 82 genomic regions with signs of selection in the domestic pig but not in the wild boar, two of which overlapped with genomic regions enriched for domestic alleles in the hybrid pool. Genes in these regions can be linked with reproductive success. Given our results, domestic introgression does not seem to be pervasive in the Sardinian wild boar. Nevertheless, we suggest monitoring the possible spread of advantageous domestic alleles in the coming years.

A Predator on the Doorstep: Kill Site Selection by a Lone Wolf in a Peri-Urban Park in a Mediterranean Area.

The aim of the study was to assess which kill site characteristics were selected by a lone wolf living in a protected Mediterranean coastal area near the city of Pisa, Italy, where both wild and domestic ungulates were available as potential prey. Between 2017 and 2019, we monitored the wolf's predatory behaviour through a combination of camera trapping and active search for kill sites and prey carcasses. The main prey found was the fallow deer (n = 82); only two wild boars and no domestic ungulates were found preyed upon. The features and habitat of kill sites were modelled to test for selection by the wolf. The habitat type of kill site was composed of meadows and pastures (89.3%), woods (7.3%), degraded coastal areas (1.9%), roads and rivers (1.1%), and marshes (0.5%). We calculated their distance from landscape features and ran a binomial generalised linear model to test the influence of such landscape variables. The distance of kill sites from landscape elements was significantly different from random control sites, and a positive selection for fences was found. In fact, the wolf pushed fallow deer towards a fence to constrain them and prevent them from escaping. We also analysed the body condition of predated fallow deer as a percentage of fat content in the bone marrow of the hind legs. Our results revealed the selection of the lone wolf for deer in good body condition. This is a possible outcome of the habitat selection shown by fallow deer in the study area, where fenced open pastures are the richest in trophic resources; therefore, our findings suggest a high efficacy for the lone wolf hunting strategy, but also the adoption of a high risk feeding strategy by deer. This study suggests that a lone predator can take advantage of human infrastructures to maximise its predatory effectiveness.

Main roads and land cover shaped the genetic structure of a Mediterranean island wild boar population.

Patterns of genetic differentiation within and among animal populations might vary due to the simple effect of distance or landscape features hindering gene flow. An assessment of how landscape connectivity affects gene flow can help guide management, especially in fragmented landscapes. Our objective was to analyze population genetic structure and landscape genetics of the native wild boar (Sus scrofa meridionalis) population inhabiting the island of Sardinia (Italy), and test for the existence of Isolation-by-Distance (IBD), Isolation-by-Barrier (IBB), and Isolation-by-Resistance (IBR). A total of 393 Sardinian wild boar samples were analyzed using a set of 16 microsatellite loci. Signals of genetic introgression from introduced non-native wild boars or from domestic pigs were revealed by a Bayesian cluster analysis including 250 reference individuals belonging to European wild populations and domestic breeds. After removal of introgressed individuals, genetic structure in the population was investigated by different statistical approaches, supporting a partition into five discrete subpopulations, corresponding to five geographic areas on the island: north-west (NW), central west (CW), south-west (SW), north-central east (NCE), and south-east (SE). To test the IBD, IBB, and IBR hypotheses, we optimized resistance surfaces using genetic algorithms and linear mixed-effects models with a maximum likelihood population effects parameterization. Landscape genetics analyses revealed that genetic discontinuities between subpopulations can be explained by landscape elements, suggesting that main roads, urban settings, and intensively cultivated areas are hampering gene flow (and thus individual movements) within the Sardinian wild boar population. Our results reveal how human-transformed landscapes can affect genetic connectivity even in a large-sized and highly mobile mammal such as the wild boar, and provide crucial information to manage the spread of pathogens, including the African Swine Fever virus, endemic in Sardinia.

Correction to: Genome-wide SNP data unveils the globalization of domesticated pigs.

An amendment to this paper has been published and can be accessed via the original article.

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 (qiw) 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.

Author Correction: Hotspots of recent hybridization between pigs and wild boars in Europe.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East-Caucasian and Central Italian populations.

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East-Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.

Hotspots of recent hybridization between pigs and wild boars in Europe.

After a strong demographic decline before World War II, wild boar populations are expanding and the species is now the second-most abundant ungulate in Europe. This increase raises concerns due to wild boar impact on crops and natural ecosystems and as potential vector of diseases. Additionally, wild boar can hybridize with domestic pigs, which could increase health risks and alter wild boar adaptive potential. We analysed 47,148 Single Nucleotide Polymorphisms in wild boar from Europe (292) and the Near East (16), and commercial (44) and local (255) pig breeds, to discern patterns of hybridization across Europe. We identified 33 wild boars with more than 10% domestic ancestry in their genome, mostly concentrated in Austria, Bosnia and Herzegovina, Bulgaria and Serbia. This difference is probably due to contrasting practices, with free-ranging vs. industrial farming but more samples would be needed to investigate larger geographic patterns. Our results suggest hybridization has occurred over a long period and is still ongoing, as we observed recent hybrids. Although wild and domestic populations have maintained their genetic distinctiveness, potential health threats raise concerns and require implementation of management actions and farming practices aimed at reducing contact between wild and domestic pigs.

Large-scale mitochondrial DNA analysis reveals new light on the phylogeography of Central and Eastern-European Brown hare (Lepus europaeus Pallas, 1778).

European brown hare, Lepus europaeus, from Central and Eastern European countries (Hungary, Poland, Serbia, Lithuania, Romania, Georgia and Italy) were sampled, and phylogenetic analyses were carried out on two datasets: 1.) 137 sequences (358 bp) of control region mtDNA; and 2.) 105 sequences of a concatenated fragment (916 bp), including the cytochrome b, tRNA-Thr, tRNA-Pro and control region mitochondrial DNA. Our sequences were aligned with additional brown hare sequences from GenBank. A total of 52 and 51 haplotypes were detected within the two datasets, respectively, and assigned to two previously described major lineages: Anatolian/Middle Eastern (AME) and European (EUR). Furthermore, the European lineage was divided into two subclades including South Eastern European (SEE) and Central European (CE). Sympatric distribution of the lineages of the brown hare in South-Eastern and Eastern Europe revealed contact zones there. BAPS analysis assigned sequences from L. europaeus to five genetic clusters, whereas CE individuals were assigned to only one cluster, and AME and SEE sequences were each assigned to two clusters. Our findings uncover numerous novel haplotypes of Anatolian/Middle Eastern brown hare outside their main range, as evidence for the combined influence of Late Pleistocene climatic fluctuations and anthropogenic activities in shaping the phylogeographic structure of the species. Our results support the hypothesis of a postglacial brown hare expansion from Anatolia and the Balkan Peninsula to Central and Eastern Europe, and suggest some slight introgression of individual haplotypes from L. timidus to L. europaeus.

Estimation of pack density in grey wolf (Canis lupus) by applying spatially explicit capture-recapture models to camera trap data supported by genetic monitoring.

BACKGROUND: Density estimation is a key issue in wildlife management but is particularly challenging and labour-intensive for elusive species. Recently developed approaches based on remotely collected data and capture-recapture models, though representing a valid alternative to more traditional methods, have found little application to species with limited morphological variation. We implemented a camera trap capture-recapture study to survey wolf packs in a 560-km2 area of Central Italy. Individual recognition of focal animals (alpha) in the packs was possible by relying on morphological and behavioural traits and was validated by non-invasive genotyping and inter-observer agreement tests. Two types (Bayesian and likelihood-based) of spatially explicit capture-recapture (SCR) models were fitted on wolf pack capture histories, thus obtaining an estimation of pack density in the area. RESULTS: In two sessions of camera trapping surveys (2014 and 2015), we detected a maximum of 12 wolf packs. A Bayesian model implementing a half-normal detection function without a trap-specific response provided the most robust result, corresponding to a density of 1.21 ± 0.27 packs/100 km2 in 2015. Average pack size varied from 3.40 (summer 2014, excluding pups and lone-transient wolves) to 4.17 (late winter-spring 2015, excluding lone-transient wolves). CONCLUSIONS: We applied for the first time a camera-based SCR approach in wolves, providing the first robust estimate of wolf pack density for an area of Italy. We showed that this method is applicable to wolves under the following conditions: i) the existence of sufficient phenotypic/behavioural variation and the recognition of focal individuals (i.e. alpha, verified by non-invasive genotyping); ii) the investigated area is sufficiently large to include a minimum number of packs (ideally 10); iii) a pilot study is carried out to pursue an adequate sampling design and to train operators on individual wolf recognition. We believe that replicating this approach in other areas can allow for an assessment of density variation across the wolf range and would provide a reliable reference parameter for ecological studies.

Genome-wide SNP data unveils the globalization of domesticated pigs.

BACKGROUND: Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection. RESULTS: A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection. CONCLUSIONS: Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

Long-lasting, kin-directed female interactions in a spatially structured wild boar social network.

Individuals can increase inclusive fitness benefits through a complex network of social interactions directed towards kin. Preferential relationships with relatives lead to the emergence of kin structures in the social system. Cohesive social groups of related individuals and female philopatry of wild boar create conditions for cooperation through kin selection and make the species a good biological model for studying kin structures. Yet, the role of kinship in shaping the social structure of wild boar populations is still poorly understood. In the present study, we investigated spatio-temporal patterns of associations and the social network structure of the wild boar Sus scrofa population in Bialowieza National Park, Poland, which offered a unique opportunity to understand wild boar social interactions away from anthropogenic factors. We used a combination of telemetry data and genetic information to examine the impact of kinship on network cohesion and the strength of social bonds. Relatedness and spatial proximity between individuals were positively related to the strength of social bond. Consequently, the social network was spatially and genetically structured with well-defined and cohesive social units. However, spatial proximity between individuals could not entirely explain the association patterns and network structure. Genuine, kin-targeted, and temporarily stable relationships of females extended beyond spatial proximity between individuals while males interactions were short-lived and not shaped by relatedness. The findings of this study confirm the matrilineal nature of wild boar social structure and show how social preferences of individuals translate into an emergent socio-genetic population structure.

Contemporary genetic structure, phylogeography and past demographic processes of wild boar Sus scrofa population in Central and Eastern Europe.

The wild boar (Sus scrofa) is one of the most widely distributed mammals in Europe. Its demography was affected by various events in the past and today populations are increasing throughout Europe. We examined genetic diversity, structure and population dynamics of wild boar in Central and Eastern Europe. MtDNA control region (664 bp) was sequenced in 254 wild boar from six countries (Poland, Hungary, Belarus, Ukraine, Moldova and the European part of Russia). We detected 16 haplotypes, all known from previous studies in Europe; 14 of them belonged to European 1 (E1) clade, including 13 haplotypes from E1-C and one from E1-A lineages. Two haplotypes belonged respectively to the East Asian and the Near Eastern clade. Both haplotypes were found in Russia and most probably originated from the documented translocations of wild boar. The studied populations showed moderate haplotype (0.714±0.023) and low nucleotide diversity (0.003±0.002). SAMOVA grouped the genetic structuring of Central and Eastern European wild boar into three subpopulations, comprising of: (1) north-eastern Belarus and the European part of Russia, (2) Poland, Ukraine, Moldova and most of Belarus, and (3) Hungary. The multimodal mismatch distribution, Fu's Fs index, Bayesian skyline plot and the high occurrence of shared haplotypes among populations did not suggest strong demographic fluctuations in wild boar numbers in the Holocene and pre-Holocene times. This study showed relatively weak genetic diversity and structure in Central and Eastern European wild boar populations and underlined gaps in our knowledge on the role of southern refugia and demographic processes shaping genetic diversity of wild boar in this part of Europe.

Domestication does not narrow MHC diversity in Sus scrofa.

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.

Mitochondrial CR-1 variation in Sardinian hares and its relationships with other Old World hares (Genus Lepus).

Among the European fauna, the Sardinian hare (Lepus sp.) is peculiar in that it differs from all other hares inhabiting the continent. Here, we report on the variation of a 461 bp sequence of hypervariable domain 1 of the mitochondrial control region, examined in 42 hares collected throughout Sardinia and compared to the corresponding sequences of different Lepus taxa. Seventeen novel haplotypes were found in the Sardinian population, resulting in a haplotype diversity of 0.840 and a nucleotide diversity of 0.012. As a result of Bayesian and principal coordinates analyses, Sardinian hares were grouped with North African hares, constituting a monophyletic clade that diverges from all other Old World hares, including Cape hares from South Africa and East Asia. Hence, our data agree that populations inhabiting North Africa and Sardinia form a distinct taxon, which could possibly be included in the L. capensis superspecies. Moreover, two corresponding lineages can be found in Sardinia and Tunisia, providing evidence of a common origin of the two populations and thus supporting the hypothesis that North African hares were introduced into the island in historical times. Our data show that the two lineages differ in their geographic distribution throughout the island and that the wild Sardinian population also shows the signature of a postintroduction demographic expansion.