Population genomics of the white-beaked dolphin (Lagenorhynchus albirostris): Implications for conservation amid climate-driven range shifts

Climate change is rapidly affecting species distributions across the globe, particularly in the North Atlantic. For highly mobile and elusive cetaceans, the genetic data needed to understand population dynamics are often scarce. Cold-water obligate species such as the white-beaked dolphin (Lagenorhynchus albirostris) face pressures from habitat shifts due to rising sea surface temperatures in addition to other direct anthropogenic threats. Unravelling the genetic connectivity between white-beaked dolphins across their range is needed to understand the extent to which climate change and anthropogenic pressures may impact species-wide genetic diversity and identify ways to protect remaining habitat. We address this by performing a population genomic assessment of white-beaked dolphins using samples from much of their contemporary range. We show that the species displays significant population structure across the North Atlantic at multiple scales. Analysis of contemporary migration rates suggests a remarkably high connectivity between populations in the western North Atlantic, Iceland and the Barents Sea, while two regional populations in the North Sea and adjacent UK and Irish waters are highly differentiated from all other clades. Our results have important implications for the conservation of white-beaked dolphins by providing guidance for the delineation of more appropriate management units and highlighting the risk that local extirpation may have on species-wide genetic diversity. In a broader context, this study highlights the importance of understanding genetic structure of all species threatened with climate change-driven range shifts to assess the risk of loss of species-wide genetic diversity.

Lion Localizer: A software tool for inferring the provenance of lions (Panthera leo) using mitochondrial DNA

The illegal poaching of lions for their body parts poses a severe threat to lion populations across Africa. Poaching accounts for 35% of all human-caused lion deaths, with 51% attributed to retaliatory killings following livestock predation. In nearly half of the retaliatory killings, lion body parts are removed, suggesting that high demand for lion body parts may fuel killings attributed to human-lion conflict. Trafficked items are often confiscated in transit or destination countries far from their country of origin. DNA from lion parts may in some cases be the only available means for examining their geographic origins. In this paper, we present the Lion Localizer, a full-stack software tool that houses a comprehensive database of lion mitochondrial DNA (mtDNA) sequences sourced from previously published studies. The database covers 146 localities from across the African continent and India, providing information on the potential provenance of seized lion body parts. Lion mtDNA sequences of 350 or 1,140 bp corresponding to the cytochrome b region can be generated from lion products and queried against the Lion Localizer database. Using the query sequence, the Lion Localizer generates a listing of exact or partial matches, which are displayed on an interactive map of Africa. This allows for the rapid identification of potential regions and localities where lions have been or are presently being targeted by poachers. By examining the potential provenance of lion samples, the Lion Localizer serves as a valuable resource in the fight against lion poaching. The software is available at https://lionlocalizer.org.

Evaluation of the Impact of Population Management on the Genetic Parameters of Selected Spiral-Horned Antelopes

The rapid loss of biodiversity and the associated reduction and fragmentation of habitats means that ex situ populations have become an important part of species conservation. These populations, which are often established from a small number of founders, require careful management to avoid the negative effects of genetic drift and inbreeding. Although the inclusion of molecular data is recommended, their availability for captive breeding management remains limited. The aim of this study was to evaluate the relationship between the levels of genetic diversity in six spiral-horned antelope taxa bred under human care and their respective management strategies, conservation status, demography, and geographic origin, using 10 nuclear DNA microsatellite loci and mitochondrial control region DNA sequences. Our findings include associations between genetic diversity and management intensity but also with the diversity and contribution of wild populations to captive founders, with some populations apparently composed of animals from divergent wild lineages elevating captive genetic diversity. When population sizes are large, the potential advantages of maximizing genetic diversity in widely outcrossed populations may need careful consideration with respect to the potential disruption of adaptive diversity. Genetic data serve as a robust tool for managing captive populations, yet their interpretation necessitates a comprehensive understanding of species biology and history.

Multidisciplinary estimates of connectivity and population structure suggest the use of multiple units for the conservation and management of meagre, Argyrosomus regius

Information on population structure and connectivity of targeted species is key for proper implementation of spatial conservation measures. We used a combination of genomics, biophysical modelling, and biotelemetry to infer the population structure and connectivity of Atlantic meagre, an important fisheries resource throughout its distribution. Genetic samples from previously identified Atlantic spawning locations (Gironde, Tejo, Guadalquivir, Banc d'Arguin) and two additional regions (Algarve and Senegal) were analysed using genome-wide SNP-genotyping and mitochondrial DNA analyses. Biophysical models were conducted to investigate larval dispersal and connectivity from the known Atlantic spawning locations. Additionally, thirteen fish were double-tagged with biotelemetry transmitters off the Algarve (Portugal) to assess movement patterns and connectivity of adult individuals. This multidisciplinary approach provided a robust overview of meagre population structure and connectivity in the Atlantic. Nuclear SNP-genotyping showed a clear differentiation between the European and African populations, with significant isolation of the few known Atlantic spawning sites. The limited level of connectivity between these subpopulations is potentially driven by adults, capable of wide-ranging movements and connecting sites 500 km apart, as evidenced by tagging studies, whilst larval dispersal inferred by modelling is much more limited (average of 52 km; 95% of connectivity events up to 174 km). Our results show sufficient evidence of population structure, particularly between Africa and Europe but also within Europe, for the meagre to be managed as separate stocks. Additionally, considering the low degree of larvae connectivity, the implementation of marine protected areas in key spawning sites could be crucial towards species sustainability.

Comparative population genomics of manta rays has global implications for management

Understanding population connectivity and genetic diversity is of fundamental importance to conservation. However, in globally threatened marine megafauna, challenges remain due to their elusive nature and wide-ranging distributions. As overexploitation continues to threaten biodiversity across the globe, such knowledge gaps compromise both the suitability and effectiveness of management actions. Here, we use a comparative framework to investigate genetic differentiation and diversity of manta rays, one of the most iconic yet vulnerable groups of elasmobranchs on the planet. Despite their recent divergence, we show how oceanic manta rays (Mobula birostris) display significantly higher heterozygosity than reef manta rays (Mobula alfredi) and that M. birostris populations display higher connectivity worldwide. Through inferring modes of colonization, we reveal how both contemporary and historical forces have likely influenced these patterns, with important implications for population management. Our findings highlight the potential for fisheries to disrupt population dynamics at both local and global scales and therefore have direct relevance for international conservation of marine species.

Phylogeography of the Sunda pangolin, Manis javanica: Implications for taxonomy, conservation management and wildlife forensics

The Sunda pangolin (Manis javanica) is the most widely distributed Asian pangolin species, occurring across much of Southeast Asia and in southern China. It is classified as Critically Endangered and is one of the most trafficked mammals in the world, which not only negatively impacts wild Sunda pangolin populations but also poses a potential disease risk to other species, including humans and livestock. Here, we aimed to investigate the species' phylogeography across its distribution to improve our understanding of the species' evolutionary history, elucidate any taxonomic uncertainties and enhance the species' conservation genetic management and potential wildlife forensics applications. We sequenced mtDNA genomes from 23 wild Sunda pangolins of known provenance originating from Malaysia to fill sampling gaps in previous studies, particularly in Borneo. To conduct phylogenetic and population genetic analyses of Sunda pangolins across their range, we integrated these newly generated mitochondrial genomes with previously generated mtDNA and nuclear DNA data sets (RAD-seq SNP data). We identified an evolutionarily distinct mtDNA lineage in north Borneo, estimated to be ~1.6 million years divergent from lineages in west/south Borneo and the mainland, comparable to the divergence time from the Palawan pangolin. There appeared to be mitonuclear discordance, with no apparent genetic structure across Borneo based on analysis of nuclear SNPs. These findings are consistent with the 'out of Borneo hypothesis', whereby Sunda pangolins diversified in Borneo before subsequently migrating throughout Sundaland, and/or a secondary contact scenario between mainland and Borneo. We have elucidated possible taxonomic issues in the Sunda/Palawan pangolin complex and highlight the critical need for additional georeferenced samples to accurately apportion its range-wide genetic variation into appropriate taxonomic and conservation units. Additionally, these data have improved forensic identification testing involving these species and permit the implementation of geographic provenance testing in some scenarios.

Conservation management strategy impacts inbreeding and mutation load in scimitar-horned oryx

In an age of habitat loss and overexploitation, small populations, both captive and wild, are increasingly facing the effects of isolation and inbreeding. Genetic management has therefore become a vital tool for ensuring population viability. However, little is known about how the type and intensity of intervention shape the genomic landscape of inbreeding and mutation load. We address this using whole-genome sequence data of the scimitar-horned oryx (Oryx dammah), an iconic antelope that has been subject to contrasting management strategies since it was declared extinct in the wild. We show that unmanaged populations are enriched for long runs of homozygosity (ROH) and have significantly higher inbreeding coefficients than managed populations. Additionally, despite the total number of deleterious alleles being similar across management strategies, the burden of homozygous deleterious genotypes was consistently higher in unmanaged groups. These findings emphasize the risks associated with deleterious mutations through multiple generations of inbreeding. As wildlife management strategies continue to diversify, our study reinforces the importance of maintaining genome-wide variation in vulnerable populations and has direct implications for one of the largest reintroduction attempts in the world.

Genetic diversity goals and targets have improved, but remain insufficient for clear implementation of the post-2020 global biodiversity framework

Genetic diversity among and within populations of all species is necessary for people and nature to survive and thrive in a changing world. Over the past three years, commitments for conserving genetic diversity have become more ambitious and specific under the Convention on Biological Diversity's (CBD) draft post-2020 global biodiversity framework (GBF). This Perspective article comments on how goals and targets of the GBF have evolved, the improvements that are still needed, lessons learned from this process, and connections between goals and targets and the actions and reporting that will be needed to maintain, protect, manage and monitor genetic diversity. It is possible and necessary that the GBF strives to maintain genetic diversity within and among populations of all species, to restore genetic connectivity, and to develop national genetic conservation strategies, and to report on these using proposed, feasible indicators.

Phylogeography and population genetic structure of the European roe deer in Switzerland following recent recolonization

In the early 1800s, the European roe deer (Capreolus capreolus) was probably extirpated from Switzerland, due to overhunting and deforestation. After a federal law was enacted in 1875 to protect lactating females and young, and limiting the hunting season, the roe deer successfully recovered and recolonized Switzerland. In this study, we use mitochondrial DNA and nuclear DNA markers to investigate the recolonization and assess contemporary genetic structure in relation to broad topographic features, in order to understand underlying ecological processes, inform future roe deer management strategies, and explore the opportunity for development of forensic traceability tools. The results concerning the recolonization origin support natural, multidirectional immigration from neighboring countries. We further demonstrate that there is evidence of weak genetic differentiation within Switzerland among topographic regions. Finally, we conclude that the genetic data support the recognition of a single roe deer management unit within Switzerland, within which there is a potential for broad‐scale geographic origin assignment using nuclear markers to support law enforcement.

Effect of diet energy level and genomic residual feed intake on bred Holstein dairy heifer growth and feed efficiency

The objective of this study was to determine growth, feed intake, and feed efficiency of postbred dairy heifers with different genomic residual feed intake (RFI) predicted as a lactating cow when offered diets differing in energy density. Postbred Holstein heifers (n = 128, ages 14-20 mo) were blocked by initial weight (high, medium-high, medium-low, and low) with 32 heifers per block. Each weight block was sorted by RFI (high or low) to obtain 2 pens of heifers with high and low genomically predicted RFI within each block (8 heifers per pen). Low RFI heifers were expected to have greater feed efficiency than high RFI heifers. Dietary treatments consisted of a higher energy control diet based on corn silage and alfalfa haylage [HE; 62.7% total digestible nutrients, 11.8% crude protein, and 45.6% neutral detergent fiber; dry matter (DM) basis], and a lower energy diet diluted with straw (LE; 57.0% total digestible nutrients, 11.7% crude protein, and 50.1% neutral detergent fiber; DM basis). Each pen within a block was randomly allocated a diet treatment to obtain a 2 × 2 factorial arrangement (2 RFI levels and 2 dietary energy levels). Diets were offered in a 120-d trial. Dry matter intake by heifers was affected by diet (11.0 vs. 10.0 kg/d for HE and LE, respectively) but not by RFI or the interaction of RFI and diet. Daily gain was affected by the interaction of RFI and diet, with low RFI heifers gaining more than high RFI heifers when fed LE (0.94 vs. 0.85 kg/d for low and high RFI, respectively), but no difference for RFI groups when fed HE (1.16 vs. 1.19 kg/d for low and high RFI, respectively). Respective feed efficiencies were improved for low RFI compared with high RFI heifers when fed LE (10.6 vs. 11.8 kg of feed DM/kg of gain), but no effect of RFI was found when fed HE (9.4 vs. 9.5 kg of DM/kg of gain for high and low RFI, respectively). No effect of RFI or diet on first-lactation performance through 150 DIM was observed. Based on these results, the feed efficiency of heifers having different genomic RFI may be dependent on diet energy level, whereby low RFI heifers utilized the LE diet more efficiently. The higher fiber straw (LE) diet controlled intake and maintained more desirable heifer weight gains. This suggests that selection for improved RFI in lactating cows may improve feed efficiency in growing heifers when fed to meet growth goals of 0.9 to 1.0 kg of gain/d.

The genome sequence of the European golden eagle, Aquila chrysaetos chrysaetos Linnaeus 1758

We present a genome assembly from an individual female Aquila chrysaetos chrysaetos (the European golden eagle; Chordata; Aves; Accipitridae). The genome sequence is 1.23 gigabases in span. The majority of the assembly is scaffolded into 28 chromosomal pseudomolecules, including the W and Z sex chromosomes.

Global Commitments to Conserving and Monitoring Genetic Diversity Are Now Necessary and Feasible

Global conservation policy and action have largely neglected protecting and monitoring genetic diversity—one of the three main pillars of biodiversity. Genetic diversity (diversity within species) underlies species’ adaptation and survival, ecosystem resilience, and societal innovation. The low priority given to genetic diversity has largely been due to knowledge gaps in key areas, including the importance of genetic diversity and the trends in genetic diversity change; the perceived high expense and low availability and the scattered nature of genetic data; and complicated concepts and information that are inaccessible to policymakers. However, numerous recent advances in knowledge, technology, databases, practice, and capacity have now set the stage for better integration of genetic diversity in policy instruments and conservation efforts. We review these developments and explore how they can support improved consideration of genetic diversity in global conservation policy commitments and enable countries to monitor, report on, and take action to maintain or restore genetic diversity.

Correction to: Authors’ Reply to Letter to the Editor: Continued improvement to genetic diversity indicator for CBD

A correction to this paper has been published: https://doi.org/10.1007/s10592-021-01376-9

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

Phylogeography of the iconic Australian pink cockatoo, Lophochroa leadbeateri

The pink cockatoo (Lophochroa leadbeateri; or Major Mitchell’s cockatoo) is one of Australia’s most iconic bird species. Two subspecies based on morphology are separated by a biogeographical divide, the Eyrean Barrier. Testing the genetic basis for this subspecies delineation, clarifying barriers to gene flow and identifying any cryptic genetic diversity will likely have important implications for conservation and management. Here, we used genome-wide single nucleotide polymorphisms (SNPs) and mitochondrial DNA data to conduct the first range-wide genetic assessment of the species. The aims were to investigate the phylogeography of the pink cockatoo, to characterize conservation units and to reassess subspecies boundaries. We found consistent but weak genetic structure between the two subspecies based on nuclear SNPs. However, phylogenetic analysis of nuclear SNPs and mitochondrial DNA sequence data did not recover reciprocally monophyletic groups, indicating incomplete evolutionary separation between the subspecies. Consequently, we have proposed that the two currently recognized subspecies be treated as separate management units rather than evolutionarily significant units. Given that poaching is suspected to be a threat to this species, we assessed the utility of our data for wildlife forensic applications. We demonstrated that a subspecies identification test could be designed using as few as 20 SNPs.

Phylogenomics and species delimitation for effective conservation of manta and devil rays

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Population history of the golden eagle inferred from whole-genome sequencing of three of its subspecies

The application of evolutionary genetic research to investigate the potential for endangered species to adapt to changing environments is important for conservation biology. Effective population size (Ne) is informative for understanding adaptive potential as it refers to the genetic variation in breeding individuals who have contributed to contemporary and historic population diversity. We reconstruct fluctuations in Ne in three golden eagle subspecies (Japanese, Scottish, North American) using the pairwise sequential Markovian coalescent (PSMC) model based on whole-genome sequence data. Our results indicate the timing of subspeciation events and suggest significant ongoing demographic reductions since the start of the Last Glacial Period. Importantly, we find evidence for gene flow from continental populations into the ancestral Japanese population resulting in a short, sharp recovery in genetic diversity. Timing agrees with the palaeogeographic estimates of land bridge connections between the Japanese archipelago and Asian continent and matches a similar Ne spike in the Scottish population, but not in the North American population. Given contemporary declines in isolated Japanese and UK island populations, our study highlights a concerning loss of local genetic diversity, but also indicates the likely response of populations to genetic reinforcement from neighbouring subspecies, increasing management options and encouraging a range-wide species conservation approach.

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.

DNA barcoding validates species labelling of certified seafood

Seafood is one of the most traded food commodities in the world with demand steadily increasing [1]. There is, however, a rising concern over the vulnerability of seafood supply chains to species mislabelling and fraud [1,2]. DNA methods have been widely used to detect species mislabelling and a recent meta-analysis of 4500 seafood product tests from 51 publications found an average of 30 percent were not the species stated on the label or menu [3]. This high rate poses a serious threat to consumer trust, reputations of seafood businesses and the sustainability of fishery resources. Seafood certification schemes may help reduce this problem. Here, we use DNA barcoding [4] to validate the species identity of 1402 certified seafood products derived from 27 species across 18 countries and find that in over 99% of cases species labelling was correct.

Museum specimens provide reliable SNP data for population genomic analysis of a widely distributed but threatened cockatoo species

Natural history museums harbour a plethora of biological specimens which are of potential use in population and conservation genetic studies. Although technical advancements in museum genomics have enabled genome-wide markers to be generated from aged museum specimens, the suitability of these data for robust biological inference is not well characterized. The aim of this study was to test the utility of museum specimens in population and conservation genomics by assessing the biological and technical validity of single nucleotide polymorphism (SNP) data derived from such samples. To achieve this, we generated thousands of SNPs from 47 red-tailed black cockatoo (Calyptorhychus banksii) traditional museum samples (i.e. samples that were not collected with the primary intent of DNA analysis) and 113 fresh tissue samples (cryopreserved liver/muscle) using a restriction site-associated DNA marker approach (DArTseq^™ ). Thousands of SNPs were successfully generated from most of the traditional museum samples (with a mean age of 44 years, ranging from 5 to 123 years), although 38% did not provide useful data. These SNPs exhibited higher error rates and contained significantly more missing data compared with SNPs from fresh tissue samples, likely due to considerable DNA fragmentation. However, based on simulation results, the level of genotyping error had a negligible effect on inference of population structure in this species. We did identify a bias towards low diversity SNPs in older samples that appears to compromise temporal inferences of genetic diversity. This study demonstrates the utility of a RADseq-based method to produce reliable genome-wide SNP data from traditional museum specimens.

How many species of giraffe are there?

In a recent paper in Current Biology, Fennessy and colleagues [1] conclude that there are four species of giraffe and that their numbers are declining in Africa. Giraffes (Giraffa camelopardalis) are presently classified as one species, with nine subspecies, which are considered 'Vulnerable' on the IUCN Red List [2]. The present consensus of one species divided into nine subspecies has previously been questioned (Supplemental information), and Fennessy and colleagues [1] provide another viewpoint on giraffe taxonomy. The fundamental reason for different taxonomic interpretations is that they are based upon different datasets that adopt different statistical techniques and follow different criteria for nomenclature.

Defining end user requirements for a field-based molecular detection system for wildlife forensic investigations

The increasing use of non-laboratory-based DNA and protein detection methods promise to provide rapid investigative intelligence and support sample prioritisation. Primarily developed for human forensic or medical applications, current systems may also show utility in the field of wildlife forensic science. However, it is currently unknown whether the requirements of the wildlife forensic community can be met by current non-laboratory based tools. Given the diverse array of stakeholders and sample types commonly encountered, it is necessary to first identify the needs of the community and then try and map their needs to current instrumentation. By using a market research style questionnaire, this study identified key requirements for a non-laboratory-based system following feedback from the wildlife forensic community. Data showed that there is strong support for field-based detection methods while highlighting concerns including contamination risks and reduced quality assurance associated with non-laboratory testing. Key species and applications were identified alongside hurdles to implementation and adoption. Broadly, the requirements align with many of the developmental drivers that have led to the rise of in-field portable detection instrumentation, specifically rapid detection within one hour, ease-of-use, and ≥95% accuracy. Several existing platforms exist that met some of the identified requirements but not all. With further collaboration between industry partners and the wildlife forensic community it is possible that new field-based systems can be developed and applied routinely.

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.

Parallel evolution and adaptation to environmental factors in a marine flatfish: Implications for fisheries and aquaculture management of the turbot (Scophthalmus maximus)

Unraveling adaptive genetic variation represents, in addition to the estimate of population demographic parameters, a cornerstone for the management of aquatic natural living resources, which, in turn, represent the raw material for breeding programs. The turbot (Scophthalmus maximus) is a marine flatfish of high commercial value living on the European continental shelf. While wild populations are declining, aquaculture is flourishing in southern Europe. We evaluated the genetic structure of turbot throughout its natural distribution range (672 individuals; 20 populations) by analyzing allele frequency data from 755 single nucleotide polymorphism discovered and genotyped by double-digest RAD sequencing. The species was structured into four main regions: Baltic Sea, Atlantic Ocean, Adriatic Sea, and Black Sea, with subtle differentiation apparent at the distribution margins of the Atlantic region. Genetic diversity and effective population size estimates were highest in the Atlantic populations, the area of greatest occurrence, while turbot from other regions showed lower levels, reflecting geographical isolation and reduced abundance. Divergent selection was detected within and between the Atlantic Ocean and Baltic Sea regions, and also when comparing these two regions with the Black Sea. Evidence of parallel evolution was detected between the two low salinity regions, the Baltic and Black seas. Correlation between genetic and environmental variation indicated that temperature and salinity were probably the main environmental drivers of selection. Mining around the four genomic regions consistently inferred to be under selection identified candidate genes related to osmoregulation, growth, and resistance to diseases. The new insights are useful for the management of turbot fisheries and aquaculture by providing the baseline for evaluating the consequences of turbot releases from restocking and farming.

Performance and precision of double digestion RAD (ddRAD) genotyping in large multiplexed datasets of marine fish species

The development of Genotyping-By-Sequencing (GBS) technologies enables cost-effective analysis of large numbers of Single Nucleotide Polymorphisms (SNPs), especially in "non-model" species. Nevertheless, as such technologies enter a mature phase, biases and errors inherent to GBS are becoming evident. Here, we evaluated the performance of double digest Restriction enzyme Associated DNA (ddRAD) sequencing in SNP genotyping studies including high number of samples. Datasets of sequence data were generated from three marine teleost species (>5500 samples, >2.5 × 10¹² bases in total), using a standardized protocol. A common bioinformatics pipeline based on STACKS was established, with and without the use of a reference genome. We performed analyses throughout the production and analysis of ddRAD data in order to explore (i) the loss of information due to heterogeneous raw read number across samples; (ii) the discrepancy between expected and observed tag length and coverage; (iii) the performances of reference based vs. de novo approaches; (iv) the sources of potential genotyping errors of the library preparation/bioinformatics protocol, by comparing technical replicates. Our results showed use of a reference genome and a posteriori genotype correction improved genotyping precision. Individual read coverage was a key variable for reproducibility; variance in sequencing depth between loci in the same individual was also identified as an important factor and found to correlate to tag length. A comparison of downstream analysis carried out with ddRAD vs single SNP allele specific assay genotypes provided information about the levels of genotyping imprecision that can have a significant impact on allele frequency estimations and population assignment. The results and insights presented here will help to select and improve approaches to the analysis of large datasets based on RAD-like methodologies.

Single-nucleotide polymorphism discovery and panel characterization in the African forest elephant

The continuing decline in forest elephant (Loxodonta cyclotis) numbers due to poaching and habitat reduction is driving the search for new tools to inform management and conservation. For dense rainforest species, basic ecological data on populations and threats can be challenging and expensive to collect, impeding conservation action in the field. As such, genetic monitoring is being increasingly implemented to complement or replace more burdensome field techniques. Single-nucleotide polymorphisms (SNPs) are particularly cost-effective and informative markers that can be used for a range of practical applications, including population census, assessment of human impact on social and genetic structure, and investigation of the illegal wildlife trade. SNP resources for elephants are scarce, but next-generation sequencing provides the opportunity for rapid, inexpensive generation of SNP markers in nonmodel species. Here, we sourced forest elephant DNA from 23 samples collected from 10 locations within Gabon, Central Africa, and applied double-digest restriction-site-associated DNA (ddRAD) sequencing to discover 31,851 tags containing SNPs that were reduced to a set of 1,365 high-quality candidate SNP markers. A subset of 115 candidate SNPs was then selected for assay design and validation using 56 additional samples. Genotyping resulted in a high conversion rate (93%) and a low per allele error rate (0.07%). This study provides the first panel of 107 validated SNP markers for forest elephants. This resource presents great potential for new genetic tools to produce reliable data and underpin a step-change in conservation policies for this elusive species.

Phylogeography of the reticulated python (Malayopython reticulatus ssp.): Conservation implications for the worlds' most traded snake species

As an important economic natural resource in Southeast Asia, reticulated pythons (Malayopython reticulatus ssp.) are primarily harvested from the wild for their skins-which are prized in the luxury leather goods industry. Trade dynamics of this CITES Appendix II listed species are complex and management approaches on the country or regional level appear obscure. Little is known about the actual geographic point-of-harvest of snakes, how genetic diversity is partitioned across the species range, how current harvest levels may affect the genetic viability of populations, and whether genetic structure could (or should) be accounted for when managing harvest quotas. As an initial survey, we use mitochondrial sequence data to define the broad-scale geographic structure of genetic diversity across a significant portion of the reticulated python's native range. Preliminary results reveal: (1) prominent phylogenetic structure across populations east and west of Huxley's modification of Wallace's line. Thirty-four haplotypes were apportioned across two geographically distinct groups, estimated to be moderately (5.2%); (2) Philippine, Bornean and Sulawesian populations appear to cluster distinctly; (3) individuals from Ambon Island suggest recent human introduction. Malayopython reticulatus is currently managed as a single taxonomic unit across Southeast Asia yet these initial results may justify special management considerations of the Philippine populations as a phylogenetically distinct unit, that warrants further examination. In Indonesia, genetic structure does not conform tightly to political boundaries and therefore we advocate the precautionary designation and use of Evolutionary Significant Units within Malayopython reticulatus, to inform and guide regional adaptive management plans.

Contrasting results from molecular and pedigree-based population diversity measures in captive zebra highlight challenges facing genetic management of zoo populations

Zoo conservation breeding programs manage the retention of population genetic diversity through analysis of pedigree records. The range of demographic and genetic indices determined through pedigree analysis programs allows the conservation of diversity to be monitored relative to the particular founder population for a species. Such approaches are based on a number of well-documented founder assumptions, however without knowledge of actual molecular genetic diversity there is a risk that pedigree-based measures will be misinterpreted and population genetic diversity misunderstood. We examined the genetic diversity of the captive populations of Grevy's zebra, Hartmann's mountain zebra and plains zebra in Japan and the United Kingdom through analysis of mitochondrial DNA sequences. Very low nucleotide variability was observed in Grevy's zebra. The results were evaluated with respect to current and historic diversity in the wild, and indicate that low genetic diversity in the captive population is likely a result of low founder diversity, which in turn suggests relatively low wild genetic diversity prior to recent population declines. Comparison of molecular genetic diversity measures with analogous diversity indices generated from the studbook data for Grevy's zebra and Hartmann's mountain zebra show contrasting patterns, with Grevy's zebra displaying markedly less molecular diversity than mountain zebra, despite studbook analysis indicating that the Grevy's zebra population has substantially more founders, greater effective population size, lower mean kinship, and has suffered less loss of gene diversity. These findings emphasize the need to validate theoretical estimates of genetic diversity in captive breeding programs with empirical molecular genetic data. Zoo Biol. 36:87-94, 2017. © 2016 Wiley Periodicals, Inc.

Genome-wide DNA markers to support genetic management for domestication and commercial production in a large rodent, the Ghanaian grasscutter (Thryonomys swinderianus)

Domestication and commercial production of the grasscutter, Thryonomys swinderianus, a large rodent, represents an important opportunity to secure sustainable animal protein for local communities in West Africa. To support production, DNA markers are required for population diversity assessment, pedigree analysis and marker-assisted selection. This study reports the application of double-digest RAD sequencing to simultaneously discover and genotype SNP markers in 24 wild and recently domesticated grasscutters. An initial panel of 1209 SNP loci was characterised from a total of more than 21 000 candidate loci containing single SNPs. This genome-wide resource represents the first application of its type to commercial production of a large rodent for food and advances the use of agricultural genomics in Ghana.

Splitting or lumping? A conservation dilemma exemplified by the critically endangered dama gazelle (Nanger dama)

Managers of threatened species often face the dilemma of whether to keep populations separate to conserve local adaptations and minimize the risk of outbreeding, or whether to manage populations jointly to reduce loss of genetic diversity and minimise inbreeding. In this study we examine genetic relatedness and diversity in three of the five last remaining wild populations of dama gazelle and a number of captive populations, using mtDNA control region and cytochrome b data. Despite the sampled populations belonging to the three putative subspecies, which are delineated according to phenotypes and geographical location, we find limited evidence for phylogeographical structure within the data and no genetic support for the putative subspecies. In the light of these data we discuss the relevance of inbreeding depression, outbreeding depression, adaptive variation, genetic drift, and phenotypic variation to the conservation of the dama gazelle and make some recommendations for its future conservation management. The genetic data suggest that the best conservation approach is to view the dama gazelle as a single species without subspecific divisions.

Nuclear and mitochondrial genetic structure in the Eurasian beaver (Castor fiber) - implications for future reintroductions

Many reintroduction projects for conservation fail, and there are a large number of factors that may contribute to failure. Genetic analysis can be used to help stack the odds of a reintroduction in favour of success, by conducting assessment of source populations to evaluate the possibility of inbreeding and outbreeding depression and by conducting postrelease monitoring. In this study, we use a panel of 306 SNP (single nucleotide polymorphism) markers and 487-489 base pairs of mitochondrial DNA control region sequence data to examine 321 individuals from possible source populations of the Eurasian beaver for a reintroduction to Scotland. We use this information to reassess the phylogenetic history of the Eurasian beavers, to examine the genetic legacy of past reintroductions on the Eurasian landmass and to assess the future power of the genetic markers to conduct ongoing monitoring via parentage analysis and individual identification. We demonstrate the capacity of medium density genetic data (hundreds of SNPs) to provide information suitable for applied conservation and discuss the difficulty of balancing the need for high genetic diversity against phylogenetic best fit when choosing source population(s) for reintroduction.

Is population structure in the European white stork determined by flyway permeability rather than translocation history?

European white stork are long considered to diverge to eastern and western migration pools as a result of independent overwintering flyways. In relatively recent times, the western and northern distribution has been subject to dramatic population declines and country-specific extirpations. A number of independent reintroduction programs were started in the mid 1950s to bring storks back to historical ranges. Founder individuals were sourced opportunistically from the Eastern and Western European distributions and Algeria, leading to significant artificial mixing between eastern and western flyways. Here we use mitochondrial and microsatellite DNA to test the contention that prior to translocation, eastern and western flyways were genetically distinct. The data show a surprising lack of structure at any spatial or temporal scale suggesting that even though birds were moved between flyways, there is evidence of natural mixing prior to the onset of translocation activities. Overall a high retention of genetic diversity, high N_ef, and an apparent absence of recent genetic bottleneck associated with early 20th century declines suggest that the species is well equipped to respond to future environmental pressures.

Cross-amplification and validation of SNPs conserved over 44 million years between seals and dogs

High-density SNP arrays developed for humans and their companion species provide a rapid and convenient tool for generating SNP data in closely-related non-model organisms, but have not yet been widely applied to phylogenetically divergent taxa. Consequently, we used the CanineHD BeadChip to genotype 24 Antarctic fur seal (Arctocephalus gazella) individuals. Despite seals and dogs having diverged around 44 million years ago, 33,324 out of 173,662 loci (19.2%) could be genotyped, of which 173 were polymorphic and clearly interpretable. Two SNPs were validated using KASP genotyping assays, with the resulting genotypes being 100% concordant with those obtained from the high-density array. Two loci were also confirmed through in silico visualisation after mapping them to the fur seal transcriptome. Polymorphic SNPs were distributed broadly throughout the dog genome and did not differ significantly in proximity to genes from either monomorphic SNPs or those that failed to cross-amplify in seals. However, the nearest genes to polymorphic SNPs were significantly enriched for functional annotations relating to energy metabolism, suggesting a possible bias towards conserved regions of the genome.

Signatures of diversifying selection in European pig breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.

The domestic pig, an important source of protein worldwide, was domesticated from the ancestral wild boar in multiple locations throughout the world. In Europe, local types were developed following domestication, but phenotypically distinct breeds only arose in the eighteenth century with the advent of systematic breeding. Recently developed molecular tools for pigs (as well as other livestock species) now allow a genetic characterisation of breed histories, including identification of regions of the genome that have been under selection in the establishment of breeds. We have applied these tools to identify genomic regions associated with breed development in a set of commercial and traditional pig breeds. We found strong evidence of genetic differentiation between breeds near genes associated with traits that are used to define breed standards, such as ear morphology and coat colour, as well as in regions of the genome that are associated with pork production traits. It is well documented that crosses with Asian pigs have been used to modify European breeds. We have found evidence of genetic influence from Asian pigs in European breeds, again in regions of the genome associated with breed standard characteristics, including ear shape and coat colour, as well as production traits.