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.

Wolf attacks on hunting dogs are negatively related to prey abundance in Finland: an analysis at the wolf territory level

Attacks by wolves (Canis lupus) on dogs (C. familiaris) presumably are motivated both by preying and elimination of potential competitors. Regardless of these alternative motivations in wolves, the risk of attacks might be higher when the density of primary prey is low. We examined how many dogs do territorial wolves in Finland kill in relation to the population density of the most abundant ungulates, moose (Alces alces), white-tailed deer (Odocoileus virginianus), and roe deer (Capreolus capreolus). Most attacks by wolves on dogs take place in hunting with dogs. The number of wolf-killed dogs was in highly significant negative relationship to the population density of white-tailed deer and to total ungulate biomass per unit area which is largely determined by the density of white-tailed deer. Our results indicate that abundant wild prey would decrease the risk at which wolves attack dogs. On the other side of the coin prevail two hard facts which wildlife managers had to take a notice. White-tailed deer, although a potential mitigator of wolf–human conflict, is an alien species and a partner in > 6000 traffic collisions annually in Finland. One factor that seemed to increase the risk of wolf attacks on dogs is the low ungulate density in regions where moose is the only remarkable ungulate prey. Higher moose densities could decrease the risk of attacks, but on the other hand, higher densities could increase the risk of serious traffic collisions and browsing damages in forests.

Assessment of the residential Finnish wolf population combines DNA captures, citizen observations and mortality data using a Bayesian state-space model

Assessment of the Finnish wolf population relies on multiple sources of information. This paper describes how Bayesian inference is used to pool the information contained in different data sets (point observations, non-invasive genetics, known mortalities) for the estimation of the number of territories occupied by family packs and pairs. The output of the assessment model is a joint probability distribution, which describes current knowledge about the number of wolves within each territory. The joint distribution can be used to derive probability distributions for the total number of wolves in all territories and for the pack status within each territory. Most of the data set comprises of both voluntary-provided point observations and DNA samples provided by volunteers and research personnel. The new method reduces the role of expert judgement in the assessment process, providing increased transparency and repeatability.

Spatial Memory Drives Foraging Strategies of Wolves, but in Highly Individual Ways

The ability of wild animals to navigate and survive in complex and dynamic environments depends on their ability to store relevant information and place it in a spatial context. Despite the centrality of spatial memory, and given our increasing ability to observe animal movements in the wild, it is perhaps surprising how difficult it is to demonstrate spatial memory empirically. We present a cognitive analysis of movements of several wolves (Canis lupus) in Finland during a summer period of intensive hunting and den-centered pup-rearing. We tracked several wolves in the field by visiting nearly all GPS locations outside the den, allowing us to identify the species, location and timing of nearly all prey killed. We then developed a model that assigns a spatially explicit value based on memory of predation success and territorial marking. The framework allows for estimation of multiple cognitive parameters, including temporal and spatial scales of memory. For most wolves, fitted memory-based models outperformed null models by 20 to 50% at predicting locations where wolves chose to forage. However, there was a high amount of individual variability among wolves in strength and even direction of responses to experiences. Some wolves tended to return to locations with recent predation success—following a strategy of foraging site fidelity—while others appeared to prefer a site switching strategy. These differences are possibly explained by variability in pack sizes, numbers of pups, and features of the territories. Our analysis points toward concrete strategies for incorporating spatial memory in the study of animal movements while providing nuanced insights into the behavioral strategies of individual predators.

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.

Genetic signature of immigrants and their effect on genetic diversity in the recently established Scandinavian wolf population

Transboundary connectivity is a key component when conserving and managing animal species that require large areas to maintain viable population sizes. Wolves Canis lupus recolonized the Scandinavian Peninsula in the early 1980s. The population is geographically isolated and relies on immigration to not lose genetic diversity and to maintain long term viability. In this study we address (1) to what extent the genetic diversity among Scandinavian wolves has recovered during 30 years since its foundation in relation to the source populations in Finland and Russia, (2) if immigration has occurred from both Finland and Russia, two countries with very different wolf management and legislative obligations to ensure long term viability of wolves, and (3) if immigrants can be assumed to be unrelated. Using 26 microsatellite loci we found that although the genetic diversity increased among Scandinavian wolves (n = 143), it has not reached the same levels found in Finland (n = 25) or in Russia (n = 19). Low genetic differentiation between Finnish and Russian wolves, complicated our ability to determine the origin of immigrant wolves (n = 20) with respect to nationality. Nevertheless, based on differences in allelic richness and private allelic richness between the two countries, results supported the occurrence of immigration from both countries. A priori assumptions that immigrants are unrelated is non-advisable, since 5.8% of the pair-wise analyzed immigrants were closely related. To maintain long term viability of wolves in Northern Europe, this study highlights the potential and need for management actions that facilitate transboundary dispersal.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07761-5.

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.

Does public information about wolf (Canis lupus) movements decrease wolf attacks on hunting dogs (C. familiaris)?

The threat that wolves (Canis lupus) pose to hunting dogs is one reason why Finnish hunters have negative attitudes towards wolves and one of the potential motivations for the illegal killing of wolves. During 2010–2017, wolves killed an average of 38 dogs (range 24–50) per year in Finland. Most of the attacks (91%) were directed at hunting dogs during the hunting season. To decrease the risk of attacks, the last seven positions (one position per hour) of GPS-collared wolves were accessible to the public with a 5 × 5 km resolution during the hunting seasons (from August 20th to February 28th) of 2013/2014 (from September 2nd onwards), 2015/2016, 2016/2017 and 2017/2018. The link was visited more than 1 million times in 3 of the 4 seasons. Fatal attacks on dogs occurred on 17% of the days during the hunting seasons of our study (n = 760 days). Both the attacks and visits peaked in September–November, which is the primary hunting season in Finland. According to the general linear model, the number of daily visits to the website was higher on days when fatal attacks occurred than on other days. Additionally, season and the number of days passed from the first day of the season were significantly related to the daily visits. Visits were temporally auto-correlated, and the parameter values in the model where the dependent variable was the number of visits on the next day were only slightly different from those in the first model. A two-way interaction between season and attack existed, and the least squares means were significantly different in 2017/2018. The change in daily visits between consecutive days was related only to the number of days from the beginning of the season. We examined whether this kind of service decreased dog attacks by wolves. Wolf attacks were recorded in 32% of the wolf territories, where at least one wolf had been collared (n = 22). However, within the territories without any GPS-collared wolves, the proportion of territories with wolf attack(s) was significantly higher than those elsewhere (50%, n = 48). Although public information decreased the risk of attacks, it did not completely protect dogs from wolf attacks and may in some cases increase the risk of illegally killing wolves. The most remarkable benefit of this kind of service to the conservation of the wolf population might be the message to the public that management is not overlooking hunters’ concerns about wolf attacks on their dogs.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

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.

Electronic supplementary material

The online version of this article (10.1186/s12898-017-0154-8) contains supplementary material, which is available to authorized users.

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 .

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.

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.

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.

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.

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 Białowieża 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.

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.

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.

Limited gene flow among brown bear populations in far Northern Europe? Genetic analysis of the east-west border population in the Pasvik Valley

Noninvasively collected genetic data can be used to analyse large-scale connectivity patterns among populations of large predators without disturbing them, which may contribute to unravel the species' roles in natural ecosystems and their requirements for long-term survival. The demographic history of brown bears (Ursus arctos) in Northern Europe indicates several extinction and recolonization events, but little is known about present gene flow between populations of the east and west. We used 12 validated microsatellite markers to analyse 1580 hair and faecal samples collected during six consecutive years (2005-2010) in the Pasvik Valley at 70°N on the border of Norway, Finland and Russia. Our results showed an overall high correlation between the annual estimates of population size (N(c) ), density (D), effective size (N(e) ) and N(e) /N(c) ratio. Furthermore, we observed a genetic heterogeneity of ∼0.8 and high N(e) /N(c) ratios of ∼0.6, which suggests gene flow from the east. Thus, we expanded the population genetic study to include Karelia (Russia, Finland), Västerbotten (Sweden) and Troms (Norway) (477 individuals in total) and detected four distinct genetic clusters with low migration rates among the regions. More specifically, we found that differentiation was relatively low from the Pasvik Valley towards the south and east, whereas, in contrast, moderately high pairwise F(ST) values (0.91-0.12) were detected between the east and the west. Our results indicate ongoing limits to gene flow towards the west, and the existence of barriers to migration between eastern and western brown bear populations in Northern Europe.

A forensic DNA profiling system for Northern European brown bears (Ursus arctos)

A set of 13 dinucleotide STR loci (G1A, G10B, G1D, G10L, MU05, MU09, MU10, MU15, MU23, MU26, MU50, MU51, MU59) were selected as candidate markers for a DNA forensic profiling system for Northern European brown bear (Ursus arctos). We present results from validation of the markers with respect to their sensitivity, species specificity and performance (precision, heterozygote balance and stutter ratios). All STRs were amplified with 0.6ng template input, and there were no false bear genotypes in the cross-species amplification tests. The validation experiments showed that stutter ratios and heterozygote balance was more pronounced than in the tetranucleotide loci used in human forensics. The elevated ratios of stutter and heterozygote balance at the loci validated indicate that these dinucleotide STRs are not well suited for interpretation of individual genotypes in mixtures. Based on the results from the experimental validations we discuss the challenges related to genotyping dinucleotide STRs in single source samples. Sequence studies of common alleles showed that, in general, the size variation of alleles corresponded with the variation in number of repeats. The samples characterized by sequence analysis may serve as standard DNA samples for inter laboratory calibration. A total of 479 individuals from eight Northern European brown bear populations were analyzed in the 13 candidate STRs. Locus MU26 was excluded as a putative forensic marker after revealing large deviations from expected heterozygosity likely to be caused by null-alleles at this locus. The remaining STRs did not reveal significant deviations from Hardy-Weinberg equilibrium expectations except for loci G10B and MU10 that showed significant deviations in one population each, respectively. There were 9 pairwise locus comparisons that showed significant deviation from linkage equilibrium in one or two out of the eight populations. Substantial genetic differentiation was detected in some of the pairwise population comparisons and the average estimate of population substructure (F(ST)) was 0.09. The average estimate of inbreeding (F(IS)) was 0.005. Accounting for population substructure and inbreeding the total average probability of identity in each of the eight populations was lower than 1.1×10(-9) and the total average probability of sibling identity was lower than 1.3×10(-4). The magnitude of these measurements indicates that if applying these twelve STRs in a DNA profiling system this would provide individual specific evidence.

Summer movements, predation and habitat use of wolves in human modified boreal forests

Grey wolves (Canis lupus), formerly extirpated in Finland, have recolonized a boreal forest environment that has been significantly altered by humans, becoming a patchwork of managed forests and clearcuts crisscrossed by roads, power lines, and railways. Little is known about how the wolves utilize this impacted ecosystem, especially during the pup-rearing summer months. We tracked two wolves instrumented with GPS collars transmitting at 30-min intervals during two summers in eastern Finland, visiting all locations in the field, identifying prey items and classifying movement behaviors. We analyzed preference and avoidance of habitat types, linear elements and habitat edges, and tested the generality of our results against lower resolution summer movements of 23 other collared wolves. Wolves tended to show a strong preference for transitional woodlands (mostly harvested clearcuts) and mixed forests over coniferous forests and to use forest roads and low use linear elements to facilitate movement. The high density of primary roads in one wolf's territory led to more constrained use of the home territory compared to the wolf with fewer roads, suggesting avoidance of humans; however, there did not appear to be large differences on the hunting success or the success of pup rearing for the two packs. In total, 90 kills were identified, almost entirely moose (Alces alces) and reindeer (Rangifer tarandus sspp.) calves of which a large proportion were killed in transitional woodlands. Generally, wolves displayed a high level of adaptability, successfully exploiting direct and indirect human-derived modifications to the boreal forest environment.