Novel efficient genome-wide SNP panels for the conservation of the highly endangered Iberian lynx

Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs' genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

First whole-genome sequence and assembly of the Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps), a critically endangered species, using Oxford Nanopore Technologies

The Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps) is currently considered one of the most endangered primates in the world and is classified as critically endangered [International union for conservation of nature (IUCN)]. It faces multiple threats, the most significant one being habitat loss due to deforestation in western Ecuador. Genomic tools are keys for the management of endangered species, but this requires a reference genome, which until now was unavailable for A. f. fusciceps. The present study reports the first whole-genome sequence and assembly of A. f. fusciceps generated using Oxford Nanopore long reads. DNA was extracted from a subadult male, and libraries were prepared for sequencing following the Ligation Sequencing Kit SQK-LSK112 workflow. Sequencing was performed using a MinION Mk1C sequencer. The sequencing reads were processed to generate a genome assembly. Two different assemblers were used to obtain draft genomes using raw reads, of which the Flye assembly was found to be superior. The final assembly has a total length of 2.63 Gb and contains 3,861 contigs, with an N50 of 7,560,531 bp. The assembly was analyzed for annotation completeness based on primate ortholog prediction using a high-resolution database, and was found to be 84.3% complete, with a low number of duplicated genes indicating a precise assembly. The annotation of the assembly predicted 31,417 protein-coding genes, comparable with other mammal assemblies. A reference genome for this critically endangered species will allow researchers to gain insight into the genetics of its populations and thus aid conservation and management efforts of this vulnerable species.

Development and evaluation of a novel single nucleotide polymorphism panel for North American bison

Genome-wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole-genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping-by-sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations.

Genomic insights into the conservation of wild and domestic animal diversity: A review

Due to environmental change and anthropogenic activities, global biodiversity has suffered an unprecedented loss, and the world is now heading toward the sixth mass extinction event. This urges the need to step up our efforts to promote the sustainable use of animal genetic resources and plan effective strategies for their conservation. Although habitat preservation and restoration are the primary means of conserving biodiversity, genomic technologies offer a variety of novel tools for identifying biodiversity hotspots and thus, support conservation efforts. Conservation genomics is a broad area of science that encompasses the application of genomic data from thousands or tens of thousands of genome-wide markers to address important conservation biology concerns. Genomic approaches have revolutionized the way we understand and manage animal populations, providing tools to identify and preserve unique genetic variants and alleles responsible for adaptive genetic variation, reducing the deleterious consequences of inbreeding, and increasing the adaptive potential of threatened species. The advancement of genomic technologies, particularly comparative genomic approaches, and the increased accessibility of genomic resources in the form of genome-enabled taxa for non-model organisms, provides a distinct advantage in defining conservation units over traditional genetics approaches. The objective of this review is to provide an exhaustive overview of the concept of conservation genomics, discuss the rationale behind the transition from conservation genetics to genomic approaches, and emphasize the potential applications of genomic techniques for conservation purposes. We also highlight interesting case studies in both livestock and wildlife species where genomic techniques have been used to accomplish conservation goals. Finally, we address some challenges and future perspectives in this field.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Applying genomic approaches to delineate conservation strategies using the freshwater mussel Margaritifera margaritifera in the Iberian Peninsula as a model

Effective conservation actions to counteract the current decline of populations and species require a deep knowledge on their genetic structure. We used Single Nucleotide Polymorphisms (SNPs) to infer the population structure of the highly threatened freshwater pearl mussel Margaritifera margaritifera in the Iberian Peninsula. A total of 130 individuals were collected from 26 locations belonging to 16 basins. We obtained 31,692 SNPs through Genotyping by Sequencing (GBS) and used this dataset to infer population structure. Genetic diversity given as observed heterozygosity was low. Pairwise F_ST comparisons revealed low levels of genetic differentiation among geographically close populations. Up to 3 major genetic lineages were determined: Atlantic, Cantabrian and Douro. This structure suggests a close co-evolutionary process with brown trout (Salmo trutta), the primordial fish host of this mussel in the studied area. Some sub-basins showed some genetic structuring, whereas in others no intrapopulation differentiation was found. Our results confirm that genetic conservation units do not match individual basins, and that knowledge about the genetic structure is necessary before planning recovery plans that may involve relocation or restocking. The same reasoning should be applied to strictly freshwater species that are sessile or have restricted dispersal abilities and are currently imperiled worldwide.

African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity

African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.

Selection of genome-wide SNPs for pooled allelotyping assays useful for population monitoring

Parasitic worms are serious pests of humans, livestock and crops worldwide. Multiple management strategies are employed in order to reduce their impact, and some of these may affect their genome and population allelic frequency distribution. The evolution of chemical resistance, ecological changes, and pest dispersal have allowed an increasing number of pests to become difficult to control with current management methods. Their lifestyle limits the use of ecological and individual-based management of populations. There is a need to develop rapid, affordable, and simple diagnostics to assess the efficacy of management strategies and delay the evolution of resistance to these strategies. This study presents a multi-locus, equal-representation, whole genome pooled SNPs selection approach as a monitoring tool for the ovine nematode parasite Haemonchus contortus. The SNP selection method used two reference genomes of different quality, then validated these SNPs against a high-quality recent genome assembly. From over 11 million high-quality SNPs identified, 334 SNPs were selected, of which 262 were species-specific, yielded similar allele frequencies when assessed as multiple individuals or as pools of individuals, and suitable to distinguish mixed nematode isolate pools from single isolate pools. As a proof-of-concept, 21 Australian H. contortus populations with various phenotypes and genotypes were screened. This analysis confirmed the overall low-level of genetic differentiation between populations collected from the field, but clearly identifying highly inbred populations, and populations showing genetic signatures associated with chemical resistance. The analysis showed that 66% of the SNPs were necessary for stability in assessing population genetic patterns, and SNP pairs did not show linkage according to allelic frequencies across the 21 populations. This method demonstrates that ongoing monitoring of parasite allelic frequencies and genetic changes can be achieved as a management assessment tool to identify drug-treatment failure, population incursions, and inbreeding signatures due to selection. The SNP selection method could also be applied to other parasite species.

Waste not, want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single-nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black-capped vireo source-sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.

A robust sequencing assay of a thousand amplicons for the high-throughput population monitoring of Alpine ibex immunogenetics

Polymorphism for immune functions can explain significant variation in health and reproductive success within species. Drastic loss in genetic diversity at such loci constitutes an extinction risk and should be monitored in species of conservation concern. However, effective implementations of genome‐wide immune polymorphism sets into high‐throughput genotyping assays are scarce. Here, we report the design and validation of a microfluidics‐based amplicon sequencing assay to comprehensively capture genetic variation in Alpine ibex (Capra ibex). This species represents one of the most successful large mammal restorations recovering from a severely depressed census size and a massive loss in diversity at the major histocompatibility complex (MHC). We analysed 65 whole‐genome sequencing sets of the Alpine ibex and related species to select the most representative markers and to prevent primer binding failures. In total, we designed ~1,000 amplicons densely covering the MHC, further immunity‐related genes as well as randomly selected genome‐wide markers for the assessment of neutral population structure. Our analysis of 158 individuals shows that the genome‐wide markers perform equally well at resolving population structure as RAD‐sequencing or low‐coverage genome sequencing data sets. Immunity‐related loci show unexpectedly high degrees of genetic differentiation within the species. Such information can now be used to define highly targeted individual translocations. Our design strategy can be realistically implemented into genetic surveys of a large range of species. In conclusion, leveraging whole‐genome sequencing data sets to design targeted amplicon assays allows the simultaneous monitoring of multiple genetic risk factors and can be translated into species conservation recommendations.

Breakthroughs and new horizons in reproductive biology of rare and endangered animal species

Because of higher extinction rates due to human and natural factors, more basic and applied research in reproductive biology is required to preserve wild species and design proper strategies leading to sustainable populations. The objective of the review is to highlight recent, inspiring breakthroughs in wildlife reproduction science that will set directions for future research and lead to more successes in conservation biology. Despite new tools and approaches allowing a better and faster understanding of key mechanisms, we still know little about reproduction in endangered species. Recently, the most striking advances have been obtained in nonmammalian species (fish, birds, amphibians, or corals) with the development of alternative solutions to preserve fertility or new information about parental nutritional influence on embryo development. A novel way has also been explored to consider the impact of environmental changes on reproduction-the allostatic load-in a vast array of species (from primates to fish). On the horizon, genomic tools are expected to considerably change the way we study wildlife reproduction and develop a concept of "precision conservation breeding." When basic studies in organismal physiology are conducted in parallel, new approaches using stem cells to create artificial gametes and gonads, innovations in germplasm storage, and more research on reproductive microbiomes will help to make a difference. Lastly, multiple challenges (for instance, poor integration of new tools in conservation programs, limited access to study animals, or few publication options) will have to be addressed if we want reproductive biology to positively impact conservation of biodiversity.

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.

From microsatellites to single nucleotide polymorphisms for the genetic monitoring of a critically endangered sturgeon

The use of genetic information is crucial in conservation programs for the establishment of breeding plans and for the evaluation of restocking success. Short tandem repeats (STRs) have been the most widely used molecular markers in such programs, but next-generation sequencing approaches have prompted the transition to genome-wide markers such as single nucleotide polymorphisms (SNPs). Until now, most sturgeon species have been monitored using STRs. The low diversity found in the critically endangered European sturgeon (Acipenser sturio), however, makes its future genetic monitoring challenging, and the current resolution needs to be increased. Here, we describe the discovery of a highly informative set of 79 SNPs using double-digest restriction-associated DNA (ddRAD) sequencing and its validation by genotyping using the MassARRAY system. Comparing with STRs, the SNP panel proved to be highly efficient and reproducible, allowing for more accurate parentage and kinship assignments' on 192 juveniles of known pedigree and 40 wild-born adults. We explore the effectiveness of both markers to estimated relatedness and inbreeding, using simulated and empirical datasets. Interestingly, we found significant correlations between STRs and SNPs at individual heterozygosity and inbreeding that give support to a reasonable representation of whole genome diversity for both markers. These results are useful for the conservation program of A. sturio in building a comprehensive studbook, which will optimize conservation strategies. This approach also proves suitable for other case studies in which highly discriminatory genetic markers are needed to assess parentage and kinship.

metapop2: Re-implementation of software for the analysis and management of subdivided populations using gene and allelic diversity

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long-term selection response than heterozygosity. Here we present a completely revised and updated re-implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross-software and cross-platform compatibility.

Genetic evaluation of the Iberian lynx ex situ conservation programme

Ex situ programmes have become critical for improving the conservation of many threatened species, as they establish backup populations and provide individuals for reintroduction and reinforcement of wild populations. The Iberian lynx was considered the most threatened felid species in the world in the wake of a dramatic decline during the second half of the 20th century that reduced its numbers to around only 100 individuals. An ex situ conservation programme was established in 2003 with individuals from the two well-differentiated, remnant populations, with great success from a demographic point of view. Here, we evaluate the genetic status of the Iberian lynx captive population based on molecular data from 36 microsatellites, including patterns of relatedness and representativeness of the two remnant genetic backgrounds among founders, the evolution of diversity and inbreeding over the years, and genetic differentiation among breeding facilities. In general terms, the ex situ population harbours most of the genetic variability found in the two wild populations and has been able to maintain reasonably low levels of inbreeding and high diversity, thus validating the applied management measures and potentially representing a model for other species in need of conservation.

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.

The conservation physiology toolbox: status and opportunities

For over a century, physiological tools and techniques have been allowing researchers to characterize how organisms respond to changes in their natural environment and how they interact with human activities or infrastructure. Over time, many of these techniques have become part of the conservation physiology toolbox, which is used to monitor, predict, conserve, and restore plant and animal populations under threat. Here, we provide a summary of the tools that currently comprise the conservation physiology toolbox. By assessing patterns in articles that have been published in 'Conservation Physiology' over the past 5 years that focus on introducing, refining and validating tools, we provide an overview of where researchers are placing emphasis in terms of taxa and physiological sub-disciplines. Although there is certainly diversity across the toolbox, metrics of stress physiology (particularly glucocorticoids) and studies focusing on mammals have garnered the greatest attention, with both comprising the majority of publications (>45%). We also summarize the types of validations that are actively being completed, including those related to logistics (sample collection, storage and processing), interpretation of variation in physiological traits and relevance for conservation science. Finally, we provide recommendations for future tool refinement, with suggestions for: (i) improving our understanding of the applicability of glucocorticoid physiology; (ii) linking multiple physiological and non-physiological tools; (iii) establishing a framework for plant conservation physiology; (iv) assessing links between environmental disturbance, physiology and fitness; (v) appreciating opportunities for validations in under-represented taxa; and (vi) emphasizing tool validation as a core component of research programmes. Overall, we are confident that conservation physiology will continue to increase its applicability to more taxa, develop more non-invasive techniques, delineate where limitations exist, and identify the contexts necessary for interpretation in captivity and the wild.

A SNP resource for studying North American moose

Background: Moose ( Alces alces) colonized the North American continent from Asia less than 15,000 years ago, and spread across the boreal forest regions of Canada and the northern United States (US).  Contemporary populations have low genetic diversity, due either to low number of individuals in the original migration (founder effect), and/or subsequent population bottlenecks in North America.  Genetic tests based on informative single nucleotide polymorphism (SNP) markers are helpful in forensic and wildlife conservation activities, but have been difficult to develop for moose, due to the lack of a reference genome assembly and whole genome sequence (WGS) data. Methods:  WGS data were generated for four individual moose from the US states of Alaska, Idaho, Wyoming, and Vermont with minimum and average genome coverage depths of 14- and 19-fold, respectively.  Cattle and sheep reference genomes were used for aligning sequence reads and identifying moose SNPs. Results:  Approximately 11% and 9% of moose WGS reads aligned to cattle and sheep genomes, respectively.  The reads clustered at genomic segments, where sequence identity between these species was greater than 95%.  In these segments, average mapped read depth was approximately 19-fold.  Sets of 46,005 and 36,934 high-confidence SNPs were identified from cattle and sheep comparisons, respectively, with 773 and 552 of those having minor allele frequency of 0.5 and conserved flanking sequences in all three species.  Among the four moose, heterozygosity and allele sharing of SNP genotypes were consistent with decreasing levels of moose genetic diversity from west to east.  A minimum set of 317 SNPs, informative across all four moose, was selected as a resource for future SNP assay design. Conclusions:  All SNPs and associated information are available, without restriction, to support development of SNP-based tests for animal identification, parentage determination, and estimating relatedness in North American moose.