A simulation test of the effectiveness of several methods for error-checking non-invasive genetic data

Estimating red deer (Cervus elaphus) population size based on non-invasive genetic sampling

Some deer species are of conservation concern; others are officially managed as a food source or for their trophies, whereas in many regions, deer are regarded as overabundant or even as a nuisance causing damages. Regardless of local management issues, in most cases, reliable data on deer population sizes and sex ratios are lacking. Non-invasive genetic approaches are promising tools for the estimation of population size and structure. We developed and tested a non-invasive genetic approach for red deer (Cervus elaphus) population size and density estimation based on faeces collected from three free-ranging red deer populations in south-western Germany. Altogether, we genotyped 2762 faecal samples, representing 1431 different individuals. We estimated population density for both sexes separately using two different approaches: spatially explicit capture-recapture (SECR) approach and a single-session urn model (CAPWIRE). The estimated densities of both approaches were similar for all three study areas, ranging between total densities of 3.3 (2.5–4.4) and 8.5 (6.4–11.3) red deer/km2. The estimated sex ratios differed significantly between the studied populations (ranging between 1:1.1 and 1:1.7), resulting in considerable consequences for management. In further research, the issues of population closure and approximation of the effectively sampled area for density estimation should be addressed. The presented approach can serve as a valuable tool for the management of deer populations, and to our knowledge, it represents the only sex-specific approach for estimation of red deer population size and density.

A pilot study-genetic diversity and population structure of snow leopards of Gilgit-Baltistan, Pakistan, using molecular techniques

Background

The Hindu Kush and Karakoram mountain ranges in Pakistan's northern areas are a natural habitat of the snow leopard (Panthera uncia syn. Uncia uncia) but the ecological studies on this animal are scarce since it is human shy by nature and lives in difficult mountainous tracts. The pilot study is conducted to exploit the genetic diversity and population structure of the snow leopard in this selected natural habitat of the member of the wildcat family in Pakistan.

Method

About 50 putative scat samples of snow leopard from five localities of Gilgit-Baltistan (Pakistan) along with a control sample of zoo maintained male snow leopard were collected for comparison. Significant quality and quantity of genomic DNA was extracted from scat samples using combined Zhang-phenol-chloroform method and successful amplification of cytochrome c oxidase I gene (190 bp) using mini-barcode primers, seven simple sequence repeats (SSR) markers and Y-linked AMELY gene (200 bp) was done.

Results

Cytochrome c oxidase I gene sequencing suggested that 33/50 (66%) scat samples were of snow leopard. AMELY primer suggested that out of 33 amplified samples, 21 (63.63%) scats were from male and 12 (36.36%) from female leopards. Through successful amplification of DNA of 25 out of 33 (75.75%) scat samples using SSR markers, a total of 68 alleles on seven SSR loci were identified, showing low heterozygosity, while high gene flow between population.

Discussion

The low gene flow rate among the population results in low genetic diversity causing decreased diversification. This affects the adaptability to climatic changes, thus ultimately resulting in decreased population size of the species.

Partitioning drivers of spatial genetic variation for a continuously distributed population of boreal caribou: Implications for management unit delineation

Isolation by distance (IBD) is a natural pattern not readily incorporated into theoretical models nor traditional metrics for differentiating populations, although clinal genetic differentiation can be characteristic of many wildlife species. Landscape features can also drive population structure additive to baseline IBD resulting in differentiation through isolation-by-resistance (IBR). We assessed the population genetic structure of boreal caribou across western Canada using nonspatial (STRUCTURE) and spatial (MEMGENE) clustering methods and investigated the relative contribution of IBD and IBR on genetic variation of 1,221 boreal caribou multilocus genotypes across western Canada. We further introduced a novel approach to compare the partitioning of individuals into management units (MU) and assessed levels of genetic connectivity under different MU scenarios. STRUCTURE delineated five genetic clusters while MEMGENE identified finer-scale differentiation across the study area. IBD was significant and did not differ for males and females both across and among detected genetic clusters. MEMGENE landscape analysis further quantified the proportion of genetic variation contributed by IBD and IBR patterns, allowing for the relative importance of spatial drivers, including roads, water bodies, and wildfires, to be assessed and incorporated into the characterization of population structure for the delineation of MUs. Local population units, as currently delineated in the boreal caribou recovery strategy, do not capture the genetic variation and connectivity of the ecotype across the study area. Here, we provide the tools to assess fine-scale spatial patterns of genetic variation, partition drivers of genetic variation, and evaluate the best management options for maintaining genetic connectivity. Our approach is highly relevant to vagile wildlife species that are of management and conservation concern and demonstrate varying degrees of IBD and IBR with clinal spatial genetic structure that challenges the delineation of discrete population boundaries.

Accurate recapture identification for genetic mark-recapture studies with error-tolerant likelihood-based match calling and sample clustering

Error-tolerant likelihood-based match calling presents a promising technique to accurately identify recapture events in genetic mark-recapture studies by combining probabilities of latent genotypes and probabilities of observed genotypes, which may contain genotyping errors. Combined with clustering algorithms to group samples into sets of recaptures based upon pairwise match calls, these tools can be used to reconstruct accurate capture histories for mark-recapture modelling. Here, we assess the performance of a recently introduced error-tolerant likelihood-based match-calling model and sample clustering algorithm for genetic mark-recapture studies. We assessed both biallelic (i.e. single nucleotide polymorphisms; SNP) and multiallelic (i.e. microsatellite; MSAT) markers using a combination of simulation analyses and case study data on Pacific walrus (Odobenus rosmarus divergens) and fishers (Pekania pennanti). A novel two-stage clustering approach is demonstrated for genetic mark-recapture applications. First, repeat captures within a sampling occasion are identified. Subsequently, recaptures across sampling occasions are identified. The likelihood-based matching protocol performed well in simulation trials, demonstrating utility for use in a wide range of genetic mark-recapture studies. Moderately sized SNP (64+) and MSAT (10-15) panels produced accurate match calls for recaptures and accurate non-match calls for samples from closely related individuals in the face of low to moderate genotyping error. Furthermore, matching performance remained stable or increased as the number of genetic markers increased, genotyping error notwithstanding.

Evaluating noninvasive genetic sampling techniques to estimate large carnivore abundance

Monitoring large carnivores is difficult because of intrinsically low densities and can be dangerous if physical capture is required. Noninvasive genetic sampling (NGS) is a safe and cost-effective alternative to physical capture. We evaluated the utility of two NGS methods (scat detection dogs and hair sampling) to obtain genetic samples for abundance estimation of coyotes, black bears and Canada lynx in three areas of Newfoundland, Canada. We calculated abundance estimates using program capwire, compared sampling costs, and the cost/sample for each method relative to species and study site, and performed simulations to determine the sampling intensity necessary to achieve abundance estimates with coefficients of variation (CV) of <10%. Scat sampling was effective for both coyotes and bears and hair snags effectively sampled bears in two of three study sites. Rub pads were ineffective in sampling coyotes and lynx. The precision of abundance estimates was dependent upon the number of captures/individual. Our simulations suggested that ~3.4 captures/individual will result in a < 10% CV for abundance estimates when populations are small (23-39), but fewer captures/individual may be sufficient for larger populations. We found scat sampling was more cost-effective for sampling multiple species, but suggest that hair sampling may be less expensive at study sites with limited road access for bears. Given the dependence of sampling scheme on species and study site, the optimal sampling scheme is likely to be study-specific warranting pilot studies in most circumstances.

Enhanced understanding of predator-prey relationships using molecular methods to identify predator species, individual and sex

Predator species identification is an important step in understanding predator-prey interactions, but predator identifications using kill site observations are often unreliable. We used molecular tools to analyse predator saliva, scat and hair from caribou calf kills in Newfoundland, Canada to identify the predator species, individual and sex. We sampled DNA from 32 carcasses using cotton swabs to collect predator saliva. We used fragment length analysis and sequencing of mitochondrial DNA to distinguish between coyote, black bear, Canada lynx and red fox and used nuclear DNA microsatellite analysis to identify individuals. We compared predator species detected using molecular tools to those assigned via field observations at each kill. We identified a predator species at 94% of carcasses using molecular methods, while observational methods assigned a predator species to 62.5% of kills. Molecular methods attributed 66.7% of kills to coyote and 33.3% to black bear, while observations assigned 40%, 45%, 10% and 5% to coyote, bear, lynx and fox, respectively. Individual identification was successful at 70% of kills where a predator species was identified. Only one individual was identified at each kill, but some individuals were found at multiple kills. Predator sex was predominantly male. We demonstrate the first large-scale evaluation of predator species, individual and sex identification using molecular techniques to extract DNA from swabs of wild prey carcasses. Our results indicate that kill site swabs (i) can be highly successful in identifying the predator species and individual responsible; and (ii) serve to inform and complement traditional methods.

Estimating Size and Trend of the North Interlake Woodland Caribou Population Using Fecal-DNA and Capture-Recapture Models

A critical step in recovery efforts for endangered and threatened species is the monitoring of population demographic parameters. As part of these efforts, we evaluated the use of fecal-DNA based capture–recapture methods to estimate population sizes and population rate of change for the North Interlake woodland caribou herd (Rangifer tarandus caribou), Manitoba, Canada. This herd is part of the boreal population of woodland caribou, listed as threatened under the federal Species at Risk Act (2003) and the provincial Manitoba Endangered Species Act (2006). Between 2004 and 2009 (9 surveys), we collected 1,080 fecal samples and identified 180 unique genotypes (102 females and 78 males). We used a robust design survey plan with 2 surveys in most years and analysed the data with Program MARK to estimate encounter rates (p), apparent survival rates (ϕ), rates of population change (λ), and population sizes (N). We estimated these demographic parameters for males and females and for 2 genetic clusters within the North Interlake. The population size estimates were larger for the Lower than the Upper North Interlake area and the proportion of males was lower in the Lower (33%) than the Upper North Interlake (49%). Population rate of change for the entire North Interlake area (2005–2009) using the robust design Pradel model was significantly <1.0 (λ = 0.90, 95% CI: 0.82–0.99) and varied between sex and area with the highest being for males in Lower North Interlake (λ = 0.98, 95% CI: 0.83–1.13) and the lowest being for females in Upper North Interlake (λ = 0.83, 95% CI: 0.69–0.97). The additivity of λ between sex and area is supported on the log scale and translates into males having a λ that is 0.09 greater than females and independent of sex, Lower North Interlake having a λ that is 0.06 greater than Upper North Interlake. Population estimates paralleled these declining trends, which correspond to trends observed in other fragmented populations of woodland caribou along the southern part of their range. The results of this study clearly demonstrate the applicability and success of non-invasive genetic sampling in monitoring populations of woodland caribou. © 2012 The Wildlife Society.

How many genetic markers to tag an individual? An empirical assessment of false matching rates among close relatives

Genetic identification of individuals is now commonplace, enabling the application of tagging methods to elusive species or species that cannot be tagged by traditional methods. A key aspect is determining the number of loci required to ensure that different individuals have non-matching multi-locus genotypes. Closely related individuals are of particular concern because of elevated matching probabilities caused by their recent co-ancestry. This issue may be addressed by increasing the number of loci to a level where full siblings (the relatedness category with the highest matching probability) are expected to have non-matching multi-locus genotypes. However, increasing the number of loci to meet this "full-sib criterion" greatly increases the laboratory effort, which in turn may increase the genotyping error rate resulting in an upward-biased mark-recapture estimate of abundance as recaptures are missed due to genotyping errors. We assessed the contribution of false matches from close relatives among 425 maternally related humpback whales, each genotyped at 20 microsatellite loci. We observed a very low (0.5-4%) contribution to falsely matching samples from pairs of first-order relatives (i.e., parent and offspring or full siblings). The main contribution to falsely matching individuals from close relatives originated from second-order relatives (e.g., half siblings), which was estimated at 9%. In our study, the total number of observed matches agreed well with expectations based upon the matching probability estimated for unrelated individuals, suggesting that the full-sib criterion is overly conservative, and would have required a 280% relative increase in effort. We suggest that, under most circumstances, the overall contribution to falsely matching samples from close relatives is likely to be low, and hence applying the full-sib criterion is unnecessary. In those cases where close relatives may present a significant issue, such as unrepresentative sampling, we propose three different genotyping strategies requiring only a modest increase in effort, which will greatly reduce the number of false matches due to the presence of related individuals.

Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears

A fundamental challenge to estimating population size with mark-recapture methods is heterogeneous capture probabilities and subsequent bias of population estimates. Confronting this problem usually requires substantial sampling effort that can be difficult to achieve for some species, such as carnivores. We developed a methodology that uses two data sources to deal with heterogeneity and applied this to DNA mark-recapture data from grizzly bears (Ursus arctos). We improved population estimates by incorporating additional DNA "captures" of grizzly bears obtained by collecting hair from unbaited bear rub trees concurrently with baited, grid-based, hair snag sampling. We consider a Lincoln-Petersen estimator with hair snag captures as the initial session and rub tree captures as the recapture session and develop an estimator in program MARK that treats hair snag and rub tree samples as successive sessions. Using empirical data from a large-scale project in the greater Glacier National Park, Montana, USA, area and simulation modeling we evaluate these methods and compare the results to hair-snag-only estimates. Empirical results indicate that, compared with hair-snag-only data, the joint hair-snag-rub-tree methods produce similar but more precise estimates if capture and recapture rates are reasonably high for both methods. Simulation results suggest that estimators are potentially affected by correlation of capture probabilities between sample types in the presence of heterogeneity. Overall, closed population Huggins-Pledger estimators showed the highest precision and were most robust to sparse data, heterogeneity, and capture probability correlation among sampling types. Results also indicate that these estimators can be used when a segment of the population has zero capture probability for one of the methods. We propose that this general methodology may be useful for other species in which mark-recapture data are available from multiple sources.

Patterns of paternity and group fission in wild multimale mountain gorilla groups

To understand variation among social systems, it is essential to know the relative reproductive success of individuals in group-living species. Particularly interesting for such studies are taxa such as mountain gorillas in which both one-male and multimale groups are common, because of the opportunity to estimate the costs and benefits to males of pursuing different reproductive strategies. We genotyped 68 individuals from two groups of multimale mountain gorilla groups in Bwindi Impenetrable National Park, Uganda to determine the distribution of paternity among the males. In both groups, the dominant male sired the majority of offspring. One group underwent a fission, and we found that the eight offspring assigned to the dominant silverback (and their mothers) remained with their father, while the two offspring of unknown paternity ended up in the small group headed by the formerly subordinate silverback. This is consistent with the proposal that the outcome of group fission in primates is not only influenced by maternal relationships among individuals, but also by patrilineal relationships. Results of this study show that subordinate males may gain reproductive benefits even while queuing for dominance status. Despite ecological differences between Bwindi and the Virunga Volcanoes, male mountain gorillas living in both populations benefit from remaining in multimale groups.