Dataset Information

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

ABSTRACT: 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.

SUBMITTER: Sethi SA

PROVIDER: S-EPMC5210676 | biostudies-literature | 2016 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

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

Sethi Suresh A SA Linden Daniel D Wenburg John J Lewis Cara C Lemons Patrick P Fuller Angela A Hare Matthew P MP

Royal Society open science 20161221 12

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 th ...[more]

PMID: 28083094

Similar Datasets

Project description:Obtaining robust estimates of population abundance is a central challenge hindering the conservation and management of many threatened and exploited species. Close-kin mark-recapture (CKMR) is a genetics-based approach that has strong potential to improve the monitoring of data-limited species by enabling estimates of abundance, survival, and other parameters for populations that are challenging to assess. However, CKMR models have received limited sensitivity testing under realistic population dynamics and sampling scenarios, impeding the application of the method in population monitoring programs and stock assessments. Here, we use individual-based simulation to examine how unmodeled population dynamics and aging uncertainty affect the accuracy and precision of CKMR parameter estimates under different sampling strategies. We then present adapted models that correct the biases that arise from model misspecification. Our results demonstrate that a simple base-case CKMR model produces robust estimates of population abundance with stable populations that breed annually; however, if a population trend or non-annual breeding dynamics are present, or if year-specific estimates of abundance are desired, a more complex CKMR model must be constructed. In addition, we show that CKMR can generate reliable abundance estimates for adults from a variety of sampling strategies, including juvenile-focused sampling where adults are never directly observed (and aging error is minimal). Finally, we apply a CKMR model that has been adapted for population growth and intermittent breeding to two decades of genetic data from juvenile lemon sharks (Negaprion brevirostris) in Bimini, Bahamas, to demonstrate how application of CKMR to samples drawn solely from juveniles can contribute to monitoring efforts for highly mobile populations. Overall, this study expands our understanding of the biological factors and sampling decisions that cause bias in CKMR models, identifies key areas for future inquiry, and provides recommendations that can aid biologists in planning and implementing an effective CKMR study, particularly for long-lived data-limited species.

Project description:TILLING (Targeting Induced Local Lesions IN Genomes) is a powerful reverse genetics method in plant functional genomics and breeding to identify mutagenized individuals with improved behavior for a trait of interest. Pooled high throughput sequencing (HTS) of the targeted genes allows efficient identification and sample assignment of variants within genes of interest in hundreds of individuals. Although TILLING has been used successfully in different crops and even applied to natural populations, one of the main issues for a successful TILLING experiment is that most currently available bioinformatics tools for variant detection are not designed to identify mutations with low frequencies in pooled samples or to perform sample identification from variants identified in overlapping pools. Our research group maintains the Next Generation Sequencing Experience Platform (NGSEP), an open source solution for analysis of HTS data. In this manuscript, we present three novel components within NGSEP to facilitate the design and analysis of TILLING experiments: a pooled variants detector, a sample identifier from variants detected in overlapping pools and a simulator of TILLING experiments. A new implementation of the NGSEP calling model for variant detection allows accurate detection of low frequency mutations within pools. The samples identifier implements the process to triangulate the mutations called within overlapping pools in order to assign mutations to single individuals whenever possible. Finally, we developed a complete simulator of TILLING experiments to enable benchmarking of different tools and to facilitate the design of experimental alternatives varying the number of pools and individuals per pool. Simulation experiments based on genes from the common bean genome indicate that NGSEP provides similar accuracy and better efficiency than other tools to perform pooled variants detection. To the best of our knowledge, NGSEP is currently the only tool that generates individual assignments of the mutations discovered from the pooled data. We expect that this development will be of great use for different groups implementing TILLING as an alternative for plant breeding and even to research groups performing pooled sequencing for other applications.

Project description:Abundance estimation of carnivore populations is difficult and has prompted the use of non-invasive detection methods, such as remotely-triggered cameras, to collect data. To analyze photo data, studies focusing on carnivores with unique pelage patterns have utilized a mark-recapture framework and studies of carnivores without unique pelage patterns have used a mark-resight framework. We compared mark-resight and mark-recapture estimation methods to estimate bobcat (Lynx rufus) population sizes, which motivated the development of a new "hybrid" mark-resight model as an alternative to traditional methods. We deployed a sampling grid of 30 cameras throughout the urban southern California study area. Additionally, we physically captured and marked a subset of the bobcat population with GPS telemetry collars. Since we could identify individual bobcats with photos of unique pelage patterns and a subset of the population was physically marked, we were able to use traditional mark-recapture and mark-resight methods, as well as the new "hybrid" mark-resight model we developed to estimate bobcat abundance. We recorded 109 bobcat photos during 4,669 camera nights and physically marked 27 bobcats with GPS telemetry collars. Abundance estimates produced by the traditional mark-recapture, traditional mark-resight, and "hybrid" mark-resight methods were similar, however precision differed depending on the models used. Traditional mark-recapture and mark-resight estimates were relatively imprecise with percent confidence interval lengths exceeding 100% of point estimates. Hybrid mark-resight models produced better precision with percent confidence intervals not exceeding 57%. The increased precision of the hybrid mark-resight method stems from utilizing the complete encounter histories of physically marked individuals (including those never detected by a camera trap) and the encounter histories of naturally marked individuals detected at camera traps. This new estimator may be particularly useful for estimating abundance of uniquely identifiable species that are difficult to sample using camera traps alone.

Project description:Accurate estimates of population abundance are a critical component of species conservation efforts in order to monitor the potential recovery of populations. Capture-mark-recapture (CMR) is a widely used approach to estimate population abundance, yet social species moving in groups violate the assumption of CMR approaches that all individuals in the population are detected independently. We developed a closed CMR model that addresses an important characteristic of group-living species-that individual-detection probability typically is conditional on group detection. Henceforth termed the Two-Step model, this approach first estimates group-detection probability and then-conditional on group detection-estimates individual-detection probability for individuals within detected groups. Overall abundance is estimated assuming that undetected groups have the same average group size as detected groups. We compared the performance of this Two-Step CMR model to a conventional (One-Step) closed CMR model that ignored group structure. We assessed model sensitivity to variation in both group- and individual-detection probability. Both models returned overall unbiased estimates of abundance, but the One-Step model returned deceptively narrow Bayesian confidence intervals (BCI) that failed to encompass the correct population abundance an average 52% of the time. Contrary, under the Two-Step model, CI coverage was on average 96%. Both models had similar root mean squared errors (RMSE), except for scenarios with low group detection probability, where the Two-Step model had much lower RMSE. For illustration with a real data set, we applied the Two-Step and regular model to non-invasive genetic capture-recapture data of mountain gorillas (Gorilla beringei beringei). As with simulations, abundance estimates under both models were similar, but the Two-Step model estimate had a wider confidence interval. Results support using the Two-Step model for species living in constant groups, particularly when group detection probability is low, to reduce risk of bias and adequately portray uncertainty in abundance estimates. Important sources of variation in detection need to be incorporated into the Two-Step model when applying it to field data.

Project description:Population parameters of poorly marked gregarious species are difficult to estimate. This is the case for common dolphins (Delphinus sp.), a genus known for its lack of distinctive marks resulting in a low mark ratio. Furthermore, the widespread nature of common dolphins results in low recaptures. We developed reliable photo-identification protocols to ensure accurate identification of individuals in the Hauraki Gulf, New Zealand. These protocols combined the use of nicks and notches and pigmentation patterns for identification and included the development of a distinctiveness threshold. The data were further stratified by the level of distinctiveness of each individual (as distinctive or highly-distinctive). Photo-identification surveys were conducted from January 2010 to December 2013. Mark-recapture techniques were implemented through a POPAN super-population approach to estimate seasonal apparent survival, capture probability and abundance of dolphins. A total of 2,083 unique adult common dolphins were identified, 51.3% were classified as D1 (highly distinctive; n = 1,069) and 48.7% as D2 (distinctive; n = 1,014). Of all individuals identified, 34.3% (n = 704) were re-sighted over subsequent years. The proportion of marked dolphins (when compared to unmarked dolphins) was 26.3% for D1 and 46.4% for D1 & D2, respectively. Apparent survival was estimated at 0.767 (CI = 0.694-0.827) for D1 animals, and 0.796 (CI = 0.729-0.850) for D1 & D2 combined. For D1 only, seasonal abundance varied from 732 (CI = 460-1,177) in autumn 2010 to 5,304 (CI = 4,745-5,930) in spring 2013. While the inclusion of D2 individuals may offer a more precise estimate of total abundance, the inability to determine additional sources of bias (for example, arising from under or overestimated mark ratios) meant that estimates for D1 individuals were deemed the least biased for this population. The photo-identification protocol, stratification of the data and steps taken to eliminate potential model violations provided a useful and novel approach to estimate population parameters for common dolphins. These approaches could be implemented for other large gregarious populations (≥500 individuals) of animals with poor natural markings.

Dataset Information

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

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets