The Threatened Species Imperative: Conservation assessments would benefit from population genomic insights

Genomic Diversity as a Key Conservation Criterion: Proof-of-Concept From Mammalian Whole-Genome Resequencing Data

Many international, national, state, and local organizations prioritize the ranking of threatened and endangered species to help direct conservation efforts. For example, the International Union for Conservation of Nature (IUCN) assesses the Green Status of species and publishes the influential Red List of threatened species. Unfortunately, such conservation yardsticks do not explicitly consider genetic or genomic diversity (GD), even though GD is positively associated with contemporary evolutionary fitness, individual viability, and with future evolutionary potential. To test whether populations of genome sequences could help improve conservation assessments, we estimated GD metrics from 82 publicly available mammalian datasets and examined their statistical association with attributes related to conservation. We also considered intrinsic biological factors, including trophic level and body mass, that could impact GD and quantified their relative influences. Our results identify key population GD metrics that are both reflective and predictive of IUCN conservation categories. Specifically, our analyses revealed that Watterson's theta (the population mutation rate) and autozygosity (a product of inbreeding) are associated with the current Red List categorization, likely because demographic declines that lead to "listing" decisions also reduce levels of standing genetic variation. We argue that by virtue of this relationship, conservation organizations like IUCN could leverage emerging genome sequence data to help categorize Red List threat rankings (especially in otherwise data-deficient species) and/or enhance Green Status assessments to establish a baseline for future population monitoring. Thus, our paper (1) outlines the theoretical and empirical justification for a new GD-based assessment criterion, (2) provides a bioinformatic pipeline for estimating GD from population genomic data, and (3) suggests an analytical framework that can be used to measure baseline GD while providing quantitative GD context for consideration by conservation authorities.

Habitat remediation followed by managed connectivity reduces unwanted changes in evolutionary trajectory of high extirpation risk populations

As we continue to convert green spaces into roadways and buildings, connectivity between populations and biodiversity will continue to decline. In threatened and endangered species, this trend is particularly concerning because the cessation of immigration can cause increased inbreeding and loss of genetic diversity, leading to lower adaptability and higher extirpation probabilities in these populations. Unfortunately, monitoring changes in genetic diversity from management actions such as assisted migration and predicting the extent of introduced genetic variation that is needed to prevent extirpation is difficult and costly in situ. Therefore, we designed an agent-based model to link population-wide genetic variability and the influx of unique alleles via immigration to population stability and extirpation outcomes. These models showed that management of connectivity can be critical in restoring at-risk populations and reducing the effects of inbreeding depression. However, the rescued populations were more similar to the migrant source population (average FST range 0.05-0.10) compared to the historical recipient population (average FST range 0.23-0.37). This means that these management actions not only recovered the populations from the effects of inbreeding depression, but they did so in a way that changed the evolutionary trajectory that was predicted and expected for these populations prior to the population crash. This change was most extreme in populations with the smallest population sizes, which are representative of critically endangered species that could reasonably be considered candidates for restored connectivity or translocation strategies. Understanding how these at-risk populations change in response to varying management interventions has broad implications for the long-term adaptability of these populations and can improve future efforts for protecting locally adapted allele complexes when connectivity is restored.

Whole-genome resequencing confirms the genetic effects of dams on an endangered fish Hemibagrus guttatus (Siluriformes: Bagridae): A case study in a tributary of the Pearl River

Dam construction in riverine ecosystems has fragmented natural aquatic habitats and has altered environmental conditions. As a result, damming has been demonstrated to threaten aquatic biodiversity by reducing species distribution ranges and hindering gene exchange, leading to the inability to adapt to environmental changes. Knowledge of the contemporary genetic diversity and genetic structure of fish populations that are separated by dams is vital to developing effective conservation strategies, particularly for endangered fish species. We chose the Lianjiang River, a tributary of the Pearl River, as a case study to assess the effects of dams on the genetic diversity and genetic structure of an endangered fish species, Hemibagrus guttatus, using whole-genome resequencing data from 63 fish samples. The results indicated low levels of genetic diversity, high levels of inbreeding and decreasing trend of effective population size in fragmented H. guttatus populations. In addition, there were significant genetic structure and genetic differentiation among populations, suggesting that the dams might have affected H. guttatus populations. Our findings may benefit management and conservation practices for this endangered species that is currently suffering from the effects of dam construction.

Conserving a threatened North American walnut: a chromosome-scale reference genome for butternut (Juglans cinerea)

With the advent of affordable and more accurate third-generation sequencing technologies, and the associated bioinformatic tools, it is now possible to sequence, assemble, and annotate more species of conservation concern than ever before. Juglans cinerea, commonly known as butternut or white walnut, is a member of the walnut family, native to the Eastern United States and Southeastern Canada. The species is currently listed as Endangered on the IUCN Red List due to decline from an invasive fungus known as Ophiognomonia clavigignenti-juglandacearum (Oc-j) that causes butternut canker. Oc-j creates visible sores on the trunks of the tree which essentially starves and slowly kills the tree. Natural resistance to this pathogen is rare. Conserving butternut is of utmost priority due to its critical ecosystem role and cultural significance. As part of an integrated undergraduate and graduate student training program in biodiversity and conservation genomics, the first reference genome for Juglans cinerea is described here. This chromosome-scale 539 Mb assembly was generated from over 100 × coverage of Oxford Nanopore long reads and scaffolded with the Juglans mandshurica genome. Scaffolding with a closely related species oriented and ordered the sequences in a manner more representative of the structure of the genome without altering the sequence. Comparisons with sequenced Juglandaceae revealed high levels of synteny and further supported J. cinerea's recent phylogenetic placement. Comparative assessment of gene family evolution revealed a significant number of contracting families, including several associated with biotic stress response.

A Chromosome-Level Genome Assembly and Annotation for the Clouded Apollo Butterfly (Parnassius mnemosyne): A Species of Global Conservation Concern

The clouded apollo (Parnassius mnemosyne) is a palearctic butterfly distributed over a large part of western Eurasia, but population declines and fragmentation have been observed in many parts of the range. The development of genomic tools can help to shed light on the genetic consequences of the decline and to make informed decisions about direct conservation actions. Here, we present a high-contiguity, chromosome-level genome assembly of a female clouded apollo butterfly and provide detailed annotations of genes and transposable elements. We find that the large genome (1.5 Gb) of the clouded apollo is extraordinarily repeat rich (73%). Despite that, the combination of sequencing techniques allowed us to assemble all chromosomes (nc = 29) to a high degree of completeness. The annotation resulted in a relatively high number of protein-coding genes (22,854) compared with other Lepidoptera, of which a large proportion (21,635) could be assigned functions based on homology with other species. A comparative analysis indicates that overall genome structure has been largely conserved, both within the genus and compared with the ancestral lepidopteran karyotype. The high-quality genome assembly and detailed annotation presented here will constitute an important tool for forthcoming efforts aimed at understanding the genetic consequences of fragmentation and decline, as well as for assessments of genetic diversity, population structure, inbreeding, and genetic load in the clouded apollo butterfly.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.