The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs.

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Calf/female ratio and population dynamics of wild forest reindeer in relation to wolf and moose abundances in a managed European ecosystem.

BACKGROUND: The alternative prey hypothesis describes the mechanism for apparent competition whereby the mortality of the secondary prey species increases (and population size decreases decreases) by the increased predation by the shared predator if the population size of the primary prey decreases. Apparent competition is a process where the abundance of two co-existing prey species are negatively associated because they share a mutual predator, which negatively affects the abundance of both prey Here, we examined whether alternative prey and/or apparent competition hypothesis can explain the population dynamics and reproductive output of the secondary prey, wild forest reindeer (Rangifer tarandus fennicus) in Finland, in a predator-prey community in which moose (Alces alces) is the primary prey and the wolf (Canis lupus) is the generalist predator. METHODS: We examined a 22-year time series (1996-2017) to determine how the population size and the calf/female ratio of wild forest reindeer in Eastern Finland were related to the abundances of wolf and moose. Only moose population size was regulated by hunting. Summer predation of wolves on reindeer focuses on calves. We used least squares regression (GLS) models (for handling autocorrelated error structures and resulting pseudo-R2s) and generalized linear mixed (GLMs) models (for avoidance of negative predictions) to determine the relationships between abundances. We performed linear and general linear models for the calf/female ratio of reindeer. RESULTS AND SYNTHESIS: The trends in reindeer population size and moose abundance were almost identical: an increase during the first years and then a decrease until the last years of our study period. Wolf population size in turn did not show long-term trends. Change in reindeer population size between consecutive winters was related positively to the calf/female ratio. The calf/female ratio was negatively related to wolf population size, but the reindeer population size was related to the wolf population only when moose abundance was entered as another independent variable. The wolf population was not related to moose abundance even though it is likely to consist the majority of the prey biomass. Because reindeer and moose populations were positively associated, our results seemed to support the alternative prey hypothesis more than the apparent competition hypothesis. However, these two hypotheses are not mutually exclusive and the primary mechanism is difficult to distinguish as the system is heavily managed by moose hunting. The recovery of wild forest reindeer in eastern Finland probably requires ecosystem management involving both habitat restoration and control of species abundances.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples.

BACKGROUND: Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples. RESULTS: Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations. CONCLUSIONS: We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Whole-genome analyses provide no evidence for dog introgression in Fennoscandian wolf populations.

Hybridization and admixture can threaten the genetic integrity of populations and be of particular concern to endangered species. Hybridization between grey wolves and dogs has been documented in many wolf populations worldwide and is a prominent example of human-mediated hybridization between a domesticated species and its wild relative. We analysed whole-genome sequences from >200 wolves and >100 dogs to study admixture in Fennoscandian wolf populations. A principal component analysis of genetic variation and admixture showed that wolves and dogs were well-separated, without evidence for introgression. Analyses of local ancestry revealed that wolves had <1% mixed ancestry, levels comparable to the degree of mixed ancestry in many dogs, and likely not resulting from recent wolf-dog hybridization. We also show that the founders of the Scandinavian wolf population were genetically inseparable from Finnish and Russian Karelian wolves, pointing at the geographical origin of contemporary Scandinavian wolves. Moreover, we found Scandinavian-born animals among wolves sampled in Finland, demonstrating bidirectional gene flow between the Scandinavian Peninsula and eastern countries. The low incidence of admixture between wolves and dogs in Fennoscandia may be explained by the fact that feral dogs are rare in this part of Europe and that careful monitoring and management act to remove hybrids before they backcross into wolf populations.

The evolutionary history of grey wolf Y chromosomes.

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Can only poorer European countries afford large carnivores?

BACKGROUND: One of the classic approaches in environmental economics is the environmental Kuznets curve, which predicts that when a national economy grows from low to medium levels, threats to biodiversity conservation increase, but they decrease when the economy moves from medium to high. We evaluated this approach by examining how population densities of the brown bear (Ursus arctos), gray wolf (Canis lupus), and Eurasian lynx (Lynx lynx) were related to the national economy in 24 European countries. METHODOLOGY/PRINCIPAL FINDINGS: We used forest proportions, the existence of a compensation system, and country group (former socialist countries, Nordic countries, other countries) as covariates in a linear model with the first- and the second-order polynomial terms of per capita gross domestic product (GDP). Country group was treated as a random factor, but remained insignificant and was ignored. All models concerning brown bear and wolf provided evidence that population densities decreased with increasing GDP, but densities of lynx were virtually independent of GDP. Models for the wolf explained >80% of the variation in densities, without a difference between the models with all independent variables and the model with only GDP. For the bear, the model with GDP alone accounted for 10%, and all three variables 33%, of the variation in densities. CONCLUSIONS: Wolves exhibit a higher capacity for dispersal and reproduction than bear or lynx, but still exists at the lowest densities in wealthy European countries. We are aware that several other factors, not available for our models, influenced large carnivore densities. Based on the pronounced differences among large carnivore species in their countrywide relationships between densities and GDP, and a strikingly high relationship for the gray wolf, we suggest that our results reflected differences in political history and public acceptance of these species among countries. The compensation paid for the damages caused by the carnivores is not a key to higher carnivore densities, but might be necessity for the presence of large carnivores to be accepted in countries with high GDP.

Non-invasive genetic monitoring involving citizen science enables reconstruction of current pack dynamics in a re-establishing wolf population.

BACKGROUND: Carnivores are re-establishing in many human-populated areas, where their presence is often contentious. Reaching consensus on management decisions is often hampered by a dispute over the size of the local carnivore population. Understanding the reproductive dynamics and individual movements of the carnivores can provide support for management decisions, but individual-level information can be difficult to obtain from elusive, wide-ranging species. Non-invasive genetic sampling can yield such information, but makes subsequent reconstruction of population history challenging due to incomplete population coverage and error-prone data. Here, we combine a collaborative, volunteer-based sampling scheme with Bayesian pedigree reconstruction to describe the pack dynamics of an establishing grey wolf (Canis lupus) population in south-west Finland, where wolf breeding was recorded in 2006 for the first time in over a century. RESULTS: Using DNA extracted mainly from faeces collected since 2008, we identified 81 individual wolves and assigned credible full parentages to 70 of these and partial parentages to a further 9, revealing 7 breeding pairs. Individuals used a range of strategies to obtain breeding opportunities, including dispersal to established or new packs, long-distance migration and inheriting breeding roles. Gene flow occurred between all packs but inbreeding events were rare. CONCLUSIONS: These findings demonstrate that characterizing ongoing pack dynamics can provide detailed, locally-relevant insight into the ecology of contentious species such as the wolf. Involving various stakeholders in data collection makes these results more likely to be accepted as unbiased and hence reliable grounds for management decisions.

Genetic substructure and admixture as important factors in linkage disequilibrium-based estimation of effective number of breeders in recovering wildlife populations.

The number of effective breeders (Nb ) and effective population size (Ne ) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium-based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993-2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on Nb and Ne , by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for Ne were larger than Nb and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of Nb and Ne . Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on Nb and Ne .

Human behaviour can trigger large carnivore attacks in developed countries.

The media and scientific literature are increasingly reporting an escalation of large carnivore attacks on humans in North America and Europe. Although rare compared to human fatalities by other wildlife, the media often overplay large carnivore attacks on humans, causing increased fear and negative attitudes towards coexisting with and conserving these species. Although large carnivore populations are generally increasing in developed countries, increased numbers are not solely responsible for the observed rise in the number of attacks by large carnivores. Here we show that an increasing number of people are involved in outdoor activities and, when doing so, some people engage in risk-enhancing behaviour that can increase the probability of a risky encounter and a potential attack. About half of the well-documented reported attacks have involved risk-enhancing human behaviours, the most common of which is leaving children unattended. Our study provides unique insight into the causes, and as a result the prevention, of large carnivore attacks on people. Prevention and information that can encourage appropriate human behaviour when sharing the landscape with large carnivores are of paramount importance to reduce both potentially fatal human-carnivore encounters and their consequences to large carnivores.

What is the animal doing? Tools for exploring behavioural structure in animal movements.

Movement data provide a window - often our only window - into the cognitive, social and biological processes that underlie the behavioural ecology of animals in the wild. Robust methods for identifying and interpreting distinct modes of movement behaviour are of great importance, but complicated by the fact that movement data are complex, multivariate and dependent. Many different approaches to exploratory analysis of movement have been developed to answer similar questions, and practitioners are often at a loss for how to choose an appropriate tool for a specific question. We apply and compare four methodological approaches: first passage time (FPT), Bayesian partitioning of Markov models (BPMM), behavioural change point analysis (BCPA) and a fitted multistate random walk (MRW) to three simulated tracks and two animal trajectories - a sea lamprey (Petromyzon marinus) tracked for 12 h and a wolf (Canis lupus) tracked for 1 year. The simulations - in which, respectively, velocity, tortuosity and spatial bias change - highlight the sensitivity of all methods to model misspecification. Methods that do not account for autocorrelation in the movement variables lead to spurious change points, while methods that do not account for spatial bias completely miss changes in orientation. When applied to the animal data, the methods broadly agree on the structure of the movement behaviours. Important discrepancies, however, reflect differences in the assumptions and nature of the outputs. Important trade-offs are between the strength of the a priori assumptions (low in BCPA, high in MRW), complexity of output (high in the BCPA, low in the BPMM and MRW) and explanatory potential (highest in the MRW). The animal track analysis suggests some general principles for the exploratory analysis of movement data, including ways to exploit the strengths of the various methods. We argue for close and detailed exploratory analysis of movement before fitting complex movement models.

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

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

Evidence of rapid change in genetic structure and diversity during range expansion in a recovering large terrestrial carnivore.

Recovery of natural populations occurs often with simultaneous or subsequent range expansions. According to population genetic theory, genetic structuring emerges at the expansion front together with decreasing genetic diversity, owing to multiple founder events. Thereupon, as the expansion proceeds and connectivity among populations is established, homogenization and a resurgence of genetic diversity are to be expected. Few studies have used a fine temporal scale combined with genetic sampling to track range expansions as they proceed in wild animal populations. As a natural experiment, the historical eradication of large terrestrial carnivores followed by their recovery and recolonization may facilitate empirical tests of these ideas. Here, using brown bear (Ursus arctos) as model species, we tested predictions from genetic theory of range expansion. Individuals from all over Finland were genotyped for every year between 1996 and 2010 using 12 validated autosomal microsatellite markers. A latitudinal shift of about 110 km was observed in the distribution and delineation of genetic clusters during this period. As the range expansion proceeded, we found, as theory predicts, that the degree of genetic structure decreased, and that both genetic variation and admixture increased. The genetic consequences of range expansions may first be detected after multiple generations, but we found major changes in genetic composition after just 1.5 generations, accompanied by population growth and increased migration. These rapid genetic changes suggest an ongoing concerted action of geographical and demographic expansion combined with substantial immigration of bears from Russia during the recovery of brown bears within the large ecosystem of northern Europe.

Y chromosome haplotype distribution of brown bears (Ursus arctos) in Northern Europe provides insight into population history and recovery.

High-resolution, male-inherited Y-chromosomal markers are a useful tool for population genetic analyses of wildlife species, but to date have only been applied in this context to relatively few species besides humans. Using nine Y-chromosomal STRs and three Y-chromosomal single nucleotide polymorphism markers (Y-SNPs), we studied whether male gene flow was important for the recent recovery of the brown bear (Ursus arctos) in Northern Europe, where the species declined dramatically in numbers and geographical distribution during the last centuries but is expanding now. We found 36 haplotypes in 443 male extant brown bears from Sweden, Norway, Finland and northwestern Russia. In 14 individuals from southern Norway from 1780 to 1920, we found two Y chromosome haplotypes present in the extant population as well as four Y chromosome haplotypes not present among the modern samples. Our results suggested major differences in genetic connectivity, diversity and structure between the eastern and the western populations in Northern Europe. In the west, our results indicated that the recovered population originated from only four male lineages, displaying pronounced spatial structuring suggestive of large-scale population size increase under limited male gene flow within the western subpopulation. In the east, we found a contrasting pattern, with high haplotype diversity and admixture. This first population genetic analysis of male brown bears shows conclusively that male gene flow was not the main force of population recovery.

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

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

Admixture and gene flow from Russia in the recovering Northern European brown bear (Ursus arctos).

Large carnivores were persecuted to near extinction during the last centuries, but have now recovered in some countries. It has been proposed earlier that the recovery of the Northern European brown bear is supported by migration from Russia. We tested this hypothesis by obtaining for the first time continuous sampling of the whole Finnish bear population, which is located centrally between the Russian and Scandinavian bear populations. The Finnish population is assumed to experience high gene flow from Russian Karelia. If so, no or a low degree of genetic differentiation between Finnish and Russian bears could be expected. We have genotyped bears extensively from all over Finland using 12 validated microsatellite markers and compared their genetic composition to bears from Russian Karelia, Sweden, and Norway. Our fine masked investigation identified two overlapping genetic clusters structured by isolation-by-distance in Finland (pairwise FST = 0.025). One cluster included Russian bears, and migration analyses showed a high number of migrants from Russia into Finland, providing evidence of eastern gene flow as an important driver during recovery. In comparison, both clusters excluded bears from Sweden and Norway, and we found no migrants from Finland in either country, indicating that eastern gene flow was probably not important for the population recovery in Scandinavia. Our analyses on different spatial scales suggest a continuous bear population in Finland and Russian Karelia, separated from Scandinavia.

Long-range gene flow and the effects of climatic and ecological factors on genetic structuring in a large, solitary carnivore: the Eurasian lynx.

Due to their high mobility, large terrestrial predators are potentially capable of maintaining high connectivity, and therefore low genetic differentiation among populations. However, previous molecular studies have provided contradictory findings in relation to this. To elucidate patterns of genetic structure in large carnivores, we studied the genetic variability of the Eurasian lynx, Lynx lynx throughout north-eastern Europe using microsatellite, mitochondrial DNA control region and Y chromosome-linked markers. Using SAMOVA we found analogous patterns of genetic structure based on both mtDNA and microsatellites, which coincided with a relatively little evidence for male-biased dispersal. No polymorphism for the cytochrome b and ATP6 mtDNA genes and Y chromosome-linked markers were found. Lynx inhabiting a large area encompassing Finland, the Baltic countries and western Russia formed a single genetic unit, while some marginal populations were clearly divergent from others. The existence of a migration corridor was suggested to correspond with distribution of continuous forest cover. The lowest variability (in both markers) was found in lynx from Norway and Bialowieza Primeval Forest (BPF), which coincided with a recent demographic bottleneck (Norway) or high habitat fragmentation (BPF). The Carpathian population, being monomorphic for the control region, showed relatively high microsatellite diversity, suggesting the effect of a past bottleneck (e.g. during Last Glacial Maximum) on its present genetic composition. Genetic structuring for the mtDNA control region was best explained by latitude and snow cover depth. Microsatellite structuring correlated with the lynx's main prey, especially the proportion of red deer (Cervus elaphus) in its diet. Eurasian lynx are capable of maintaining panmictic populations across eastern Europe unless they are severely limited by habitat continuity or a reduction in numbers. Different correlations of mtDNA and microsatellite population divergence patterns with climatic and ecological factors may suggest separate selective pressures acting on males and females in this solitary carnivore.

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

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

Free-ranging Eurasian lynx (Lynx lynx) as host of Toxoplasma gondii in Finland.

We investigated the presence of Toxoplasma gondii infections in Eurasian lynx (Lynx lynx) in Finland by analyzing samples from 337 lynx that were legally hunted during the 2010-2011 season and by performing a retrospective nationwide database search of postmortem toxoplasmosis diagnoses in this species. We detected specific anti-T. gondii IgG antibodies in 290 (86.1%) of the 337 lynx. The method used was a direct agglutination test, and samples positive at the used dilution 1:40 were defined as antibody positive. Older lynx had 14.3 times higher odds of being antibody-positive than did lynx of the presumed age of 7-10 mo, and lynx weighing >15 kg had 16.7 times higher odds of being antibody positive than did those ≤ 15 kg. Lynx from the southwest were more often antibody positive, with an odds ratio 6.3, than lynx from the northeast. None of the 332 fecal samples available was positive for the presence of T. gondii-like oocysts with a quantitative MgSO4 flotation technique, and none of the 167 free-ranging Eurasian lynx examined postmortem by veterinary pathologists from January 2000 to May 2010 had died from toxoplasmosis. Although Finnish lynx were confirmed to commonly encounter T. gondii, we found no evidence of an ongoing contribution to the environmental oocyst burden nor of the lynx dying from the infection.

Endoparasites of the Eurasian lynx (Lynx lynx) in Finland.

We sampled 339 fecal samples, 296 intestines, and 82 lungs from 371 lynx hunted during the 2010-2011 season in Finland. The fecal samples were analyzed for endoparasites by a quantitative flotation method, and helminths from intestines were studied morphologically, while lungs were investigated for pulmonary parasites. From fecal samples, eggs and oocysts of at least 6 different endoparasite species were identified, with a mean of 1.5 (range 0-4) parasite species per host. In the intestines, at least 4 different helminth species were found, with the mean of 2.0 (range 1-4) species per infected host. The prevalence of eggs in feces and the prevalence of worms in intestines were 71% and 93% for Toxocara cati , 29% and 68% for Taenia spp., and 5% and 2% for Diphyllobothrium sp., respectively. Only eggs were detected for Capillaria sp. (46%) and Uncinaria sp. (0.6%) nematodes, and only adults were detected for Mesocestoides sp. cestodes (0.3%). Significant positive correlations were evident between the number of T. cati (r = 0.664; P = 0.01) and Diphyllobothrium sp. (r = 0.645; P = 0.01) eggs per gram of feces and adult worms detected in intestine. In addition to the metazoan parasites, protozoan Isospora sp. oocysts were also found (0.6%). Pulmonary samples were all negative for parasites. These data demonstrate that lynx commonly harbor various endoparasites, some of which are zoonotic.