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.

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 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

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.

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.

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.

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), 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.

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.

Brown and polar bears have become prominent examples in phylogeography, but previous phylogeographic studies relied largely on maternally inherited mitochondrial DNA (mtDNA) or were geographically restricted. The male-specific Y chromosome, a natural counterpart to mtDNA, has remained underexplored. Although this paternally inherited chromosome is indispensable for comprehensive analyses of phylogeographic patterns, technical difficulties and low variability have hampered its application in most mammals. We developed 13 novel Y-chromosomal sequence and microsatellite markers from the polar bear genome and screened these in a broad geographic sample of 130 brown and polar bears. We also analyzed a 390-kb-long Y-chromosomal scaffold using sequencing data from published male ursine genomes. Y chromosome evidence support the emerging understanding that brown and polar bears started to diverge no later than the Middle Pleistocene. Contrary to mtDNA patterns, we found 1) brown and polar bears to be reciprocally monophyletic sister (or rather brother) lineages, without signals of introgression, 2) male-biased gene flow across continents and on phylogeographic time scales, and 3) male dispersal that links the Alaskan ABC islands population to mainland brown bears. Due to female philopatry, mtDNA provides a highly structured estimate of population differentiation, while male-biased gene flow is a homogenizing force for nuclear genetic variation. Our findings highlight the importance of analyzing both maternally and paternally inherited loci for a comprehensive view of phylogeographic history, and that mtDNA-based phylogeographic studies of many mammals should be reevaluated. Recent advances in sequencing technology render the analysis of Y-chromosomal variation feasible, even in nonmodel organisms.

SUMMARYEchinococcus granulosusis the causative agent of cystic echinococcosis. The disease is a significant global public health concern and human infections are most commonly associated withE. granulosussensu stricto (s. s.) genotype G1. The objectives of this study were to: (i) analyse the genetic variation and phylogeography ofE. granulosuss. s. G1 in part of its main distribution range in Europe using 8274 bp of mtDNA; (ii) compare the results with those derived from previously used shorter mtDNA sequences and highlight the major differences. We sequenced a total of 91E. granulosuss. s. G1 isolates from six different intermediate host species, including humans. The isolates originated from seven countries representing primarily Turkey, Italy and Spain. Few samples were also from Albania, Greece, Romania and from a patient originating from Algeria, but diagnosed in Finland. The analysed 91 sequences were divided into 83 haplotypes, revealing complex phylogeography and high genetic variation ofE. granulosuss. s. G1 in Europe, particularly in the high-diversity domestication centre of western Asia. Comparisons with shorter mtDNA datasets revealed that 8274 bp sequences provided significantly higher phylogenetic resolution and thus more power to reveal the genetic relations between different haplotypes.

SUMMARYRed fox (Vulpes vulpes) is the most widely distributed canid in the world and an important source of multiple zoonotic pathogens capable of causing life-threatening diseases, such as rabies and alveolar echinococcosis. Informing general public of potential risks related to foxes is becoming more important since the fox densities have increased in many countries and the species is colonizing urban areas in Europe and around the world with increasing pace, bringing zoonotic pathogens to the immediate neighbourhood of humans and their companion animals. The aim of this study was to examine the parasite fauna of red foxes in Estonia. We found in Estonian foxes a total of 17 endoparasite taxa, including ten zoonotic species. All the analysed individuals were infected and the average parasite species richness was 6·37. However, the infection rates varied to a very large extent for different parasite species, ranging from 0·9 to 91·5%. Of zoonotic species, the highest infection rate was observed forAlaria alata(90·7%),Eucoleus aerophilus(87·6%) andUncinaria stenocephala(84·3%). The prevalence of tapewormEchinococcus multilocularis, a causative agent for alveolar echinococcosis, was also relatively high (31·5%), presenting a potential risk to human health.

Echinococcus granulosus is a taeniid cestode and the etiological agent of an infectious zoonotic disease known as cystic echinococcosis (CE) or hydatid disease. CE is a serious public health concern in many parts of the world, including the Americas, where it is highly endemic in many regions. Echinococcus granulosus displays high intraspecific genetic variability and is divided into multiple genotypes (G1-G8, G10) with differences in their biology and etiology. Of these, genotype G1 is responsible for the majority of human and livestock infections and has the broadest host spectrum. However, despite the high significance to the public and livestock health, the data on genetic variability and regional genetic differences of genotype G1 in America are scarce. The aim of this study was to evaluate the genetic variability and phylogeography of G1 in several countries in America by sequencing a large portion of the mitochondrial genome. We analysed 8279bp of mtDNA for 52 E. granulosus G1 samples from sheep, cattle and pigs collected in Argentina, Brazil, Chile and Mexico, covering majority of countries in the Americas where G1 has been reported. The phylogenetic network revealed 29 haplotypes and a high haplotype diversity (Hd=0.903). The absence of phylogeographic segregation between different regions in America suggests the importance of animal transportation in shaping the genetic structure of E. granulosus G1. In addition, our study revealed many highly divergent haplotypes, indicating a long and complex evolutionary history of E. granulosus G1 in the Americas.

Outdoor farming has been traditionally implemented in some European rural areas, but was mostly limited to backyard or small farming systems, however we are now observing an increment in this practice, as a result of both increased attention to animal welfare and consumers interest in organic food. At the same time, Europe is showing interest in rewilding projects. This combination might expand the wildlife-livestock interface, with enhanced sanitary risk and hybridization.Hybridization events between the domestic pig (Sus scrofa domestica; hereafter DP) and the wild boar (Sus scrofa; hereafter WB) have been reported in several European countries and might provide insight on human practices that can lead to increased, although spatially delimited, risk of contagion. The markers most commonly used to investigate this aspect are mitochondrial DNA (mtDNA), MC1R and microsatellite (STR). All of them present some advantages and some disadvantages. MtDNA and MC1R can be easily compared among laboratories, while STR require calibration. Furthermore, STR can only detect recent hybridization events, while MC1R and mtDNA can also identify more ancient episodes. On the other hand STR and MC1R are biparentally inherited, thus representing the whole population, while mtDNA only shows the maternal lineage. In recent years the development and decrease in price of genomics techniques have allowed the implementation of genomic markers (like whole genome sequencing or Single Nucleotide Polymorphisms - SNP) for the study of genomic ancestry and hybridization. SNP markers are easily comparable between laboratories, biparentally inherited and allow the simultaneous study of loci neutral and under selection.Here we present a comparison based on data from several European countries of genetic hybridization levels, detected with SNP markers, with data on hybridization perception, collected through an online survey. The survey covers 29 countries, with a major contribution from academic researchers (70%) but not negligible answer rate from the non-academic component (28%). Genetic data are based on 235 WB from 22 areas 149 DP (49 from 5 commercial lines and 100 from 9 local breeds) analysed with the Porcine SNP60 Beadchip (64,232 SNP).Although data are not perfectly comparable, due to incomplete overlap of sampled areas and differences of methods, results show some interesting aspects. Genetic data show a clear separation between WB and DP, with a limited number of hybrids in both populations. The introgression level varies considerably among populations, from non-detectable to very high. Perceived presence of hybrids, based on phenotypic characteristics or historic data, is usually higher and widespread than that observed with SNP data, especially in the WB population. Participants that stated hybridization was a past issue, often still report the presence of introgressed phenotypes in WB (50%) or both WB and DP (29%).These results show the importance of considering both available information and hard data as they might help interpretation of outcomes or complement each other. At the same time, they point out the need to investigate further the genomes of hybrids individuals, as they might show the presence of regions under selection introgressed from one of the two forms into the other that might help better understand selection, introgression and genetic basis of WB invasiveness and that would be useful in the development of contingency or management plans.

Aim Using sequences of complete mitochondrial genomes, our aims were: (1) to investigate the matrilineal phylogeographical structure, migration patterns and lineage coalescence times in a large, continuous population of brown bears (Ursus arctos); and (2) to develop a novel spatial genetic method to identify migration corridors and barriers. Location North-western Eurasia: from eastern European Russia to the Baltic Sea. Methods We sequenced the complete mitochondrial genomes of 95 brown bears. The phylogeographical resolution of complete genomes was compared to that derived from subsets of the genome, including the most commonly used shorter sequences. We conducted network and Bayesian phylogeographical analyses and developed a novel, spatially explicit, individual-based approach (called DResD) for identifying migration corridors and barriers. Results Analysis of mitogenome sequences revealed five haplogroups, specific to particular geographical areas, exhibiting far greater resolving power than shorter sequences. Estimated coalescence times for the haplogroups ranged from 7.7 to 15.2 ka, suggesting that their divergence took place after the last glaciation. We found several migration trends, including a large westward migration from eastern European Russia towards Finland. We also found evidence of a potential barrier and a migration corridor in the south-west of the study area. Main conclusions The use of complete mitochondrial genomes from a brown bear population in north-western Eurasia allowed us to identify phylogeographical structure, signatures of demographic history and spatial processes that had not previously been detected using shorter sequences. These findings have implications for studies on other species and populations, especially those exhibiting low mtDNA diversity. The relatively recent divergence estimates for haplogroups highlight the significance not only of the last glaciation but also of climatic fluctuations during the post-glacial period for the divergence of mammal populations in Europe. Our spatial genetic method represents a new tool for the analysis of genetic data in a geographical context and is applicable to any data that yield genetic distance matrices, including microsatellites, amplified fragment length polymorphisms (AFLPs) and single-nucleotide polymorphisms (SNPs).

The ignita species group within the genus Chrysis includes over 100 cuckoo wasp species, which all lead a parasitic lifestyle and exhibit very similar morphology. The lack of robust, diagnostic morphological characters has hindered phylogenetic reconstructions and contributed to frequent misidentification and inconsistent interpretations of species in this group. Therefore, molecular phylogenetic analysis is the most suitable approach for resolving the phylogeny and taxonomy of this group. We present a well-resolved phylogeny of the Chrysis ignita species group based on mitochondrial sequence data from 41 ingroup and six outgroup taxa. Although our emphasis was on European taxa, we included samples from most of the distribution range of the C. ignita species group to test for monophyly. We used a continuous mitochondrial DNA sequence consisting of 16S rRNA, tRNA(Val), 12S rRNA and ND4. The location of the ND4 gene at the 3' end of this continuous sequence, following 12S rRNA, represents a novel mitochondrial gene arrangement for insects. Due to difficulties in aligning rRNA genes, two different Bayesian approaches were employed to reconstruct phylogeny: (1) using a reduced data matrix including only those positions that could be aligned with confidence; or (2) using the full sequence dataset while estimating alignment and phylogeny simultaneously. In addition maximum-parsimony and maximum-likelihood analyses were performed to test the robustness of the Bayesian approaches. Although all approaches yielded trees with similar topology, considerably more nodes were resolved with analyses using the full data matrix. Phylogenetic analysis supported the monophyly of the C. ignita species group and divided its species into well-supported clades. The resultant phylogeny was only partly in accordance with published subgroupings based on morphology. Our results suggest that several taxa currently treated as subspecies or names treated as synonyms may in fact constitute separate species. Our study provides a solid basis for further systematic investigations of this enigmatic insect group.

This review provides an up-to-date synthesis of the matrilineal phylogeography of a uniquely wellstudied Holarctic mammal, the brown bear.We extend current knowledge by presenting a DNA sequence derived from one of the earliest known fossils of a polar bear (dated to 115 000 years before present), a species that shares a paraphyletic mitochondrial association with brown bears. A molecular clock analysis of 140 mitochondrial DNA sequences, including our new polar bear sequence, provides novel insights into the times of origin for different brown bear clades. We propose a number of regional biogeographic scenarios based on genetic data, divergence time estimates and paleontological records. The case of the brown bear provides an example for researchers working with less well-studied taxa: it shows clearly that phylogeographic models based on patterns of modern genetic variation alone can be substantially improved by including data on historical patterns of genetic diversity in the form of ancient DNA sequences derived from accurately dated samples and by using an approach to divergence-time estimation that suits the data under analysis. Using such approaches it has been possible to (i) establish that the processes shaping modern genetic diversity in brown bears acted recently, within the last three glacial cycles; (ii) distinguish among hypotheses concerning species' responses to climatic oscillations in accordance with the lack of phylogeographic structure that existed in brown bears prior to the last glacial maximum (LGM); (iii) reassess theories linking monophyletic brown bear populations to particular LGM refuge areas; and (iv) identify vicariance events and track analogous patterns of migration by brown bears out of Eurasia to North America and Japan. Refereed/Peer-reviewed

Knowledge of population structure and genetic diversity and the spatio-temporal demographic processes affecting populations is crucial for effective wildlife preservation, yet these factors are still poorly understood for organisms with large continuous ranges. Available population genetic data reveal that widespread mammals have for the most part only been carefully studied at the local population scale, which is insufficient for understanding population processes at larger scales. Here, we provide data on population structure, genetic diversity and gene flow in a brown bear population inhabiting the large territory of northwestern Eurasia. Analysis of 17 microsatellite loci indicated significant population substructure, consisting of four genetic groups. While three genetic clusters were confined to small geographical areas-located in Estonia, southern Finland and Leningrad oblast, Russia-the fourth cluster spanned a very large area broadly falling between northern Finland and the Arkhangelsk and Kirov oblasts of Russia. Thus, the data indicate a complex pattern where a fraction of the population exhibits large-scale gene flow that is unparalleled by other wild mammals studied to date, while the remainder of the population appears to have been structured by a combination of demographic history and landscape barriers. These results based on nuclear data are generally in good agreement with evidence previously derived using mitochondrial markers, and taken together, these markers provide complementary information about female-specific and population-level processes. Moreover, this study conveys information about spatial processes occurring over multiple generations that cannot be readily gained using other approaches, e.g. telemetry.

Omnivores exploit numerous sources of protein and other nutrients throughout the year, and meat is generally considered a high-quality resource. However, it is unknown if there is any general association between latitude and carnivorous behavior in omnivorous mammals. We examined the relative importance of meat and other dietary components, including anthropogenic food items, in the diet of brown bears (Ursus arctos) in Estonia using conventional scat- and stomach-content analyses as well as stable-isotope (δ15N, δ13C) analyses. When food habits of brown bears in Estonia were compared with those of other populations in central and northern Europe, the proportion of animal prey in the diet was positively correlated with latitude. Further comparison with the data on the diet of two other omnivorous mammals, the European badger (Meles meles) and the European pine marten (Martes martes), provides evidence that increased carnivory towards northern latitudes may be a general adaptation in omnivorous mammals.

The brown bear has proved a useful model for studying Late Quaternary mammalian phylogeography. However, information is lacking from northern continental Eurasia, which constitutes a large part of the species' current distribution. We analysed mitochondrial DNA sequences (totalling 1943 bp) from 205 bears from northeast Europe and Russia in order to characterize the maternal phylogeography of bears in this region. We also estimated the formation times of the sampled brown bear lineages and those of its extinct relative, the cave bear. Four closely related haplogroups belonging to a single mitochondrial subclade were identified in northern continental Eurasia. Several haplotypes were found throughout the whole study area, while one haplogroup was restricted to Kamchatka. The haplotype network, estimated divergence times and various statistical tests indicated that bears in northern continental Eurasia recently underwent a sudden expansion, preceded by a severe bottleneck. This brown bear population was therefore most likely founded by a small number of bears that were restricted to a single refuge area during the last glacial maximum. This pattern has been described previously for other mammal species and as such may represent one general model for the phylogeography of Eurasian mammals. Bayesian divergence time estimates are presented for different brown and cave bear clades. Moreover, our results demonstrate the extent of substitution rate variation occurring throughout the phylogenetic tree, highlighting the need for appropriate calibration when estimating divergence times.

The tribe Lythriini is a small group of diurnally active geometrid moths consisting of a single Palaearctic genus Lythria with five species. The systematic placement of Lythriini has remained controversial: though traditionally it has been placed into the subfamily Larentiinae, a number of morphological characters link this tribe with the subfamily Sterrhinae. A molecular phylogenetic study was conducted to verify the systematic position of Lythriini, using sequences of both mitochondrial and nuclear genes: elongation factor 1α (EF-1α), wingless (wgl), 28S rRNA expansion segment D2 (28S D2), cytochrome oxidase subunit 1 (COI) and NADH dehydrogenase subunit 1 (ND1) (a total of 3784 bp). Phylogenetic analysis reliably demonstrated that Lythriini belong to the subfamily Sterrhinae. Therefore, we propose to remove tribe Lythriini from Larentiinae and unite it with Sterrhinae. Moreover, our analysis supports the monophyly of both Sterrhinae and Larentiinae. However, although both morphological data and interspecific genetic distances insinuated that Lythria cruentaria and L. sanguinaria are sister species, the latter formed a clade of sister taxa together with L. purpuraria.

We analyzed phylogenetic relationships of brown bear (Ursus arctos) maternal lineages in Finland using nucleotide sequences of the mitochondrial control region of 135 individuals. A total of 14 haplotypes, which belong to 2 haplogroups, were characterized. Haplogroup I is spread across the entire bear range, whereas haplogroup D is restricted to the southeastern part of the country. Most unconventional was the concentration of haplotype C in a small territory in central Finland, which is explained by the natal philopatry of bears that descended from one or more translocated females. Analysis of Finnish brown bear haplotypes together with 15 European eastern lineage haplotypes revealed that the Romanian–Slovakian region was probably an important refuge (or route) for recolonization of Finland by brown bears after the last ice age.

We estimated the phylogenetic relationships of brown bear maternal haplotypes from countries of northeastern Europe (Estonia, Finland and European Russia), using sequences of mitochondrial DNA (mtDNA) control region of 231 bears. Twenty-five mtDNA haplotypes were identified. The brown bear population in northeastern Europe can be divided into three haplogroups: one with bears from all three countries, one with bears from Finland and Russia, and the third composed almost exclusively of bears from European Russia. Four haplotypes from Finland and European Russia matched exactly with haplotypes from Slovakia, suggesting the significance of the current territory of Slovakia in ancient demographic processes of brown bears. Based on the results of this study and those from the recent literature, we hypothesize that the West Carpathian Mountains have served either as one of the northernmost refuge areas or as an important movement corridor for brown bears of the Eastern lineage towards northern Europe during or after the last ice age. Bayesian analyses were performed to investigate the temporal framework of brown bear lineages in Europe. The molecular clock was calibrated using Beringian brown bear sequences derived from radiocarbon-dated ancient samples, and the estimated mutation rate was 29.8% (13.3%− − − 47.6%) per million years. The whole European population and Western and Eastern lineages formed about 175 000, 70 000 and 25 000 years before present, respectively. Our approach to estimating the time frame of brown bear evolution demonstrates the importance of using an appropriate mutation rate, and this has implications for other studies of Pleistocene populations.