Towards high–throughput analyses of fecal samples from wildlife

High–throughput sequencing offers new possibilities in molecular ecology and conservation studies. However, its potential has not yet become fully exploited for noninvasive studies of free–ranging animals, such as those based on feces. High–throughput sequencing allows sequencing of short DNA fragments and could allow simultaneous genotyping of a very large number of samples and markers at a low cost. The application of high throughput genotyping to fecal samples from wildlife has been hindered by several labor–intensive steps. We evaluate alternative protocols which could allow higher throughput for two of these steps: sample collection and DNA extraction. Two different field sampling and seven different DNA extraction methods are tested here on grey wolf (Canis lupus) feces. There was high variation in genotyping success rates. The field sampling method based on surface swabbing performed much worse than the extraction from a fecal fragment. In addition, there is a lot of room for improvement in the DNA extraction step. Optimization of protocols can lead to very much more efficient, cheaper and higher throughput noninvasive monitoring. Selection of appropriate markers is still of paramount importance to increase genotyping success.

Correlates of species richness in the largest Neotropical amphibian radiation

Although tropical environments are often considered biodiversity hotspots, it is precisely in such environments where least is known about the factors that drive species richness. Here, we use phylogenetic comparative analyses to study correlates of species richness for the largest Neotropical amphibian radiation: New World direct-developing frogs. Clade-age and species richness were nonsignificantly, negatively correlated, suggesting that clade age alone does not explain among-clade variation in species richness. A combination of ecological and morphological traits explained 65% of the variance in species richness. A more vascularized ventral skin, the ability to colonize high-altitude ranges, encompassing a large variety of vegetation types, correlated significantly with species richness, whereas larger body size was marginally correlated with species richness. Hence, whereas high-altitude ranges play a role in shaping clade diversity in the Neotropics, intrinsic factors, such as skin structures and possibly body size, might ultimately determine which clades are more speciose than others.

Hybridization between white-headed ducks and introduced ruddy ducks in Spain

The ruddy duck, Oxyura jamaicensis, was introduced to Great Britain in the mid-20th century and has recently spread to other Western European countries. In Spain, ruddy ducks hybridize with the globally endangered white-headed duck, Oxyura leucocephala. We assessed the effects of hybridization on the Spanish white-headed ducks, which constitute 25% of the global population of this species, using a panel of eight nuclear intron markers, 10 microsatellite loci, and mtDNA control region sequences. These data allowed parental individuals, F(1) hybrids, and the progeny of backcrossing to be reliably distinguished. We show that hybrids between the two species are fertile and produce viable offspring in backcrosses with both parental species. To date, however, we found no extensive introgression of ruddy duck genes into the Spanish white-headed duck population, probably due to the early implementation of an effective ruddy duck and hybrid control programme. We also show that genetic diversity in the expanding European ruddy duck population, which was founded by just seven individuals, exceeds that of the native Spanish white-headed duck population, which recently recovered from a severe bottleneck. Unless effective control of ruddy ducks is continued, genetic introgression will compromise the unique behavioural and ecological adaptations of white-headed ducks and consequently their survival as a genetically and evolutionary distinct species.

Genetic diversity, population structure, effective population size and demographic history of the Finnish wolf population

The Finnish wolf population (Canis lupus) was sampled during three different periods (1996-1998, 1999-2001 and 2002-2004), and 118 individuals were genotyped with 10 microsatellite markers. Large genetic variation was found in the population despite a recent demographic bottleneck. No spatial population subdivision was found even though a significant negative relationship between genetic relatedness and geographic distance suggested isolation by distance. Very few individuals did not belong to the local wolf population as determined by assignment analyses, suggesting a low level of immigration in the population. We used the temporal approach and several statistical methods to estimate the variance effective size of the population. All methods gave similar estimates of effective population size, approximately 40 wolves. These estimates were slightly larger than the estimated census size of breeding individuals. A Bayesian model based on Markov chain Monte Carlo simulations indicated strong evidence for a long-term population decline. These results suggest that the contemporary wolf population size is roughly 8% of its historical size, and that the population decline dates back to late 19th century or early 20th century. Despite an increase of over 50% in the census size of the population during the whole study period, there was only weak evidence that the effective population size during the last period was higher than during the first. This may be caused by increased inbreeding, diminished dispersal within the population, and decreased immigration to the population during the last study period.

Unequal contribution of sexes in the origin of dog breeds

Dogs (Canis familiaris) were domesticated from the gray wolf (Canis lupus) at least 14,000 years ago, and there is evidence of dogs with phenotypes similar to those in modern breeds 4000 years ago. However, recent genetic analyses have suggested that modern dog breeds have a much more recent origin, probably <200 years ago. To study the origin of contemporaneous breeds we combined the analysis of paternally inherited Y chromosome markers with maternally inherited mitochondrial DNA and biparentally inherited autosomal microsatellite markers in both domestic dogs and their wild ancestor, the gray wolf. Our results show a sex bias in the origin of breeds, with fewer males than females contributing genetically, which clearly differs from the breeding patterns in wild gray wolf populations where both sexes have similar contributions. Furthermore, a comparison of mitochondrial DNA and Y chromosome diversity in dog groups recognized by the World Canine Organization, as well as in groups defined by the breeds' genetic composition, shows that paternal lineages are more differentiated among groups than maternal lineages. This demonstrates a lower exchange of males than of females between breeds belonging to different groups, which illustrates how breed founders may have been chosen.

Two centuries of the Scandinavian wolf population: patterns of genetic variability and migration during an era of dramatic decline

The grey wolf (Canis lupus) was numerous on the Scandinavian peninsula in the early 19th century. However, as a result of intense persecution, the population declined dramatically and was virtually extinct from the peninsula by the 1960s. We examined historical patterns of genetic variability throughout the period of decline, from 1829 to 1979. Contemporary Finnish wolves, considered to be representative of a large eastern wolf population, were used for comparison. Mitochondrial DNA (mtDNA) variability among historical Scandinavian wolves was significantly lower than in Finland while Y chromosome variability was comparable between the two populations. This may suggest that long-distance migration from the east has been male-biased. Importantly though, as the historical population was significantly differentiated from contemporary Finnish wolves, the overall immigration rate to the Scandinavian peninsula appears to have been low. Levels of variability at autosomal microsatellite loci were high by the early 1800s but declined considerably towards the mid-20th century. At this time, approximately 40% of the allelic diversity and 30% of the heterozygosity had been lost. After 1940, however, there is evidence of several immigration events, coinciding with episodes of marked population increase in Russian Karelia and subsequent westwards migration.

Combined use of maternal, paternal and bi-parental genetic markers for the identification of wolf-dog hybrids

The identification of hybrids is often a subject of primary concern for the development of conservation and management strategies, but can be difficult when the hybridizing species are closely related and do not possess diagnostic genetic markers. However, the combined use of mitochondrial DNA (mtDNA), autosomal and Y chromosome genetic markers may allow the identification of hybrids and of the direction of hybridization. We used these three types of markers to genetically characterize one possible wolf-dog hybrid in the endangered Scandinavian wolf population. We first characterized the variability of mtDNA and Y chromosome markers in Scandinavian wolves as well as in neighboring wolf populations and in dogs. While the mtDNA data suggested that the target sample could correspond to a wolf, its Y chromosome type had not been observed before in Scandinavian wolves. We compared the genotype of the target sample at 18 autosomal microsatellite markers with those expected in pure specimens and in hybrids using assignment tests. The combined results led to the conclusion that the animal was a hybrid between a Scandinavian female wolf and a male dog. This finding confirms that inter-specific hybridization between wolves and dogs can occur in natural wolf populations. A possible correlation between hybridization and wolf population density and disturbance deserves further research.

Y chromosome haplotyping in Scandinavian wolves (Canis lupus) based on microsatellite markers

The analysis of mitochondrial DNA sequences has for a long time been the most extensively used genetic tool for phylogenetic, phylogeographic and population genetic studies. Since this approach only considers female lineages, it tends to give a biased picture of the population history. The use of protein polymorphisms and microsatellites has helped to obtain a more unbiased view, but complementing population genetic studies with Y chromosome markers could clarify the role of each sex in natural processes. In this study we analysed genetic variability at four microsatellite loci on the canid Y chromosome. With these four microsatellites we constructed haplotypes and used them to study the genetic status of the Scandinavian wolf population, a population that now contains 60-70 animals but was thought to have been extinct in the 1970s. In a sample of 100 male wolves from northern Europe we found 17 different Y chromosome haplotypes. Only two of these were found in the current Scandinavian population. This indicates that there should have been at least two males involved in the founding of the Scandinavian wolf population after the bottleneck in the 1970s. The two Scandinavian Y chromosome haplotypes were not found elsewhere in northern Europe, which indicates low male gene flow between Scandinavia and the neighbouring countries.

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.

Widespread origins of domestic horse lineages

Domestication entails control of wild species and is generally regarded as a complex process confined to a restricted area and culture. Previous DNA sequence analyses of several domestic species have suggested only a limited number of origination events. We analyzed mitochondrial DNA (mtDNA) control region sequences of 191 domestic horses and found a high diversity of matrilines. Sequence analysis of equids from archaeological sites and late Pleistocene deposits showed that this diversity was not due to an accelerated mutation rate or an ancient domestication event. Consequently, high mtDNA sequence diversity of horses implies an unprecedented and widespread integration of matrilines and an extensive utilization and taming of wild horses. However, genetic variation at nuclear markers is partitioned among horse breeds and may reflect sex-biased dispersal and breeding.

Tripartite genetic subdivisions in the ornate shrew (Sorex ornatus)

We examined cytochrome b sequence variation in 251 ornate shrews (Sorex ornatus) from 20 localities distributed throughout their geographical range. Additionally, vagrant (S. vagrans) and montane (S. monticolus) shrews from four localities were used as outgroups. We found 24 haplotypes in ornate shrews from California (USA) and Baja California (Mexico) that differed by 1-31 substitutions in 392 bp of mitochondrial DNA (mtDNA) sequence. In a subset of individuals, we sequenced 699 bp of cytochrome b to better resolve the phylogeographic relationships of populations. The ornate shrew is phylogeographically structured into three haplotype clades representing southern, central and northern localities. Analysis of allozyme variation reveals a similar pattern of variation. Several other small California vertebrates have a similar tripartite pattern of genetic subdivision. We suggest that topographic barriers and expansion and contraction of wetland habitats in the central valley during Pleistocene glacial cycles account for these patterns of genetic variation. Remarkably, the northern ornate shrew clade is phylogenetically clustered with another species of shrew suggesting that it may be a unique lowland form of the vagrant shrew that evolved in parallel to their southern California counterparts.

Genetic variation and population structure in Scandinavian wolverine (Gulo gulo) populations

Wolverine (Gulo gulo) numbers in Scandinavia were significantly reduced during the early part of the century as a result of predator removal programmes and hunting. Protective legislation in both Sweden and Norway in the 1960s and 1970s has now resulted in increased wolverine densities in Scandinavia. We report here the development of 15 polymorphic microsatellite markers in wolverine and their use to examine the population sub-structure and genetic variability in free-ranging Scandinavian wolverine populations as well as in a sample of individuals collected before 1970. Significant subdivision between extant populations was discovered, in particular for the small and isolated population of southern Norway, which represents a recent recolonization. Overall genetic variability was found to be lower than previously reported for other mustelids, with only two to five alleles per locus and observed heterozygosities (H(O)) ranging from 0.269 to 0.376 across the examined populations, being lowest in southern Norway. Analysis of the mitochondrial DNA control region revealed no variation throughout the surveyed populations. As the historical sample did not show higher levels of genetic variability, our results are consistent with a reduction in the genetic variation in Scandinavian wolverines that pre-dates the demographic bottleneck observed during the last century. The observed subdivision between populations calls for management caution when issuing harvest quotas, especially for the geographically isolated south Norwegian population.

Phylogenetic relationships, evolution, and genetic diversity of the domestic dog

The spectacular diversity in size, conformation, and pelage that characterizes the domestic dog reflects not only the intensity of artificial selection but ultimately the genetic variability of founding populations. Here we review past molecular genetic data that are relevant to understanding the origin and phylogenetic relationships of the dog. DNA-DNA hybridization data show that the dog family Canidae diverged about 50 million years ago from other carnivore families. In contrast, the extant canids are very closely related and diverged from a common ancestor about 10 million years ago. The evidence supporting a close relationship of dogs with gray wolves is overwhelming. However, dogs are remarkably diverse in mitochondrial and nuclear genes. Mitochondrial DNA analysis suggests a more ancient origin of dogs than has been indicated by the fossil record. In addition, dogs have originated from or interbred with wolves throughout their history at different times and different places. We test the possibility of an independent domestication event in North America by analysis of mtDNA variation in the Xoloitzcuintli. This unusual breed is believed to have been kept isolated for thousands of years and may be one of the most ancient breeds in North America. Our results do not support a New World domestication of dogs nor a close association of the Xoloitzcuintli with other hair-less breeds of dogs. Despite their phenotypic uniformity, the Xoloitzcuintli has a surprisingly high level of mtDNA sequence variation. Other breeds are also genetically diverse, suggesting that dog breeds were often founded with a large number of dogs from outbred populations.

Conservation genetics of the endangered Pampas deer (Ozotoceros bezoarticus)

The Pampas deer (Ozotoceros bezoarticus L. 1758) is the most endangered neotropical cervid, and in the past occupied a wide range of open habitats including grassland, pampas, savanna, and cerrado (Brazil) from 5 degrees to 41 degrees S. To better understand the effect of habitat fragmentation on gene flow and genetic variation, and to uncover genetic units for conservation, we examined DNA sequences from the mitochondrial control region of 54 individuals from six localities distributed throughout the present geographical range of the Pampas deer. Our results suggest that the control region of the Pampas deer is one of the most polymorphic of any mammal. This remarkably high variability probably reflects large historic population sizes of millions of individuals in contrast to numbers of fewer than 80,000 today. Gene flow between populations is generally close to one migrant per generation and, with the exception of two populations from Argentina, all populations are significantly differentiated. The degree of gene flow was correlated with geographical distance between populations, a result consistent with limited dispersal being the primary determinant of genetic differentiation between populations. The molecular genetic results provide a mandate for habitat restoration and reintroduction of Pampas deer so that levels of genetic variation can be preserved and historic patterns of abundance can be reconstructed. However, the source of individuals for reintroduction generally should be from populations geographically closest to those now in danger of extinction.

Multiple and ancient origins of the domestic dog

Mitochondrial DNA control region sequences were analyzed from 162 wolves at 27 localities worldwide and from 140 domestic dogs representing 67 breeds. Sequences from both dogs and wolves showed considerable diversity and supported the hypothesis that wolves were the ancestors of dogs. Most dog sequences belonged to a divergent monophyletic clade sharing no sequences with wolves. The sequence divergence within this clade suggested that dogs originated more than 100,000 years before the present. Associations of dog haplotypes with other wolf lineages indicated episodes of admixture between wolves and dogs. Repeated genetic exchange between dog and wolf populations may have been an important source of variation for artificial selection.

Use of faeces for scent marking in Iberian wolves (Canis lupus)

The study of the characteristics of 232 wolf (Canis lupus) scats from northwestern Spain reveals that faeces were not left at random. They accumulated near trail junctions, where the probability of interception by other traveling wolves is maximal, and were commonly associated with evident reference points (i.e., ashes, decaying bark, isolated shrubs, carrion, promontories, etc.) that may visually or olfactorily increase the signal's effectiveness. The center of the trail was clearly preferred to the borders, facilitating location of the scat by other individuals. If one side of the trail was more exposed to the wind, scats accumulated more on that side than on the other, and their deposition on the trail margin may increase the radius of their influence even more. These results suggest that faeces may be more important for spacing among Iberian wolves than expected, considering the minor role assigned to the scats in most previous studies of wolf scent marking.

Use of tiletamine and zolazepam to immobilize captive Iberian wolves (Canis lupus)

A mixture of tiletamine and zolazepam (Zoletil) was used to immobilize 29 captive born Iberian wolves. Based on their excitability during handling procedures the animals were categorized as excited (n = 15) and unexcited (n = 14). We observed differences in the responses of these groups to the drugs. Although immobilized with higher doses (mean +/- SD of 6.94 +/- 2.13 versus 5.04 +/- 1.74 mg/kg for the unexcited) the excited individuals had an irregular and less predictable response, with five individuals needing additional dosages in the excited group compared to one animal in the unexcited group. Arousal time and duration of immobilization of excited wolves was not correlated to initial drug doses, but was in unexcited animals; the excited group had a poorer thermal regulation. Differences in arousal time and duration could be the a result of the different doses used. Excited wolves were older than unexcited (5.4 +/- 3.07 versus 2.86 +/- 2.11 years, respectively). For captive wolves, doses of about 5 mg/kg are recommended for non-excited and 10 mg/kg for excited individuals.

Tooth losses and anomalies in the wolf (Canis lupus)

We observed tooth losses (and subsequent alveolar resorptions) and other cases of anodonty and polydonty and anomalies in the shape and occlusion of the teeth in skulls of wolves from the former Soviet Union. Resorptions increase in number with age, but do not seem to have much of an effect on individual survival (they probably do not affect hunting), and are found in half of individuals over 8 years of age. The highest numbers of resorptions, as well as cases of polydonty and dental anomalies, were in the premolars. This suggests that the premolars are not greatly subject to selective pressure. At least 6.8% of the skulls have malformed teeth or malocclusions.