Importance of a pilot study for non-invasive genetic sampling: genotyping errors and population size estimation in red deer

Population size information is critical for managing endangered or harvested populations. Population size can now be estimated from non-invasive genetic sampling. However, pitfalls remain such as genotyping errors (allele dropout and false alleles at microsatellite loci). To evaluate the feasibility of non-invasive sampling (e.g., for population size estimation), a pilot study is required. Here, we present a pilot study consisting of (i) a genetic step to test loci amplification and to estimate allele frequencies and genotyping error rates when using faecal DNA, and (ii) a simulation step to quantify and minimise the effects of errors on estimates of population size. The pilot study was conducted on a population of red deer in a fenced natural area of 5440 ha, in France. Twelve microsatellite loci were tested for amplification and genotyping errors. The genotyping error rates for microsatellite loci were 0–0.83 (mean=0.2) for allele dropout rates and 0–0.14 (mean=0.02) for false allele rates, comparable to rates encountered in other non-invasive studies. Simulation results suggest we must conduct 6 PCR amplifications per sample (per locus) to achieve approximately 97% correct genotypes. The 3% error rate appears to have little influence on the accuracy and precision of population size estimation. This paper illustrates the importance of conducting a pilot study (including genotyping and simulations) when using non-invasive sampling to study threatened or managed populations.     

megasat: automated inference of microsatellite genotypes from sequence data

megasat is software that enables genotyping of microsatellite loci using next‐generation sequencing data. Microsatellites are amplified in large multiplexes, and then sequenced in pooled amplicons. megasat reads sequence files and automatically scores microsatellite genotypes. It uses fuzzy matches to allow for sequencing errors and applies decision rules to account for amplification artefacts, including nontarget amplification products, replication slippage during PCR (amplification stutter) and differential amplification of alleles. An important feature of megasat is the generation of histograms of the length–frequency distributions of amplification products for each locus and each individual. These histograms, analogous to electropherograms traditionally used to score microsatellite genotypes, enable rapid evaluation and editing of automatically scored genotypes. megasat is written in Perl, runs on Windows, Mac OS X and Linux systems, and includes a simple graphical user interface. We demonstrate megasat using data from guppy, Poecilia reticulata. We genotype 1024 guppies at 43 microsatellites per run on an Illumina MiSeq sequencer. We evaluated the accuracy of automatically called genotypes using two methods, based on pedigree and repeat genotyping data, and obtained estimates of mean genotyping error rates of 0.021 and 0.012. In both estimates, three loci accounted for a disproportionate fraction of genotyping errors; conversely, 26 loci were scored with 0–1 detected error (error rate ≤0.007). Our results show that with appropriate selection of loci, automated genotyping of microsatellite loci can be achieved with very high throughput, low genotyping error and very low genotyping costs.     

Quantitative PCR assessment of microsatellite and SNP genotyping with variable quality DNA extracts

In this study we developed eight quantitative PCR (qPCR) assays to evaluate the starting copy number of nuclear and mitochondrial DNA fragments ranging from 75 to 350 base-pairs in DNA extracts from Chinook salmon tissues with varying quality. Samples were genotyped with 13 microsatellite and 29 SNP assays and average genotyping success for good, intermediate, and poor quality samples was 96%, 24%, and 24% for microsatellite loci, and 98%, 97%, and 79% for SNPs, respectively. As measured by qPCR, good quality samples had a consistently high number of starting copies across all fragment sizes with little change between the smallest and largest size. In contrast, the intermediate and poor quality samples displayed decreases in starting copy number as fragment size increased, and was most pronounced with poor samples. Logistic regression of genotyping success by starting copy number indicated that in order to achieve at least 90% genotyping success, approximately 1,000 starting copies of nuclear DNA are necessary for microsatellite loci, and as few as 14 starting copies for SNP assays (but we recommend at least 50 copies to reduce genotyping error). While these guidelines apply specifically to Chinook salmon and the genetic markers included in this study, the principles are transferable to other species and markers due to the underlying process associated with template quantity and PCR amplification.     

Multiplex pre-amplification for noninvasive genetic sampling: is the extra effort worth it?

Microsatellite genotyping of hair and faeces using standard polymerase chain reaction (PCR) resulted in low success rates and high error rates in a 2003–2004 pilot study using noninvasive genetic sampling for the brown bear (Ursus arctos) in the Italian Alps. Thus, we evaluated the performance of multiplex pre-amplification for improving microsatellite genotyping results. Brown bear faecal DNA extracts of varying quality (n = 33) and hair DNA extracts of poor (n = 32) and good (n = 34) quality were used to compare standard PCR and pre-amplification. In contrast to previous studies, there was no significant difference between methods for individual locus amplification success, genotyping error and genotyping success rates for scat and hair samples. The use of pre-amplification requires an additional investment of time and resources, and our results raise questions about the universal value of pre-amplification approaches. We suggest that researchers carefully evaluate the performance of pre-amplification compared to standard PCR using field-collected samples from the study area of interest before engaging in large-scale noninvasive genetic analyses.     

Application of automated DNA sizing technology for genotyping microsatellite loci.

Highly polymorphic microsatellite loci offer great promise for gene mapping studies, but fulfillment of this potential will require substantial improvements in methods for accurate and efficient genotyping. Here, we report a genotyping method based on fluorescently labeled PCR primers and size characterization of PCR products using an automated DNA fragment analyzer. We capitalize on the availability of three distinct fluorescent dyes to label uniquely loci that overlap in size, and this innovation increases by threefold the number of loci that can be analyzed simultaneously. We label size standards with a fourth dye and combine these with the microsatellite PCR products in each gel lane. Computer programs provide very rapid and accurate sizing of microsatellite alleles and efficient data management. In addition, fluorescence signals are linear over a much greater range of intensity than conventional autoradiography. This facilitates multiplexing of loci (since signal intensities often vary greatly) and helps distinguish major peaks from artifacts, thereby improving genotyping accuracy.     

Low genotyping error rates in wild ungulate faeces sampled in winter

We show that Alpine ibex (Capra ibex) and Corsican mouflon (Ovis musimon) faeces yield useful DNA for microsatellite analysis, however, we detected higher genotyping error rates for spring faeces than for winter faeces. We quantified the genotyping error rate by repeatedly genotyping four microsatellites. Respectively, 99 and 95% of mouflon and ibex genotyping repetitions provided a correct genotype using winter samples, whereas spring samples provided only 52 and 59% correct genotypes. Thus, before starting a noninvasive study, we recommend that researchers conduct a pilot study to quantify genotyping error rates for each season, population and species to be studied.     

A test of the multiplex pre-amplification approach in microsatellite genotyping of wolverine faecal DNA

Recently, a two-step PCR approach, referred to as multiplex pre-amplification, was proposed to improve microsatellite amplification from non-invasive samples such as faecal DNA. Here, we compare this new approach to standard PCR with respect to amplification success and genotyping error rates in microsatellite analysis (18 markers) of wolverine faecal DNA (48 extracts initially shown to contain amplifiable DNA). The multiplex pre-amplification approach was clearly advantageous both in terms of successful PCR amplifications (91% vs. 80%) and allelic dropout rate (2.4% vs. 12.5%). However, dropouts were to a high extent repeated in all second-step amplifications following multiplex pre-amplification, indicative of being generated during the initial PCR. Analysing more than one PCR from the initial multiplex PCR product may thus be of limited value. We instead suggest to perform two initial multiplex PCRs and to analyse a single second-step PCR from each of them. This was tested for 22 extracts at 18 loci and proved to be an effective way to obtaining a correct genotype.     

Multiplex preamplification PCR and microsatellite validation enables accurate single nucleotide polymorphism genotyping of historical fish scales

Incorporating historical tissues into the study of ecological, conservation and management questions can broaden the scope of population genetic research by enhancing our understanding of evolutionary processes and anthropogenic influences on natural populations. Genotyping historical and low-quality samples has been plagued by challenges associated with low amounts of template DNA and the potential for pre-existing DNA contamination among samples. We describe a two-step process designed to (i) accurately genotype large numbers of historical low-quality scale samples in a high-throughput format and (ii) screen samples for pre-existing DNA contamination. First, we describe how an efficient multiplex preamplification PCR of 45 single nucleotide polymorphisms (SNPs) can generate highly accurate genotypes with low failure and error rates in subsequent SNP genotyping reactions of individual historical scales from sockeye salmon (Oncorhynchus nerka). Second, we demonstrate how the method can be modified for the amplification of microsatellite loci to detect pre-existing DNA contamination. A total of 760 individual historical scale and 182 contemporary fin clip samples were genotyped and screened for contamination. Genotyping failure and error rates were exceedingly low and similar for both historical and contemporary samples. Pre-existing contamination in 21% of the historical samples was successfully identified by screening the amplified microsatellite loci. The advantages of automation, low failure and error rates, and ability to multiplex both the preamplification and subsequent genotyping reactions combine to make the protocol ideally suited for efficiently genotyping large numbers of potentially contaminated low-quality sources of DNA.     

Effect of genotyping error in model-free linkage analysis using microsatellite or single-nucleotide polymorphism marker maps

Errors while genotyping are inevitable and can reduce the power to detect linkage. However, does genotyping error have the same impact on linkage results for single-nucleotide polymorphism (SNP) and microsatellite (MS) marker maps? To evaluate this question we detected genotyping errors that are consistent with Mendelian inheritance using large changes in multipoint identity-by-descent sharing in neighboring markers. Only a small fraction of Mendelian consistent errors were detectable (e.g., 18% of MS and 2.4% of SNP genotyping errors). More SNP genotyping errors are Mendelian consistent compared to MS genotyping errors, so genotyping error may have a greater impact on linkage results using SNP marker maps. We also evaluated the effect of genotyping error on the power and type I error rate using simulated nuclear families with missing parents under 0, 0.14, and 2.8% genotyping error rates. In the presence of genotyping error, we found that the power to detect a true linkage signal was greater for SNP (75%) than MS (67%) marker maps, although there were also slightly more false-positive signals using SNP marker maps (5 compared with 3 for MS). Finally, we evaluated the usefulness of accounting for genotyping error in the SNP data using a likelihood-based approach, which restores some of the power that is lost when genotyping error is introduced.     

A test of the efficacy of whole‐genome amplification on DNA obtained from low‐yield samples

Conservation and population genetic studies are sometimes hampered by insufficient quantities of high quality DNA. One potential way to overcome this problem is through the use of whole genome amplification (WGA) kits. We performed rolling circle WGA on DNA obtained from matched hair and tissue samples of North American red squirrels (Tamiasciurus hudsonicus). Following polymerase chain reaction (PCR) at four microsatellite loci, we compared genotyping success for DNA from different source tissues, both pre‐ and post‐WGA. Genotypes obtained with tissue were robust, whether or not DNA had been subjected to WGA. DNA extracted from hair produced results that were largely concordant with matched tissue samples, although amplification success was reduced and some allelic dropout was observed. WGA of hair samples resulted in a low genotyping success rate and an unacceptably high rate of allelic dropout and genotyping error. The problem was not rectified by conducting PCR of WGA hair samples in triplicate. Therefore, we conclude that WGA is only an effective method of enhancing template DNA quantity when the initial sample is from high‐yield material.     

DNA degradation in avian faecal samples and feasibility of non-invasive genetic studies of threatened capercaillie populations

We evaluated the feasibility of using faeces as a non-invasively collected DNA source for the genetic study of an endangered bird population (capercaillie; Tetrao urogallus). We used a multitube approach, and for our panel of 11 microsatellites genotyping reliability was estimated at 98% with five repetitions. Experiments showed that free DNases in faecal material were the major cause of DNA degradation. Our results demonstrate that using avian faeces as a source of DNA, reliable microsatellite genotyping can be obtained with a reasonable number of PCR replicates.     

Population size estimation in Yellowstone wolves with error-prone noninvasive microsatellite genotypes

Determining population sizes can be difficult, but is essential for conservation. By counting distinct microsatellite genotypes, DNA from noninvasive samples (hair, faeces) allows estimation of population size. Problems arise because genotypes from noninvasive samples are error-prone, but genotyping errors can be reduced by multiple polymerase chain reaction (PCR). For faecal genotypes from wolves in Yellowstone National Park, error rates varied substantially among samples, often above the 'worst-case threshold' suggested by simulation. Consequently, a substantial proportion of multilocus genotypes held one or more errors, despite multiple PCR. These genotyping errors created several genotypes per individual and caused overestimation (up to 5.5-fold) of population size. We propose a 'matching approach' to eliminate this overestimation bias.     

Improved noninvasive genotyping method: application to brown bear (Ursus arctos) faeces

We redesigned new microsatellite primers and one sex‐specific primer for amplification of faecal DNA from brown bears (Ursus arctos). We also combined a semi‐nested polymerase chain reaction (PCR) with a newly developed multiplex preamplification method in order to increase the quality of the amplified DNA fragments. In comparison with a conventional PCR approach, the genotyping error rate was substantially reduced and the amplification rate was increased. This new approach could be transposed to other species where conventional PCR methods experience low success due to limited DNA concentration and/or quality.     

Two-step multiplex polymerase chain reaction improves the speed and accuracy of genotyping using DNA from noninvasive and museum samples

Many studies in molecular ecology rely upon the genotyping of large numbers of low‐quantity DNA extracts derived from noninvasive or museum specimens. To overcome low amplification success rates and avoid genotyping errors such as allelic dropout and false alleles, multiple polymerase chain reaction (PCR) replicates for each sample are typically used. Recently, two‐step multiplex procedures have been introduced which drastically increase the success rate and efficiency of genotyping. However, controversy still exists concerning the amount of replication needed for suitable control of error. Here we describe the use of a two‐step multiplex PCR procedure that allows rapid genotyping using at least 19 different microsatellite loci. We applied this approach to quantified amounts of noninvasive DNAs from western chimpanzee, western gorilla, mountain gorilla and black and white colobus faecal samples, as well as to DNA from ~100‐year‐old gorilla teeth from museums. Analysis of over 45 000 PCRs revealed average success rates of > 90% using faecal DNAs and 74% using museum specimen DNAs. Average allelic dropout rates were substantially reduced compared to those obtained using conventional singleplex PCR protocols, and reliable genotyping using low (< 25 pg) amounts of template DNA was possible. However, four to five replicates of apparently homozygous results are needed to avoid allelic dropout when using the lowest concentration DNAs (< 50 pg/reaction), suggesting that use of protocols allowing routine acceptance of homozygous genotypes after as few as three replicates may lead to unanticipated errors when applied to low‐concentration DNAs.     

Microsatellite markers for population-genetic studies of shiitake (Lentinula edodes) strains

In order to analyze the genotyping of shiitake (Lentinula edodes, Berk.), one of commercially and widely grown edible mushrooms, we examined group of total 89 strains that are registered in Korea, Japan, and China respectively, using five microsatellite markers (Led A2, Led A8, Led B2, Led B6, and Led D6) registered in NCBI (National Center for Biotechnology Information). After we prepared synthesis of primers modified with 5′-FAM fluorescent dye and conducted PCR program, we obtained reasonable products. And then we performed microsatellite genotyping analysis using an ABI 3730xl Genetic Analyzer. According to genotyping analysis, the number of alleles of each microsatellite marker ranged from 5 to 14 with the average value of 8.2. The expected and the observed heterozygosity over all microsatellite ranged from 0.27–0.83 and 0.14–0.61. Among 5 microsatellite markers PIC (Polymorphism Information Content) values of Led A8, Led D6, and Led B6 resulted 0.82, 0.60, and 0.57 respectively. We observed that these values showed relatively high values discriminating from other values of microsatellite marker. The average of total values of Led A8, Led D6, and Led B6 came out 0.53 as result, which is higher than 0.5. As a result, this average can be subject to significant value to be used as marker. By using microsatellite marker we can analogically establish population relationships among shiitake (L. edodes) strains grown in East Asian region, develop new varieties, and propose discriminative criteria for different breeding and variety classification. Moreover, microsatellite markers enable us to obtain the genetic inheritance and information that can be used for protection of intellectual property rights.     

Fast and cost‐effective single nucleotide polymorphism (SNP) detection in the absence of a reference genome using semideep next‐generation Random Amplicon Sequencing (RAMseq)

Biodiversity has suffered a dramatic global decline during the past decades, and monitoring tools are urgently needed providing data for the development and evaluation of conservation efforts both on a species and on a genetic level. However, in wild species, the assessment of genetic diversity is often hampered by the lack of suitable genetic markers. In this article, we present Random Amplicon Sequencing (RAMseq), a novel approach for fast and cost‐effective detection of single nucleotide polymorphisms (SNPs) in nonmodel species by semideep sequencing of random amplicons. By applying RAMseq to the Eurasian otter (Lutra lutra), we identified 238 putative SNPs after quality filtering of all candidate loci and were able to validate 32 of 77 loci tested. In a second step, we evaluated the genotyping performance of these SNP loci in noninvasive samples, one of the most challenging genotyping applications, by comparing it with genotyping results of the same faecal samples at microsatellite markers. We compared (i) polymerase chain reaction (PCR) success rate, (ii) genotyping errors and (iii) Mendelian inheritance (population parameters). SNPs produced a significantly higher PCR success rate (75.5% vs. 65.1%) and lower mean allelic error rate (8.8% vs. 13.3%) than microsatellites, but showed a higher allelic dropout rate (29.7% vs. 19.8%). Genotyping results showed no deviations from Mendelian inheritance in any of the SNP loci. Hence, RAMseq appears to be a valuable tool for the detection of genetic markers in nonmodel species, which is a common challenge in conservation genetic studies.     

A nondestructive,rapid, reliable and inexpensive method to sample,store and extract high‐quality DNA from fish body mucus and buccal cells

A rapid, nondestructive, reproducible and cheap DNA extraction method from body mucus and buccal cells of northern pike and brown trout is described. Buccal cells and body mucus were sampled on FTA Cards; the captured DNA was used directly for microsatellite and polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analyses. A complete concordance with control DNA was found. The genotyping error rate for microsatellite ranged from 1.9% to 3.3% for the northern pike and brown trout, respectively. This methodology, using for the first time these materials as a fish DNA source, combines speed of sampling and processing, with a twofold to a threefold time and costs saving.     

An efficient method for screening faecal DNA genotypes and detecting new individuals and hybrids in the red wolf (<Emphasis Type="Italic">Canis rufus</Emphasis>) experimental population area

Previously, sequencing of mitochondrial DNA (mtDNA) from non-invasively collected faecal material (scat) has been used to help manage hybridization in the wild red wolf (Canis rufus) population. This method is limited by the maternal inheritance of mtDNA and the inability to obtain individual identification. Here, we optimize the use of nuclear DNA microsatellite markers on red wolf scat DNA to distinguish between individuals and detect hybrids. We develop a data filtering method in which scat genotypes are compared to known blood genotypes to reduce the number of PCR amplifications needed. We apply our data filtering method and the more conservative maximum likelihood ratio method (MLR) of Miller et al. (2002 Genetics 160:357–366) to a scat dataset previously screened for hybrids by sequencing of mtDNA. Using seven microsatellite loci, we obtained genotypes for 105 scats, which were matched to 17 individuals. The PCR amplification success rate was 50% and genotyping error rates ranged from 6.6% to 52.1% per locus. Our data filtering method produced comparable results to the MLR method, and decreased the time and cost of analysis by 25%. Analysis of this dataset using our data filtering method verified that no hybrid individuals were present in the Alligator River National Wildlife Refuge, North Carolina in 2000. Our results demonstrate that nuclear DNA microsatellite analysis of red wolf scats provides an efficient and accurate approach to screen for new individuals and hybrids.     

Undetected genotyping errors cause apparent overtransmission of common alleles in the transmission/disequilibrium test

The transmission/disequilibrium test (TDT), a family-based test of linkage and association, is a popular and intuitive statistical test for studies of complex inheritance, as it is nonparametric and robust to population stratification. We carried out a literature search and located 79 significant TDT-derived associations between a microsatellite marker allele and a disease. Among these, there were 31 (39%) in which the most common allele was found to exhibit distorted transmission to affected offspring, implying that the allele may be associated with either susceptibility to or protection from a disease. In 27 of these 31 studies (87%), the most common allele appeared to be overtransmitted to affected offspring (a risk factor), and, in the remaining 4 studies, the most common allele appeared to be undertransmitted (a protective factor). In a second literature search, we identified 92 case-control studies in which a microsatellite marker allele was found to have significantly different frequencies in case and control groups. Of these, there were 37 instances (40%) in which the most common allele was involved. In 12 of these 37 studies (32%), the most common allele was enriched in cases relative to controls (a risk factor), and, in the remaining 25 studies, the most common allele was enriched in controls (a protective factor). Thus, the most common allele appears to be a risk factor when identified through the TDT, and it appears to be protective when identified through case-control analysis. To understand this phenomenon, we incorporated an error model into the calculation of the TDT statistic. We show that undetected genotyping error can cause apparent transmission distortion at markers with alleles of unequal frequency. We demonstrate that this distortion is in the direction of overtransmission for common alleles. Therefore, we conclude that undetected genotyping errors may be contributing to an inflated false-positive rate among reported TDT-derived associations and that genotyping fidelity must be increased.     

Comparison of random and SSR-enriched shotgun pyrosequencing for microsatellite discovery and single multiplex PCR optimization in Acacia harpophylla F. Muell. Ex Benth

Streamlining the development and genotyping of microsatellites in species for which no genetic information is available represents an important technical challenge to overcome in order to enable mainstream application of state-of-the-art population genetic analysis techniques in nonmodel organisms. Using the example of Acacia harpophylla, an acacia tree endemic of north-eastern Australia, we show that high-throughput shotgun pyrosequencing technology, so-called second-generation sequencing, reduces time and cost of microsatellite marker discovery in nonmodel organisms and of their large-scale typing in natural populations. We found that 0.5% of short sequence reads generated on 454 Genome Sequencer FLX Titanium from random genome sampling and 2.2% of reads generated with prior microsatellite enrichment yielded microsatellite markers with designed polymerase chain reaction (PCR) primers, suggesting that enrichment increases efficiency of pyrosequencing when microsatellite discovery is the primary goal. Using stringent selection criteria to facilitate downstream PCR multiplex design, we identified 1435 microsatellite loci with designed primers from a total of 200,908 short sequence reads. From a subset of 96 loci tested for amplification, 38 were validated for population genetics applications, leading to the optimization of a cost-effective multiplex PCR protocol for the simultaneous typing of nine microsatellites in natural populations of A. harpophylla.