Non-invasive genetic monitoring involving citizen science enables reconstruction of current pack dynamics in a re-establishing wolf population
© The Author(s) 2017
Received: 6 November 2017
Accepted: 7 December 2017
Published: 19 December 2017
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.
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.
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.
Large carnivores are important for ecosystem sustainability and diversity, but in many places they have declined drastically due to anthropogenic pressure [1, 2]. Carnivores frequently come into conflict with human interests as they prey on domestic and game animals and are often perceived as a threat to human safety [2–4]. As a consequence, many large carnivores have been eradicated from wide geographic ranges to date. However, strict conservation measures have recently enabled their return in some areas, including Europe  and North America . Re-establishment of such emotive species to urban and agricultural landscapes often triggers public concern as local people are no longer accustomed to the presence of predators. Indeed, a better understanding of the re-establishing species’ biology has been shown to be essential to minimizing conflict and fostering sustainable coexistence of humans and wild carnivores [7, 8].
One of the most controversial and conflict-prone carnivores is the grey wolf (Canis lupus). Following near extirpation from most of Europe and North America, the species is now resettling in some parts of its previous home ranges [5, 9]. One such area is southwestern Finland, which wolves have recently recolonized as breeders after an absence of ca. 120 years [10, 11]. This area has a tragic history with wolves that affects the social environment for their conservation and management: in 1879–1882, > 20 children were killed in putative wolf attacks and local wolves were subsequently eradicated by professional hunters . Rebounding started in 2006 when the first family pack was confirmed . These historical events have contributed to a pronounced conflict over the re-establishment of wolves between enforced, top–down conservation aims (e.g. strict protection under the European Union’s Habitats Directive) and local inhabitants’ concerns for the safety of their families, livestock and pets [8, 11, 13]. This cultural context, together with incidences of wolf sightings close to human settlements , has provoked strong calls for the complete eradication of wolves from the area [8, 11].
A key element in conflicts over the management of carnivores is typically a lack of information or disagreement regarding basic biological facts related to the carnivores. In the case of wolves in SW Finland, there have been several debates about the number of wolf packs and individuals, in which the opinions of different stakeholders vary and some local residents mistrust the official estimates . Reliable identification of individual wolves in this area could help to resolve the dispute over numbers of wolves, providing an unbiased basis for conservation and management, but this is not possible with the common census approach of track-based counts. To this end, genetic methods using non-invasively collected samples that allow individual identification but do not require capturing, handling or direct physical contact with the animals have recently provided a useful tool for studying the population size and structure of rare and elusive species [15–18] including various wolf populations [19–21]. However, few genetic monitoring studies have thus far taken advantage of the insight that such data can provide on management-relevant reproductive dynamics [e.g. 21, 22]. In a re-establishing population, robust characterization of relatedness patterns allows (1) reconstruction of the colonisation history at the individual level, (2) identifying breeding pairs, and (3) inferring likely future breeding events. Understanding these dynamics provides a basis for targeted management actions , both immediate—e.g. identifying problematic groups or individuals for control—and anticipatory, such as predicting changes in pack configurations based on past responses to similar pressures. Investigating these patterns of population change is particularly valuable in a recently-established and currently expanding population such as the SW Finland wolves, where pack development can be followed in almost real time and predictions of change tested in the context of e.g. harvesting of identified breeding individuals.
Particularly when data collection is based on citizen science—in which interested non-experts contribute to data collection—non-invasive genetic sampling provides an efficient alternative to deliver detailed monitoring of local individuals that does not rely on field observations of these elusive, wide-ranging animals and at the same time enables inclusion of different stakeholders in knowledge generation. However, non-invasive sampling also poses challenges for conventional genetic pedigree reconstruction techniques [23–26]. The genetic material can be of poor quality and prone to genotyping errors, and many parents of identified individuals may not be sampled but without a priori knowledge of which or how many individuals. Further, individuals may be both parents and/or offspring, but there is little prior cohort information to constrain these relationships. Pedigree reconstruction in a Bayesian framework provides a versatile tool to overcome these constraints in a single analytical framework, reducing potential bias caused by assumptions around unknown variables and allowing phenotypic data to be considered together with genetic data in estimating parentages [23, 25, 27].
We used a combination of non-invasive and opportunistic samples to obtain genetic data from which to characterize the relatedness structure of the re-establishing wolf population in SW Finland in a Bayesian framework. To enhance the social acceptance of our outcomes in this contentious debate, we employed a collaborative sampling framework that is primarily based on a joint effort between local people, wildlife administration and researchers. Thus, stakeholders representing a range of interests contributed to knowledge generation, with a view to decreasing conflict over knowledge ownership and bias. Using an established genotyping panel of 17 microsatellite markers , we first identified 81 genetically unique individuals represented in the samples. We then reconstructed a multi-generational pedigree to determine reproductive pairs, identify family packs and characterize their habitat areas. Lastly, supplementing the family structure with temporally explicit observational data, we constructed a population history of wolf re-colonisation in this human-populated area. These data are already being used as a basis for local population management and the entire monitoring framework is being expanded over a larger area.
This study used genetic data from non-invasive samples collected in 2013–2016 supplemented with genetic information available from wolves that were killed or trapped in the area in 2008–2015. Our main source material for DNA extraction was non-invasively collected fecal samples, which are relatively easy to find either by following wolf tracks, especially during winter when snow-cover makes tracks more visible, or near wolf kill sites or alongside forest roads, which wolves use when moving across the territory [31, 32]. Sample collection was initially organized and performed in collaboration with local wildlife agencies in the study area (Varsinais-Suomen and Satakunnan Riistakeskus, the Wildlife Agencies in Varsinais-Suomi and Satakunta), the Regional Wildlife Council (alueellinen riistaneuvosto) of Varsinais-Suomi, and the Natural Resources Institute Finland (Luke). Potential collaborators were invited to two meetings, and a voluntary scheme was established to opportunistically collect fecal samples. Sample collectors were mainly voluntary large carnivore personnel (hunters) that have been trained to identify wolf tracks and signs [see 10], with additional samples collected by other local hunters and nature enthusiasts that make recreational use of the area.
Sample collectors were provided with instructions and collection kits to ensure consistency in the sampling protocol and associated data. Whole fecal lobes were collected using a plastic bag and the following information was recorded: name and e-mail address of the collector, date, coordinates of the sampling location or exact written description about the sampling site, an estimate of how many days the sample may have been in the field and how many tracks were seen around the sample (if applicable). Additional information was often added by sample collectors, for example to note uncertainty that the sample was produced by a wolf. To minimize DNA degradation, samples were frozen upon collection and subsequently transported in cool bags to the laboratory, where they were stored at − 20 °C. In order to motivate sample providers, each participant received frequent feedback about the progress of the respective samples. Additionally, two progress reports were sent during and at the end of the sample collection season. From summer 2014 to spring 2016, new samples were collected by Natural Resources Institute Finland (Luke) assisted by the volunteers; the sampling scheme was not however advertised as widely as it was in 2013–2014 due to funding uncertainty. In total, this collection protocol yielded 516 scat samples, 368 of which yielded at least partial genotypes, covering adults and litters produced in the breeding seasons 2013–2015.
Sample sizes at different stages of the study protocol. DNA amplification was attempted on all samples collected during the study period (2013–2015); the number of genotypes obtained and used in the analysis is lower as many samples did not yield a sufficient number of markers to assign the sample to an individual
Scats, study period
Tissue/secretions, study period
Scats, before study period
Tissue secretions, before study period
Samples from a range of sources (tissue, feces, hair, urine and secretion samples) were processed as a part of ongoing monitoring procedures over several years, and genotyping was conducted by several different researchers. Therefore, the DNA extractions and genotyping protocols were modified during the course of the research. Below, a summary of these procedures is provided; full details can be found in Additional file 1: Appendix 1.
DNA was extracted using commercially available kits with minor modifications of the manufacturer’s protocols. The kit and protocol varied depending on the DNA source material (details in Additional file 1: Appendix 1). Negative controls were added in all extraction batches where samples were expected to have low quality or quantity DNA. The quantity and quality of all extracted DNA were quantified using a NanoDrop ND-1000 spectrophotometer.
All samples were genotyped at 17 established canid microsatellite loci (Additional file 1: Table S1) [33, 34] that have been previously used in wolf population genetic studies in Finland . Separate amplification protocols were used for tissue and other sample types.
Amplification was conducted in 3–5 multiplex PCRs depending on the DNA source (Additional file 1: Appendix 1) in 10–12 µl volumes containing 1–3 µl template DNA, 1 × Qiagen multiplex PCR master mix (QMP) and 0.1–0.4 µM of each primer (Additional file 1: Table S1). The thermal cycling program was: initial denaturation at 95 °C for 15 min, 35–40 amplification cycles of denaturation at 94 °C for 30 s, annealing at 60 °C for 90 s and extension at 72 °C for 60 s, and a final extension step at 60 °C for 10–30 min.
Individual multi-locus genotypes were generated using fragment analysis of the pooled fluorescently labeled PCR products. The allele sizes of the microsatellite markers were determined by capillary electrophoresis using an ABI PRISM 3130xl Genetic Analyzer-sequencer (Applied Biosystems) and individual genotypes were called using the genotyping program GeneMarker v.2.2.0 (Softgenetics, Inc.). After an automated initial allele call, all alleles were manually checked and confirmed before a final genotype was assigned.
Genotyping was carried out a number of times on each sample, as recommended to ensure reliability of genotypes when analyzing very dilute and/or degraded DNA such as that extracted from feces [35–37]. Each sample was first amplified in three independent replicates; at each locus, heterozygous genotypes were considered reliable if both alleles were observed in at least two replicates and for homozygous genotypes if the same genotype was observed in all three [38, 39]. If any locus did not meet these requirements, a second series of three replicates was performed. If any locus yielded no reliable genotype, that locus was treated in the analysis as missing data.
The sex of individuals identified through microsatellite genotyping was determined genetically by amplifying sex specific fragments from DBY intron 7 and DBX intron 6. Primers used for amplification are canid specific and amplify two differently sized fragments in males and one in females . PCR amplification consisted of one multiplexed 10 µl reaction with primer concentrations of 0.1–0.2 uM, 3 µl template DNA and 1 × Qiagen multiplex PCR master mix (QMP) primer (Additional file 1: Table S1). The PCR profile had an initial denaturation at 95 °C for 15 min followed by a touchdown profile with 20 cycles of denaturation at 94 °C for 30 s, annealing at 60 °C for 90 s with a 0.2 °C decrease in temperature after each cycle and extension at 72 °C for 60 s, followed by 20 more cycles with annealing at 55 °C and a final extension step at 60 °C for 30 min. Amplified fragments were visualized on a 2% agarose gel.
Individual identification analysis was done with the software Gimlet v.1.3.3 and Cervus 3.0.7 . Matching multilocus genotypes were assumed to represent the same individual. Due to the likelihood of genotyping errors due to the low quality and quantity of the non-invasively collected DNA, we allowed maximum two mismatching alleles (one at each of two loci) between genotypes for them to be considered to have originated from the same individual. All assembled consensus individual genotypes were successfully genotyped at a minimum of 11 loci (median 17 loci).
Molecular checks for potential dog samples were necessary to ascertain that the collected fecal samples were indeed from wolves as the droppings and tracks of wolves and domestic dogs can look very similar in the field. Dog fecal samples were identified against the dog reference samples using either factorial correspondence analysis (FCA) in the software Genetix  or NewHybrids  and excluded from further analysis.
The pedigree was estimated using the MasterBayes package in R, which implements a Monte Carlo Markov chain (MCMC) sampling approach to simultaneously estimate the pedigree alongside other population-level parameters in a Bayesian framework . Opportunistically or non-invasively collected genetic data can be challenging for constructing pedigrees [23–25]. Several precise pedigree reconstruction techniques are able to overcome these data limitations by making assumptions on unknown variables, but this can be problematic as it risks propagating errors into the parentage assignments [24, 26]. Moreover, resolving the direction of parent–offspring relationships in a multigenerational pedigree is an unsolved challenge in for exact pedigree reconstruction methods . The Bayesian approach provides a robust solution to improve reliability of parentage assignment in incompletely sampled populations using low-quality DNA sources [23–25]. Unknown parameters such as the unsampled population size and genotyping error rate are modelled including uncertainty, avoiding assumptions and thus reducing bias in parentage assignments, and the iterative sampling of pedigree configurations enables identification of the most likely directions of relationships between individuals without demographic information [23–25].
Allele frequencies were taken from a previously published data-set covering the entire range of the Finnish wolf population  (Additional file 1, n = 185 individuals) to avoid skewing the likelihoods of allele sharing by descent in the small, potentially closely related population of the current study. Alleles in the SW Finland sub-population that were not present in the Finland-wide sample were added manually at low frequency. All pedigrees treated genotyping error rate as two distinct components, drop-out rate (non-amplification, E1) and stochastic error (genotype wrongly scored, E2) . All pedigrees used the same priors on unsampled population size, broad lognormal distributions with means in the likely real parameter range (for males µ = log(7), σ = 0.5; for females, µ = log(2), σ = 0.5), and all MCMC chains were thinned every 2 iterations. Individuals that mismatched offspring at more than 4 loci were excluded as potential parents; the main pedigree was robust to a range of mismatch thresholds. For each pedigree, the MCMC sampling was tuned to achieve the recommended acceptance rate for efficient mixing (30–50%). Individuals known to have died at younger than 2 years of age (the earliest known reproduction of wild wolves ) were excluded as potential parents and the pedigree was checked for legality (no loops) at each iteration. The data was ordered to make later-observed individuals more likely to be fitted as offspring, but no explicit temporal information was available to inform age or cohort structure. To optimize our use of the available data and improve the mixing of the Markov chains, the final pedigree was obtained in several stages: (1) deriving a preliminary pedigree to inform estimation of genotyping errors; (2) estimating error rates; (3) restricting candidate parents to likely individuals; (4) fitting pedigrees using these estimated parameters; (5) averaging out stochasticity in the sampling to obtain a robust main pedigree; (6) assessing the biological sensibility of all assignments.
First, a preliminary pedigree (step 1) was estimated using the consensus individual genotypes with fixed error rates (E1 = 4%, ; E2 = 2.2%, ), run for 20,000 iterations with a burn-in of 2000, and the modal pedigree extracted from the sampled posterior distribution of pedigrees. Parentages with over 90% confidence from this preliminary pedigree were used as the input pedigree for the error rate estimation (step 2) from the per-sample dataset, which thus used information both from repeated samples of individuals and family relationships. Error rates were estimated for each locus. The pedigree was kept fixed, priors on the error rates were broad with means as above (shape parameters of beta distribution: E1, α = 4, β = 100; E2, α = 2, β = 100) and the model was initiated at a high but realistic value for both error types (0.1) to fully explore parameter space. In pedigree estimation (steps 3 and 4), the estimated error rates were used for individuals represented only by scat samples, which made up the majority of both repeat samples and samples from individuals with preliminary pedigree information. For tissue samples, error rates could not be reliably estimated from this data (only 6 individuals had more than one tissue sample), so a constant E1 and E2 from previous work for a range of tissue types were used across the panel . To further reduce the dimensionality of the main pedigree estimation problem and hence increase confidence in distinguishing between parents with numerically similar likelihoods, before fitting the final pedigree, we first identified potential breeding individuals (step 3). This model used fixed, derived error rates from step 2 and was run 20 times, each run using 25,000 iterations after a short burn-in of 1000 to ensure sufficient exploration of probability space. From across all 20 modal pedigrees, all individuals assigned in ≥ 5% of parentages were taken as “possible breeders” in the main pedigree estimation.
For the main pedigree fit (step 4), we restricted potential parents to be only the identified possible breeders and used a chain of 30,000 iterations with a burn-in of 3000. To make the ultimate output pedigree robust to the stochasticity inherent in a random sampling scheme, this main model was run 20 times and the modal pedigree extracted from each run including likelihoods for each assignment. For each offspring, the average likelihood of each parent pair was calculated across these 20 pedigrees (step 5) and the most likely overall pair assigned as parents in the main pedigree. These assignments were robust to lability among poorly resolved relationships. However, due to a lack of accurate temporal information and low diversity at many of the microsatellite loci, some assignments required manual correction post hoc (step 6) according to biological understanding (“biologist corrections” ). These corrections were based on likely birth years inferred from the obtained pedigree, a minimum age of 2 years at reproduction, and wolves’ mating system of only one pair in each territory producing a litter of several cubs . Because reliable age data was only available for very few individuals and could thus not be used to constrain candidate parents, some assignments were found to be temporally impossible (insufficient time for given number of generations) based on birth years (cohort) inferred from sample dates for individuals throughout the pedigree structure. Such anachronisms were manually removed (mothers of 2 offspring, details in Additional file 1). In addition, certain medium-confidence (70–80%) parent pairs were biologically unlikely, appearing in the pedigree as parents to only one individual, whereas all reliably assigned pairs (with at least one offspring assigned to both parents with over 95% likelihood in the final pedigree) were represented by several offspring. In all such cases of uncertain assignments, parentages could be confirmed or offspring could manually be reassigned to a reliably identified breeding pair, in the same territory as the offspring, whose genotypes were consistent with the offspring’s.
Simulation tests of assignment accuracy
To explore constraints on the reliability of assignments, pedigrees were fitted to simulated genotypes using different numbers of microsatellite markers. Genotypes were simulated using the simgenotypes function in MasterBayes for the sampled population, using the configuration of the final pedigree as the “true” pedigree and using the same base allele frequencies, error rates and set of possible parents as above. Genotypes were simulated for panels of 5, 10, 17, 20, 25, 30 and 34 loci, randomly sampled without replacement (within each group of up to 17 loci) from the real markers [following 26]. For each panel size, 10 sets of genotypes were simulated, and for each simulated genotype set, a pedigree was fitted according to steps 3–5 above. This output pedigree was compared to the true (input) pedigree.
Analysis and pedigree interpretation
All analysis was carried out in R . Expected heterozygosity (HE) was calculated using the library pegas , population-level inbreeding based on allele frequencies using the library adegenet , and inbreeding based on pedigree links using the library pedigree . The simulated pedigrees were assessed in terms of the number of confident assignments made, and the number of confident assignments that were correct, i.e. matched the input pedigree. Unless otherwise stated, we term assignments “confident” at over 80% likelihood. All descriptive statistics are given as mean ± standard error.
From the pedigree, we interpreted further details of breeding timing and dispersal routes using knowledge of the biology of the system and the place and time of collection of each sample. Likely breeding years, and hence birth years of new cubs, were inferred from the simultaneous appearance of sibling groups, supported by the estimated ages of killed young wolves (where available). Individuals assigned to two unsampled parents were assumed to be immigrants. All these interpretations are henceforth presented as “likely” scenarios, as we cannot rule out that individuals were present at other times or places but simply avoided detection, which could affect interpretation.
In total, 516 scat samples were collected in the study area by at least 48 named collectors. Of these samples, 10 were identified as non-wolf (dog or other canid). Of the wolf samples, 368 were successfully genotyped at ≥ 10 loci and could be assigned to a known or new individual. Among all samples, including those from other sources, we identified 81 distinct individuals, represented by 1–19 samples each (median 4 samples per individual). Of these, 23% of samples required a second (78 samples) or third (6 samples) triplicate PCR run to obtain a sufficiently complete genotype.
Summary of the number of confident assignments made in the pedigree fit (main pedigree) and following knowledge-based confirmations of low-confidence assignments or adjustments of temporally impossible or biologically unlikely assignments (final pedigree). Complete individual assignments from the main and final pedigrees are given in Additional file 1
Nr. of individuals
Panel-wide HE was 0.58 ± 0.04 (range per locus 0.20–0.78). Although inbreeding appeared high at the population level on the basis of allele frequencies (across all individuals, mean F = 0.15 ± 0.01, range 0.05–0.47), the pedigree structure gave evidence of only one recent inbreeding event between a grandfather and granddaughter (which gives a pedigree-based F = 0.007).
In contrast to Mynämäki, offspring in the Köyliö pack have been produced by only one male, a likely immigrant, paired with two different females (a likely immigrant for 2011–14, followed in 2015 by the male’s granddaughter dispersing from the Mynämäki pack). Three daughters of the first Köyliö pair dispersed over several years to breed in neighbouring territories with inheriting (Mynämäki pack), likely immigrant (Raasepori pack) or adult disperser (Pöytyä pack) males. Thus, in 2014, the three reproductive females in the core area were a mother and two of her daughters (Fig. 2).
The pedigree revealed two occasions of mate switches; in both cases a male’s initial female disappeared from the record but we do not have direct evidence of mortality. It is therefore unclear whether these switches were active mate choice by the male or a response to losing an established partner. There were no instances of new breeding pairs being established between existing pack members. Full details of the interpretation of pack dynamics are given in Additional file 1.
To examine the statistical robustness of our parentage assignments, we repeated the pedigree fitting on a population of simulated genotypes based on the data. This test showed that 17 markers, i.e. the size of the current panel, gave a mean of 131.9 correct and confident assignments (81.4% of 162 possible assignments). Increasing the panel size to 20 slightly improved accuracy to a mean of 141.7 (87.5%) correct, confident assignment, but further increases up to 34 loci made no significant improvement (linear models: change to number of correct assignments with 20 compared to 17 markers, 3.27 ± 1.05, t = 3.125, p = 0.006; effect of number of markers from 20 to 34 on number of correct assignments − 0.06 ± 0.20, t = − 0.32, p = 0.753) (Fig. 4). Hence, while our pedigree-fitting pipeline is not able to completely resolve all pedigree uncertainties, the confidently-resolved relationships in the empirical data are most likely reliable: at all simulated panel sizes, very few confident assignments were made that were not correct (Fig. 4). At the current panel size of 17 loci, a mean of 132.4 parents were confidently assigned, of which 131.9 were correct, giving a false positive rate of only 0.3% among confident assignments in the simulation.
Here, we have described the pack dynamics of a newly-established and expanding wolf population in south-west Finland based on a collaborative genetic sampling framework in which stakeholders representing a range of interests contributed to knowledge generation. We subsequently used this non-invasively collected genetic data for multi-generational pedigree reconstruction, allowing us to characterize reproductive dynamics and individual dispersal patterns in and between the local wolf packs. This ongoing study provides an excellent opportunity to closely follow current changes in pack structure and reproductive decisions in a dynamic population. The reconstructed pedigree shows a closely interrelated population consisting of sibling groups from 7 breeding (alpha) pairs across four established packs in 2009–2015, with individuals dispersing/moving between all the packs. Despite this, there was little direct inbreeding among the range of strategies utilized to obtain and maintain reproductive opportunities. These results illustrate the value of molecular data for elucidating reproductive dynamics and movements of packs and individuals in detail that not achievable through on-the-ground monitoring alone. Our analytical routine maximized confidence in the final pedigree despite the potentially low quality of the primarily faecal DNA used, the collection of field data mainly by volunteers, and the inherent challenges of multigenerational pedigree reconstruction in incompletely sampled populations . While non-invasive genetic monitoring of individuals is a well-established management tool [15, 20, 37, 53, 54], our extension of previous approaches allows us to detail ongoing changes in the size and structure of packs and the population, providing accessible, unbiased information that is much-needed in the fierce debate concerning management of the grey wolf.
Characterizing pack dynamics
We found that all packs in this small recovering population are related through several different pathways. While calculations based on allele frequencies suggested substantial inbreeding, the pedigree structure indicated that the high interrelatedness is largely due to the population being mainly composed of related sibling groups, rather than a result of recent inbreeding events. Despite the newly-established population expanding to four simultaneously reproducing pairs in recent years, the first instance of direct inbreeding was only detected in 2015. Here, the female was the male’s granddaughter, a level of relatedness (r = 0.25) at which wolves have previously been shown to avoid inbreeding [between half-sibs; 55]. The rarity of inbreeding during the study period suggests that inbreeding avoidance occurs in this population, particularly as most potential mates are related and many detected immigrants have become breeders; inbreeding avoidance occurs in other grey wolf populations as well as in other canids [21, 55, 56]. However, it is important to note that we have insufficient information to rule out the possibility that the first breeding individuals in the data (in the first pairs in Mynämäki and Köyliö) could have been relatives. Relatedness among the founders would for example contribute to the high apparent inbreeding suggested by allele frequencies. In addition, in this small population in which only a small proportion of individuals are likely to breed, genetic drift is expected to decrease diversity over time, and opportunities for inbreeding may become more frequent as the population expands . Despite these influences tending to reduce diversity, our results suggest that natural mechanisms such as inbreeding avoidance can at the same time contribute to limiting the potential negative effects of associated inbreeding depression [as in 56], supporting the short-term viability of the SW Finland wolf population. Continued monitoring of this expanding population affords a rare opportunity to track ongoing changes in inbreeding dynamics, which could enable the design of management actions to provide opportunities for outbreeding, such as dispersal corridors, and hence maintain diversity and support reproduction [56, 57].
Potential and established breeders used a broad range of strategies to obtain and maintain reproductive opportunities, with inbreeding apparently the least preferred option after dispersal to existing or new packs, inheriting breeding positions in the natal pack, or (for established breeding males) mate switching. The only instance of inbreeding occurred when a dispersing female joined an established male, whose previous mate disappeared from our records the year before the new litter was recorded. This diversity of breeding strategies and strong inbreeding avoidance mirrors findings from a reintroduced North American wolf population . Overall, females dispersed between known packs more frequently than males, although our pedigree shows that this may not be sufficient to ensure inbreeding avoidance in a densely occupied habitat. In contrast, males that became established as breeders during the study tended to come from unknown natal packs, i.e. likely long-distance dispersers. This pattern of shorter female and longer male dispersal is also seen in other canids [African wild dogs Lycaon pictus, 58]. Together with the observation that there were more novel breeding opportunities for females in our population, this suggests that while dispersal at the local scale may be sufficient to secure a breeding position for females, males must disperse over longer distances. Such immigration is extremely valuable in maintaining genetic diversity and population viability, particularly in small and interrelated populations [56, 57]. Dispersal is common in other parts of the expanding Finnish wolf population: in a 6-year study, over half of tracked non-alpha individuals dispersed from their natal pack, but with only a minority of dispersers (33%) reproducing in a new territory . Moreover, our female-biased local dispersal but male-biased immigration suggests that harvesting breeding males may limit population expansion more than removing females, which could be more easily replaced locally. In any case, continued population expansion does appear to be possible: two new packs were established (Pöytyä and Raasepori) during the study, and suitable breeding territories may yet remain unfilled. Such insights into current reproductive dynamics illustrates the value of applying genetic monitoring approaches to a newly-establishing population, in contrast to previous work focusing on established populations and historical patterns.
These insights into the behaviour of SW Finland wolves are underpinned by the robustness of our pedigree: 74% of all parentage assignments in the final pedigree were made with high confidence, even in this interrelated population with relatively few founders and consequently limited genetic diversity (HE of 0.58, compared to 3-year temporal estimates across the Finnish population during 1995–2009, HE > 0.67, ). Of the low-confidence assignments, 15 could be confirmed based on biological understanding of the system, but conversely, such interpretation suggested that 4 confident assignments and 16 low-confidence assignments were biologically unlikely. The majority of these manual changes were made to parents that were most likely too young or too many generations back, once likely birth years had been interpreted from the pedigree configuration of sibling groups and their detection dates. Reliable age data would contribute substantially to increasing confidence in our pedigree structure, but this is a challenge to obtain from non-invasive samples. Due to the inherent opportunism in production and detection of faecal material, an individual may be in the population long before or after its first or last detection, even when the sample collection protocol is systematic. Therefore, we did not have sufficiently accurate chronological data for most individuals (e.g. cohort or death date) with which to explicitly constrain potential parent–offspring relationships. It would be valuable to investigate whether faecal lobe size could be used to indicate the age of individual wolves, as has been done in non-invasive genetic studies of western gorillas .
Despite these practical challenges, a non-invasive sampling scheme involving citizen scientists provides a productive tool to underpin management decisions for contentious populations such as the SW Finland wolves. Enabling local people from different backgrounds to contribute concretely to data collection can mitigate the perceived conflict over who holds ownership of information [8, 11, 13]. Moreover, direct involvement with monitoring can increase the public’s understanding of wolf biology, which contributes to mitigating conflicts with the predator [7, 8].
Pedigree reconstruction from low-quality genetic material
In our study, the drawbacks of a non-invasive sampling scheme were dealt with effectively by combining it with Bayesian estimation of the pedigree. The lack of temporal resolution, individuals missing from the samples, and potential for high error rates associated with non-invasive sampling of genetic material  are also inherent constraints on exact pedigree reconstruction techniques [23–25, 27]. The MasterBayes estimation framework currently provides the most utilitarian solution to this problem, with the possibility of using chronological or other phenotypic data (where available), to improve relationship resolution while retaining the ability to assign unsampled parents  and thus being robust to incomplete sampling of the population. However, even in this flexible pedigree-fitting routine, it remains a challenge to separate between parent–offspring and sibling relationships (both r = 0.5) where one parent is unsampled and there is no temporal information to inform the direction of the relationship [23, 25]. For example, if the true father is missing from the data, the brother of an individual is likely to be assigned as that individual’s father . This may have contributed to the difficulty in confidently resolving relationships among the historical members of the Mynämäki pack. The earliest litters of individual AUL-001, a long-standing breeding female in that pack, were produced before our systematic monitoring began (see SI), likely with an unsampled male, reducing our confidence in distinguishing sibships from parentages. Similarly, individuals SAL-001 and REN-001 were assigned as father and son with an unsampled mother, but may in fact be brothers with both parents unsampled. This uncertainty must be borne in mind when interpreting this type of genetic monitoring data from incompletely sampled populations, as in the majority of conservation applications.
Nonetheless, our simulations showed that at present, our panel of 17 established microsatellites [28, 33, 34] is sufficient to make over 80% of assignments confidently and correctly, and that this could be increased to a maximum of almost 90% by adding 3 extra loci. Broadly, the main limitation on inference from the pedigree was confidence rather than accuracy of assignments: 10–20% of the simulation assignments were made with low confidence, whereas for panels larger than 5 loci only an average of 2.4% of assignments were incorrect. Confidence is likely limited in part by population structure and in part by assumptions in the analysis. Firstly, the close relatedness structure and limited genetic diversity of this population make it difficult to choose between genetically similar potential parents , often with descendants or siblings identified as parents. For example, AUL-001 carried common alleles at many loci, giving several other individuals, particularly her offspring, similar likelihoods of having produced her offspring. Secondly, we estimated genotyping error rates conservatively—per sample, not accounting for consensus genotypes based on multiple samples—and thus they are likely too high. This relaxation of the assumption that genotypes are true reduces the chance of erroneous assignments, but also reduces confidence as mismatches could be due to error either in the genotype or the assignment. Using a different marker type, single-nucleotide polymorphisms (SNPs), has the potential to improve confidences by reducing genotyping error due to the shorter sequences amplified [61, 62]. SNP-based pedigree reconstruction has already generated management-relevant information on reproductive behaviour in brown bears (Ursus arctos)  and a 96-locus multiplex SNP panel suitable for non-invasively collected genetic material has recently been developed for individual identification of wolves . However, the limited variation per SNP locus compared to multiallelic microsatellites, possibly exacerbated by relatively high rates of missing data (up to 10% of samples per locus; ), means that pedigree reconstruction requires a large number of markers to capture sufficient variation to resolve relationships. This sort of study would likely require a larger number of markers than currently established for wild wolves: using non-invasive samples from a wide-ranging brown bear population, 96 SNPs were sufficient to identify 433 individuals but to assign at least one confident parent to only 82 of these (19%; ).
The pipeline we have developed provides a reliable framework for pedigree reconstruction in this newly-established, incompletely and non-invasively sampled, multigenerational population. While some uncertainty remained in the computational outcome, this could mostly be satisfactorily resolved through manual corrections based on knowledge gained from the output pedigree. From the final pedigree, many aspects of population change of direct relevance to management decisions can be interpreted, such as dispersal and reproductive strategies, reconfigurations of mating pairs following mortality and average breeding tenure. For example, understanding dispersal patterns allows inference of likely responses of packs to the removal of breeding individuals, and hence the potential to manage gene flow between packs to some extent by targeting certain individuals during population-control hunting. Ongoing data collection will provide an opportunity to test such measures by examining changes in pack structure as a consequence of several breeding individuals being legally shot in 2016. These successes and opportunities support the continued use and planned expansion of this monitoring technique in this wolf population. If sampling and analysis can be conducted frequently enough, pack dynamics could be tracked in close to real time, perhaps even enabling prediction of certain reproductive events. Understanding how pack structure changes in the early stages of recolonisation illustrates how natural establishments and management interventions might proceed in other conservation contexts, such as reintroductions [16, 56]. As well as providing management-relevant information, the pedigree approach to characterizing population dynamics provides a qualitative description of the history and current structure of the population that may be more likely to be understood and accepted as unbiased by different stakeholders in the conflict. Unbiased information and public involvement are crucial to help move the fierce debate away from generalizations and emotion-driven claims towards a more evidence-based discussion, facilitating effective decision-making around the management of the wolves in this area and other re-establishing carnivores in densely-populated areas globally.
The non-invasive sampling scheme producing the main data of this study was initiated by TL, OT and JP, and continued and lead by IK who also contributed other types of samples. JP and OT produced and contributed one part of the genotype data, JH and JA another one, and ML and CP a third one. HG-W, CP and TL planned, and HG-W performed the pedigree analysis and other analyses. HG-W and TL wrote a draft manuscript, which all authors commented and substantially improved. All authors read and approved the final manuscript.
We thank all volunteer sample collectors, the wildlife agencies of Varsinais-Suomi and Satakunta and the wildlife council of Varsinais-Suomi for facilitating sample collection, and Antti Härkälä from Natural Resources Institute (Luke) for great help in logistics and in resolving the spatial configurations of the wolf territories. We are grateful to Dr. Sari Granholm for the dog samples and to Katja Salminen for skilled laboratory assistance.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset supporting the conclusions of this article is available in the Dryad repository (https://doi.org/10.5061/dryad.76tr6).
Consent for publication
Ethics approval and consent to participate
Funding for the initial stages of the laboratory work was awarded by Varsinais-Suomen maakuntaliitto and for the pedigree analysis from the University of Turku strategic research funds. OT was funded by a Marie Curie Intra European Fellowship within the 7th European Community Framework Program. We thank the Finnish Ministry for Agriculture and Forestry and the Natural Resources Institute Finland (Luke) for their support.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Estes JA, Terborgh J, Brashares JS, Power ME, Berger J, Bond WJ, et al. Trophic downgrading of planet Earth. Science. 2011;333:301–6.View ArticlePubMedGoogle Scholar
- Ripple WJ, Estes JA, Beschta RL, Wilmers CC, Ritchie EG, Hebblewhite M, et al. Status and ecological effects of the world’s largest carnivores. Science. 2014;343:151–62. https://doi.org/10.1126/science.1241484.View ArticleGoogle Scholar
- Loe J, Roskaft E. Large carnivores and human safety: a review. Ambio. 2004;33:283–8.View ArticlePubMedGoogle Scholar
- Penteriani V, Delgado MM, Pinchera F, Naves J, Fernandez-Gil A, Kojola I, et al. Human behaviour can trigger large carnivore attacks in developed countries. Sci Rep. 2016;6:20552. https://doi.org/10.1038/srep20552.View ArticlePubMedPubMed CentralGoogle Scholar
- Chapron G, Kaczensky P, Linnell JDC, von Arx M, Huber D, Andrén H, et al. Recovery of large carnivores in Europe’s modern human-dominated landscapes. Science. 2014;346:1517–9. https://doi.org/10.1126/science.1257553.View ArticlePubMedGoogle Scholar
- Smith DW, Peterson RO, Houston DB. Yellowstone after wolves. Bioscience. 2003;53:330–43.View ArticleGoogle Scholar
- Glikman JA, Vaske JJ, Bath AJ, Ciucci P, Boitani L. Residents’ support for wolf and bear conservation: the moderating influence of knowledge. Eur J Wildl Res. 2011;58:295–302.View ArticleGoogle Scholar
- Hiedanpää J. Institutional misfits: law and habits in Finnish wolf policy. Ecol Soc. 2013;18:24.View ArticleGoogle Scholar
- Boitani L. Wolf conservation and recovery. In: Mech LD, Boitani L, editors. Wolves behavior, ecology, and conservation. Chicago: University of Chicago Press; 2003. p. 317–40.Google Scholar
- Kojola I, Helle P, Heikkinen S, Lindén H, Paasivaara A, Wikman M. Tracks in snow and population size estimation: the wolf Canis lupus in Finland. Wildl Biol. 2014;20:279–84.View ArticleGoogle Scholar
- Hiedanpää J, Salo M, Kotilainen J. Teleodynamics and institutional change: the hardship of protecting the Amur tiger, big-leaf mahogany, and gray wolf. J Nat Cons. 2015;26:36–44.View ArticleGoogle Scholar
- Linnell JDC, Andersen R, Andersone Z, Balciauskas L, Blanco JC, Boitani L, et al. The fear of wolves: a review of wolf attacks on humans. NINA Oppdragsmeld. 2002;731:1–65.Google Scholar
- Pohja-Mykrä M, Kurki S. Strong community support for illegal killing challenges wolf management. Eur J Wildl Res. 2014;60:759–70.View ArticleGoogle Scholar
- Kojola I, Hallikainen V, Mikkola K, Gurarie E, Heikkinen S, Kaartinen S, Nikula A, Nivala V. Wolf visitations close to human residences in Finland: the role of age, residence density, and time of day. Biol Cons. 2016;198:9–14.View ArticleGoogle Scholar
- Waits LP, Paetkau D. Noninvasive genetic sampling tools for wildlife biologists: a review of applications and recommendations for accurate data collection. J Wildl Man. 2005;69:1419–33.View ArticleGoogle Scholar
- Bergl RA, Vigilant L. Genetic analysis reveals population structure and recent migration within the highly fragmented range of the Cross River gorilla (Gorilla gorilla diehli). Mol Ecol. 2007;16:501–16.View ArticlePubMedGoogle Scholar
- Schwartz MK, Monfort SL. Genetic and endocrine tools for carnivore surveys. In: Long RA, MacKay P, Zielinski WJ, Ray JC, editors. Non-invasive survey methods for carnivores. Washington: Island Press; 2008.Google Scholar
- Smith O, Wang J. When can noninvasive samples provide sufficient information in conservation genetics studies? Mol Ecol Res. 2014;14:1011–23.Google Scholar
- Lucchini V, Fabbri E, Marucco F, Ricci S, Boitani L, Randi E. Noninvasive molecular tracking of colonizing wolf (Canis lupus) packs in the western Italian Alps. Mol Ecol. 2002;11:857–68.View ArticlePubMedGoogle Scholar
- Stenglein JL, Waits LP, Ausband DE, Zager P, Mack CM. Estimating gray wolf pack size and family relationships using noninvasive genetic sampling at rendezvous sites. J Mammal. 2011;92:784–95.View ArticleGoogle Scholar
- Caniglia R, Fabbri E, Galaverni M, Milanesi P, Randi E. Noninvasive sampling and genetic variability, pack structure, and dynamics in an expanding wolf population. J Mammal. 2014;95:41–59.View ArticleGoogle Scholar
- Norman AJ, Spong G. Single nucleotide polymorphism-based dispersal estimates using noninvasive sampling. Ecol Evol. 2015;5:3056–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Hadfield JD, Richardson DS, Burke T. Towards unbiased parentage assignment: combining genetic, behavioural and spatial data in a Bayesian framework. Mol Ecol. 2006;15:3715–30.View ArticlePubMedGoogle Scholar
- Koch M, Hadfield JD, Sefc KM, Sturmbauer C. Pedigree reconstruction in wild cichlid fish populations. Mol Ecol. 2008;17:4500–11.View ArticlePubMedGoogle Scholar
- Jones A, Small C, Paczolt K, Ratterman N. A practical guide to methods of parentage analysis. Mol Ecol Res. 2010;10:6–30.View ArticleGoogle Scholar
- Aykanat T, Johnston SE, Cotter D, Cross TF, Poole R, Prodőhl PA, Reed T, Rogan G, McGinnity P, Primmer CR. Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis. BMC Evol Biol. 2014;14:1–17.View ArticleGoogle Scholar
- Walling CA, Pemberton JM, Hadfield JD, Kruuk LEB. Comparing parentage inference software: reanalysis of a red deer pedigree. Mol Ecol. 2010;19:1914–28.View ArticlePubMedGoogle Scholar
- Jansson E, Ruokonen M, Kojola I, Aspi J. Rise and fall of a wolf population: genetic diversity and structure during recovery, rapid expansion and drastic decline. Mol Ecol. 2012;21:5178–93.View ArticlePubMedGoogle Scholar
- Kojola I, Aspi J, Hakala A, Heikkinen S, Ilmoni C, Ronkainen S. Dispersal in an expanding wolf population in Finland. J Mammal. 2006;87:281–6.View ArticleGoogle Scholar
- Kojola I, Kaartinen S, Hakala A, Heikkinen S, Voipio H-M. Dispersal behavior and the connectivity between wolf populations in Northern Europe. J Wildl Man. 2009;72:309–13.View ArticleGoogle Scholar
- Barja I, de Miguel FJ, Barcena F. The importance of crossroads in faecal marking behaviour of the wolves (Canis lupus). Naturwissenschaften. 2004;91:489–92.View ArticlePubMedGoogle Scholar
- Gurarie E, Suutarinen J, Kojola I, Ovaskainen O. Summer movements, predation and habitat use of wolves in human modified boreal forests. Oecologia. 2011;165:891–903.View ArticlePubMedGoogle Scholar
- Fredholm M, Wintero AK. Variation of short tandem repeats within and between species belonging to the canidae family. Mamm Gen. 1995;6:11–8.View ArticleGoogle Scholar
- Francisco LV, Langston AA, Mellersh CS, Neal CL, Ostrander EA. A class of highly polymorphic tetranucleotide repeats for canine genetic mapping. Mamm Gen. 1996;7:359–62.View ArticleGoogle Scholar
- Morin PA, Chambers KE, Boesch C, Vigilant L. Quantitative polymerase chain reaction analysis of DNA from noninvasive samples for accurate microsatellite genotyping of wild chimpanzees (Pan troglodytes verus). Mol Ecol. 2001;10:1835–44.View ArticlePubMedGoogle Scholar
- Taberlet P, Luikart G, Waits LP. Noninvasive genetic sampling: look before you leap. Trends Ecol Evol. 1999;14:323–7.View ArticlePubMedGoogle Scholar
- Arandjelovic M, Head J, Kuhl H, Boesch C, Robbins MM, Maisels F, Vigilant L. Effective non-invasive genetic monitoring of multiple wild western gorilla groups. Biol Cons. 2010;143:1780–91.View ArticleGoogle Scholar
- Jansson E, Harmoinen J, Ruokonen M, Aspi J. Living on the edge: reconstructing the genetic history of the Finnish wolf population. BMC Evol Biol. 2014;14:64.View ArticlePubMedPubMed CentralGoogle Scholar
- Stansbury CR, Ausband DE, Zager P. A long-term population monitoring approach for a wide-ranging carnivore: noninvasive genetic sampling of gray wolf rendezvous sites in Idaho, USA. J Wildl Man. 2014;78:1040–9.View ArticleGoogle Scholar
- Seddon JM. Canid-specific primers for molecular sexing using tissue or non-invasive samples. Cons Gen. 2005;6:147–9.View ArticleGoogle Scholar
- Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16:1099–106.View ArticlePubMedGoogle Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Université de Montpellier II, Montpellier (France). 1996–2004.Google Scholar
- Anderson EC, Thompson EA. A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002;160:1217–29.PubMedPubMed CentralGoogle Scholar
- Wang J. Sibship reconstruction from genetic data with typing errors. Genetics. 2004;166:1963–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Mech LD, Boitani L. Wolves: behavior, ecology, and conservation. Chicago: University of Chicago Press; 2003.View ArticleGoogle Scholar
- Kuismin J. Assessing wolf’s (Canis lupus) population size in south-western Finland using genetic methods. Master’s thesis, University of Turku, Finland; 2015.Google Scholar
- Scandura M, Capitani C, Iacolina L, Marco A. An empirical approach for reliable microsatellite genotyping of wolf DNA from multiple noninvasive sources. Cons Gen. 2006;7:813–23.View ArticleGoogle Scholar
- Knapp SM, Craig BA, Waits LP. Incorporating genotyping error into non-invasive DNA-based mark—recapture population estimates. J Wildl Man. 2009;73:598–604.View ArticleGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna; 2013.Google Scholar
- Paradis E. Pegas: an R package for population genetics with an integrated—modular approach. Bioinformatics. 2010;26:419–20.View ArticlePubMedGoogle Scholar
- Jombart T. Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5.View ArticlePubMedGoogle Scholar
- Coster A. Pedigree: pedigree functions. R package version 1.4. https://CRAN.R-project.org/package=pedigree; 2013.
- Schwartz MK, Luikart G, Waples RS. Genetic monitoring as a promising tool for conservation and management. Trends Ecol Evol. 2007;22:25–33.View ArticlePubMedGoogle Scholar
- De Barba M, Waits LP, Garton EO, Genovesi P, Randi E, Mustoni A, Groff C. The power of genetic monitoring for studying demography, ecology and genetics of a reintroduced brown bear population. Mol Ecol. 2010;19:3938–51.View ArticlePubMedGoogle Scholar
- Geffen E, Kam M, Hefner R, Hersteinsson P, Angerbjorn A, Dalen L, Fuglei E, Noren K, Adams JR, Vucetich J, Meier TJ, Mech LD, vonHoldt BM, Stahler DR, Wayne RK. Kin encounter rate and inbreeding avoidance in canids. Mol Ecol. 2011;20:5348–58.View ArticlePubMedGoogle Scholar
- von Holdt BM, Stahler DR, Smith DW, Earl DA, Pollinger JP, Wayne RK. The genealogy and genetic viability of reintroduced Yellowstone grey wolves. Mol Ecol. 2008;17:252–74.View ArticleGoogle Scholar
- Åkesson M, Liberg O, Sand H, Wabakken P, Bensch S, Flagstad O. Genetic rescue in a severely inbred wolf population. Mol Ecol. 2016;25:4745–56.View ArticlePubMedPubMed CentralGoogle Scholar
- McNutt JW. Sex-biased dispersal in African wild dogs, Lycaon pictus. Anim Behav. 1996;52:1067–77.View ArticleGoogle Scholar
- Bradley BJ, Doran-Sheehy DM, Lukas D, Boesch C, Vigilant L. Dispersed male networks in western gorillas. Curr Biol. 2004;14:510–3.View ArticlePubMedGoogle Scholar
- Marshall TC, Slate J, Kruuk LEB, Pemberton J. Statistical confidence for likelihood-based paternity inference in natural populations. Mol Ecol. 1998;7:639–55.View ArticlePubMedGoogle Scholar
- Kraus RHS, vonHoldt B, Cocchiararo B, Harms V, Bayerl H, Kuhn R, Forster DW, Fickel J, Roos C, Nowak C. A single-nucleotide polymorphism-based approach for rapid and cost-effective genetic wolf monitoring in Europe based on noninvasively collected samples. Mol Ecol Res. 2015;15:295–305.View ArticleGoogle Scholar
- Fabbri E, Caniglia R, Mucci N, Thomsen HP, Krag K, Pertoldi C, Loeschke V, Randi E. Comparison of single nucleotide polymorphisms and microsatellites in non-invasive genetic monitoring of a wolf population. Arch Biol Sci. 2012;64:321–35.View ArticleGoogle Scholar