Noninvasive genetics based on microsatellite markers has become an indispensable tool for wildlife monitoring and conservation research over the past decades. However, microsatellites have several drawbacks, such as the lack of standardisation between laboratories and high error rates. Here, we propose an alternative single-nucleotide polymorphism (SNP)-based marker system for noninvasively collected samples, which promises to solve these problems. Using nanofluidic SNP genotyping technology (Fluidigm), we genotyped 158 wolf samples (tissue, scats, hairs, urine) for 192 SNP loci selected from the Affymetrix v2 Canine SNP Array. We carefully selected an optimised final set of 96 SNPs (and discarded the worse half), based on assay performance and reliability. We found rates of missing data in this SNP set of <10% and genotyping error of~1%, which improves genotyping accuracy by nearly an order of magnitude when compared to published data for other marker types. Our approach provides a tool for rapid and cost-effective genotyping of noninvasively collected wildlife samples. The ability to standardise genotype scoring combined with low error rates promises to constitute a major technological advancement and could establish SNPs as a standard marker for future wildlife monitoring.

By combining high-throughput sequencing with target enrichment ('hybridization capture'), researchers are able to obtain molecular data from genomic regions of interest for projects that are otherwise constrained by sample quality (e.g. degraded and contamination-rich samples) or a lack of a priori sequence information (e.g. studies on nonmodel species). Despite the use of hybridization capture in various fields of research for many years, the impact of enrichment conditions on capture success is not yet thoroughly understood. We evaluated the impact of a key parameter--hybridization temperature--on the capture success of mitochondrial genomes across the carnivoran family Felidae. Capture was carried out for a range of sample types (fresh, archival, ancient) with varying levels of sequence divergence between bait and target (i.e. across a range of species) using pools of individually indexed libraries on Agilent SureSelect(™) arrays. Our results suggest that hybridization capture protocols require specific optimization for the sample type that is being investigated. Hybridization temperature affected the proportion of on-target sequences following capture: for degraded samples, we obtained the best results with a hybridization temperature of 65 °C, while a touchdown approach (65 °C down to 50 °C) yielded the best results for fresh samples. Evaluation of capture performance at a regional scale (sliding window approach) revealed no significant improvement in the recovery of DNA fragments with high sequence divergence from the bait at any of the tested hybridization temperatures, suggesting that hybridization temperature may not be the critical parameter for the enrichment of divergent fragments.

The colonization history of Madeiran house mice was investigated by analysing the complete mitochondrial (mt) D-loop sequences of 156 mice from the island of Madeira and mainland Portugal, extending on previous studies. The numbers of mtDNA haplotypes from Madeira and mainland Portugal were substantially increased (17 and 14 new haplotypes respectively), and phylogenetic analysis confirmed the previously reported link between the Madeiran archipelago and northern Europe. Sequence analysis revealed the presence of four mtDNA lineages in mainland Portugal, of which one was particularly common and widespread (termed the 'Portugal Main Clade'). There was no support for population bottlenecks during the formation of the six Robertsonian chromosome races on the island of Madeira, and D-loop sequence variation was not found to be structured according to karyotype. The colonization time of the Madeiran archipelago by Mus musculus domesticus was approached using two molecular dating methods (mismatch distribution and Bayesian skyline plot). Time estimates based on D-loop sequence variation at mainland sites (including previously published data from France and Turkey) were evaluated in the context of the zooarchaeological record of M. m. domesticus. A range of values for mutation rate (mu) and number of mouse generations per year was considered in these analyses because of the uncertainty surrounding these two parameters. The colonization of Portugal and Madeira by house mice is discussed in the context of the best-supported parameter values. In keeping with recent studies, our results suggest that mutation rate estimates based on interspecific divergence lead to gross overestimates concerning the timing of recent within-species events.

No abstract

Determining species distributions can be extremely challenging but is crucial to ecological and conservation research. Environmental DNA (eDNA) approaches have shown particular promise in aquatic systems for several vertebrate and invertebrate species. For terrestrial animals, however, eDNA‐based surveys are considerably more difficult due to the lack of or difficulty in obtaining appropriate sampling substrate. In water‐limited ecosystem where terrestrial mammals are often forced to congregate at waterholes, water and sediment from shared water sources may be a suitable substrate for noninvasive eDNA approaches. We characterized mitochondrial DNA sequences from a broad range of terrestrial mammal species in two different African ecosystems (in Namibia and Tanzania) using eDNA isolated from native water, sediment and water filtered through glass fibre filters. A hybridization capture enrichment with RNA probes targeting the mitochondrial genomes of 38 mammal species representing the genera/families expected at the respective ecosystems was employed, and 16 species were identified, with a maximum mitogenome coverage of 99.8%. Conventional genus‐specific PCRs were tested on environmental samples for two genera producing fewer positive results than hybridization capture enrichment. An experiment with mock samples using DNA from non‐African mammals showed that baits covering 30% of nontarget mitogenomes produced 91% mitogenome coverage after capture. In the mock samples, over‐representation of DNA of one species still allowed for the detection of DNA of other species that was at a 100‐fold lower concentration. Hybridization capture enrichment of eDNA is therefore an effective method for monitoring terrestrial mammal species from shared water sources.

The phylogeographic architecture of the common vole, Microtus arvalis, has been well-studied using mitochondrial DNA and used to test hypotheses relating to glacial refugia. The distribution of the five described cytochrome b (cyt b) lineages in Europe west of Russia has been interpreted as a consequence of postglacial expansion from both southern and central European refugia. A recently proposed competing model suggests that the 'cradle' of the M. arvalis lineages is in western central Europe from where they dispersed in different directions after the Last Glacial Maximum. In the present study, we report a new cyt b lineage of the common vole from the Balkans that is not closely related to any other lineage and whose presence might help resolve these issues of glacial refugia. The Balkan phylogroup occurs along the southern distributional border of M. arvalis in central and eastern Bosnia and Herzegovina, Montenegro, and eastern Serbia. Further north and west in Slovenia, Bosnia and Herzegovina, and Serbia, common voles belong to the previously-described Eastern lineage, whereas both lineages are sympatric in one site in Bosnia and Herzegovina. The Balkan phylogroup most reasonably occupied a glacial refugium already known for various Balkan endemic species, in contrast to the recently proposed model. South-east Europe is an absolutely crucial area for understanding the postglacial colonization history of small mammals in Europe and the present study adds to the very few previous detailed phylogeographic studies of this region.

The Leopard cat Prionailurus bengalensis is a habitat generalist that is widely distributed across Southeast Asia. Based on morphological traits, this species has been subdivided into 12 subspecies. Thus far, there have been few molecular studies investigating intraspecific variation, and those had been limited in geographic scope. For this reason, we aimed to study the genetic structure and evolutionary history of this species across its very large distribution range in Asia. We employed both PCR-based (short mtDNA fragments, 94 samples) and high throughput sequencing based methods (whole mitochondrial genomes, 52 samples) on archival, noninvasively collected and fresh samples to investigate the distribution of intraspecific genetic variation. Our comprehensive sampling coupled with the improved resolution of a mitochondrial genome analyses provided strong support for a deep split between Mainland and Sundaic Leopard cats. Although we identified multiple haplogroups within the species' distribution, we found no matrilineal evidence for the distinction of 12 subspecies. In the context of Leopard cat biogeography, we cautiously recommend a revision of the Prionailurus bengalensis subspecific taxonomy: namely, a reduction to 4 subspecies (2 mainland and 2 Sundaic forms).

BackgroundResolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts?ResultsIn this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these.ConclusionsThe phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.Electronic supplementary materialThe online version of this article (10.1186/s12862-018-1268-0) contains supplementary material, which is available to authorized users.