Rehabilitation and release of orphaned Eurasian lynx (Lynx lynx) in Europe: Implications for management and conservation.

Rehabilitation of injured or immature individuals has become an increasingly used conservation and management tool. However, scientific evaluation of rehabilitations is rare, raising concern about post-release welfare as well as the cost-effectiveness of spending scarce financial resources. Over the past 20 years, events of juvenile Eurasian lynx presumably orphaned have been observed in many European lynx populations. To guide the management of orphaned lynx, we documented survival, rehabilitation and fate after the release and evaluated the potential relevance of lynx orphan rehabilitation for population management and conservation implications. Data on 320 orphaned lynx was collected from 1975 to 2022 from 13 countries and nine populations. The majority of orphaned lynx (55%) were taken to rehabilitation centres or other enclosures. A total of 66 orphans were released back to nature. The portion of rehabilitated lynx who survived at least one year after release was 0.66. Release location was the best predictor for their survival. Of the 66 released lynx, ten have reproduced at least once (8 females and 2 males). Conservation implications of rehabilitation programmes include managing genetic diversity in small, isolated populations and reintroducing species to historical habitats. The lynx is a perfect model species as most reintroduced populations in Central Europe show significantly lower observed heterozygosity than most of the autochthonous populations, indicating that reintroduction bottlenecks, isolation and post-release management have long-term consequences on the genetic composition of populations. The release of translocated orphans could be a valuable contribution to Eurasian lynx conservation in Europe. It is recommended to release orphans at the distribution edge or in the frame of reintroduction projects instead of a release in the core area of a population where it is not necessary from a demographic and genetic point of view. Rehabilitation programmes can have conservation implications that extend far beyond individual welfare benefits.

Wolf genetic diversity compared across Europe using the yardstick method.

Integrating data across studies with traditional microsatellite genetic markers requires careful calibration and represents an obstacle for investigation of wide-ranging species where populations require transboundary management. We used the "yardstick" method to compare results published across Europe since 2002 and new wolf (Canis lupus) genetic profiles from the Carpathian Mountains in Central Europe and the Dinaric Mountains in Southeastern Europe, with the latter as our reference population. We compared each population with Dinaric wolves, considering only shared markers (range 4-17). For each population, we calculated standard genetic diversity indices plus calibrated heterozygosity (Hec) and allelic richness (Ac). Hec and Ac in Dinaric (0.704 and 9.394) and Carpathian wolves (0.695 and 7.023) were comparable to those observed in other large and mid-sized European populations, but smaller than those of northeastern Europe. Major discrepancies in marker choices among some studies made comparisons more difficult. However, the yardstick method, including the new measures of Hec and Ac, provided a direct comparison of genetic diversity values among wolf populations and an intuitive interpretation of the results. The yardstick method thus permitted the integration of diverse sources of publicly available microsatellite data for spatiotemporal genetic monitoring of evolutionary potential.

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

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ 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.

Lead and Other Trace Element Levels in Brains of Croatian Large Terrestrial Carnivores: Influence of Biological and Ecological Factors.

Trace element pollution can adversely affect the brains of individuals and thus impact the entire population of apex predators, such as large European carnivores. We assessed exposure to prominent neurotoxicants As, Cd, Hg and Pb by measuring their brain stem levels in brown bears (n = 114), grey wolves (n = 8), Eurasian lynx (n = 3), and golden jackals (n = 2) sampled in 2015-2022 in Croatia. The highest of the non-essential elements was the Pb level in the bears' brains (median, Q1-Q3; 11.1, 7.13-24.1 µg/kg wet mass), with 4% of animals, all subadults, exceeding the established normal bovine levels (100 µg/kg wet mass). Species-specific differences were noted for Ca, Cd, Cu, Fe, Pb and Se brain levels. Female brown bears had higher As brain levels than males. Cubs and yearlings had lower brain Cd, but higher Zn, while subadults had higher Cu than adult bears. Hepatic As, Cd, Cu and Hg levels were shown to be a moderate proxy for estimating brain levels in bears (rS = 0.30-0.69). Multiple associations of As, Cd, Hg and Pb with essential elements pointed to a possible interaction and disturbance of brain Ca, Cu, Fe, Se and Zn homeostasis. Non-essential element levels in the brains of four studied species were lower than reported earlier for terrestrial meso-carnivores and humans. The age and sex of animals were highlighted as essential factors in interpreting brain element levels in ecotoxicological studies of large carnivores.

Population genetic structure of wolves in the northwestern Dinaric-Balkan region.

The Balkan Peninsula and the Dinaric Mountains possess extraordinary biodiversity and support one of the largest and most diverse wolf (Canis lupus) populations in Europe. Results obtained with diverse genetic markers show west-east substructure, also seen in various other species, despite the absence of obvious barriers to movement. However, the spatial extent of the genetic clusters remains unresolved, and our aim was to combine fine-scale sampling with population and spatial genetic analyses to improve resolution of wolf genetic clusters. We analyzed 16 autosomal microsatellites from 255 wolves sampled in Slovenia, Croatia, Bosnia and Herzegovina (BIH), and Serbia and documented three genetic clusters. These comprised (1) Slovenia and the regions of Gorski kotar and Lika in Croatia, (2) the region of Dalmatia in southern Croatia and BIH, and (3) Serbia. When we mapped the clusters geographically, we observed west-east genetic structure across the study area, together with some specific structure in BIH-Dalmatia. We observed that cluster 1 had a smaller effective population size, consistent with earlier reports of population recovery since the 1980s. Our results provide foundation for future genomic studies that would further resolve the observed west-east population structure and its evolutionary history in wolves and other taxa in the region and identify focal areas for habitat conservation. They also have immediate importance for conservation planning for the wolves in one of the most important parts of the species' European range.

Fear, foraging and olfaction: how mesopredators avoid costly interactions with apex predators.

Where direct killing is rare and niche overlap low, sympatric carnivores may appear to coexist without conflict. Interference interactions, harassment and injury from larger carnivores may still pose a risk to smaller mesopredators. Foraging theory suggests that animals should adjust their behaviour accordingly to optimise foraging efficiency and overall fitness, trading off harvest rate with costs to fitness. The foraging behaviour of red foxes, Vulpes vulpes, was studied with automated cameras and a repeated measures giving-up density (GUD) experiment where olfactory risk cues were manipulated. In Plitvice Lakes National Park, Croatia, red foxes increased GUDs by 34% and quitting harvest rates by 29% in response to wolf urine. In addition to leaving more food behind, foxes also responded to wolf urine by spending less time visiting food patches each day and altering their behaviour in order to compensate for the increased risk when foraging from patches. Thus, red foxes utilised olfaction to assess risk and experienced foraging costs due to the presence of a cue from gray wolves, Canis lupus. This study identifies behavioural mechanisms which may enable competing predators to coexist, and highlights the potential for additional ecosystem service pathways arising from the behaviour of large carnivores. Given the vulnerability of large carnivores to anthropogenic disturbance, a growing human population and intensifying resource consumption, it becomes increasingly important to understand ecological processes so that land can be managed appropriately.

Apex predatory mammals as bioindicator species in environmental monitoring of elements in Dinaric Alps (Croatia).

Tissue element investigations of apex terrestrial mammals are very scarce in Europe. We quantified 16 essential and nonessential elements in the kidney cortex, liver, and muscle tissue of 467 brown bears (Ursus arctos), 125 gray wolves (Canis lupus), one Eurasian lynx (Lynx lynx), and three golden jackals (Canis aureus) from Croatia by inductively coupled plasma mass spectrometry (ICP-MS). Renal cadmium (0.6% of animals) and lead (1%) and hepatic lead (5%) were found in toxicologically relevant levels for mammals only in bears, while the other elements were within normal range. The association of age, sex, season, and region with measured tissue elements in bear and wolf was estimated by multiple regression analyses. Age-related accumulation of cadmium was observed in bears and wolves. Lead tissue content increased with the age of bears but declined in wolves. Female bears and wolves had higher arsenic, iron, and thallium than males in some tissues. Also, cadmium, mercury, copper, zinc, selenium, molybdenum, and uranium were more abundant only in female bears. Male bears had higher potassium, zinc, and magnesium, while male wolves had higher calcium in some tissues compared to female wolves. Seasonal differences were mainly observed for bears' tissues and region-specific differences only in wolves. The bear kidneys had the highest levels of cobalt, copper, molybdenum, cadmium, and lead among the four studied species. The element levels reported for bears and wolves represent baseline values for the Dinaric population.

The prevalence and impact of Babesia canis and Theileria sp. in free-ranging grey wolf (Canis lupus) populations in Croatia.

BACKGROUND: Babesia spp. and Theileria spp. are important emerging causes of disease in dogs. Alongside these domesticated hosts, there is increasing recognition that these piroplasms can also be found in a range of wild animals with isolated reports describing the presence of these pathogen in foxes (Vulpes vulpes) and captive grey wolves (Canis lupus). The prevalence and impact of these infections in free-ranging populations of canids are unknown. To gain a better insight into the epidemiology and pathogenesis of piroplasm infections in free-ranging grey wolves, pathological and molecular investigations into captive and free-ranging grey wolves in Croatia were performed. RESULTS: The carcasses of 107 free-ranging wolves and one captive wolf were the subjects of post-mortem investigations and sampling for molecular studies. A blood sample from one live captured wolf for telemetric tracking was also used for molecular analysis. PCR amplification targeting the 18S RNA gene revealed that 21 of 108 free-ranging wolves and one captive animal were positive for Theileria/Babesia DNA. Subsequent sequencing of a fragment of the 18S RNA gene revealed that 7/22 animals were positive for Babesia canis while the other amplified sequence were found to be identical with corresponding 18S rDNA sequences of Theileria capreoli isolated from wild deer (15/22). Haematological and cytological analysis revealed the presence of signet-ring shaped or pear-shaped piroplasms in several animals with the overall parasite burden in all positive animals assessed to be very low. Pathological investigation of the captive animal revealed fatal septicemia as a likely outcome of hemolytic anaemia. There was little or no evidence of hemolytic disease consistent with babesiosis in other animals. CONCLUSION: Importantly, the presence of B. canis in free-ranging grey wolves has not been described before but has been reported in a single fox and domestic dogs only. That B. canis infections cause disease in dogs but have little impact on wolf health possibly suggests that the wolf is the natural and the domestic dog is a secondary host. Surprisingly, the frequent finding of Theileria capreoli in wolves suggests that this Theileria species is not restricted to ungulates (cervids) but commonly infects also this carnivore species. Nevertheless, the potential role that these asymptomatically infected animals may play in the dispersal of these pathogens to susceptible sympatric species such as domesticated dogs requires further investigation.

Protozoan and helminth parasite fauna of free-living Croatian wild wolves (Canis lupus) analyzed by scat collection.

The European wolf (Canis lupus) is a large carnivore species present in limited areas of Europe with several small populations still being considered as endangered. Wolves can be infected by a wide range of protozoan and metazoan parasites with some of them affecting free-living wolf health condition. On this account, an epidemiological survey was conducted to analyze the actual parasite fauna in Croatian wild wolves. In total, 400 individual faecal samples were collected during field studies on wolf ecology in the years 2002-2011. Parasite stages were identified by the sodium acetate acetic acid formalin (SAF)-technique, carbolfuchsin-stained faecal smears and Giardia/Cryptosporidium coproantigen-ELISAs. A subset of taeniid eggs-positive wolf samples was additionally analyzed by PCR and subsequent sequencing to identify eggs on Echinococcus granulosus/E. multilocularis species level. In total 18 taxa of parasites were here detected. Sarcocystis spp. (19.1%) occurred most frequently in faecal samples, being followed by Capillaria spp. (16%), ancylostomatids (13.1%), Crenosoma vulpis (4.6%), Angiostrongylus vasorum (3.1%), Toxocara canis (2.8%), Hammondia/Neospora spp. (2.6 %), Cystoisospora ohioensis (2.1%), Giardia spp. (2.1%), Cystoisospora canis (1.8%), Cryptosporidium spp. (1.8%), Trichuris vulpis (1.5%), Taenia spp. (1.5%), Diphyllobothrium latum (1.5%), Strongyloides spp. (0.5%), Opisthorchis felineus (0.5%), Toxascaris leonina (0.3%), Mesocestoides litteratus (0.3%) and Alaria alata (0.3%). Some of the here identified parasites represent relevant pathogens for wolves, circulating between these carnivorous definitive hosts and a variety of mammalian intermediate hosts, e. g. Taenia spp. and Sarcocystis spp., while others are considered exclusively pathogenic for canids (e.g. A. vasorum, C. vulpis, T. vulpis, Cystoisospora spp.). This study provides first records on the occurrence of the two relevant anthropozoonotic parasites, Giardia spp. and Cryptosporidium spp., in wild wolves from Croatia.

Border Security Fencing and Wildlife: The End of the Transboundary Paradigm in Eurasia?

The ongoing refugee crisis in Europe has seen many countries rush to construct border security fencing to divert or control the flow of people. This follows a trend of border fence construction across Eurasia during the post-9/11 era. This development has gone largely unnoticed by conservation biologists during an era in which, ironically, transboundary cooperation has emerged as a conservation paradigm. These fences represent a major threat to wildlife because they can cause mortality, obstruct access to seasonally important resources, and reduce effective population size. We summarise the extent of the issue and propose concrete mitigation measures.

Demographically-Based Evaluation of Genomic Regions under Selection in Domestic Dogs.

Controlling for background demographic effects is important for accurately identifying loci that have recently undergone positive selection. To date, the effects of demography have not yet been explicitly considered when identifying loci under selection during dog domestication. To investigate positive selection on the dog lineage early in the domestication, we examined patterns of polymorphism in six canid genomes that were previously used to infer a demographic model of dog domestication. Using an inferred demographic model, we computed false discovery rates (FDR) and identified 349 outlier regions consistent with positive selection at a low FDR. The signals in the top 100 regions were frequently centered on candidate genes related to brain function and behavior, including LHFPL3, CADM2, GRIK3, SH3GL2, MBP, PDE7B, NTAN1, and GLRA1. These regions contained significant enrichments in behavioral ontology categories. The 3rd top hit, CCRN4L, plays a major role in lipid metabolism, that is supported by additional metabolism related candidates revealed in our scan, including SCP2D1 and PDXC1. Comparing our method to an empirical outlier approach that does not directly account for demography, we found only modest overlaps between the two methods, with 60% of empirical outliers having no overlap with our demography-based outlier detection approach. Demography-aware approaches have lower-rates of false discovery. Our top candidates for selection, in addition to expanding the set of neurobehavioral candidate genes, include genes related to lipid metabolism, suggesting a dietary target of selection that was important during the period when proto-dogs hunted and fed alongside hunter-gatherers.

Dietary items as possible sources of (137)Cs in large carnivores in the Gorski Kotar forest ecosystem, Western Croatia.

The mountain forest ecosystem of Gorski Kotar is distant from any significant sources of environmental pollution, though recent findings have revealed that this region is among the most intense (137)Cs contaminated area in Croatia. Therefore, the aim of this study was to investigate (137)Cs and (40)K load in three large predator species in the mountain forest ecosystem. Radionuclides mass activities were determined by the gamma-spectrometric method in the muscle tissue of brown bear (47), wolf (7), lynx (1) and golden jackal (2). The highest (137)Cs mass activity was found in lynx (153 Bq kg(-1)), followed by brown bear (132 Bq kg(-1)), wolf (22.2 Bq kg(-1)), and golden jackal (2.48 Bq kg(-1)). Analysis of 63 samples of dietary items suggests that they are not all potentially dominant sources of (137)Cs for wildlife. The most important source of radionuclides for the higher parts of the food-chain from the study area were found to be the mushroom species wood hedgehog (Hydnum repandum), with a transfer factor TF of 5.166, and blueberry (Vaccinium myrtillus) as a plant species (TF=2.096). Food items of animal origin indicated higher mass activity of radionuclides and therefore are possible moderate bioindicators of environmental pollution. The results also revealed that possible unknown wild animal food sources are a caesium source in the study region, and further study is required to illuminate this issue.

Genome-wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore.

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north-south population differentiation. We investigated more than 67-K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric-Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome-wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long-term separation among population clusters.

Organochlorine pesticides and polychlorinated biphenyl congeners in wild terrestrial mammals from Croatia: Interspecies comparison of residue levels and compositions.

In this pilot study, we investigated levels of organochlorine pesticides (OCPs) and polychlorinated biphenyls (PCBs) in the adipose tissues of two free-ranging terrestrial carnivores from Croatia sampled in 2010 and 2011: the brown bear (Ursus arctos; N=32) and the grey wolf (Canis lupus; N=29). Concentrations of ∑OCPs and ΣPCBs ranged from 0.45 to 4.09 ng g(-1) lipid mass (lm) and from 0.93 to 8.52 ng g(-1) lm in brown bear, and from 1.18 to 5.67 ng g(-1) lm and 2.68 to 48.9 ng g(-1) lm in grey wolf adipose tissues, respectively. PCBs were dominant accounting for over 72% of total analyzed persistent organic pollutants (POPs). The sum of six indicator non-dioxin-like PCBs (Σ6IndNDL PCBs) made up 60-93% and 58-85% of the total congener concentrations in brown bears and wolves, respectively. Although the levels of the measured parameters were significantly higher in grey wolves than in bears, the contaminant profiles of the two species were similar, with γ-HCH, HCB, ß-HCH and DDE as major OCP contaminants, and PCB-153>PCB-180≈PCB-170>PCB-138 as the dominant congeners. The sum of hexachlorocyclohexanes (ΣHCHs) and 8 toxicologically relevant dioxin-like PCBs (Σ8ToxDL PCBs) was higher in the males than in the females of the brown bear. Concentrations of ΣDDTs, HCB, ΣOCP, ΣPCBs, Σ6IndNDL PCBs, and toxicologically relevant non-dioxin-like PCBs (ΣToxNDL PCBs) were significantly positively correlated with lipid content in the grey wolf. Concentrations of OCPs and PCBs in brown bears and wolves from Croatia were low and normal for large terrestrial mammals.

Recovery of large carnivores in Europe's modern human-dominated landscapes.

The conservation of large carnivores is a formidable challenge for biodiversity conservation. Using a data set on the past and current status of brown bears (Ursus arctos), Eurasian lynx (Lynx lynx), gray wolves (Canis lupus), and wolverines (Gulo gulo) in European countries, we show that roughly one-third of mainland Europe hosts at least one large carnivore species, with stable or increasing abundance in most cases in 21st-century records. The reasons for this overall conservation success include protective legislation, supportive public opinion, and a variety of practices making coexistence between large carnivores and people possible. The European situation reveals that large carnivores and people can share the same landscape.

Cadmium, lead, and mercury exposure assessment among croatian consumers of free-living game.

Free-living game can be an important source of dietary cadmium and lead; the question is whether exposure to these two elements is such that it might cause adverse health effects in the consumers. The aim of this study was to estimate dietary exposure to cadmium, lead, and mercury from free-living big game (fallow deer, roe deer, red deer, wild boar, and brown bear), and to mercury from small game (pheasant and hare), hunted in Croatia from 1990 to 2012. The exposure assessment was based on available literature data and our own measurements of metal levels in the tissues of the game, by taking into account different consumption frequencies (four times a year, once a month and once a week). Exposure was expressed as percentage of (provisional) tolerable weekly intake [(P)TWI] values set by the European Food Safety Authority (EFSA). Consumption of game meat (0.002-0.5 % PTWI) and liver (0.005-6 % PTWI) assumed for the general population (four times a year) does not pose a health risk to consumers from the general population, nor does monthly (0.02-6 % PTWI) and weekly (0.1-24 % PTWI) consumption of game meat. However, because of the high percentage of free-living game liver and kidney samples exceeding the legislative limits for cadmium (2-99 %) and lead (1-82 %), people should keep the consumption of certain game species' offal as low as possible. Children and pregnant and lactating women should avoid eating game offal altogether. Free-living game liver could be an important source of cadmium if consumed on a monthly basis (3-74 % TWI), and if consumed weekly (11-297 % TWI), it could even give rise to toxicological concern.

Selenium in brown bears (Ursus arctos) from Croatia: Relation to cadmium and mercury.

Muscle (n = 111), liver (n = 111), and kidney cortex (n = 101) samples from brown bears (Ursus arctos) were collected in the 2009 and 2010 hunting seasons in Croatia and analysed for selenium (Se), cadmium (Cd), and total mercury (Hg). The aim was to assess the levels of these elements according to age, sex, and season of collection, and to investigate possible Se/Cd and Se/Hg interactions. Median Se concentrations were 0.139 µg/g in muscle, 0.409 µg/g in liver and 1.75 µg/g wet mass in kidney cortex. Median Cd and Hg were 0.0078 and 0.0018 µg/g in muscle, 1.09 and 0.031 µg/g in liver, and 16.5 and 0.206 µg/g wet mass in the renal cortex, respectively. Se/Cd molar ratios were less than 1 in the kidney cortex, and close to or above 1 in liver and muscle, respectively. Toxic Cd and Hg correlated with Se in all of the studied tissues. Sex differences were found for all three elements (except Se in liver), with females having higher tissue concentration than males. Only Cd showed age-dependence. Bear samples collected in fall had higher Se in muscles, and Hg in muscles and liver compared to samples collected in spring. Element concentrations in brown bear tissues were within the range of previously reported studies. Bear meat is considered a rich source of Se, safe for consumption with regard to its Cd and Hg content. According to the molar ratio and correlation results, we assume that Se binding is not the primary detoxification pathway for Cd and Hg in brown bears.

Multilocus detection of wolf x dog hybridization in italy, and guidelines for marker selection.

Hybridization and introgression can impact the evolution of natural populations. Several wild canid species hybridize in nature, sometimes originating new taxa. However, hybridization with free-ranging dogs is threatening the genetic integrity of grey wolf populations (Canis lupus), or even the survival of endangered species (e.g., the Ethiopian wolf C. simensis). Efficient molecular tools to assess hybridization rates are essential in wolf conservation strategies. We evaluated the power of biparental and uniparental markers (39 autosomal and 4 Y-linked microsatellites, a melanistic deletion at the ß-defensin CBD103 gene, the hypervariable domain of the mtDNA control-region) to identify the multilocus admixture patterns in wolf x dog hybrids. We used empirical data from 2 hybrid groups with different histories: 30 presumptive natural hybrids from Italy and 73 Czechoslovakian wolfdogs of known hybrid origin, as well as simulated data. We assessed the efficiency of various marker combinations and reference samples in admixture analyses using 69 dogs of different breeds and 99 wolves from Italy, Balkans and Carpathian Mountains. Results confirmed the occurrence of hybrids in Italy, some of them showing anomalous phenotypic traits and exogenous mtDNA or Y-chromosome introgression. Hybridization was mostly attributable to village dogs and not strictly patrilineal. The melanistic ß-defensin deletion was found only in Italian dogs and in putative hybrids. The 24 most divergent microsatellites (largest wolf-dog FST values) were equally or more informative than the entire panel of 39 loci. A smaller panel of 12 microsatellites increased risks to identify false admixed individuals. The frequency of F1 and F2 was lower than backcrosses or introgressed individuals, suggesting hybridization already occurred some generations in the past, during early phases of wolf expansion from their historical core areas. Empirical and simulated data indicated the identification of the past generation backcrosses is always uncertain, and a larger number of ancestry-informative markers is needed.

Genome sequencing highlights the dynamic early history of dogs.

To identify genetic changes underlying dog domestication and reconstruct their early evolutionary history, we generated high-quality genome sequences from three gray wolves, one from each of the three putative centers of dog domestication, two basal dog lineages (Basenji and Dingo) and a golden jackal as an outgroup. Analysis of these sequences supports a demographic model in which dogs and wolves diverged through a dynamic process involving population bottlenecks in both lineages and post-divergence gene flow. In dogs, the domestication bottleneck involved at least a 16-fold reduction in population size, a much more severe bottleneck than estimated previously. A sharp bottleneck in wolves occurred soon after their divergence from dogs, implying that the pool of diversity from which dogs arose was substantially larger than represented by modern wolf populations. We narrow the plausible range for the date of initial dog domestication to an interval spanning 11-16 thousand years ago, predating the rise of agriculture. In light of this finding, we expand upon previous work regarding the increase in copy number of the amylase gene (AMY2B) in dogs, which is believed to have aided digestion of starch in agricultural refuse. We find standing variation for amylase copy number variation in wolves and little or no copy number increase in the Dingo and Husky lineages. In conjunction with the estimated timing of dog origins, these results provide additional support to archaeological finds, suggesting the earliest dogs arose alongside hunter-gathers rather than agriculturists. Regarding the geographic origin of dogs, we find that, surprisingly, none of the extant wolf lineages from putative domestication centers is more closely related to dogs, and, instead, the sampled wolves form a sister monophyletic clade. This result, in combination with dog-wolf admixture during the process of domestication, suggests that a re-evaluation of past hypotheses regarding dog origins is necessary.

North-South differentiation and a region of high diversity in European wolves (Canis lupus).

European wolves (Canis lupus) show population genetic structure in the absence of geographic barriers, and across relatively short distances for this highly mobile species. Additional information on the location of and divergence between population clusters is required, particularly because wolves are currently recolonizing parts of Europe. We evaluated genetic structure in 177 wolves from 11 countries using over 67K single nucleotide polymorphism (SNP) loci. The results supported previous findings of an isolated Italian population with lower genetic diversity than that observed across other areas of Europe. Wolves from the remaining countries were primarily structured in a north-south axis, with Croatia, Bulgaria, and Greece (Dinaric-Balkan) differentiated from northcentral wolves that included individuals from Finland, Latvia, Belarus, Poland and Russia. Carpathian Mountain wolves in central Europe had genotypes intermediate between those identified in northcentral Europe and the Dinaric-Balkan cluster. Overall, individual genotypes from northcentral Europe suggested high levels of admixture. We observed high diversity within Belarus, with wolves from western and northern Belarus representing the two most differentiated groups within northcentral Europe. Our results support the presence of at least three major clusters (Italy, Carpathians, Dinaric-Balkan) in southern and central Europe. Individuals from Croatia also appeared differentiated from wolves in Greece and Bulgaria. Expansion from glacial refugia, adaptation to local environments, and human-related factors such as landscape fragmentation and frequent killing of wolves in some areas may have contributed to the observed patterns. Our findings can help inform conservation management of these apex predators and the ecosystems of which they are part.