Search

Cystic echinococcosis (CE) is considered the most severe parasitic disease that ever affected the human population in Iceland. Before the start of eradication campaign in the 1860s, Iceland was a country with very high prevalence of human CE, with approximately every fifth person infected. Eradication of CE from Iceland by 1979 was a huge success story and served as a leading example for other countries on how to combat such a severe One Health problem. However, there is no genetic information on Echinococcus parasites before eradication. Here, we reveal the genetic identity for one of the last Echinococcus isolates in Iceland, obtained from a sheep 46 years ago (1977). We sequenced a large portion of the mitochondrial genome (8141 bp) and identified the isolate as Echinococcus granulosus sensu stricto genotype G1. As G1 is known to be highly infective genotype to humans, it may partly explain why such a large proportion of human population in Iceland was infected at a time . The study demonstrates that decades-old samples hold significant potential to uncover genetic identities of parasites in the past.

Echinococcus granulosus sensu lato is a group of tapeworm species known to cause cystic echinococcosis. Within this group, the Echinococcus canadensis cluster includes genotypes G8 and G10 that have a predominantly sylvatic life cycle – transmission occurs between wild cervids and wolves. Relatively few studies have explored the genetic variation of the elusive G8 and G10, and their extent of genetic variation is yet to be investigated at the complete mitochondrial (mt) genome level. The aim was to explore the genetic variation of these 2 genotypes in Europe using complete mtDNA sequences and provide a high-quality reference dataset for future studies. Sequences of complete mt genomes were produced for 29 samples of genotype G8 and G10 from wolves, moose, reindeer and roe deer, originating from Finland, Sweden, Russia, Poland, Latvia and Estonia. Genetic variation was explored based on phylogenetic network analysis, revealing marked differences between G8 and G10 (over 400 mutations), and more detailed patterns of variability within the 2 genotypes than previously observed. Understanding the mt genetic composition of a species provides a baseline for future studies aiming to understand whether this mt distinctiveness is mirrored in the nuclear genome and whether it has any impact on any phenotypic traits or parasite transmission.

Biotic interactions are expected to influence species’ responses to climate change, but they are usually not included when predicting future range shifts. We assessed the importance of biotic interactions to understand future consequences of climate and land use change for biodiversity using as a model system the brown bear (Ursus arctos) in Europe. By including biotic interactions using the spatial variation of energy contribution and habitat models of each food species, we showed that the use of biotic factors considerably improves our understanding of the distribution of brown bears. Predicted future range shifts, which included changes in the distribution of food species, varied greatly when considering various scenarios of change in biotic factors, warning about future indirect climate change effects. Our study confirmed that advancing our understanding of ecological networks of species interactions will improve future scenarios of biodiversity change, which is key for conserving biodiversity and ecosystem services.

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

Countries in the northern Baltic Sea region have been historically linked for thousands of years, and animal husbandry is one of the many information sources that enables the study of historical connections. Although goat husbandry in this part of Europe has been evidenced by scant archaeological materials, zooarchaeological and historical evidence has revealed its continuity, at least since the Late Iron Age. To explore the historical relationship between goat lineages and investigate affinities between the past and present-day populations in the Baltic Sea region, we analysed a 476-bp fragment of the mitochondrial DNA control region in 14 ancient goats from Finland and Estonia and 10 extant goats from Finland. The results revealed high mitochondrial diversity among the ancient goats. Two maternal lineages were shared between the Late Iron Age and medieval individuals from Finland and Estonia. Moreover, ancient Finnish and Estonian goats showed maternal affinity to extant Finngoat and Swedish Landrace breeds. Overall, the analysis of maternal goat lineages confirmed tight historical connections in the region.

In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (Ne) was estimated to be 199 individuals, and total combined population size (NC) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe's low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans.

The concept of One Health emphasizes the interdependence of human, animal and environmental health and is of growing significance, in part owing to the problems related to emerging infectious diseases of wildlife origin. Wild mammalian predators are a potential risk factor for transmission of zoonotic pathogens to domesticated animals and humans. This is especially relevant in rural areas, where transmission of zoonotic pathogens can occur particularly efficiently when free-ranging dogs are present. The main aim of this study was to determine helminth infections among wild mammalian predators and evaluate the overlap between helminth faunas of wild mammals and dogs. Scat samples of predators were collected in coastal areas of Western Estonia and genetic methodology applied for the correct identification of predator species from their scat. Parasitic helminths of mammalian predators in the scat samples were analysed and compared with dog data from a previous study. High helminth prevalence (~90%) was found in dominant predator species in the area, namely the red fox (Vulpes vulpes) and golden jackal (Canis aureus). Moreover, the helminth fauna of both wild species, including potentially zoonotic helminths, overlapped largely with that of rural dogs in the same area. The results, together with the ones from earlier parasitological studies among humans in Estonia, emphasize the potential risk of pathogen transmission from wild mammalian predators to dogs and from dogs to humans, making parasitic diseases of wildlife a One Health concern.

Cystic echinococcosis (CE), a zoonotic disease caused by the species complex of Echinococcus granulosus sensu lato (s.l.), is endemic in Middle East and the Mediterranean basin, where pastoral activity is widespread. Despite the chronic endemicity of the disease in Lebanon and neighbouring countries, recent data are scant. The objectives of this survey were to evaluate the current epidemiology of CE in Lebanon, investigate the prevalence and determine the infecting genotypes in locally raised sheep and goats. A multidimensional approach combining post-mortem inspection of slaughtered animals and molecular diagnosis of the parasite was conducted to this end. From 2018 to 2020, 62.9% of the sheep and 20.9% of the goats were found positive for CE. The presence of hydatids varied between organs, showing higher prevalence in the liver of sheep vs the lungs of goats, however, a higher fertility rate of hydatid cyst was observed in lungs for both animals. Molecular diagnosis based on partial mitochondrial cox1 gene (795 bp) showed that the majority of isolates were identified as E. granulosus sensu stricto (98.7%) and only one isolate from goat was identified as Echinococcus canadensis (genotype G7; 1.3%). Echinococcus granulosus s.s. population among the sheep and goats was represented by 22 haplotypes having very little genetic differentiation and relatively moderate haplotype diversity. Population demographics explored through neutrality indices suggested expanding population within the intermediate hosts. These results document the high prevalence of CE in the livestock of Lebanon and reveal for the first time the presence of three different genotypes G1, G3 and G7.

Regardless of the highly efficient anthelmintics available and the control measures taken by dog owners and veterinarians, gastrointestinal parasites, especially zoonotic helminths, are still abundant in dogs and pose a health risk to humans. Free-ranging dogs in rural areas can be an important source of helminth infection. The aims of the present work were to collect scats of rural dogs, determine the environmental contamination caused by helminth infections among rural dogs of Western Estonia, analyse how diet affects helminth infection rate and compare the findings to a previous study focusing on dog helminths in urban areas of Estonia. To differentiate the scats of dogs from other sympatric canids, a genetic method was applied. Of 328 samples, genetic analysis identified 84 scats belonging to dogs, of which 87.0% were infected with helminths. A high proportion of rural dog scats harboured eggs of Taeniidae (65.5%), followed by Trichuris spp./Eucoleus spp. (15.5%), Uncinaria stenocephala (14.7%) and Toxocara canis (4.3%). Coinfections occurred in 34.5% of the samples, being the most common between Taeniidae and U. stenocephala (41.4%). The intensity model indicated higher helminth infection rate in rural dogs preying on rodents and game. In comparison to urban dogs, rural dogs were nine times more likely to be infected with intestinal parasites. These results emphasize the need to implement measures to reduce helminth infections in dogs living in rural areas of Western Estonia. Among a complex of measures to be taken, we suggest that it is also important to diagnose which gastrointestinal parasite species infect dogs to determine specific anthelmintic treatment against these parasites.

Cystic echinococcosis is a socioeconomically important parasitic disease caused by the larval stage of the canid tapeworm Echinococcus granulosus, afflicting millions of humans and animals worldwide. The development of a vaccine (called EG95) has been the most notable translational advance in the fight against this disease in animals. However, almost nothing is known about the genomic organisation/location of the family of genes encoding EG95 and related molecules, the extent of their conservation or their functions. The lack of a complete reference genome for E. granulosus genotype G1 has been a major obstacle to addressing these areas. Here, we assembled a chromosomal-scale genome for this genotype by scaffolding to a high quality genome for the congener E. multilocularis, localised Eg95 gene family members in this genome, and evaluated the conservation of the EG95 vaccine molecule. These results have marked implications for future explorations of aspects such as developmentally-regulated gene transcription/expression (using replicate samples) for all E. granulosus stages; structural and functional roles of non-coding genome regions; molecular ‘cross-talk’ between oncosphere and the immune system; and defining the precise function(s) of EG95. Applied aspects should include developing improved tools for the diagnosis and chemotherapy of cystic echinococcosis of humans.

Aim: The Expansion-Contraction model has been used to explain the responses of species to climatic changes. During periods of unfavourable climatic conditions, species retreat to refugia from where they may later expand. This paper focuses on the palaeoecology of red deer over the past 54 ka across Europe and the Urals, to reveal patterns of change in their range and explore the role of environmental conditions in determining their distribution. Location: Europe and western Asia to 63°E. Taxon: Red deer (Cervus elaphus). Methods: We collected 984 records of radiocarbon-dated red deer subfossils from the Late Pleistocene and the Holocene, including 93 original dates. For each deer sample we compiled climatic and biome type data for the corresponding time intervals. Results: During the last 54 ka changes in red deer range in Europe and the Urals were asynchronous and differed between western and eastern Europe and western Asia due to different environmental conditions in those regions. The range of suitable areas for deer during the Last Glacial Maximum (LGM) was larger than previously thought and covered vast regions not only in southern but also in western and eastern Europe. Throughout the period investigated the majority of specimens inhabited forests in the temperate climatic zone. The contribution of forests in deer localities significantly decreased during the last 4 ka, due to deforestation of Europe caused by humans. Mean January temperature was the main limiting factor for species distribution. Over 90% of the samples were found in areas where mean January temperature was above −10°C. Main conclusions: Red deer response to climatic oscillations are in agreement with the Expansion-Contraction model but in contradiction to the statement of only the southernmost LGM refugia of the species. During the last 54 ka red deer occurred mostly in forests of the temperate climatic zone. © 2020 John Wiley & Sons Ltd. European Social Fund, Grant/Award Number: UDA-POKL.04.01.01-00-072/09-00; University of Wroclaw, Grant/Award Number: 0410/2990/18; Institute of Environmental Biology, University of Wrocław, Grant/Award Number: 0410/2990/18; Mammal Research Institute Polish Academy of Sciences; Narodowe Centrum Nauki , Grant/Award Number: DEC-2013/11/B/NZ8/00888 and UMO-2016/23/B/HS3/00387; Romanian National Authority for Scientific Research, UEFISCDI, Grant/Award Number: PN-IIIP4-ID-PCE-2016-0676; National Centre for Atmospheric Science and the Centre for Environmental Data Analysis, UK; Faunal Database of the Stage Three Project; Leverhulme Trust, Grant/Award Number: F00568W.

BackgroundReduction of genetic diversity can lead to reduced fitness of species, such as the loss of adaptability to changing environments. The native Eurasian beaver (Castor fiber) was hunted to extinction from Finland and many other countries in Europe in the nineteenth century. In Finland, the species was re-introduced in the 1930s with only a few individuals from Norway. Re-introductions were performed also in other countries of northern Europe and as a result, Eurasian beaver populations have undergone population bottlenecks leading to low levels of genetic diversity.Materials and MethodsHere, 200 Eurasian beaver samples from Finland, Estonia, Lithuania, and Russian Karelia were investigated using 12 microsatellite markers to examine the level of genetic diversity and relationship between the populations.ResultsWhile Russian and Estonian populations were genetically the closest, the Finnish population was clearly distinct from all others and had the lowest genetic variability among the study populations. This may be deleterious to the population especially in a changing environment.ConclusionsGenetic rescue could be the best solution to increase the genetic diversity and improve the future prospects of the population, although more studies are required to resolve the optimal source population.

The brown bear (Ursus arctos) is an iconic carnivoran species of the Northern Hemisphere. Its population history has been studied extensively using mitochondrial markers, which demonstrated signatures of multiple waves of migration, arguably connected with glaciation periods. Among Eurasian brown bears, Siberian populations remain understudied. We have sequenced complete mitochondrial genomes of four ancient (~4.5–40 kya) bears from South Siberia and 19 modern bears from South Siberia and the Russian Far East. Reconstruction of phylogenetic relationships between haplotypes and evaluation of modern population structure have demonstrated that all the studied samples belong to the most widespread Eurasian clade 3. One of the ancient haplotypes takes a basal position relative to the whole of clade 3; the second is basal to the haplogroup 3a (the most common subclade), and two others belong to clades 3a1 and 3b. Modern Siberian bears retain at least some of this diversity; apart from the most common haplogroup 3a, we demonstrate the presence of clade 3b, which was previously found mainly in mainland Eurasia and Northern Japan. Our findings highlight the importance of South Siberia as a refugium for northern Eurasian brown bears and further corroborate the hypothesis of several waves of migration in the Pleistocene.

The present phylogeographic pattern of red deer in Eurasia is not only a result of the contraction of their distribution range into glacial refugia and postglacial expansion, but probably also an effect of replacement of some red deer s.l. mtDNA lineages by others during the last 50 000 years. To better recognize this process, we analysed 501 sequences of mtDNA cytochrome b, including 194 ancient and 75 contemporary samples newly obtained for this study. The inclusion of 161 radiocarbon-dated samples enabled us to study the phylogeny in a temporal context and conduct divergence-time estimation and molecular dating. Depending on methodology, our estimate of divergence between Cervus elaphus and Cervus canadensis varied considerably (370 000 or 1.37 million years BP, respectively). The divergence times of genetic lineages and haplogroups corresponded to large environmental changes associated with stadials and interstadials of the Late Pleistocene. Due to the climatic oscillations, the distribution of C. elaphus and C. canadensis fluctuated in north–south and east–west directions. Some haplotypes dated to pre-Last Glacial Maximum periods were not detected afterwards, representing possibly extinct populations. We indicated with a high probability the presence of red deer sensu lato in south-eastern Europe and western Asia during the Last Glacial Maximum.

SummaryThe Sicilian wolf represented the only population of wolves living on a Mediterranean island until the first half of the twentieth century (1930s-1960s) 1–7. Previous studies hypothesised that they remained isolated from mainland wolves from the end of the Last Glacial Maximum (LGM) 8,9, until human persecutions led them to extinction 1–7.There are only seven known Sicilian wolf specimens from the 19th and 20th century preserved in museums in Italy and recent morphometric analyses assigned them to the new subspecies Canis lupus cristaldii10. To better understand the origins of the Sicilian wolf, and its relationship to other wolf populations, we sequenced four whole genomes (3.8×-11.6×) and five mitogenomes. We investigated the relationship between Sicilian wolves and other modern breeds to identify potential admixture. Furthermore, considering that the last land-bridge between Sicily and Italy disappeared after the LGM 11, around 17 kya, we explored the possibility that the Sicilian wolf retained ancestry from ancient wolf and dog lineages. Additionally, we explored whether the long-term isolation might have affected the genomic diversity, inbreeding levels and genetic load of the Sicilian wolf.Our findings show that the Sicilian wolves shared most ancestry with the modern Italian wolf population but are better modelled as admixed with European dog breeds, and shared traces of Eneolithic and Bronze age European dogs. We also find signatures of severe inbreeding and low genomic diversity at population and individual levels due to long-term isolation and drift, suggesting also low effective population size.

Review The neglected zoonosis cystic echinococcosis affects mainly pastoral and rural communities in both low-income and upper-middle-income countries. In Europe, it should be regarded as an orphan and rare disease. Although human cystic echinococcosis is a notifiable parasitic infectious disease in most European countries, in practice it is largely under-reported by national health systems. To fill this gap, we extracted data on the number, incidence, and trend of human cases in Europe through a systematic review approach, using both the scientific and grey literature and accounting for the period of publication from 1997 to 2021. The highest number of possible human cases at the national level was calculated from various data sources to generate a descriptive model of human cystic echinococcosis in Europe. We identified 64 745 human cystic echinococcosis cases from 40 European countries. The mean annual incidence from 1997 to 2020 throughout Europe was 0·64 cases per 100 000 people and in EU member states was 0·50 cases per 100 000 people. Based on incidence rates and trends detected in this study, the current epicentre of cystic echinococcosis in Europe is in the southeastern European countries, whereas historical endemic European Mediterranean countries have recorded a decrease in the number of cases over the time. Key messages: - This study aims to shed light on the unrecognised incidence of cystic echinococcosis in Europe, unveiling its epidemiological effect by providing a quantitative measure of number, incidence, and trends of human cases documented within the period 1997–2021; - Since human cystic echinococcosis cases are generally under-reported and data have uncertainty (partly due to misdiagnosis), data provided in this study should be considered as a conservative estimate of the real impact of this zoonotic infection historically occurring in Europe; - For the years 2017–19, we identified a total number of cystic echinococcosis cases four-fold higher than for The European Surveillance System (TESSy) data; - Decreasing trends have been recorded in most southern Mediterranean and some eastern European countries, where cystic echinococcosis has traditionally been highly prevalent; - Increasing trends have been identified in some eastern and southeastern European countries but, unexpectedly, also in most non-endemic countries of northern and western Europe; - Based on incidence and trends from 2017–19, the current epicentre of cystic echinococcosis in Europe is represented by the Balkan Peninsula; - Cystic echinococcosis in Europe remains a relevant public health issue and findings from this study should be used to support the planning of surveillance and control programmes in Europe according to the WHO 2021–2030 roadmap for neglected tropical diseases. This research was funded by the MEmE project from the EU’s Horizon 2020 Research and Innovation programme under grant agreement number 773830: One Health European Joint Programme. This work was also partially supported by research funding (grant PRG1209) from the Estonian Ministry of Education and Research info:eu-repo/semantics/publishedVersion

International audience; The cestode Echinococcus multilocularis is the causative agent of alveolar echinococcosis, a severe helminthic zoonotic disease distributed in the Northern Hemisphere. The lifecycle of the parasite is mainly sylvatic, involving canid and rodent hosts. The absence of genetic data from most eastern European countries is a major knowledge gap, affecting the study of associations with parasite populations in Western Europe. In this study, EmsB microsatellite genotyping of E. multilocularis was performed to describe the genetic diversity and relatedness of 785 E. multilocularis isolates from four western and nine eastern European countries, as well as from Armenia and the Asian parts of Russia and Turkey. The presence of the same E. multilocularis populations in the Benelux resulting from expansion from the historical Alpine focus can be deduced from the main profiles shared between these countries. All 33 EmsB profiles obtained from 528 samples from the nine eastern European countries belonged to the European clade, except one Asian profile form Ryazan Oblast, Russia. The expansion of E. multilocularis seems to have progressed from the historical Alpine focus through Hungary, Slovakia, the Czech Republic and southern Poland towards Latvia and Estonia. Most of the samples from Asia belong to the Asian clade, with one EmsB profile shared between Armenia and Turkey, and two between Turkey and Russia. However, two European profiles were described from two foxes in Turkey, including one harboring worms from both European and Asian clades. Three EmsB profiles from three Russian samples were associated with the Arctic clade. Two E. multilocularis profiles from rodents from Lake Baikal belonged to the Mongolian clade, described for the first time here using EmsB. Further worldwide studies on the genetic diversity of E. multilocularis using both mitochondrial sequencing and EmsB genotyping are needed to understand the distribution and expansion of the various clades.

By the beginning of the twentieth century, many brown bear populations in Europe were on the brink of extinction due to relentless hunting pressure and habitat loss. The situation was critical also in Estonia, where in the 1920s the population went through a severe demographic bottleneck. Thanks to the protective measures implemented in the 1930s, the population started to recover. However, the process has been slow, especially in the western and southern areas of the country. To study the effects of the demographic bottleneck, we analysed 216 brown bear samples from throughout their main range in Estonia. In combination with widely used methods of population genetics, a recently developed spatially explicit analysis (distribution of residual dissimilarity, DResD) was also applied. Three genetic clusters were revealed, of which two were most likely founded by the survivors of the bottleneck. The DResD analysis revealed several contact zones near the Estonian-Russian border, suggesting that the third cluster was influenced by gene flow from the neighbouring population in Russia. The DResD analysis revealed also a male biased movement corridor along the forested south-north axis in the central part of Estonia. In comparison to other European populations, the genetic diversity of the Estonian population is relatively low and is comparable with other populations that have gone through a severe bottleneck. This work has important implications for brown bear conservation and highlights once again the dangers of excessive hunting.

Background Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method. Methods We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow. Results We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441–445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat. Conclusions We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such ‘new’ repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

Domestication of horses fundamentally transformed long-range mobility and warfare. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling at Botai, Central Asia around 3500 bc3. Other longstanding candidate regions for horse domestication, such as Iberia and Anatolia, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 bc, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 bc driving the spread of Indo-European languages. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium bc Sintashta culture., The work by G. Boeskorov is done on state assignment of DPMGI SB RAS. This project was supported by the University Paul Sabatier IDEX Chaire d’Excellence (OURASI); Villum Funden miGENEPI research programme; the CNRS ‘Programme de Recherche Conjoint’ (PRC); the CNRS International Research Project (IRP AMADEUS); the France Génomique Appel à Grand Projet (ANR-10-INBS-09-08, BUCEPHALE project); IB10131 and IB18060, both funded by Junta de Extremadura (Spain) and European Regional Development Fund; Czech Academy of Sciences (RVO:67985912); the Zoological Institute ZIN RAS (АААА-А19-119032590102-7); and King Saud University Researchers Supporting Project (NSRSP–2020/2). The research was carried out with the financial support of the Russian Foundation for Basic Research (19-59-15001 and 20-04-00213), the Russian Science Foundation (16-18-10265, 20-78-10151, and 21-18-00457), the Government of the Russian Federation (FENU-2020-0021), the Estonian Research Council (PRG29), the Estonian Ministry of Education and Research (PRG1209), the Hungarian Scientific Research Fund (Project NF 104792), the Hungarian Academy of Sciences (Momentum Mobility Research Project of the Institute of Archaeology, Research Centre for the Humanities); and the Polish National Science Centre (2013/11/B/HS3/03822). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie (grant agreement 797449). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreements 681605, 716732 and 834616).

Pakistan is a neglected endemic focus for Echinococcus granulosus sensu lato, a zoonotic parasite species complex with the ability to infect wide spectrum of hosts. Wide gaps exist in literature for etiological agents of cystic echinococcosis (CE) in Pakistan due to a very low number of studies on identifying the exact genotypes involved in epidemiological manifestation of this disease. Focusing on transmission patterns and epidemiological dynamics, this study aimed at investigating infective genotypes among the cattle population of south Punjab, Pakistan, employing a mitochondrial marker nad5 (680 bp). Nucleotide sequences retrieved from 28 hydatid cyst isolates displayed considerable intraspecific variation revealing the existence of G3 and G1 strains of Echinococcus granulosus sensu stricto. The G3 genotype emerged as the predominant cause (78.57%) of hydatidosis in cattle. Apart from this, to understand phylogeographical relations, homologous nucleotide sequences of the partial nad5 gene from six major regions of the world were employed in the population genetics analysis to have an insight into genetic variability and demographics of G3 genotype in particular. Diversification of G3 and its haplotypes in Pakistan (n = 11) and other regions of the world (India, Iran, Turkey, Italy and France) was demonstrated. It was further demonstrated that the South Asian population (Pakistan and India) was highly differentiated from the other regions. It could, therefore, be speculated that G3 is diverging and expanding its population with South Asia as the main focal point.

Cats are important hosts for different zoonotic parasites that can be hazardous to human health. To date, few studies have attempted to identify the factors affecting parasitic infections in shelter animals. This study aims to analyse the presence of endoparasites in shelter cats in Tartu, Estonia, and identify factors affecting endoparasite prevalence and intensity. The risk factors considered were age, location (urban vs rural cats) and time spent in shelter. In total, 290 faecal samples were collected from cats at an animal shelter in 2015-2016 and investigated for endoparasites using the concentration flotation technique. In total, 138 shelter cats (47.6%) were infected with endoparasites and their overall prevalence was: Toxocara cati (36.6%), Cystoisospora spp. (12.4%), Taeniidae gen. sp. (4.1%), Toxoplasma gondii/Hammondia hammondi (3.4%), Eucoleus aerophilus (2.1%), Cryptosporidium spp. (2.1%), Ancylostoma sp. (0.7%) and Giardia sp. (0.7%). Coinfections occurred in 38 cats (13.1%) most frequently of T. cati and Cystoisospora spp. (4.5%), Cystoisospora spp. and T. gondii/H. hammondi (2.1%). Where species identification of cestode and nematode samples was not possible according to morphology, genetic analysis of the mitochondrial cox1 gene was carried out. DNA was successfully analysed for 6 out of 13 samples that required genetic identification, revealing Ancylostoma tubaeforme in one nematode sample and Hydatigera taeniaeformis in five cestode samples. Cats from rural areas had significantly higher endoparasite prevalence than cats from urban areas. Helminth prevalence decreased to some extent due to anthelmintic treatment in cats available for adoption (held ≥15 days in the shelter), whereas the prevalence of infection with protists increased significantly in these animals. It is important to note that the analysis revealed lower infection intensity for quarantine cats (held 1-14 days in the shelter) compared with cats available for adoption. The relatively high prevalence of endoparasites (including zoonotic) in shelter cats ready for adoption suggests that current anthelminthic procedures require improvements.

Native animal breeds constitute an invaluable pool of genetic resources in a changing environment. Discovering native breeds and safeguarding their genetic diversity through specific conservation programs is therefore of high importance. Endogenous retroviruses have proved to be a reliable genetic marker for studying the demographic history of sheep (Ovis aries). Previous research has revealed two migratory episodes of domesticated sheep from the Middle East to Europe. The first episode included predominantly ‘primitive populations’, while the second and most recent is hypothesised to have included sheep with markedly improved wool production. To examine whether the recently discovered Kihnu native sheep in Estonia have historically been part of the first migratory episode and to what extent they have preserved primitive genetic characters, we analysed retroviral insertions in 80 modern Kihnu sheep and 83 ancient sheep from the Bronze Age to Modern Period (850 BCE–1950 CE). We identified that the Kihnu sheep have preserved ‘primitive’, ‘Nordic’, and other ‘ancient’ retrotypes that were present both in archaeological and modern samples, confirming their shared ancestry and suggesting that contemporary Kihnu native sheep originate from the first migratory episode. However, over the course of history, there has been a gradual decrease in the frequency of primitive retrotypes. Furthermore, Kihnu sheep possessed several ‘novel’ retrotypes that were absent in archaeological individuals, but were shared with improvement breeds, suggesting recent crossing within the last two centuries. To preserve these ancient lineages, our results are being applied in the conservation program of the Kihnu Native Sheep Society.

Oxen – castrated male cattle (Bos taurus) – were used as draught animals in Finland, Estonia and Sweden for centuries. From the medieval period onwards, the wider phenomenon of oxen culture is evidenced by zooarchaeological material and historical documents, which reflect, among other things, the castration age. When males were castrated at a young age – a method used in Sweden, according to the ethnological sources – it affected the growth of the metacarpal bones; thus, the metacarpals of these oxen should be distinct from those of cows and bulls. However, when males were castrated as adults – a method used in Estonia – it had less effect on the bones, so the metacarpals of these oxen are more likely to resemble those of bulls. Documents from the early modern period indicate that in Finland, cattle were castrated as calves (as they were in Sweden), but some medieval bone finds indicate that they were castrated in adulthood (as they were in Estonia). To resolve this apparent contradiction between the zooarchaeological data and the written evidence, and to understand the cultural transitions in farming practices relating to draught oxen in the northern Baltic Sea region, we studied morphometric changes in 621 metacarpals in Finland and Estonia dating from the medieval and post-medieval periods (13th to 20th century) and compared them with previously published contemporaneous data from Sweden. To investigate the size and shape of the individuals and to distinguish oxen from bulls and cows, we applied widely used osteometric criteria, including the maximum breadth of the distal epiphysis and Mennerich indices. For a set of uncertain samples from medieval Turku, Finland, we identified the sex by amplifying the TSPY gene and X and Y chromosome-specific variations in the ancient DNA. Our results confirm that throughout the medieval and post-medieval periods, late castration was practised in Estonia and early castration was practised in Sweden. In Finland, however, we detected a temporal change in the castration age: cattle were castrated as adults in the medieval period and at a young age in the post-medieval period.

Although many large mammal species went extinct at the end of the Pleistocene epoch, their DNA may persist due to past episodes of interspecies admixture. However, direct empirical evidence of the persistence of ancient alleles remains scarce. Here, we present multifold coverage genomic data from four Late Pleistocene cave bears (Ursus spelaeus complex) and show that cave bears hybridized with brown bears (Ursus arctos) during the Pleistocene. We develop an approach to assess both the directionality and relative timing of gene flow. We find that segments of cave bear DNA still persist in the genomes of living brown bears, with cave bears contributing 0.9 to 2.4% of the genomes of all brown bears investigated. Our results show that even though extinction is typically considered as absolute, following admixture, fragments of the gene pool of extinct species can survive for tens of thousands of years in the genomes of extant recipient species.

Cystic echinococcosis (CE) is a severe parasitic disease caused by the species complexEchinococcus granulosussensu lato. Human infections are most commonly associated withE. granulosussensu stricto (s.s.), comprising genotypes G1 and G3. The objective of the current study was to provide first insight into the genetic diversity and phylogeography of genotype G3. Despite the epidemiological importance of the genotype, it has remained poorly explored due to the ambiguity in the definition of the genotype. However, it was recently demonstrated that long sequences of mitochondrial DNA (mtDNA) provide a reliable method to discriminate G1 and G3 from each other. Therefore, we sequenced near-complete mtDNA of 39 G3 samples, covering most of the known distribution range and host spectra of the genotype. The phylogenetic network revealed high genetic variation withinE. granulosus s.s.G3 and while G3 is significantly less prevalent worldwide than G1, the genetic diversity of both of the genotypes is equally high. We also present the results of the Bayesian phylogeographic analysis, which yielded several well-supported diffusion routes of genotype G3 originating from Turkey and Iran, suggesting the Middle East as the origin of the genotype.

Predation on livestock is a cause of serious and long-lasting conflict between farmers and wildlife, promoting negative public attitudes and endangering conservation of large carnivores. However, while large carnivores, especially the grey wolf (Canis lupus), are often blamed for killing sheep and other farm animals, free-ranging dogs may also act as predators. To develop appropriate measures for livestock protection, reliable methods for identifying predator species are critical. Identification of predators from visual examination of livestock wounds can be ambiguous and genetic analysis is strongly preferable for accurate species determination. To estimate the proportion of wolves and dogs implicated in sheep predation, we developed a sensitive genetic assay to distinguish between wolves and domestic dogs. A total of 183 predator saliva samples collected from killed sheep in Estonia were analysed. The assay identified the predator species in 143 cases (78%). Sheep were most often killed by wolves (81%); however, predation by dogs was substantial (15%). We compared the molecular results with field observations conducted by local environmental officials and recorded some disagreement, with the latter underestimating the role of dogs. As predator saliva samples collected from prey are often of poor quality, we suggest using mitochondrial DNA as a primary tool to maximise the number of successfully analysed samples. We also suggest adopting forensic DNA analysis more widely in livestock predation assessments as a legislative measure since misidentification that is biased against wolves can be counterproductive for conservation by enhancing conflict with society and leading to increased culling and poaching.

Understanding how the diet of an ecologically important mesopredator varies throughout a large region can provide an insight into the trophic responses of generalist carnivores to global change. We review the diet of the red fox Vulpes vulpes in Europe. The main aims of this study were: (1) to identify biogeographical patterns in red fox diet throughout Europe; (2) to investigate the significance of the human footprint index (a measure of anthropogenic influence on habitats) for the Europe-wide food habits of red foxes; and (3) to discuss the ecosystem role of red foxes in the face of global change. We hypothesised that in the red fox: (1) diet composition is likely to vary along a latitudinal gradient, between seasons (pooled into warm and cold periods of the year), and in relation to human impact on habitats; and (2) dietary breadth will increase with latitude and be greater during the warm period than during the cold period. We conducted a Europe-wide analysis of red fox diet, using data from 66 studies in 17 European countries. Multivariate analysis and generalised linear models were used to evaluate variation in diet composition and niche breadth in relation to latitude, season, and human footprint index. The analysis revealed that red fox diet composition varies considerably in response to latitude, human influence, and warm and cold periods. Dietary breadth decreases at higher latitudes in the cold period, but not in the warm period, and increases in areas with high human impact. Among individual dietary items, the consumption of lagomorphs is more frequent in areas with higher human footprint index. The occurrence of rodents and birds in fox diet increases uniformly with latitude, while that of plants and invertebrates responds differently to latitude in the warm and cold periods. With climate warming, red foxes may incorporate a wider spectrum of food categories into their diet in northern areas in the future. Widening of fox dietary niche breadth could have important consequences for ecosystems and for the spread of zoonotic diseases.

International audience; The larval stages of tapeworms in the species complex Echinococcus granulosus sensu lato cause a zoonotic disease known as cystic echinococcosis (CE). Within this species complex, genotypes G6 and G7 are among the most common genotypes associated with human CE cases worldwide. However, our understanding of ecology, biology and epidemiology of G6 and G7 is still limited. An essential first step towards this goal is correct genotype identification, but distinguishing genotypes G6 and G7 has been challenging. A recent analysis based on complete mitogenome data revealed that the conventional sequencing of the cox1 (366 bp) gene fragment mistakenly classified a subset of G7 samples as G6. On the other hand, sequencing complete mitogenomes is not practical if only genotype or haplogroup identification is needed. Therefore, a simpler and less costly method is required to distinguish genotypes G6 and G7. We compared 93 complete mitogenomes of G6 and G7 from a wide geographical range and demonstrate that a combination of nad2 (714 bp) and nad5 (680 bp) gene fragments would be the best option to distinguish G6 and G7. Moreover, this method allows assignment of G7 samples into haplogroups G7a and G7b. However, due to very high genetic variability of G6 and G7, we suggest to construct a phylogenetic network based on the nad2 and nad5 sequences in order to be absolutely sure in genotype assignment. For this we provide a reference dataset of 93 concatenated nad2 and nad5 sequences (1394 bp in total) containing representatives of G6 and G7 (and haplogroups G7a and G7b), which can be used for the reconstruction of phylogenetic networks.

The larval stage of the species complex Echinococcus granulosus sensu lato (s.l.) is the cause of a widespread zoonotic disease known as cystic echinococcosis (CE). The disease is highly prevalent in southern Italy and represents a serious public health issue. The main aim of this study was to characterize E. granulosus s.l. genotypes from wild boar on a continental area of Italy (Campania region), using recently developed mtDNA markers of nad2 and nad5 for reliable identification of different genotypes. Here, nad5 (680 bp) allowed for a clear identification of G1 and G3, whereas a combination of nad2 (714 bp) and nad5 (1394 bp in total) did the same for genotype G7 and its haplogroups G7a and G7b. The results of this study revealed for the first time the presence of genotype G7 in continental Italy. While haplogroup G7b was previously shown to be restricted to the islands of Corsica and Sardinia, here we demonstrate that haplogroup G7b is also present on the mainland of Italy. This work has implications in designing future strategies to reduce CE in Italy.

Anthropogenic hybridization is widely perceived as a threat to the conservation of biodiversity. Nevertheless, to date, relevant policy and management interventions are unresolved and highly convoluted. While this is due to the inherent complexity of the issue, we hereby hypothesize that a lack of agreement concerning management goals and approaches, within the scientific community, may explain the lack of social awareness on this phenomenon, and the absence of effective pressure on decision-makers. By focusing on wolf x dog hybridization in Europe, we hereby (a) assess the state of the art of issues on wolf x dog hybridization within the scientific community, (b) assess the conceptual bases for different viewpoints, and (c) provide a conceptual framework aiming at reducing the disagreements. We adopted the Delphi technique, involving a three-round iterative survey addressed to a selected sample of experts who published at Web of Science listed journals, in the last 10 years on wolf x dog hybridization and related topics. Consensus was reached that admixed individuals should always be defined according to their genetic profile, and that a reference threshold for admixture (i.e., q-value in assignment tests) should be formally adopted for their identification. To mitigate hybridization, experts agreed on adopting preventive, proactive and, when concerning small and recovering wolf populations, reactive interventions. Overall, experts' consensus waned as the issues addressed became increasingly practical, including the adoption of lethal removal. We suggest three non-mutually exclusive explanations for this trend: (i) value-laden viewpoints increasingly emerge when addressing practical issues, and are particularly diverging between experts with different disciplinary backgrounds (e.g., ecologists, geneticists); (ii) some experts prefer avoiding the risk of potentially giving carte blanche to wolf opponents to (illegally) remove wolves, based on the wolf x dog hybridization issue; (iii) room for subjective interpretation and opinions result from the paucity of data on the effectiveness of different management interventions. These results have management implications and reveal gaps in the knowledge on a wide spectrum of issues related not only to the management of anthropogenic hybridization, but also to the role of ethical values and real-world management concerns in the scientific debate.

Tapeworms of the species complex of Echinococcus granulosus sensu lato (s. l.) are the cause of a severe zoonotic disease – cystic echinococcosis, which is listed among the most severe parasitic diseases in humans and is prioritized by the World Health Organization. A stable taxonomy of E. granulosus s. l. is essential to the medical and veterinary communities for accurate and effective communication of the role of different species in this complex on human and animal health. E. granulosus s. l. displays high genetic diversity and has been divided into different species and genotypes. Despite several decades of research, the taxonomy of E. granulosus s. l. has remained controversial, especially the species status of genotypes G6–G10. Here the Bayesian phylogeny based on six nuclear loci (7387 bp in total) demonstrated, with very high support, the clustering of G6/G7 and G8/G10 into two separate clades. According to the evolutionary species concept, G6/G7 and G8/G10 can be regarded as two distinct species. Species differentiation can be attributed to the association with distinct host species, largely separate geographical distribution and low level of cross-fertilization. These factors have limited the gene flow between genotypic groups G6/G7 and G8/G10, resulting in the formation of distinct species. We discuss ecological and epidemiological differences that support the validity of these species.

The grey wolf (Canis lupus) is an iconic large carnivore that has increasingly been recognized as an apex predator with intrinsic value and a keystone species. However, wolves have also long represented a primary source of human–carnivore conflict, which has led to long-term persecution of wolves, resulting in a significant decrease in their numbers, genetic diversity and gene flow between populations. For more effective protection and management of wolf populations in Europe, robust scientific evidence is crucial. This review serves as an analytical summary of the main findings from wolf population genetic studies in Europe, covering major studies from the ‘pre-genomic era’ and the first insights of the ‘genomics era’. We analyse, summarize and discuss findings derived from analyses of three compartments of the mammalian genome with different inheritance modes: maternal (mitochondrial DNA), paternal (Y chromosome) and biparental [autosomal microsatellites and single nucleotide polymorphisms (SNPs)]. To describe large-scale trends and patterns of genetic variation in European wolf populations, we conducted a meta-analysis based on the results of previous microsatellite studies and also included new data, covering all 19 European countries for which wolf genetic information is available: Norway, Sweden, Finland, Estonia, Latvia, Lithuania, Poland, Czech Republic, Slovakia, Germany, Belarus, Russia, Italy, Croatia, Bulgaria, Bosnia and Herzegovina, Greece, Spain and Portugal. We compared different indices of genetic diversity in wolf populations and found a significant spatial trend in heterozygosity across Europe from south-west (lowest genetic diversity) to north-east (highest). The range of spatial autocorrelation calculated on the basis of three characteristics of genetic diversity was 650−850 km, suggesting that the genetic diversity of a given wolf population can be influenced by populations up to 850 km away. As an important outcome of this synthesis, we discuss the most pressing issues threatening wolf populations in Europe, highlight important gaps in current knowledge, suggest solutions to overcome these limitations, and provide recommendations for science-based wolf conservation and management at regional and Europe-wide scales.

The raccoon dog (Nyctereutes procyonoides) is a canid that is indigenous in East Asia and alien in Europe, where it was introduced more than half a century ago. The aim of this study was to compare the parasite faunas associated with raccoon dogs in their native and introduced ranges, and to identify zoonotic parasite species. We examined 255 carcasses of hunted raccoon dogs from Estonia and recorded a total of 17 helminth species: 4 trematodes, 4 cestodes and 9 nematodes. The most prevalent parasite species were Uncinaria stenocephala (97.6%) and Alaria alata (68.3%). Average parasite species richness was 2.86 (the highest was 9) and only two animals were not parasitized at all. Although the infection intensity was determined by weight and not by sex, all animals infected with more than five helminth species were males. We also found that animals infected with higher numbers of helminth species fed significantly more on natural plants. Intentional consumption of grass may represent a self-medicating behaviour among raccoon dogs. We included the Estonian data into a wider comparison of raccoon dog parasite faunas and found a total of 54 helminth taxa, including 28 of zoonotic potential. In Europe, raccoon dogs are infected with a minimum of 32 helminth species of which 19 are zoonotic; in the native range they are infected with 26 species of which 17 are zoonotic. Most species were nematodes or trematodes, with fewer cestodes described. The recent increase in the number and range of raccoon dogs in Europe and the relatively high number of zoonotic parasite taxa that it harbours suggests that this species should be considered an important source of environmental contamination with zoonotic agents in Europe.

Summary Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D-loop of mitochondrial DNA (mtDNA; 213 bp) and the Y-chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G-oY1 of the Y-chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.

International audience; Cystic echinococcosis (CE) is a zoonotic disease caused by the larval stage of the species complex Echinococcus granulosus sensu lato. Within this complex, genotypes G6 and G7 have been frequently associated with human CE worldwide. Previous studies exploring the genetic variability and phylogeography of genotypes G6 and G7 have been based on relatively short mtDNA sequences, and the resolution of these studies has often been low. Moreover, using short sequences, the distinction between G6 and G7 has in some cases remained challenging. The aim here was to sequence complete mitochondrial genomes (mitogenomes) to obtain deeper insight into the genetic diversity, phylogeny and population structure of genotypes G6 and G7. We sequenced complete mitogenomes of 94 samples collected from 15 different countries worldwide. The results demonstrated that (i) genotypes G6 and G7 can be clearly distinguished when mitogenome sequences are used; (ii) G7 is represented by two major haplogroups, G7a and G7b, the latter being specific to islands of Corsica and Sardinia; (iii) intensive animal trade, but also geographical isolation, have likely had the largest impact on shaping the genetic structure and distribution of genotypes G6 and G7. In addition, we found phylogenetically highly divergent haplotype from Mongolia (Gmon), which had a higher affinity to G6.

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.

International audience; Echinococcus granulosus sensu stricto (s.s.) is the major cause of human cystic echinococcosis worldwide and is listed among the most severe parasitic diseases of humans. To date, numerous studies have investigated the genetic diversity and population structure of E. granulosus s.s. in various geographic regions. However, there has been no global study. Recently, using mitochondrial DNA, it was shown that E. granulosus s.s. G1 and G3 are distinct genotypes, but a larger dataset is required to confirm the distinction of these genotypes. The objectives of this study were to: (i) investigate the distinction of genotypes G1 and G3 using a large global dataset; and (ii) analyse the genetic diversity and phylogeography of genotype G1 on a global scale using near-complete mitogenome sequences. For this study, 222 globally distributed E. granulosus s.s. samples were used, of which 212 belonged to genotype G1 and 10 to G3. Using a total sequence length of 11,682 bp, we inferred phylogenetic networks for three datasets: E. granulosus s.s. (n = 222), G1 (n = 212) and human G1 samples (n = 41). In addition, the Bayesian phylogenetic and phylogeographic analyses were performed. The latter yielded several strongly supported diffusion routes of genotype G1 originating from Turkey, Tunisia and Argentina. We conclude that: (i) using a considerably larger dataset than employed previously, E. granulosus s.s. G1 and G3 are indeed distinct mitochondrial genotypes; (ii) the genetic diversity of E. granulosus s.s. G1 is high globally, with lower values in South America; and (iii) the complex phylogeographic patterns emerging from the phylogenetic and geographic analyses suggest that the current distribution of genotype G1 has been shaped by intensive animal trade.

International audience; Cystic echinococcosis (CE), a zoonotic disease caused by tapeworms of the species complex Echinococcus granulosus sensu lato, represents a substantial global health and economic burden. Within this complex, E. granulosus sensu stricto (genotypes G1 and G3) is the most frequent causative agent of human CE. Currently, there is no fully reliable method for assigning samples to genotypes G1 and G3, as the commonly used mitochondrial cox1 and nad1 genes are not sufficiently consistent for the identification and differentiation of these genotypes. Thus, a new genetic assay is required for the accurate assignment of G1 and G3. Here we use a large dataset of near-complete mtDNA sequences (n = 303) to reveal the extent of genetic variation of G1 and G3 on a broad geographical scale and to identify reliable informative positions for G1 and G3. Based on extensive sampling and sequencing data, we developed a new method, that is simple and cost-effective, to designate samples to genotypes G1 and G3. We found that the nad5 is the best gene in mtDNA to differentiate between G1 and G3, and developed new primers for the analysis. Our results also highlight problems related to the commonly used cox1 and nad1. To guarantee consistent identification of G1 and G3, we suggest using the sequencing of the nad5 gene region (680 bp). This region contains six informative positions within a relatively short fragment of the mtDNA, allowing the differentiation of G1 and G3 with confidence. Our method offers clear advantages over the previous ones, providing a significantly more consistent means to distinguish G1 and G3 than the commonly used cox1 and nad1.

In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (N e) was estimated to be 199 individuals, and total combined population size (N C) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe's low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans., Using microsatellites and mtDNA, we reveal the genetic structure, diversity, population size of the brown bear in Greece, and its detailed phylogeography for the first time, revealing a new matrilinear contact zone

Protecting animal aggregation sites is intuitively an efficient conservation approach, particularly for threatened species in fragmented landscapes. However, the appropriate scale of protection depends on accompanying threats, such as predator attraction to the same aggregation sites, to which the target species may need to respond by moving to alternative habitat patches. We performed experiments using artificial ground nests around protected capercaillie leks in commercially managed forest landscape in Estonia. We considered two scales: the landscape up to 3 km from the lek centre, and among 10–30-ha forest compartments, which differed in predator abundance. We found that nest predation significantly declined with distance from the lek. Nests were depredated by multiple mammalian and avian predator species and total predator abundance explained most of the between-forest variation in predation rate. Our results indicate that, in this hemiboreal area, (i) ground nest predation is largely determined by landscape-scale distribution of predators and (ii) predators can aggregate at capercaillie leks. The implication is that lek-centred habitat protection used for the capercaillie and other grouse may be ineffective unless the peripheral habitat quality and predator abundance are specifically addressed. More generally, predation pressure can be a serious problem in small set-asides within hostile landscapes for threatened ground-nesting birds.

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.

Cystic echinococcosis, a zoonotic disease caused by Echinococcus granulosus sensu lato (s.l.), is a significant global public health concern. Echinococcus granulosus s. l. is currently divided into numerous genotypes (G1-G8 and G10) of which G1-G3 are the most frequently implicated genotypes in human infections. Although it has been suggested that G1-G3 could be regarded as a distinct species E. granulosus sensu stricto (s. s.), the evidence to support this is inconclusive. Most importantly, data from nuclear DNA that provide means to investigate the exchange of genetic material between G1-G3 is lacking as none of the published nuclear DNA studies have explicitly included G2 or G3. Moreover, the commonly used relatively short mtDNA sequences, including the complete coxl gene, have not allowed unequivocal differentiation of genotypes G1-G3. Therefore, significantly longer mtDNA sequences are required to distinguish these genotypes with confidence. The main aim of this study was to evaluate the phylogenetic relations and taxonomy of genotypes G1-G3 using sequences of nearly complete mitogenomes (11,443 bp) and three nuclear loci (2984 bp). A total of 23 G1-G3 samples were analysed, originating from 5 intermediate host species in 10 countries. The mtDNA data demonstrate that genotypes G1 and G3 are distinct mitochondrial genotypes (separated by 37 mutations), whereas G2 is not a separate genotype or even a monophyletic cluster, but belongs to G3. Nuclear data revealed no genetic separation of G1 and G3, suggesting that these genotypes form a single species due to ongoing gene flow. We conclude that: (a) in the taxonomic sense, genotypes G1 and G3 can be treated as a single species E. granulosus s. s.; (b) genotypes G1 and G3 should be regarded as distinct genotypes only in the context of mitochondrial data; (c) we recommend excluding G2 from the genotype list. (C) 2017 Elsevier B.V. All rights reserved.

Sarcoptic mange is a highly contagious zoonotic skin disease that can have severe effect on population dynamics of many wild mammals. However, very little is known about its effect on the activity and space use of infected animals. In this study we equipped two raccoon dogs (Nyctereutes procyonoides) in a mange outbreak area with radio-collars and observed the spatial and temporal effects of sarcoptic mange to these individuals. Initially, the raccoon dogs had no external symptoms of mange infection, but developed these during the study period. One of the raccoon dogs died just 32 and the other 52 days after collaring. During a relatively short period before their death, there was an abrupt and drastic decline in their home range size. For one of the animals it started about 1 month before its death and the home range size reduced >1000×, whereas for the other raccoon dog it took place within the last week, decreasing >200×. The daily covered distances also declined considerably. These results indicate that at the later stage sarcoptic mange affected the physiological state of the diseased animals so forcefully that made them almost immobile. Our results show for the first time how rapid and severe could the impact of sarcoptic mange be for one of its most important reservoir and vector species.

Capercaillie (Tetrao urogallus) and other grouse species represent conservation concerns across Europe due to their negative abundance trends. In addition to habitat deterioration, predation is considered a major factor contributing to population declines. While the role of generalist predators on grouse predation is relatively well known, the impact of the omnivorous wild boar has remained elusive. We hypothesize that wild boar is an important predator of ground-nesting birds, but has been neglected as a bird predator because traditional morphological methods underestimate the proportion of birds in wild boar diet. To distinguish between different mammalian predator species, as well as different grouse prey species, we developed a molecular method based on the analysis of mitochondrial DNA that allows accurate species identification. We collected 109 wild boar faeces at protected capercaillie leks and surrounding areas and analysed bird consumption using genetic methods and classical morphological examination. Genetic analysis revealed that the proportion of birds in wild boar faeces was significantly higher (17.3%; 4.5×) than indicated by morphological examination (3.8%). Moreover, the genetic method allowed considerably more precise taxonomic identification of consumed birds compared to morphological analysis. Our results demonstrate: (i) the value of using genetic approaches in faecal dietary analysis due to their higher sensitivity, and (ii) that wild boar is an important predator of ground-nesting birds, deserving serious consideration in conservation planning for capercaillie and other grouse.

The survival of isolated small populations is threatened by both demographic and genetic factors. Large carnivores declined for centuries in most of Europe due to habitat changes, overhunting of their natural prey and direct persecution. However, the current rewilding trends are driving many carnivore populations to expand again, possibly reverting the erosion of their genetic diversity. In this study we reassessed the extent and origin of the genetic variation of the Italian Wolf population, which is expanding after centuries of decline and isolation. We genotyped wolves from Italy and other nine populations at four mtDNA regions (control-region, ATP6, COIII and ND4) and 39 autosomal microsatellites. Results of phylo-genetic analyses and assignment procedures confirmed in the Italian wolves a second private mtDNA haplotype, which belongs to a haplogroup distributed mostly in southern Europe. Coalescent analyses showed that the unique mtDNA haplotypes in the Italian wolves likely originated during the late Pleistocene. ABC simulations concordantly showed that the extant Wolf populations in Italy and in south-western Europe started to be isolated and declined right after the last glacial maximum. Thus, the standing genetic variation in the Italian wolves principally results from the historical isolation south of the Alps.