Generation of genome-wide SNP markers from minimally invasive sampling in endangered animals and applications in species ecology and conservation

High-density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high-quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA.