Genome sequencing and conservation genomics in the Scandinavian wolverine population

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Prioritizing Endangered Species in Genome Sequencing: Conservation Genomics in Action with the First Platinum-Standard Reference-Quality Genome of the Critically Endangered European Mink Mustela lutreola L., 1761

The European mink Mustela lutreola (Mustelidae) ranks among the most endangered mammalian species globally, experiencing a rapid and severe decline in population size, density, and distribution. Given the critical need for effective conservation strategies, understanding its genomic characteristics becomes paramount. To address this challenge, the platinum-quality, chromosome-level reference genome assembly for the European mink was successfully generated under the project of the European Mink Centre consortium. Leveraging PacBio HiFi long reads, we obtained a 2586.3 Mbp genome comprising 25 scaffolds, with an N50 length of 154.1 Mbp. Through Hi-C data, we clustered and ordered the majority of the assembly (>99.9%) into 20 chromosomal pseudomolecules, including heterosomes, ranging from 6.8 to 290.1 Mbp. The newly sequenced genome displays a GC base content of 41.9%. Additionally, we successfully assembled the complete mitochondrial genome, spanning 16.6 kbp in length. The assembly achieved a BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness score of 98.2%. This high-quality reference genome serves as a valuable genomic resource for future population genomics studies concerning the European mink and related taxa. Furthermore, the newly assembled genome holds significant potential in addressing key conservation challenges faced by M. lutreola. Its applications encompass potential revision of management units, assessment of captive breeding impacts, resolution of phylogeographic questions, and facilitation of monitoring and evaluating the efficiency and effectiveness of dedicated conservation strategies for the European mink. This species serves as an example that highlights the paramount importance of prioritizing endangered species in genome sequencing projects due to the race against time, which necessitates the comprehensive exploration and characterization of their genomic resources before their populations face extinction.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes)

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

A chromosome-level genome assembly of the yellow-throated marten (Martes flavigula)

The yellow-throated marten (Martes flavigula) is a medium-sized carnivore that is widely distributed across much of Asia and occupies an extensive variety of habitats. We reported a high-quality genome assembly of this organism that was generated using Oxford Nanopore and Hi-C technologies. The final genome sequences contained 215 contigs with a total size of 2,449.15 Mb and a contig N50 length of 68.60 Mb. Using Hi-C analysis, 2,419.20 Mb (98.78%) of the assembled sequences were anchored onto 21 linkage groups. Merqury evaluation suggested that the genome was 94.95% complete with a QV value of 43.75. Additionally, the genome was found to comprise approximately 39.74% repeat sequences, of which long interspersed elements (LINE) that accounted for 26.13% of the entire genome, were the most abundant. Of the 20,464 protein-coding genes, prediction and functional annotation was successfully performed for 20,322 (99.31%) genes. The high-quality, chromosome-level genome of the marten reported in this study will serve as a reference for future studies on genetic diversity, evolution, and conservation biology.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

A chromosome-level genome assembly reveals genomic characteristics of the American mink (Neogale vison)

Availability of a contiguous chromosome-level genome assembly is the foundational step to develop genome-based studies in American mink (Neogale vison). The main objective of this study was to provide a high quality chromosome-level genome assembly for American mink. An initial draft of the genome assembly was generated using 2,884,047 PacBio long reads. Integration of Hi-C data into the initial draft led to an assembly with 183 scaffolds and scaffold N50 of 220 Mb. This gap-free genome assembly of American mink (ASM_NN_V1) had a length of 2.68 Gb in which about 98.6% of the whole genome was covered by 15 chromosomes. In total, 25,377 genes were predicted across the American mink genome using the NCBI Eukaryotic Genome Annotation Pipeline. In addition, gene orthology, demographic history, synteny blocks, and phylogenetic relationships were studied in connection with the genomes of other related Carnivora. Furthermore, population-based statistics of 100 sequenced mink were presented using the newly assembled genome. Remarkable improvements were observed in genome contiguity, the number of scaffolds, and annotation compared to the first draft of mink genome assembly (NNQGG.v01). This high-quality genome assembly will support the development of efficient breeding strategies as well as conservation programs for American mink.

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.

Chromosomal-level reference genome assembly of the North American wolverine (Gulo gulo luscus): a resource for conservation genomics

We report a chromosomal-level genome assembly of a male North American wolverine (Gulo gulo luscus) from the Kugluktuk region of Nunavut, Canada. The genome was assembled directly from long-reads, comprising: 758 contigs with a contig N50 of 36.6 Mb; contig L50 of 20; base count of 2.39 Gb; and a near complete representation (99.98%) of the BUSCO 5.2.2 set of 9,226 genes. A presumptive chromosomal-level assembly was generated by scaffolding against two chromosomal-level Mustelidae reference genomes, the ermine and the Eurasian river otter, to derive a final scaffold N50 of 144.0 Mb and a scaffold L50 of 7. We annotated a comprehensive set of genes that have been associated with models of aggressive behavior, a trait which the wolverine is purported to have in the popular literature. To support an integrated, genomics-based wildlife management strategy at a time of environmental disruption from climate change, we annotated the principal genes of the innate immune system to provide a resource to study the wolverine’s susceptibility to new infectious and parasitic diseases. As a resource, we annotated genes involved in the modality of infection by the coronaviruses, an important class of viral pathogens of growing concern as shown by the recent spillover infections by severe acute respiratory syndrome coronavirus-2 to naïve wildlife. Tabulation of heterozygous single nucleotide variants in our specimen revealed a heterozygosity level of 0.065%, indicating a relatively diverse genetic pool that would serve as a baseline for the genomics-based conservation of the wolverine, a rare cold-adapted carnivore now under threat.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes, as well as searching for species-specific structural variants (SVs). Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observe species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This strongly argues for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

A reduced SNP panel to trace gene flow across southern European wolf populations and detect hybridization with other Canis taxa

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ Q_i ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Q_i ≤ 0.827 to parental groups). An assay design and test with invasive and non-invasive canid samples performed successfully for 178 SNPs. By separating natural population admixture from inter-specific hybridization, our reduced panel can help advance evolutionary research, monitoring, and timely conservation management.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Genome-Wide Evidence for Complex Hybridization and Demographic History in a Group of Cycas From China

Cycads represent one of the most ancestral living seed plants as well as one of the most threatened plant groups in the world. South China is a major center and potential origin of Cycas, the most rapidly diversified lineage of cycads. However, genomic-wide diversity of Cycas remains poorly understood due to the challenge of generating genomic markers associated with their inherent large genomes. Here, we perform a comprehensive conservation genomic study based on restriction-site associated DNA sequencing (RADseq) data in six representative species of Cycas in South China. Consistently low genetic diversity and strong genetic differentiation were detected across species. Both phylogenetic inference and genetic structure analysis via several methods revealed generally congruent groups among the six Cycas species. The analysis with ADMIXTURE showed low mixing of genetic composition among species, while individuals of C. dolichophylla exhibited substantial genetic admixture with C. bifida, C. changjiangensis, and C. balansae. Furthermore, the results from Treemix, f₄-statistic, and ABBA-BABA test were generally consistent and revealed the complex patterns of interspecific gene flow. Relatively strong signals of hybridization were detected between C. dolichophylla and C. szechuanensis, and the ancestor of C. taiwaniana and C. changjiangensis. Distinct patterns of demographic history were inferred for these species by Stairway Plot, and our results suggested that both climate fluctuation and frequent geological activities during the late Pleistocene exerted deep impacts on the population dynamics of these species in South China. Finally, we explore the practical implications of our findings for the development of conservation strategies in Cycas. The present study demonstrates the efficiency of RADseq for conservation genomic studies on non-model species with large and complex genomes. Given the great significance of cycads as a radical transition in the evolution of plant biodiversity, our study provides important insights into the mechanisms of diversification in such recently radiated living fossil taxa.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum

Pleistocene glacial cycles influenced the diversification of high‐latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger‐scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20–30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre‐LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high‐latitude species.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.

Limits and convergence properties of the sequentially Markovian coalescent

Several methods based on the sequentially Markovian coalescent (SMC) make use of full genome sequence data from samples to infer population demographic history including past changes in population size, admixture, migration events and population structure. More recently, the original theoretical framework has been extended to allow the simultaneous estimation of population size changes along with other life history traits such as selfing or seed banking. The latter developments enhance the applicability of SMC methods to nonmodel species. Although convergence proofs have been given using simulated data in a few specific cases, an in-depth investigation of the limitations of SMC methods is lacking. In order to explore such limits, we first develop a tool inferring the best case convergence of SMC methods assuming the true underlying coalescent genealogies are known. This tool can be used to quantify the amount and type of information that can be confidently retrieved from given data sets prior to the analysis of the real data. Second, we assess the inference accuracy when the assumptions of SMC approaches are violated due to departures from the model, namely the presence of transposable elements, variable recombination and mutation rates along the sequence, and SNP calling errors. Third, we deliver a new interpretation of SMC methods by highlighting the importance of the transition matrix, which we argue can be used as a set of summary statistics in other statistical inference methods, uncoupling the SMC from hidden Markov models (HMMs). We finally offer recommendations to better apply SMC methods and build adequate data sets under budget constraints.

Dynamic landscapes in northwestern North America structured populations of wolverines (Gulo gulo)

Cyclic climatic and glacial fluctuations of the Late Quaternary produced a dynamic biogeographic history for high latitudes. To refine our understanding of this history in northwestern North America, we explored geographic structure in a wide-ranging carnivore, the wolverine (Gulo gulo). We examined genetic variation in populations across mainland Alaska, coastal Southeast Alaska, and mainland western Canada using nuclear microsatellite genotypes and sequence data from the mitochondrial DNA (mtDNA) control region and Cytochrome b (Cytb) gene. Data from maternally inherited mtDNA reflect stable populations in Northwest Alaska, suggesting the region harbored wolverine populations since at least the Last Glacial Maximum (LGM; 21 Kya), consistent with their persistence in the fossil record of Beringia. Populations in Southeast Alaska are characterized by minimal divergence, with no genetic signature of long-term refugial persistence (consistent with the lack of pre-Holocene fossil records there). The Kenai Peninsula population exhibits mixed signatures depending on marker type: mtDNA data indicate stability (i.e., historical persistence) and include a private haplotype, whereas biparentally inherited microsatellites exhibit relatively low variation and a lack of private alleles consistent with a more recent Holocene colonization of the peninsula. Our genetic work is largely consistent with the early 20th century taxonomic hypothesis that wolverines on the Kenai Peninsula belong to a distinct subspecies. Our finding of significant genetic differentiation of wolverines inhabiting the Kenai Peninsula, coupled with the peninsula’s burgeoning human population and the wolverine’s known sensitivity to anthropogenic impacts, provides valuable foundational data that can be used to inform conservation and management prescriptions for wolverines inhabiting these landscapes.

Sample identification and pedigree reconstruction in Wolverine (Gulo gulo) using SNP genotyping of non-invasive samples

For conservation genetic studies using non-invasively collected samples, genome-wide data may be hard to acquire. Until now, such studies have instead mostly relied on analyses of traditional genetic markers such as microsatellites (SSRs). Recently, high throughput genotyping of single nucleotide polymorphisms (SNPs) has become available, expanding the use of genomic methods to include non-model species of conservation concern. We have developed a 96-marker SNP array for use in applied conservation monitoring of the Scandinavian wolverine (Gulo gulo) population. By genotyping more than a thousand non-invasively collected samples, we were able to obtain precise estimates of different types of genotyping errors and sample dropout rates. The SNP panel significantly outperforms the SSR markers (and DBY intron markers for sexing) both in terms of precision in genotyping, sex assignment and individual identification, as well as in the proportion of samples successfully genotyped. Furthermore, SNP genotyping offers a simplified laboratory and analysis pipeline with fewer samples needed to be repeatedly genotyped in order to obtain reliable consensus data. In addition, we utilised a unique opportunity to successfully demonstrate the application of SNP genotype data for reconstructing pedigrees in wild populations, by validating the method with samples from wild individuals with known relatedness. By offering a simplified workflow with improved performance, we anticipate this methodology will facilitate the use of non-invasive samples to improve genetic management of many different types of populations that have previously been challenging to survey.

High diversity, inbreeding and a dynamic Pleistocene demographic history revealed by African buffalo genomes

Genomes retain records of demographic changes and evolutionary forces that shape species and populations. Remnant populations of African buffalo (Syncerus caffer) in South Africa, with varied histories, provide an opportunity to investigate signatures left in their genomes by past events, both recent and ancient. Here, we produce 40 low coverage (7.14×) genome sequences of Cape buffalo (S. c. caffer) from four protected areas in South Africa. Genome-wide heterozygosity was the highest for any mammal for which these data are available, while differences in individual inbreeding coefficients reflected the severity of historical bottlenecks and current census sizes in each population. PSMC analysis revealed multiple changes in N_e between approximately one million and 20 thousand years ago, corresponding to paleoclimatic changes and Cape buffalo colonisation of southern Africa. The results of this study have implications for buffalo management and conservation, particularly in the context of the predicted increase in aridity and temperature in southern Africa over the next century as a result of climate change.

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDI) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

Evolvability and evolutionary rescue

The survival prospects of threatened species or populations can sometimes be improved by adaptive change. Such evolutionary rescue is particularly relevant when the threat comes from changing environments, or when long-term population persistence requires range expansion into new habitats. Conservation biologists are therefore often interested in whether or not populations or lineages show a disposition for adaptive evolution, that is, if they are evolvable. Here, we discuss four alternative perspectives that target different causes of evolvability and outline some of the key challenges those perspectives are designed to address. Standing genetic variation provides one familiar estimate of evolvability. Yet, the mere presence of genetic variation is often insufficient to predict if a population will adapt, or how it will adapt. The reason is that adaptive change not only depends on genetic variation, but also on the extent to which this genetic variation can be realized as adaptive phenotypic variation. This requires attention to developmental systems and how plasticity influences evolutionary potential. Finally, we discuss how a better understanding of the different factors that contribute to evolvability can be exploited in conservation practice.

On-site genetic analysis for species identification using lab-on-a-chip

This paper presents a microfluidic device capable of performing genetic analysis on dung samples to identify White Rhinoceros (Ceratotherium simum). The development of a microfluidic device, which can be used in the field, offers a portable and cost-effective solution for DNA analysis and species identification to aid conservation efforts. Optimization of the DNA extraction processes produced equivalent yields compared to conventional kit-based methods within just 5 minutes. The use of a color-changing loop-mediated isothermal amplification reaction for simultaneous detection of the cytochrome B sequence of C. simum enabled positive results to be obtained within as little as 30 minutes. Field testing was performed at Knowsley Safari to demonstrate real-world applicability of the microfluidic device for testing of biological samples.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Chromosomal-level genome assembly of the scimitar-horned oryx: Insights into diversity and demography of a species extinct in the wild

Captive populations provide a valuable insurance against extinctions in the wild. However, they are also vulnerable to the negative impacts of inbreeding, selection and drift. Genetic information is therefore considered a critical aspect of conservation management. Recent developments in sequencing technologies have the potential to improve the outcomes of management programmes; however, the transfer of these approaches to applied conservation has been slow. The scimitar-horned oryx (Oryx dammah) is a North African antelope that has been extinct in the wild since the early 1980s and is the focus of a large-scale and long-term reintroduction project. To enable the selection of suitable founder individuals, facilitate post-release monitoring and improve captive breeding management, comprehensive genomic resources are required. Here, we used 10X Chromium sequencing together with Hi-C contact mapping to develop a chromosomal-level genome assembly for the species. The resulting assembly contained 29 chromosomes with a scaffold N50 of 100.4 Mb, and displayed strong chromosomal synteny with the cattle genome. Using resequencing data from six additional individuals, we demonstrated relatively high genetic diversity in the scimitar-horned oryx compared to other mammals, despite it having experienced a strong founding event in captivity. Additionally, the level of diversity across populations varied according to management strategy. Finally, we uncovered a dynamic demographic history that coincided with periods of climate variation during the Pleistocene. Overall, our study provides a clear example of how genomic data can uncover valuable insights into captive populations and contributes important resources to guide future management decisions of an endangered species.

Sequencing smart: De novo sequencing and assembly approaches for a non-model mammal

Background

Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses.

Results

Generally, assemblies containing more data types achieved better scores in our ranking system. However, when accounting for misassemblies, this was not always the case for Bionano and low-coverage 10x Genomics (for scaffolding only). We also find that the extra cost associated with combining multiple data types is not necessarily associated with better genome assemblies.

Conclusions

The high degree of variability between each de novo assembly method (assessed from the 7 key metrics) highlights the importance of carefully devising the sequencing strategy to be able to carry out the desired analysis. Adding more data to genome assemblies does not always result in better assemblies, so it is important to understand the nuances of genomic data integration explained here, in order to obtain cost-effective value for money when sequencing genomes.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

A practical introduction to sequentially Markovian coalescent methods for estimating demographic history from genomic data

A common goal of population genomics and molecular ecology is to reconstruct the demographic history of a species of interest. A pair of powerful tools based on the sequentially Markovian coalescent have been developed to infer past population sizes using genome sequences. These methods are most useful when sequences are available for only a limited number of genomes and when the aim is to study ancient demographic events. The results of these analyses can be difficult to interpret accurately, because doing so requires some understanding of their theoretical basis and of their sensitivity to confounding factors. In this practical review, we explain some of the key concepts underpinning the pairwise and multiple sequentially Markovian coalescent methods (PSMC and MSMC, respectively). We relate these concepts to the use and interpretation of these methods, and we explain how the choice of different parameter values by the user can affect the accuracy and precision of the inferences. Based on our survey of 100 PSMC studies and 30 MSMC studies, we describe how the two methods are used in practice. Readers of this article will become familiar with the principles, practice, and interpretation of the sequentially Markovian coalescent for inferring demographic history.

From reference genomes to population genomics: comparing three reference-aligned reduced-representation sequencing pipelines in two wildlife species

Background

Recent advances in genomics have greatly increased research opportunities for non-model species. For wildlife, a growing availability of reference genomes means that population genetics is no longer restricted to a small set of anonymous loci. When used in conjunction with a reference genome, reduced-representation sequencing (RRS) provides a cost-effective method for obtaining reliable diversity information for population genetics. Many software tools have been developed to process RRS data, though few studies of non-model species incorporate genome alignment in calling loci. A commonly-used RRS analysis pipeline, Stacks, has this capacity and so it is timely to compare its utility with existing software originally designed for alignment and analysis of whole genome sequencing data. Here we examine population genetic inferences from two species for which reference-aligned reduced-representation data have been collected. Our two study species are a threatened Australian marsupial (Tasmanian devil Sarcophilus harrisii; declining population) and an Arctic-circle migrant bird (pink-footed goose Anser brachyrhynchus; expanding population). Analyses of these data are compared using Stacks versus two widely-used genomics packages, SAMtools and GATK. We also introduce a custom R script to improve the reliability of single nucleotide polymorphism (SNP) calls in all pipelines and conduct population genetic inferences for non-model species with reference genomes.

Results

Although we identified orders of magnitude fewer SNPs in our devil dataset than for goose, we found remarkable symmetry between the two species in our assessment of software performance. For both datasets, all three methods were able to delineate population structure, even with varying numbers of loci. For both species, population structure inferences were influenced by the percent of missing data.

Conclusions

For studies of non-model species with a reference genome, we recommend combining Stacks output with further filtering (as included in our R pipeline) for population genetic studies, paying particular attention to potential impact of missing data thresholds. We recognise SAMtools as a viable alternative for researchers more familiar with this software. We caution against the use of GATK in studies with limited computational resources or time.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5806-y) contains supplementary material, which is available to authorized users.

Genetic tagging in the Anthropocene: scaling ecology from alleles to ecosystems

The Anthropocene is an era of marked human impact on the world. Quantifying these impacts has become central to understanding the dynamics of coupled human-natural systems, resource-dependent livelihoods, and biodiversity conservation. Ecologists are facing growing pressure to quantify the size, distribution, and trajectory of wild populations in a cost-effective and socially acceptable manner. Genetic tagging, combined with modern computational and genetic analyses, is an under-utilized tool to meet this demand, especially for wide-ranging, elusive, sensitive, and low-density species. Genetic tagging studies are now revealing unprecedented insight into the mechanisms that control the density, trajectory, connectivity, and patterns of human-wildlife interaction for populations over vast spatial extents. Here, we outline the application of, and ecological inferences from, new analytical techniques applied to genetically tagged individuals, contrast this approach with conventional methods, and describe how genetic tagging can be better applied to address outstanding questions in ecology. We provide example analyses using a long-term genetic tagging dataset of grizzly bears in the Canadian Rockies. The genetic tagging toolbox is a powerful and overlooked ensemble that ecologists and conservation biologists can leverage to generate evidence and meet the challenges of the Anthropocene.