A Pipeline and Recommendations for Population and Individual Diagnostic SNP Selection in Non-Model Species

Despite substantial reductions in the cost of sequencing over the last decade, genetic panels remain relevant due to their cost-effectiveness and flexibility across a variety of sample types. In particular, single nucleotide polymorphism (SNP) panels are increasingly favoured for conservation applications. SNP panels are often used because of their adaptability, effectiveness with low-quality samples, and cost-efficiency for population monitoring and forensics. However, the selection of diagnostic SNPs for population assignment and individual identification can be challenging. The consequences of poor SNP selection are under-powered panels, inaccurate results, and monetary loss. Here, we develop a novel and user-friendly SNP selection pipeline (mPCRselect) that can be used to select SNPs for population assignment and/or individual identification. mPCRselect allows any researcher, who has sufficient SNP-level data, to design a successful and cost-effective SNP panel for a diploid species of conservation concern.

Advancements in noninvasive koala monitoring through combining Chlamydia detection with a targeted koala genotyping assay

Wildlife diseases are major players in local and global extinctions. Effective disease surveillance, management and conservation strategies require accurate estimates of pathogen prevalence. Yet pathogen detection in wild animals remains challenging. Current gold standards often require samples collected through veterinary examination, but this method is costly, intensive, invasive, and requires specialised staff and equipment. Collection of non-invasive samples, such as scats, is an effective monitoring tool which can be deployed at large scale, as scats contain DNA of both host and pathogens. The koala (Phascolarctos cinereus) is listed as 'endangered' under the EPBC Act 1999, with chlamydial disease representing a major threat. Here, we present a new approach that combines restriction-enzyme associated sequencing and targeted-sequence-capture genotyping, namely DArTcap, to detect Chlamydia pecorum in koala scats. We found this method has similar accuracy to current gold standards (qPCR of swab samples), with a sensitivity of 91.7% and a specificity of 100%. This method can be incorporated into existing koala genetic studies using marker panels, where population attributes can be estimated alongside C. pecorum presence, using the same scat samples, with the option to add further markers of interest. Such a one-stop-shop panel would considerably reduce processing times and cost.

Species-specific SNP arrays for non-invasive genetic monitoring of a vulnerable bat

Genetic tagging from scats is one of the minimally invasive sampling (MIS) monitoring approaches commonly used to guide management decisions and evaluate conservation efforts. Microsatellite markers have traditionally been used but are prone to genotyping errors. Here, we present a novel method for individual identification in the Threatened ghost bat Macroderma gigas using custom-designed Single Nucleotide Polymorphism (SNP) arrays on the MassARRAY system. We identified 611 informative SNPs from DArTseq data from which three SNP panels (44-50 SNPs per panel) were designed. We applied SNP genotyping and molecular sexing to 209 M. gigas scats collected from seven caves in the Pilbara, Western Australia, employing a two-step genotyping protocol and identifying unique genotypes using a custom-made R package, ScatMatch. Following data cleaning, the average amplification rate was 0.90 ± 0.01 and SNP genotyping errors were low (allelic dropout 0.003 ± 0.000) allowing clustering of scats based on one or fewer allelic mismatches. We identified 19 unique bats (9 confirmed/likely males and 10 confirmed/likely females) from a maternity and multiple transitory roosts, with two male bats detected using roosts, 9 km and 47 m apart. The accuracy of our SNP panels enabled a high level of confidence in the identification of individual bats. Targeted SNP genotyping is a valuable tool for monitoring and tracking of non-model species through a minimally invasive sampling approach.

Evaluating the use of non-invasive hair sampling and ddRAD to characterize populations of endangered species: Application to a peripheral population of the European mink

The application of next-generation sequencing (NGS) to non-invasive samples is one of the most promising methods in conservation genomics, but these types of samples present significant challenges for NGS. The European mink (Mustela lutreola) is critically endangered throughout its range. However, important aspects such as census size and inbreeding remain still unknown in many populations, so it is crucial to develop new methods to monitor this species. In this work, we placed hair tubes along riverbanks in a border area of the Iberian population, which allowed the genetic identification of 76 European mink hair samples. We then applied a reduced representation genomic sequencing (ddRAD) technique to a subset of these samples to test whether we could extract sufficient genomic information from them. We show that several problems with the DNA, including contamination, fragmentation, oxidation, and possibly sample mixing, affected the samples. Using various bioinformatic techniques to reduce these problems, we were able to unambiguously genotype 19 hair samples belonging to six individuals. This small number of individuals showed that the demographic status of the species in this peripheral population is worse than expected. The data obtained also allowed us to perform preliminary analyses of relatedness and inbreeding. Although further improvements in sampling and analysis are needed, the application of the ddRAD technique to non-invasively obtained hairs represents a significant advance in the genomic study of endangered species.

Individual Identification of Large Felids in Field Studies: Common Methods, Challenges, and Implications for Conservation Science

Large felids represent some of the most threatened large mammals on Earth, critical for both tourism economies and ecosystem function. Most populations are in a state of decline, and their monitoring and enumeration is therefore critical for conservation. This typically rests on the accurate identification of individuals within their populations. We review the most common and current survey methods used in individual identification studies of large felid ecology (body mass > 25 kg). Remote camera trap photography is the most extensively used method to identify leopards, snow leopards, jaguars, tigers, and cheetahs which feature conspicuous and easily identifiable coat patterning. Direct photographic surveys and genetic sampling are commonly used for species that do not feature easily identifiable coat patterning such as lions. We also discuss the accompanying challenges encountered in several field studies, best practices that can help increase the precision and accuracy of identification and provide generalised ratings for the common survey methods used for individual identification.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

Variable outcomes of hybridization between declining Alosa alosa and Alosa fallax

Hybridization dynamics between co-occurring species in environments where human-mediated changes take place are important to quantify for furthering our understanding of human impacts on species evolution and for informing management. The allis shad Alosa alosa (Linnaeus, 1758) and twaite shad Alosa fallax (Lacépède, 1803), two clupeids sister species, have been severely impacted by human activities across Europe. The shrinkage of A. alosa distribution range along with the decline of the remaining populations' abundance threatens its persistence. The main objective was to evaluate the extent of hybridization and introgression between those interacting species. We developed a set of 77 species-specific SNP loci that allowed a better resolution than morphological traits as they enabled the detection of hybrids up to the third generation. Variable rates of contemporary hybridization and introgression patterns were detected in 12 studied sites across the French Atlantic coast. Mitochondrial markers revealed a cyto-nuclear discordance almost invariably involving A. alosa individuals with an A. fallax mitochondrial DNA and provided evidence of historical asymmetric introgression. Overall, contemporary and historical introgression revealed by nuclear and mitochondrial markers strongly suggests that a transfer of genes occurs from A. fallax toward A. alosa genome since at least four generations. Moreover, the outcomes of introgression greatly depend on the catchments where local processes are thought to occur. Undoubtedly, interspecific interaction and gene flow should not be overlooked when considering the management of those species.

Canadian polar bear population structure using genome-wide markers

Predicting the consequences of environmental changes, including human-mediated climate change on species, requires that we quantify range-wide patterns of genetic diversity and identify the ecological, environmental, and historical factors that have contributed to it. Here, we generate baseline data on polar bear population structure across most Canadian subpopulations (n = 358) using 13,488 genome-wide single nucleotide polymorphisms (SNPs) identified with double-digest restriction site-associated DNA sequencing (ddRAD). Our ddRAD dataset showed three genetic clusters in the sampled Canadian range, congruent with previous studies based on microsatellites across the same regions; however, due to a lack of sampling in Norwegian Bay, we were unable to confirm the existence of a unique cluster in that subpopulation. These data on the genetic structure of polar bears using SNPs provide a detailed baseline against which future shifts in population structure can be assessed, and opportunities to develop new noninvasive tools for monitoring polar bears across their range.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

De Novo Assembly and Annotation from Parental and F1 Puma Genomes of the Florida Panther Genetic Restoration Program

In the mid-1990s, the population size of Florida panthers became so small that many individuals manifested traits associated with inbreeding depression (e.g., heart defects, cryptorchidism, high pathogen-parasite load). To mitigate these effects, pumas from Texas were introduced into South Florida to augment genetic variation in Florida panthers. In this study, we report a de novo puma genome assembly and annotation after resequencing 10 individual genomes from partial Florida-Texas-F₁ trios. The final genome assembly consisted of ∼2.6 Gb and 20,561 functionally annotated protein-coding genes. Foremost, expanded gene families were associated with neuronal and embryological development, whereas contracted gene families were associated with olfactory receptors. Despite the latter, we characterized 17 positively selected genes related to the refinement of multiple sensory perceptions, most notably to visual capabilities. Furthermore, genes under positive selection were enriched for the targeting of proteins to the endoplasmic reticulum, degradation of mRNAs, and transcription of viral genomes. Nearly half (48.5%) of ∼6.2 million SNPs analyzed in the total sample set contained putative unique Texas alleles. Most of these alleles were likely inherited to subsequent F₁ Florida panthers, as these individuals manifested a threefold increase in observed heterozygosity with respect to their immediate, canonical Florida panther predecessors. Demographic simulations were consistent with a recent colonization event in North America by a small number of founders from South America during the last glacial period. In conclusion, we provide an extensive set of genomic resources for pumas and elucidate the genomic effects of genetic rescue on this iconic conservation success story.

Improving cost-efficiency of faecal genotyping: New tools for elephant species

Despite the critical need for non-invasive tools to improve monitoring of wildlife populations, especially for endangered and elusive species, faecal genetic sampling has not been adopted as regular practice, largely because of the associated technical challenges and cost. Substantial work needs to be undertaken to refine sample collection and preparation methods in order to improve sample set quality and provide cost-efficient tools that can effectively support wildlife management. In this study, we collected an extensive set of forest elephant (Loxodonta cyclotis) faecal samples throughout Gabon, Central Africa, and prepared them for genotyping using 107 single-nucleotide polymorphism assays. We developed a new quantitative polymerase chain reaction (PCR) assay targeting a 130-bp nuclear DNA fragment and demonstrated its suitability for degraded samples in all three elephant species. Using this assay to compare the efficacy of two sampling methods for faecal DNA recovery, we found that sampling the whole surface of a dung pile with a swab stored in a small tube of lysis buffer was a convenient method producing high extraction success and DNA yield. We modelled the influence of faecal quality and storage time on DNA concentration in order to provide recommendations for optimized collection and storage. The maximum storage time to ensure 75% success was two months for samples collected within 24 hours after defecation and extended to four months for samples collected within one hour. Lastly, the real-time quantitative PCR assay allowed us to predict genotyping success and pre-screen DNA samples, thus further increasing the cost-efficiency of our approach. We recommend combining the validation of an efficient sampling method, the build of in-country DNA extraction capacity for reduced storage time and the development of species-specific quantitative PCR assays in order to increase the cost-efficiency of routine non-invasive DNA analyses and expand the use of next-generation markers to non-invasive samples.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.