Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod.

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.

Using temporal genomics to understand contemporary climate change responses in wildlife.

Monitoring the evolutionary responses of species to ongoing global climate change is critical for informing conservation. Population genomic studies that use samples from multiple time points ("temporal genomics") are uniquely able to make direct observations of change over time. Consequently, only temporal studies can show genetic erosion or spatiotemporal changes in population structure. Temporal genomic studies directly examining climate change effects are currently rare but will likely increase in the coming years due to their high conservation value. Here, we highlight four key genetic indicators that can be monitored using temporal genomics to understand how species are responding to climate change. All indicators crucially rely on having a suitable baseline that accurately represents the past condition of the population, and we discuss aspects of study design that must be considered to achieve this.

Temporal Monitoring of the Floreana Island Galapagos Giant Tortoise Captive Breeding Program.

Captive breeding programs benefit from genetic analyses that identify relatedness between individuals, assign parentage to offspring, and track levels of genetic diversity. Monitoring these parameters across breeding cycles is critical to the success of a captive breeding program as it allows conservation managers to iteratively evaluate and adjust program structure. However, in practice, genetic tracking of breeding outcomes is rarely conducted. Here, we examined the first three offspring cohorts (2017-2020) of the genetically informed captive breeding program for the Floreana Island Galapagos giant tortoise, Chelonoidis niger. This captive breeding program is unique as the Floreana tortoise has been extinct since the 1800s, but its genome persists, in part, in the form of living hybrids with the extant Volcano Wolf tortoise, C. becki. Breeding over the study period took place at the Galapagos National Park Directorate breeding facility in four corrals, each containing three females and two males. Using 17 microsatellite markers, we were able to assign parentage to 94 of the 98 offspring produced over the study period. We observe that despite the addition of more founders since the pilot breeding program, the effective population size remains low, and changes to the arrangements of breeding corrals may be necessary to encourage more equal reproductive output from the males. This study demonstrates the value of hybrids for species restoration and the importance of continually reassessing the outcomes of captive breeding.

The Galapagos giant tortoise Chelonoidis phantasticus is not extinct.

The status of the Fernandina Island Galapagos giant tortoise (Chelonoidis phantasticus) has been a mystery, with the species known from a single specimen collected in 1906. The discovery in 2019 of a female tortoise living on the island provided the opportunity to determine if the species lives on. By sequencing the genomes of both individuals and comparing them to all living species of Galapagos giant tortoises, here we show that the two known Fernandina tortoises are from the same lineage and distinct from all others. The whole genome phylogeny groups the Fernandina individuals within a monophyletic group containing all species with a saddleback carapace morphology and one semi-saddleback species. This grouping of the saddleback species is contrary to mitochondrial DNA phylogenies, which place the saddleback species across several clades. These results imply the continued existence of lineage long considered extinct, with a current known population size of a single individual.

A new lineage of Galapagos giant tortoises identified from museum samples.

The Galapagos Archipelago is recognized as a natural laboratory for studying evolutionary processes. San Cristóbal was one of the first islands colonized by tortoises, which radiated from there across the archipelago to inhabit 10 islands. Here, we sequenced the mitochondrial control region from six historical giant tortoises from San Cristóbal (five long deceased individuals found in a cave and one found alive during an expedition in 1906) and discovered that the five from the cave are from a clade that is distinct among known Galapagos giant tortoises but closely related to the species from Española and Pinta Islands. The haplotype of the individual collected alive in 1906 is in the same clade as the haplotype in the contemporary population. To search for traces of a second lineage in the contemporary population on San Cristóbal, we closely examined the population by sequencing the mitochondrial control region for 129 individuals and genotyping 70 of these for both 21 microsatellite loci and >12,000 genome-wide single nucleotide polymorphisms [SNPs]. Only a single mitochondrial haplotype was found, with no evidence to suggest substructure based on the nuclear markers. Given the geographic and temporal proximity of the two deeply divergent mitochondrial lineages in the historical samples, they were likely sympatric, raising the possibility that the lineages coexisted. Without the museum samples, this important discovery of an additional lineage of Galapagos giant tortoise would not have been possible, underscoring the value of such collections and providing insights into the early evolution of this iconic radiation.

Genotyping-in-thousands by sequencing (GT-seq) of noninvasive faecal and degraded samples: A new panel to enable ongoing monitoring of Canadian polar bear populations.

Genetic monitoring using noninvasive samples provides a complement or alternative to traditional population monitoring methods. However, next-generation sequencing approaches to monitoring typically require high quality DNA and the use of noninvasive samples (e.g., scat) is often challenged by poor DNA quality and contamination by nontarget species. One promising solution is a highly multiplexed sequencing approach called genotyping-in-thousands by sequencing (GT-seq), which can enable cost-efficient genomics-based monitoring for populations based on noninvasively collected samples. Here, we develop and validate a GT-seq panel of 324 single nucleotide polymorphisms (SNPs) optimized for genotyping of polar bears based on DNA from noninvasively collected faecal samples. We demonstrate (1) successful GT-seq genotyping of DNA from a range of sample sources, including successful genotyping (>50% loci) of 62.9% of noninvasively collected faecal samples determined to contain polar bear DNA; and (2) that we can reliably differentiate individuals, ascertain sex, assess relatedness, and resolve population structure of Canadian polar bear subpopulations based on a GT-seq panel of 324 SNPs. Our GT-seq data reveal spatial-genetic patterns similar to previous polar bear studies but at lesser cost per sample and through use of noninvasively collected samples, indicating the potential of this approach for population monitoring. This GT-seq panel provides the foundation for a noninvasive toolkit for polar bear monitoring and can contribute to community-based programmes - a framework which may serve as a model for wildlife conservation and management for species worldwide.

Ancient and historical DNA in conservation policy.

Although genetic diversity has been recognized as a key component of biodiversity since the first Convention on Biological Diversity (CBD) in 1993, it has rarely been included in conservation policies and regulations. Even less appreciated is the role that ancient and historical DNA (aDNA and hDNA, respectively) could play in unlocking the temporal dimension of genetic diversity, allowing key conservation issues to be resolved, including setting baselines for intraspecies genetic diversity, estimating changes in effective population size (Ne), and identifying the genealogical continuity of populations. Here, we discuss how genetic information from ancient and historical specimens can play a central role in preserving biodiversity and highlight specific conservation policies that could incorporate such data to help countries meet their CBD obligations.

Demographic history and patterns of molecular evolution from whole genome sequencing in the radiation of Galapagos giant tortoises.

Whole genome sequencing provides deep insights into the evolutionary history of a species, including patterns of diversity, signals of selection, and historical demography. When applied to closely related taxa with a wealth of background knowledge, population genomics provides a comparative context for interpreting population genetic summary statistics and comparing empirical results with the expectations of population genetic theory. The Galapagos giant tortoises (Chelonoidis spp.), an iconic rapid and recent radiation, offer such an opportunity. Here, we sequenced whole genomes from three individuals of the 12 extant lineages of Galapagos giant tortoise and estimate diversity measures and reconstruct changes in coalescent rate over time. We also compare the number of derived alleles in each lineage to infer how synonymous and nonsynonymous mutation accumulation rates correlate with population size and life history traits. Remarkably, we find that patterns of molecular evolution are similar within individuals of the same lineage, but can differ significantly among lineages, reinforcing the evolutionary distinctiveness of the Galapagos giant tortoise species. Notably, differences in mutation accumulation among lineages do not align with simple population genetic predictions, suggesting that the drivers of purifying selection are more complex than is currently appreciated. By integrating results from earlier population genetic and phylogeographic studies with new findings from the analysis of whole genomes, we provide the most in-depth insights to date on the evolution of Galapagos giant tortoises, and identify discrepancies between expectation from population genetic theory and empirical data that warrant further scrutiny.

Opportunities and challenges of macrogenetic studies.

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

Macrogenetic studies must not ignore limitations of genetic markers and scale.

Millette et al. (Ecology Letters, 2020, 23:55-67) reported no consistent worldwide anthropogenic effects on animal genetic diversity using repurposed mitochondrial DNA sequences. We reexamine data from this study, describe genetic marker and scale limitations which might lead to misinterpretations with conservation implications, and provide advice to improve future macrogenetic studies.

Canadian polar bear population structure using genome-wide markers.

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

A real-time PCR assay to accurately quantify polar bear DNA in fecal extracts.

DNA extracted from fecal samples contains DNA from the focal species, food, bacteria and pathogens. Most DNA quantification methods measure total DNA and cannot differentiate among sources. Despite the desirability of noninvasive fecal sampling for studying wildlife populations, low amounts of focal species DNA make it difficult to use for next-generation sequencing (NGS), where accurate DNA quantification is critical for normalization. Two factors are required prior to using fecal samples in NGS libraries: (1) an accurate quantification method for the amount of target DNA and (2) a determination of the relative amount of target DNA needed for successful single nucleotide polymorphism genotyping assays. Here, we address these needs by developing primers to amplify a 101 bp region of the nuclear F2 gene and a quantitative PCR (qPCR) assay that allows the accurate quantification of the amount of polar bear (Ursus maritimus) DNA in fecal extracts. We test the assay on pure polar bear DNA extracted from muscle tissue and find a high correlation between fluorometric and qPCR quantifications. The qPCR assay was also successfully used to quantify the amount of DNA derived from polar bears in fecal extractions. Orthologs of the F2 gene have been identified across vertebrates; thus, similar qPCR assays could be developed for other species to enable noninvasive studies.

Genotyping on the ark: A synthesis of genetic resources available for species in zoos.

Using molecular genetic information to guide population management can improve the sustainability of species in captivity. However, empirical population genetics has not been commonly applied to species management programs in zoos. One limitation may be the availability of genetic resources (e.g., markers, primers, etc.) for species held in zoos. To assess the extent to which species held in zoos have been studied using population genetics in the wild, we conducted a systematic literature review of close to 8,000 papers. We synthesized information on the availability and scale of population genetics studies across amphibian, bird, mammal, and reptile species held in zoos, and discussed their potential for informing ex situ management. We found that more than half of the species in zoos (52%) already have some genetic markers described in the literature specific for them, or a congeneric species, that could be further developed to aid the management of zoo populations, and the accumulation of these resources has been steady over the past decades. Furthermore, the proportion of species with genetic resources is even higher (62%) for species that are being managed through a formal breeding program in zoos. Our study provides encouraging results for captive program managers interested in integrating population genetics into ex situ management strategies.

Population genomics through time provides insights into the consequences of decline and rapid demographic recovery through head-starting in a Galapagos giant tortoise.

Population genetic theory related to the consequences of rapid population decline is well-developed, but there are very few empirical studies where sampling was conducted before and after a known bottleneck event. Such knowledge is of particular importance for species restoration, given links between genetic diversity and the probability of long-term persistence. To directly evaluate the relationship between current genetic diversity and past demographic events, we collected genome-wide single nucleotide polymorphism data from prebottleneck historical (c.1906) and postbottleneck contemporary (c.2014) samples of Pinzón giant tortoises (Chelonoidis duncanensis; n = 25 and 149 individuals, respectively) endemic to a single island in the Galapagos. Pinzón giant tortoises had a historically large population size that was reduced to just 150-200 individuals in the mid 20th century. Since then, Pinzón's tortoise population has recovered through an ex situ head-start programme in which eggs or pre-emergent individuals were collected from natural nests on the island, reared ex situ in captivity until they were 4-5 years old and subsequently repatriated. We found that the extent and distribution of genetic variation in the historical and contemporary samples were very similar, with the latter group not exhibiting the characteristic genetic patterns of recent population decline. No population structure was detected either spatially or temporally. We estimated an effective population size (N e) of 58 (95% CI = 50-69) for the postbottleneck population; no prebottleneck N e point estimate was attainable (95% CI = 39-infinity) likely due to the sample size being lower than the true N e. Overall, the historical sample provided a valuable benchmark for evaluating the head-start captive breeding programme, revealing high retention of genetic variation and no skew in representation despite the documented bottleneck event. Moreover, this work demonstrates the effectiveness of head-starting in rescuing the Pinzón giant tortoise from almost certain extinction.

Genome-Wide Assessment of Diversity and Divergence Among Extant Galapagos Giant Tortoise Species.

Genome-wide assessments allow for fuller characterization of genetic diversity, finer-scale population delineation, and better detection of demographically significant units to guide conservation compared with those based on "traditional" markers. Galapagos giant tortoises (Chelonoidis spp.) have long provided a case study for how evolutionary genetics may be applied to advance species conservation. Ongoing efforts to bolster tortoise populations, which have declined by 90%, have been informed by analyses of mitochondrial DNA sequence and microsatellite genotypic data, but could benefit from genome-wide markers. Taking this next step, we used double-digest restriction-site associated DNA sequencing to collect genotypic data at >26000 single nucleotide polymorphisms (SNPs) for 117 individuals representing all recognized extant Galapagos giant tortoise species. We then quantified genetic diversity, population structure, and compared results to estimates from mitochondrial DNA and microsatellite loci. Our analyses detected 12 genetic lineages concordant with the 11 named species as well as previously described structure within one species, C. becki. Furthermore, the SNPs provided increased resolution, detecting admixture in 4 individuals. SNP-based estimates of diversity and differentiation were significantly correlated with those derived from nuclear microsatellite loci and mitochondrial DNA sequences. The SNP toolkit presented here will serve as a resource for advancing efforts to understand tortoise evolution, species radiations, and aid conservation of the Galapagos tortoise species complex.

Temporal Mitogenomics of the Galapagos Giant Tortoise from Pinzón Reveals Potential Biases in Population Genetic Inference.

Empirical population genetic studies generally rely on sampling subsets of the population(s) of interest and of the nuclear or organellar genome targeted, assuming each is representative of the whole. Violations of these assumptions may impact population-level parameter estimation and lead to spurious inferences. Here, we used targeted capture to sequence the full mitochondrial genome from 123 individuals of the Galapagos giant tortoise endemic to Pinzón Island (Chelonoidis duncanensis) sampled at 2 time points pre- and postbottleneck (circa 1906 and 2014) to explicitly assess differences in diversity estimates and demographic reconstructions based on subsets of the mitochondrial genome versus the full sequences and to evaluate potential biases associated with diversity estimates and demographic reconstructions from postbottlenecked samples alone. Haplotypic diversities were equal between the temporal samples based on the full mitochondrial genome, but single gene estimates suggested either decreases or increases in diversity depending upon the region. Demographic reconstructions based on the full sequence were more similar between the temporal samples than those based on the control region alone, or a subset of 3 regions, where the trends in population size changes shifted in magnitude and direction between the temporal samples. In all cases, the estimated coalescent point was more distant for the historical than contemporary sample. In summary, our results empirically demonstrate the influence of sampling bias when interpreting population genetic patterns and punctuate the need for careful consideration of potentially conflicting evolutionary signal across the mitochondrial genome.

I-HEDGE: determining the optimum complementary sets of taxa for conservation using evolutionary isolation.

In the midst of the current biodiversity crisis, conservation efforts might profitably be directed towards ensuring that extinctions do not result in inordinate losses of evolutionary history. Numerous methods have been developed to evaluate the importance of species based on their contribution to total phylogenetic diversity on trees and networks, but existing methods fail to take complementarity into account, and thus cannot identify the best order or subset of taxa to protect. Here, we develop a novel iterative calculation of the heightened evolutionary distinctiveness and globally endangered metric (I-HEDGE) that produces the optimal ranked list for conservation prioritization, taking into account complementarity and based on both phylogenetic diversity and extinction probability. We applied this metric to a phylogenetic network based on mitochondrial control region data from extant and recently extinct giant Galápagos tortoises, a highly endangered group of closely related species. We found that the restoration of two extinct species (a project currently underway) will contribute the greatest gain in phylogenetic diversity, and present an ordered list of rankings that is the optimum complementarity set for conservation prioritization.

Founded: Genetic Reconstruction of Lineage Diversity and Kinship Informs Ex situ Conservation of Cuban Amazon Parrots (Amazona leucocephala).

Captive breeding is a widespread conservation strategy, yet such programs rarely include empirical genetic data for assessing management assumptions and meeting conservation goals. Cuban Amazon parrots (Amazona leucocephala) are considered vulnerable, and multiple on-island captive populations have been established from wild-caught and confiscated individuals of unknown ancestry. Here, we used mitochondrial haplotypic and nuclear genotypic data at 9 microsatellite loci to quantify the extent and distribution of genetic variation within and among captive populations in Zapata Swamp and Managua, Cuba, and to estimate kinship among breeders (n = 88). Using Bayesian clustering analysis, we detected 2 distinct clusters within the Zapata population, one of which was shared with Managua. Individuals from the cluster unique to Zapata possessed mitochondrial haplotypes with affinities to Cuban subspecies (A. l. leucocephala, A. l. palmarum); the shared cluster was similar, but also included haplotypes closely related to the subspecies restricted to Cayman Brac (A. l. hesterna). Overall mean kinship was low within each captive population (-0.026 to -0.012), with 19 and 11 recommended breeding pairs in Zapata and Managua, respectively, ranked according to mean kinship and informed by molecular sexing. Our results highlight the importance of understanding population history within ex situ management programs, while providing genetic information to directly inform Cuban parrot conservation.