Design of SNP markers for Aldabra giant tortoises using low coverage ddRAD-seq

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two remaining giant tortoise species worldwide. Captive-bred A. gigantea are being used in rewilding projects in the Western Indian Ocean to functionally replace the extinct endemic giant tortoise species and restore degraded island ecosystems. Furthermore, large-scale translocations may become necessary as rising sea levels threaten the only wild population on the low-lying Aldabra Atoll. Critical management decisions would be greatly facilitated by insights on the genetic structure of breeding populations. We used a double-digest restriction-associated DNA sequencing (ddRAD-seq) approach to identify single nucleotide polymorphisms (SNP) among the wild population and two additional captive populations of A. gigantea. A total of 1674 unlinked, putatively neutral genome-wide SNPs were identified. The values of expected heterozygosity ranged from 0.33 to 0.5, whereas the minor allele frequency ranged from 0.20 to 0.5. These novel SNP markers will serve as useful tools for informing the conservation of A. gigantea

Microsatellite-based genotyping of MHC class II DRB1 gene in Iberian and Alpine ibex

In an analysis of a microsatellite locus (OLADRB1) linked to the MHC DRB1 gene of Iberian and Alpine ibex (Capra pyrenaica and Capra ibex), we detected strong linkage disequilibrium between both loci. The allele length polymorphism at OLADRB1 was unambiguously linked to a particular DRB1 allele. This allowed us to develop a DRB-STR matching method for both ibex species. Validation of the DRB-STR matching method was performed in 160 Iberian ibex from Spain and 98 Alpine ibex from Switzerland and Italy. This simple and relatively inexpensive protocol may find wide applications in a variety of research areas (e.g., mate choice, pathogen-driven selection) and in the biological conservation and management of the Western European ibex population

Purging and accumulation of genetic load in conservation

Our ability to assess the threat posed by the genetic load to small and declining populations has been greatly improved by advances in genome sequencing and computational approaches. Yet, considerable confusion remains around the definitions of the genetic load and its dynamics, and how they impact individual fitness and population viability. We illustrate how both selective purging and drift affect the distribution of deleterious mutations during population size decline and recovery. We show how this impacts the composition of the genetic load, and how this affects the extinction risk and recovery potential of populations. We propose a framework to examine load dynamics and advocate for the introduction of load estimates in the management of endangered populations

Disentangling adaptation from drift in bottlenecked and reintroduced populations of Alpine ibex

Identifying local adaptation in bottlenecked species is essential for conservation management. Selection detection methods have an important role in species management plans, assessments of adaptive capacity, and looking for responses to climate change. Yet, the allele frequency changes exploited in selection detection methods are similar to those caused by the strong neutral genetic drift expected during a bottleneck. Consequently, it is often unclear what accuracy selection detection methods have across bottlenecked populations. In this study, simulations were used to explore if signals of selection could be confidently distinguished from genetic drift across 23 bottlenecked and reintroduced populations of Alpine ibex (Capra ibex). The meticulously recorded demographic history of the Alpine ibex was used to generate comprehensive simulated SNP data. The simulated SNPs were then used to benchmark the confidence we could place in outliers identified in empirical Alpine ibex RADseq derived SNP data. Within the simulated data set, the false positive rates were high for all selection detection methods (FST outlier scans and Genetic-Environment Association analyses) but fell substantially when two or more methods were combined. True positive rates were consistently low and became negligible with increased stringency. Despite finding many outlier loci in the empirical Alpine ibex SNPs, none could be distinguished from genetic drift-driven false positives. Unfortunately, the low true positive rate also prevents the exclusion of recent local adaptation within the Alpine ibex. The baselines and stringent approach outlined here should be applied to other bottlenecked species to ensure the risk of false positive, or negative, signals of selection are accounted for in conservation management plans

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is considered vulnerable due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts focused on both wild and ex-situ populations. A high-quality genome would also open avenues to investigate the genetic basis of the exceptionally long lifespan. Here, we produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture (Hi-C). We produced a 2.37 Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNAseq-assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. We also performed a low-coverage re-sequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. The high-quality chromosome-level reference genome for A. gigantea is one of the most complete turtle genomes available. It is a powerful tool to assess the population structure in the wild population and reveal the geographic origins of ex-situ individuals relevant for genetic diversity management and rewilding efforts

Ancient mitochondrial and modern whole genomes unravel massive genetic diversity loss during near extinction of Alpine ibex

Population bottlenecks can have dramatic consequences for the health and long-term survival of a species. Understanding of historic population size and standing genetic variation prior to a contraction allows estimating the impact of a bottleneck on the species' genetic diversity. Although historic population sizes can be modelled based on extant genomics, uncertainty is high for the last 10–20 millenia. Hence, integrating ancient genomes provides a powerful complement to retrace the evolution of genetic diversity through population fluctuations. Here, we recover 15 high-quality mitogenomes of the once nearly extinct Alpine ibex spanning 8601 BP to 1919 CE and combine these with 60 published modern whole genomes. Coalescent demography simulations based on modern whole genomes indicate population fluctuations coinciding with the last major glaciation period. Using our ancient and historic mitogenomes, we investigate the more recent demographic history of the species and show that mitochondrial haplotype diversity was reduced to a fifth of the prebottleneck diversity with several highly differentiated mitochondrial lineages having coexisted historically. The main collapse of mitochondrial diversity coincides with elevated human population growth during the last 1–2 kya. After recovery, one lineage was spread and nearly fixed across the Alps due to recolonization efforts. Our study highlights that a combined approach integrating genomic data of ancient, historic and extant populations unravels major long-term population fluctuations from the emergence of a species through its near extinction up to the recent past

Purging and accumulation of genetic load in conservation

Our ability to assess the threat posed by the genetic load to small and declining populations has been greatly improved by advances in genome sequencing and computational approaches. Yet, considerable confusion remains around the definitions of the genetic load and its dynamics, and how they impact individual fitness and population viability. We illustrate how both selective purging and drift affect the distribution of deleterious mutations during population size decline and recovery. We show how this impacts the composition of the genetic load, and how this affects the extinction risk and recovery potential of populations. We propose a framework to examine load dynamics and advocate for the introduction of load estimates in the management of endangered populations

Genomic Evidence for Cryptic Speciation in Tree Frogs From the Apennine Peninsula, With Description of Hyla perrini sp. nov

Despite increasing appreciation of the speciation continuum, delimiting and describing new species is a major yet necessary challenge of modern phylogeography to help optimize conservation efforts. In amphibians, the lack of phenotypic differences between closely-related taxa, their complex, sometimes unresolved phylogenetic relationships, and their potential to hybridize all act to blur taxonomic boundaries. Here we implement a multi-disciplinary approach to evaluate the nature of two deeply-diverged mitochondrial lineages previously documented in Italian tree frogs (Hyla intermedia s. l.), distributed north and south of the Northern Apennine Mountains. Based on evidence from mitochondrial phylogenetics, nuclear phylogenomics, hybrid zone population genomics, niche modeling analyses, and biometric assessments, we propose that these lineages be considered distinct, cryptic species. Both mitochondrial and nuclear data affirm that they belong to two monophyletic clades of Pliocene divergence (~3.5 My), only admixing over a relatively narrow contact zone restricted to the southeast of the Po Plain (50–100 km). These characteristics are comparable to similarly-studied parapatric amphibians bearing a specific status. Inferred from their current geographic distribution, the two Italian tree frogs feature distinct ecological niches (<15% of niche overlap), raising questions regarding potential adaptive components contributing to their incipient speciation. However, we found no diagnostic morphological and bioacoustic differences between them. This system illustrates the speciation continuum of Western-Palearctic tree frogs and identifies additional cryptic lineages of similar divergence to be treated as separate species (H. cf. meridionalis). We recommend combined approaches using genomic data as applied here for the future taxonomic assessment of cryptic diversity in alloparapatric radiations of terrestrial vertebrates, especially in controversial taxa. Finally, we formally described the northern Italian tree frogs as a new species, Hyla perrini sp. nov

raw SNP file

This is the raw SNP file used for all analyses. It was produced using GATK and pre-filtered using vcftools (i.e. minimal genotyping rate of 0.45 and a maximal normalized coverage of 2). This includes SNPs which are not within 90bp of cut sites (hence SNPs located on R2 reads are included here). For analyses presented in the manuscript, further filters were applied. See methods part in the manuscript