Evaluating the interaction of faecal pellet deposition rates and DNA degradation rates to optimize sampling design for DNA-based mark-recapture analysis of Sonoran pronghorn.

Knowledge of population demographics is important for species management but can be challenging in low-density, wide-ranging species. Population monitoring of the endangered Sonoran pronghorn (Antilocapra americana sonoriensis) is critical for assessing the success of recovery efforts, and noninvasive DNA sampling (NDS) could be more cost-effective and less intrusive than traditional methods. We evaluated faecal pellet deposition rates and faecal DNA degradation rates to maximize sampling efficiency for DNA-based mark-recapture analyses. Deposition data were collected at five watering holes using sampling intervals of 1-7 days and averaged one pellet pile per pronghorn per day. To evaluate nuclear DNA (nDNA) degradation, 20 faecal samples were exposed to local environmental conditions and sampled at eight time points from one to 124 days. Average amplification success rates for six nDNA microsatellite loci were 81% for samples on day one, 63% by day seven, 2% by day 14 and 0% by day 60. We evaluated the efficiency of different sampling intervals (1-10 days) by estimating the number of successful samples, success rate of individual identification and laboratory costs per successful sample. Cost per successful sample increased and success and efficiency declined as the sampling interval increased. Results indicate NDS of faecal pellets is a feasible method for individual identification, population estimation and demographic monitoring of Sonoran pronghorn. We recommend collecting samples <7 days old and estimate that a sampling interval of four to seven days in summer conditions (i.e., extreme heat and exposure to UV light) will achieve desired sample sizes for mark-recapture analysis while also maximizing efficiency [Corrected].

A method for estimating population sex ratio for sage-grouse using noninvasive genetic samples.

Population sex ratio is an important metric for wildlife management and conservation, but estimates can be difficult to obtain, particularly for sexually monomorphic species or for species that differ in detection probability between the sexes. Noninvasive genetic sampling (NGS) using polymerase chain reaction (PCR) has become a common method for identifying sex from sources such as hair, feathers or faeces, and is a potential source for estimating sex ratio. If, however, PCR success is sex-biased, naively using NGS could lead to a biased sex ratio estimator. We measured PCR success rates and error rates for amplifying the W and Z chromosomes from greater sage-grouse (Centrocercus urophasianus) faecal samples, examined how success and error rates for sex identification changed in response to faecal sample exposure time, and used simulation models to evaluate precision and bias of three sex assignment criteria for estimating population sex ratio with variable sample sizes and levels of PCR replication. We found PCR success rates were higher for females than males and that choice of sex assignment criteria influenced the bias and precision of corresponding sex ratio estimates. Our simulations demonstrate the importance of considering the interplay between the sex bias of PCR success, number of genotyping replicates, sample size, true population sex ratio and accuracy of assignment rules for designing future studies. Our results suggest that using faecal DNA for estimating the sex ratio of sage-grouse populations has great potential and, with minor adaptations and similar marker evaluations, should be applicable to numerous species.

Using a reference population yardstick to calibrate and compare genetic diversity reported in different studies: an example from the brown bear.

In species with large geographic ranges, genetic diversity of different populations may be well studied, but differences in loci and sample sizes can make the results of different studies difficult to compare. Yet, such comparisons are important for assessing the status of populations of conservation concern. We propose a simple approach of using a single well-studied reference population as a 'yardstick' to calibrate results of different studies to the same scale, enabling comparisons. We use a well-studied large carnivore, the brown bear (Ursus arctos), as a case study to demonstrate the approach. As a reference population, we genotyped 513 brown bears from Slovenia using 20 polymorphic microsatellite loci. We used this data set to calibrate and compare heterozygosity and allelic richness for 30 brown bear populations from 10 different studies across the global distribution of the species. The simplicity of the reference population approach makes it useful for other species, enabling comparisons of genetic diversity estimates between previously incompatible studies and improving our understanding of how genetic diversity is distributed throughout a species range.

The power of genetic monitoring for studying demography, ecology and genetics of a reintroduced brown bear population.

Genetic monitoring has rarely been used for wildlife translocations despite the potential benefits this approach offers, compared to traditional field-based methods. We applied genetic monitoring to the reintroduced brown bear population in northern Italy. From 2002 to 2008, 2781 hair and faecal samples collected noninvasively plus 12 samples obtained from captured or dead bears were used to follow the demographic and geographical expansion and changes in genetic composition. Individual genotypes were used to reconstruct the wild pedigree and revealed that the population increased rapidly, from nine founders to >27 individuals in 2008 (lambda=1.17-1.19). Spatial mapping of bear samples indicated that most bears were distributed in the region surrounding the translocation site; however, individual bears were found up to 163 km away. Genetic diversity in the population was high, with expected heterozygosity of 0.74-0.79 and allelic richness of 4.55-5.41. However, multi-year genetic monitoring data showed that mortality rates were elevated, immigration did not occur, one dominant male sired all cubs born from 2002 to 2005, genetic diversity declined, relatedness increased, inbreeding occurred, and the effective population size was extremely small (Ne=3.03, ecological method). The comprehensive information collected through genetic monitoring is critical for implementing future conservation plans for the brown bear population in the Italian Alps. This study provides a model for other reintroduction programmes by demonstrating how genetic monitoring can be implemented to uncover aspects of the demography, ecology and genetics of small and reintroduced populations that will advance our understanding of the processes influencing their viability, evolution, and successful restoration.

Comparative landscape genetics of two pond-breeding amphibian species in a highly modified agricultural landscape.

Evaluating fine-scale population structure of multiple species in the same landscape increases our ability to identify common patterns as well as discern ecological differences among species' landscape genetic relationships. In the Palouse bioregion of northern Idaho, USA, 99% of the native prairie has been converted to nonirrigated agriculture and exotic grasslands. Columbia spotted frogs (Rana luteiventris) and long-toed salamanders (Ambystoma macrodactylum) in this area breed almost entirely in artificial ponds on private land. We used genetic distances (F(ST) and D(c)) derived from eight microsatellite loci in 783 samples to evaluate the relationships among sympatric breeding populations (N = 20 and 26) of these species in a 213-km(2) landscape. Both species showed a pattern of isolation by distance that was not improved when distance was measured along drainages instead of topographically corrected straight lines (P < 0.01). After testing for autocorrelation among genetic distances, we used an information theoretic approach to model landscape resistance based on slope, soil type, solar insolation, and land cover, and multi-model inference to rank the resistance of landscape surfaces to dispersal (represented by genetic distance). For both species, urban and rural developed land cover provided the highest landscape resistances. Resistance values for long-toed salamanders followed a moisture gradient where forest provided the least resistance, while agriculture and shrub/clearcut provided the least resistance for Columbia spotted frogs. Comparative landscape genetics can be a powerful tool for detecting similarities and differences between codistributed species, and resulting models can be used to predict species-specific responses to landscape change.

Impacts of sampling location within a faeces on DNA quality in two carnivore species.

We investigated the influence of sampling location within a faeces on DNA quality by sampling from both the outside and inside of 25 brown bear (Ursus arctos) scats and the side and the tip of 30 grey wolf (Canis lupus) scats. The outside of the bear scat and side of the wolf scat had significantly lower nuclear DNA microsatellite allelic dropout error rates (U. arctos: P = 0.017; C. lupus: P = 0.025) and significantly higher finalized genotyping success rates (U. arctos: P = 0.017; C. lupus: P = 0.012) than the tip and inside of the scat. A review of the faecal DNA literature indicated that <45% of studies report the sampling location within a faeces indicating that this methodological consideration is currently underappreciated. Based on our results, we recommend sampling from the side of canid scats and the outside portion of ursid scats to obtain higher quality DNA samples. The sampling location within a faeces should be carefully considered and reported as it can directly influence laboratory costs and efficiency, as well as the ability to obtain reliable genotypes.

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.

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.

Putting the "landscape" in landscape genetics.

Landscape genetics has emerged as a new research area that integrates population genetics, landscape ecology and spatial statistics. Researchers in this field can combine the high resolution of genetic markers with spatial data and a variety of statistical methods to evaluate the role that landscape variables play in shaping genetic diversity and population structure. While interest in this research area is growing rapidly, our ability to fully utilize landscape data, test explicit hypotheses and truly integrate these diverse disciplines has lagged behind. Part of the current challenge in the development of the field of landscape genetics is bridging the communication and knowledge gap between these highly specific and technical disciplines. The goal of this review is to help bridge this gap by exposing geneticists to terminology, sampling methods and analysis techniques widely used in landscape ecology and spatial statistics but rarely addressed in the genetics literature. We offer a definition for the term "landscape genetics", provide an overview of the landscape genetics literature, give guidelines for appropriate sampling design and useful analysis techniques, and discuss future directions in the field. We hope, this review will stimulate increased dialog and enhance interdisciplinary collaborations advancing this exciting new field.

Phylogeography and mitochondrial diversity of extirpated brown bear (Ursus arctos) populations in the contiguous United States and Mexico.

The fossil record indicates that the brown bear (Ursus arctos) colonized North America from Asia over 50 000 years ago. The species historically occupied the western United States and northern Mexico but has been extirpated from over 99% of this range in the last two centuries. To evaluate colonization hypotheses, subspecific classifications, and historical patterns and levels of genetic diversity in this region, we sequenced 229 nucleotides of the mitochondrial DNA control region in 108 museum specimens. The work was set in a global context by synthesizing all previous brown bear control region sequences from around the world. In mid-latitude North America a single moderately diverse clade is observed, represented by 23 haplotypes with up to 3.5% divergence. Only eight of 23 haplotypes (35%) are observed in the extensively sampled extant populations suggesting a substantial loss of genetic variability. The restriction of all haplotypes from mid-latitude North America to a single clade suggests that this region was founded by bears with a similar maternal ancestry. However, the levels and distributions of diversity also suggest that the colonizing population was not a small founder event, and that expansion occurred long enough ago for local mutations to accrue. Our data are consistent with recent genetic evidence that brown bears were south of the ice prior to the last glacial maximum. There is no support for previous subspecies designations, although bears of the southwestern United States may have had a distinctive, but recent, pattern of ancestry.

Assessing population structure and gene flow in Montana wolverines (Gulo gulo) using assignment-based approaches.

In North America, wolverines once occupied a continuous range from Alaska southward to New Mexico. In the lower 48 states, small remnant populations remain only in the northwestern United States. Among these remnant populations, the Montana population has the highest probability of long-term persistence given its size and proximity to healthy populations in Canada. In this study, we evaluate population genetic structure and gene flow among Montana wolverines using 10 polymorphic microsatellite loci. Bayesian and frequency-based assignment tests revealed significant population substructure and provide support for at least three subpopulations in Montana. FST values between subpopulations ranged from 0.08 to 0.10 and provide evidence for male-biased dispersal. The high degree of population substructure and low levels of gene flow contrast results from wolverine population genetic studies in less fragmented landscapes of Alaska and Canada. This study provides additional support for the hypothesis that large carnivore populations of Montana are becoming increasingly fragmented due to human development and disturbance.

Using faecal DNA sampling and GIS to monitor hybridization between red wolves (Canis rufus) and coyotes (Canis latrans).

The US Fish and Wildlife Service's (USFWS) Red Wolf Recovery Program recognizes hybridization with coyotes as the primary threat to red wolf recovery. Efforts to curb or stop hybridization are hampered in two ways. First, hybrid individuals are difficult to identify based solely on morphology. Second, managers need to effectively search 6000 km(2) for the presence of coyotes and hybrids. We develop a noninvasive method to screen large geographical areas for coyotes and hybrids with maternal coyote ancestry by combining mitochondrial DNA sequence analysis of faeces (scat) and geographic information system (GIS) technology. This method was implemented on the Alligator River National Wildlife Refuge (1000 km(2)) in northeastern North Carolina. A total of 956 scats were collected in the spring of 2000 and 2001 and global positioning system (GPS) coordinates were recorded. Seventy-five percent of the scats were assigned to species and five coyote/hybrid scats were detected. Placement of scat location coordinates on a map of the experimental population area revealed that four of the coyote/hybrid scats were detected within the home ranges of sterilized hybrids. The other coyote/hybrid scat indicated the presence of a previously unknown individual. We suggest this method be expanded to include more of the experimental population area and be optimized for use with nuclear markers to improve detection of hybrid and back-crossed individuals.

Widespread occurrence of a domestic dog mitochondrial DNA haplotype in southeastern US coyotes.

Sequence analysis of the mitochondrial DNA control region from 112 southeastern US coyotes (Canis latrans) revealed 12 individuals with a haplotype closely related to those in domestic dogs. Phylogenetic analyses grouped this new haplotype in the dog/grey wolf (Canis familiaris/Canis lupus) clade with 98% bootstrap support. These results demonstrate that a male coyote hybridized with a female dog, and female hybrid offspring successfully integrated into the coyote population. The widespread distribution of this haplotype from Florida to West Virginia suggests that the hybridization event occurred long ago before the southeastern USA was colonized by coyotes. However, it could have occurred in the southeastern USA before the main front of coyotes arrived in the area between male coyotes released for sport and a local domestic dog. The introgression of domestic dog genes into the southeastern coyote population does not appear to have substantially affected the coyote's genetic, morphological, or behavioural integrity. However, our results suggest that, contrary to previous reports, hybridization can occur between domestic and wild canids, even when the latter is relatively abundant. Therefore, hybridization may be a greater threat to the persistence of wild canid populations than previously thought.

Microsatellite analysis of genetic variation among and within Alpine marmot populations in the French Alps.

The genetic structure of the Alpine marmot, Marmota marmota, was studied by an analysis of five polymorphic microsatellite loci. Eight locations were sampled in the French Alps, one from Les Ecrins valley (n = 160), another from La Sassière valley (n = 289) and the six others from the Maurienne valley (n = 139). Information on social group structure was available for both Les Ecrins and La Sassière but not for the other samples. The high levels of genetic diversity observed are at odds with the results obtained using microsatellites, minisatellites and allozymes on Alpine marmots from Germany, Austria and Switzerland. Strong deficits in heterozygotes were found in Les Ecrins and La Sassière. They are caused by a Wahlund effect due to the family structure (i.e. differentiation between the family groups). The family groups exhibit excess of heterozygotes rather than deficits. This may be caused by outbreeding and this is compatible with recent results from the genetics of related social species when information on the social structure is taken into account. The observed outbreeding could be the result of females mating with transient males or males coming from neighbouring colonies. Both indicate that the species may not be as monogamous as is usually believed. The results are also compatible with a male-biased dispersal but do not allow us to exclude some female migration. We also found a significant correlation between geographical and genetic distance indicating that isolation by distance could be an issue in marmots. This study is the first that analysed populations of marmots taking into account the social structure within populations and assessing inbreeding at different levels (region, valley, population, and family groups). Our study clearly demonstrated that the sampling strategy and behavioural information can have dramatic effects on both the results and interpretation of the genetic data.

Estimating the probability of identity among genotypes in natural populations: cautions and guidelines.

Individual identification using DNA fingerprinting methods is emerging as a critical tool in conservation genetics and molecular ecology. Statistical methods that estimate the probability of sampling identical genotypes using theoretical equations generally assume random associations between alleles within and among loci. These calculations are probably inaccurate for many animal and plant populations due to population substructure. We evaluated the accuracy of a probability of identity (P(ID)) estimation by comparing the observed and expected P(ID), using large nuclear DNA microsatellite data sets from three endangered species: the grey wolf (Canis lupus), the brown bear (Ursus arctos), and the Australian northern hairy-nosed wombat (Lasiorinyus krefftii). The theoretical estimates of P(ID) were consistently lower than the observed P(ID), and can differ by as much as three orders of magnitude. To help researchers and managers avoid potential problems associated with this bias, we introduce an equation for P(ID) between sibs. This equation provides an estimator that can be used as a conservative upper bound for the probability of observing identical multilocus genotypes between two individuals sampled from a population. We suggest computing the actual observed P(ID) when possible and give general guidelines for the number of codominant and dominant marker loci required to achieve a reasonably low P(ID) (e.g. 0.01-0.0001).

Rapid radiation events in the family Ursidae indicated by likelihood phylogenetic estimation from multiple fragments of mtDNA.

The bear family (Ursidae) presents a number of phylogenetic ambiguities as the evolutionary relationships of the six youngest members (ursine bears) are largely unresolved. Recent mitochondrial DNA analyses have produced conflicting results with respect to the phylogeny of ursine bears. In an attempt to resolve these issues, we obtained 1916 nucleotides of mitochondrial DNA sequence data from six gene segments for all eight bear species and conducted maximum likelihood and maximum parsimony analyses on all fragments separately and combined. All six single-region gene trees gave different phylogenetic estimates; however, only for control region data was this significantly incongruent with the results from the combined data. The optimal phylogeny for the combined data set suggests that the giant panda is most basal followed by the spectacled bear. The sloth bear is the basal ursine bear, and there is weak support for a sister taxon relationship of the American and Asiatic black bears. The sun bear is sister taxon to the youngest clade containing brown bears and polar bears. Statistical analyses of alternate hypotheses revealed a lack of strong support for many of the relationships. We suggest that the difficulties surrounding the resolution of the evolutionary relationships of the Ursidae are linked to the existence of sequential rapid radiation events in bear evolution. Thus, unresolved branching orders during these time periods may represent an accurate representation of the evolutionary history of bear species.

Plucked hair samples as a source of DNA: reliability of dinucleotide microsatellite genotyping.

To test whether plucked hairs are a reliable source of DNA for genotyping microsatellite loci, we carried out experiments using one, three, or 10 hairs per extract for 50 alpine marmots. For each extract, seven independent genotypings were performed for the same locus (multiple-tubes approach). Two types of genotyping errors were recorded: a false homozygote defined as the detection of only one allele of a true heterozygote, and a false allele defined as a PCR-generated allele that was not one of the alleles of the true genotype. Using DNA extracted from one, three, or 10 hairs the overall error rate was 14.00%, 4.86%, and 0.29%, respectively. Based on our results, we conclude that 10 hairs should be used to obtain consistently reliable genotypings using the single-tube approach, and that a single plucked hair could represent a reliable source of DNA if the multiple-tubes approach is used. For future studies of dinucleotide repeat diversity using DNA extracted from one to three shed or plucked hairs, we strongly recommend initiating an appropriate pilot study to quantify the error rate and to determine the reliability of the single-tube approach.

An empirical evaluation of genetic distance statistics using microsatellite data from bear (Ursidae) populations.

A large microsatellite data set from three species of bear (Ursidae) was used to empirically test the performance of six genetic distance measures in resolving relationships at a variety of scales ranging from adjacent areas in a continuous distribution to species that diverged several million years ago. At the finest scale, while some distance measures performed extremely well, statistics developed specifically to accommodate the mutational processes of microsatellites performed relatively poorly, presumably because of the relatively higher variance of these statistics. At the other extreme, no statistic was able to resolve the close sister relationship of polar bears and brown bears from more distantly related pairs of species. This failure is most likely due to constraints on allele distributions at microsatellite loci. At intermediate scales, both within continuous distributions and in comparisons to insular populations of late Pleistocene origin, it was not possible to define the point where linearity was lost for each of the statistics, except that it is clearly lost after relatively short periods of independent evolution. All of the statistics were affected by the amount of genetic diversity within the populations being compared, significantly complicating the interpretation of genetic distance data.

Noninvasive genetic tracking of the endangered Pyrenean brown bear population.

Pyrenean brown bears Ursus arctos are threatened with extinction. Management efforts to preserve this population require a comprehensive knowledge of the number and sex of the remaining individuals and their respective home ranges. This goal has been achieved using a combination of noninvasive genetic sampling of hair and faeces collected in the field and corresponding track size data. Genotypic data were collected at 24 microsatellite loci using a rigorous multiple-tubes approach to avoid genotyping errors associated with low quantities of DNA. Based on field and genetic data, the Pyrenean population was shown to be composed at least of one yearling, three adult males, and one adult female. These data indicate that extinction of the Pyrenean brown bear population is imminent without population augmentation. To preserve the remaining Pyrenean gene pool and increase genetic diversity, we suggest that managers consider population augmentation using only females. This study demonstrates that comprehensive knowledge of endangered small populations of mammals can be obtained using noninvasive genetic sampling.