Eva Bellemain

How to track and assess genotyping errors in population genetics studies

Genotyping errors occur when the genotype determined after molecular analysis does not correspond to the real genotype of the individual under consideration. Virtually every genetic data set includes some erroneous genotypes, but genotyping errors remain a taboo subject in population genetics, even though they might greatly bias the final conclusions, especially for studies based on individual identification. Here, we consider four case studies representing a large variety of population genetics investigations differing in their sampling strategies (noninvasive or traditional), in the type of organism studied (plant or animal) and the molecular markers used [microsatellites or amplified fragment length polymorphisms (AFLPs)]. In these data sets, the estimated genotyping error rate ranges from 0.8% for microsatellite loci from bear tissues to 2.6% for AFLP loci from dwarf birch leaves. Main sources of errors were allelic dropouts for microsatellites and differences in peak intensities for AFLPs, but in both cases human factors were non‐negligible error generators. Therefore, tracking genotyping errors and identifying their causes are necessary to clean up the data sets and validate the final results according to the precision required. In addition, we propose the outline of a protocol designed to limit and quantify genotyping errors at each step of the genotyping process. In particular, we recommend (i) several efficient precautions to prevent contaminations and technical artefacts; (ii) systematic use of blind samples and automation; (iii) experience and rigor for laboratory work and scoring; and (iv) systematic reporting of the error rate in population genetics studies.

Genotyping errors: causes, consequences and solutions

Although genotyping errors affect most data and can markedly influence the biological conclusions of a study, they are too often neglected. Errors have various causes, but their occurrence and effect can be limited by considering these causes in the production and analysis of the data. Procedures that have been developed for dealing with errors in linkage studies, forensic analyses and non-invasive genotyping should be applied more broadly to any genetic study. We propose a protocol for estimating error rates and recommend that these measures be systemically reported to attest the reliability of published genotyping studies.

New perspectives in diet analysis based on DNA barcoding and parallel pyrosequencing: the trnL approach

The development of DNA barcoding (species identification using a standardized DNA sequence), and the availability of recent DNA sequencing techniques offer new possibilities in diet analysis. DNA fragments shorter than 100–150 bp remain in a much higher proportion in degraded DNA samples and can be recovered from faeces. As a consequence, by using universal primers that amplify a very short but informative DNA fragment, it is possible to reliably identify the plant taxon that has been eaten. According to our experience and using this identification system, about 50% of the taxa can be identified to species using the trnL approach, that is, using the P6 loop of the chloroplast trnL (UAA) intron. We demonstrated that this new method is fast, simple to implement, and very robust. It can be applied for diet analyses of a wide range of phytophagous species at large scales. We also demonstrated that our approach is efficient for mammals, birds, insects and molluscs. This method opens new perspectives in ecology, not only by allowing large‐scale studies on diet, but also by enhancing studies on resource partitioning among competing species, and describing food webs in ecosystems.

Fifty thousand years of Arctic vegetation and megafaunal diet

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25–15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding.

Global biodiversity in freshwater and the oceans is declining at high rates. Reliable tools for assessing and monitoring aquatic biodiversity, especially for rare and secretive species, are important for efficient and timely management. Recent advances in DNA sequencing have provided a new tool for species detection from DNA present in the environment. In this study, we tested whether an environmental DNA (eDNA) metabarcoding approach, using water samples, can be used for addressing significant questions in ecology and conservation. Two key aquatic vertebrate groups were targeted: amphibians and bony fish. The reliability of this method was cautiously validated in silico, in vitro and in situ. When compared with traditional surveys or historical data, eDNA metabarcoding showed a much better detection probability overall. For amphibians, the detection probability with eDNA metabarcoding was 0.97 (CI = 0.90–0.99) vs. 0.58 (CI = 0.50–0.63) for traditional surveys. For fish, in 89% of the studied sites, the number of taxa detected using the eDNA metabarcoding approach was higher or identical to the number detected using traditional methods. We argue that the proposed DNA‐based approach has the potential to become the next‐generation tool for ecological studies and standardized biodiversity monitoring in a wide range of aquatic ecosystems.

A Multiplex Pre-Amplification Method that Significantly Improves Microsatellite Amplification and Error Rates for Faecal DNA in Limiting Conditions

Maxine P. Piggott*, Eva Bellemain, Pierre Taberlet & Andrea C. Taylor School of Biological Sciences, Monash University, Victoria, 3800, Australia; Laboratoire d’Ecologie Alpine, Génomique des Populations et Biodiversité, CNRS UMR 5553, Université Joseph Fourier, B.P. 53, 38041, Grenoble, Cedex 9, France and Department of Ecology and Natural Resource Management, Agricultural University of Norway, Postbox 5003, NO-1432 AS, Norway; Laboratoire d’Ecologie Alpine, Génomique des Populations et Biodiversité, CNRS UMR 5553, Université Joseph Fourier, B.P. 53, 38041, Grenoble, Cedex 9, France (*Author for correspondence: fax: +61-3-99055613; e-mail: Maxine.Piggott @sci.monash.edu.au)

Kin-related spatial structure in brown bears Ursus arctos

Kin-related social structure may influence reproductive success and survival and, hence, the dynamics of populations. It has been documented in many gregarious animal populations, but few solitary species. Using molecular methods and field data we tested: (1) whether kin-related spatial structure exists in the brown bear (Ursus arctos), which is a solitary carnivore, (2) whether home ranges of adult female kin overlap more than those of nonkin, and (3) whether multigenerational matrilinear assemblages, i.e., aggregated related females, are formed. Pairwise genetic relatedness between adult (5 years and older) female dyads declined significantly with geographic distance, whereas this was not the case for male–male dyads or opposite sex dyads. The amount of overlap of multiannual home ranges was positively associated with relatedness among adult females. This structure within matrilines is probably due to kin recognition. Plotting of multiannual home-range centers of adult females revealed formation of two types of matrilines, matrilinear assemblages exclusively using an area and dispersed matrilines spread over larger geographic areas. The variation in matrilinear structure might be due to differences in competitive abilities among females and habitat limitations. The influence of kin-related spatial structure on inclusive fitness needs to be clarified in solitary mammals.

New environmental metabarcodes for analysing soil DNA: potential for studying past and present ecosystems

Metabarcoding approaches use total and typically degraded DNA from environmental samples to analyse biotic assemblages and can potentially be carried out for any kinds of organisms in an ecosystem. These analyses rely on specific markers, here called metabarcodes, which should be optimized for taxonomic resolution, minimal bias in amplification of the target organism group and short sequence length. Using bioinformatic tools, we developed metabarcodes for several groups of organisms: fungi, bryophytes, enchytraeids, beetles and birds. The ability of these metabarcodes to amplify the target groups was systematically evaluated by (i) in silico PCRs using all standard sequences in the EMBL public database as templates, (ii) in vitro PCRs of DNA extracts from surface soil samples from a site in Varanger, northern Norway and (iii) in vitro PCRs of DNA extracts from permanently frozen sediment samples of late‐Pleistocene age (∼16 000–50 000 years bp) from two Siberian sites, Duvanny Yar and Main River. Comparison of the results from the in silico PCR with those obtained in vitro showed that the in silico approach offered a reliable estimate of the suitability of a marker. All target groups were detected in the environmental DNA, but we found large variation in the level of detection among the groups and between modern and ancient samples. Success rates for the Pleistocene samples were highest for fungal DNA, whereas bryophyte, beetle and bird sequences could also be retrieved, but to a much lesser degree. The metabarcoding approach has considerable potential for biodiversity screening of modern samples and also as a palaeoecological tool.

Estimating population size and trends of the Swedish brown bear Ursus arctos population

Abstract Estimating population size and trends are key issues in the conservation and management of large carnivores. The rebounding brown bear Ursus arctos population in Sweden is monitored by two different systems, both relying on voluntary resources. Population estimates have been calculated using Capture-Mark-Recapture methods, based on DNA-based scat surveys in five of the six Swedish counties with established bear populations. A total of 1,358 genotypes were identified using DNA extracted from collected scats. An independent ongoing programme, the Large Carnivore Observation Index (LCOI), was initiated in 1998. The LCOI uses effort-corrected observations of bears by moose Alces alces hunters during the moose hunt (> 2 million observation hours/year) and has shown a good correlation with relative population density of bears using the DNA-based method. From this, we have calculated population trends during the period 1998-2007. Using an exponential model, we estimated the yearly increase in the bear population to be 4.5% at the national level, varying between 0 and 10.2% in different counties. We used the regional population estimates and the trends from the LCOI, taking the variation from both systems into account using parametric bootstrapping, to calculate the regional as well as the national population size in Sweden in fall 2008. In one case (the northernmost county; Norrbotten) a DNA-scat survey was lacking, so we used assumptions based on data from the neighbouring county to estimate population size. We estimated the Swedish brown bear population to be 3,298 individuals (2,968-3,667; 95% confidence intervals) in 2008. Our results suggest that reliable information, necessary for the management of the brown bear population can be obtained from volunteers using standardised methods.

A new individual-based spatial approach for identifying genetic discontinuities in natural populations

The population concept is central in evolutionary and conservation biology, but identifying the boundaries of natural populations is often challenging. Here, we present a new approach for assessing spatial genetic structure without the a priori assumptions on the locations of populations made by adopting an individual‐centred approach. Our method is based on assignment tests applied in a moving window over an extensively sampled study area. For each individual, a spatially explicit probability surface is constructed, showing the estimated probability of finding its multilocus genotype across the landscape, and identifying putative migrants. Population boundaries are localized by estimating the mean slope of these probability surfaces over all individuals to identify areas with genetic discontinuities. The significance of the genetic discontinuities is assessed by permutation tests. This new approach has the potential to reveal cryptic population structure and to improve our ability to understand gene flow dynamics across landscapes. We illustrate our approach by simulations and by analysing two empirical datasets: microsatellite data of Ursus arctos in Scandinavia, and amplified fragment length polymorphism (AFLP) data of Rhododendron ferrugineum in the Alps.