Genetic monitoring and complex population dynamics: insights from a 12-year study of the Rio Grande silvery minnow

The positive impact of conservation action

Governments recently adopted new global targets to halt and reverse the loss of biodiversity. It is therefore crucial to understand the outcomes of conservation actions. We conducted a global meta-analysis of 186 studies (including 665 trials) that measured biodiversity over time and compared outcomes under conservation action with a suitable counterfactual of no action. We find that in two-thirds of cases, conservation either improved the state of biodiversity or at least slowed declines. Specifically, we find that interventions targeted at species and ecosystems, such as invasive species control, habitat loss reduction and restoration, protected areas, and sustainable management, are highly effective and have large effect sizes. This provides the strongest evidence to date that conservation actions are successful but require transformational scaling up to meet global targets.

Genetic erosion in an endangered desert fish during a megadrought despite long-term supportive breeding

Human water use combined with a recent megadrought have reduced river and stream flow through the southwest United States and led to periodic drying of formerly perennial river segments. Reductions in snowmelt runoff and increased extent of drying collectively threaten short-lived, obligate aquatic species, including the endangered Rio Grande silvery minnow (Hybognathus amarus). This species is subject to boom-and-bust population dynamics, under which large fluctuations in abundance are expected to lower estimates of effective population size and erode genetic diversity over time. Rates of diversity loss are also affected by additions of hatchery-origin fish used to supplement the wild population. We used demographic and genetic data from wild and hatchery individuals to examine the relationship of genetic diversity and effective population size to abundance over the last two decades. Genetic diversity was low during the early 2000s, but diversity and demographic metrics stabilized after the hatchery program was initiated and environmental conditions improved. Yet, from 2017 onward, allelic diversity declined (Cohen's d = 1.34) and remained low despite hatchery stocking and brief wild population recovery. Across the time series, single-sample estimates of effective population size based on linkage disequilibrium (LD Ne ) were positively associated (r = 0.53) with wild abundance and total abundance, but as the proportion of hatchery-origin spawners increased, LD Ne declined (r = -0.55). Megadrought limited wild spawner abundance and precluded refreshment of hatchery brood stocks with wild fish; hence, we predict a riverine population increasingly dominated by hatchery-origin individuals and accelerated loss of genetic diversity despite supplementation. We recommend an adaptive and accelerated management plan that integrates river flow management and hatchery operations to slow the pace of genetic diversity loss exacerbated by megadrought.

Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae)

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

Population genomic monitoring provides insight into conservation status but no correlation with demographic estimates of extinction risk in a threatened trout

The current extinction crisis requires effective assessment and monitoring tools. Genetic approaches are appealing given the relative ease of field sampling required to estimate genetic diversity characteristics assumed related to population size, evolutionary potential, and extinction risk, and to evaluate hybridization with non-native species simultaneously. However, linkages between population genetic metrics of diversity from survey-style field collections and demographic estimates of population size and extinction risk are still in need of empirical examples, especially for remotely distributed species of conservation concern where the approach might be most beneficial. We capitalized on an exceptional opportunity to evaluate congruence between genetic diversity metrics and demographic-based estimates of abundance and extinction risk from a comprehensive Multiple Population Viability Analysis (MPVA) in a threatened fish, the Lahontan cutthroat trout (LCT). We sequenced non-native trout reference samples and recently collected and archived tissue samples of most remaining populations of LCT (N = 60) and estimated common genetic assessment metrics, predicting minimal hybridization with non-native trout, low diversity, and declining diversity over time. We further hypothesized genetic metrics would correlate positively with MPVA-estimated abundance and negatively with extinction probability. We uncovered several instances of hybridization that pointed to immediate management needs. After removing hybridized individuals, cautious interpretation of low effective population sizes (2-63) suggested reduced evolutionary potential for many LCT populations. Other genetic metrics did not decline over time nor correlate with MPVA-based estimates of harmonic mean abundance or 30-year extinction probability. Our results demonstrate benefits of genetic monitoring for efficiently detecting hybridization and, though genetic results were disconnected from demographic assessment of conservation status, they suggest reduced evolutionary potential and likely a higher conservation risk than currently recognized for this threatened fish. We emphasize that genetic information provides essential complementary insight, in addition to demographic information, for evaluating species status.

Longitudinal monitoring of neutral and adaptive genomic diversity in a reintroduction

Restoration programs in the form of ex-situ breeding combined with reintroductions are becoming critical to counteract demographic declines and species losses. Such programs are increasingly using genetic management to improve conservation outcomes. However, the lack of long-term monitoring of genetic indicators following reintroduction prevents assessments of the trajectory and persistence of reintroduced populations. We carried out an extensive monitoring program in the wild for a threatened small-bodied fish (southern pygmy perch, Nannoperca australis) to assess the long-term genomic effects of its captive breeding and reintroduction. The species was rescued prior to its extirpation from the terminal lakes of Australia's Murray-Darling Basin, and then used for genetically informed captive breeding and reintroductions. Subsequent annual or biannual monitoring of abundance, fitness, and occupancy over a period of 11 years, combined with postreintroduction genetic sampling, revealed survival and recruitment of reintroduced fish. Genomic analyses based on data from the original wild rescued, captive born, and reintroduced cohorts revealed low inbreeding and strong maintenance of neutral and candidate adaptive genomic diversity across multiple generations. An increasing trend in the effective population size of the reintroduced population was consistent with field monitoring data in demonstrating successful re-establishment of the species. This provides a rare empirical example that the adaptive potential of a locally extinct population can be maintained during genetically informed ex-situ conservation breeding and reintroduction into the wild. Strategies to improve biodiversity restoration via ex-situ conservation should include genetic-based captive breeding and longitudinal monitoring of standing genomic variation in reintroduced populations.

A river runs through it: The causes, consequences, and management of intraspecific diversity in river networks

Rivers are fascinating ecosystems in which the eco-evolutionary dynamics of organisms are constrained by particular features, and biologists have developed a wealth of knowledge about freshwater biodiversity patterns. Over the last 10 years, our group used a holistic approach to contribute to this knowledge by focusing on the causes and consequences of intraspecific diversity in rivers. We conducted empirical works on temperate permanent rivers from southern France, and we broadened the scope of our findings using experiments, meta-analyses, and simulations. We demonstrated that intraspecific (genetic) diversity follows a spatial pattern (downstream increase in diversity) that is repeatable across taxa (from plants to vertebrates) and river systems. This pattern can result from interactive processes that we teased apart using appropriate simulation approaches. We further experimentally showed that intraspecific diversity matters for the functioning of river ecosystems. It indeed affects not only community dynamics, but also key ecosystem functions such as litter degradation. This means that losing intraspecific diversity in rivers can yield major ecological effects. Our work on the impact of multiple human stressors on intraspecific diversity revealed that-in the studied river systems-stocking of domestic (fish) strains strongly and consistently alters natural spatial patterns of diversity. It also highlighted the need for specific analytical tools to tease apart spurious from actual relationships in the wild. Finally, we developed original conservation strategies at the basin scale based on the systematic conservation planning framework that appeared pertinent for preserving intraspecific diversity in rivers. We identified several important research avenues that should further facilitate our understanding of patterns of local adaptation in rivers, the identification of processes sustaining intraspecific biodiversity-ecosystem function relationships, and the setting of reliable conservation plans.

A Conservation Hatchery Population of Delta Smelt Shows Evidence of Genetic Adaptation to Captivity After 9 Generations

Genetic adaptation to captivity is a concern for threatened and endangered species held in conservation hatcheries. Here, we present evidence of genetic adaptation to captivity in a conservation hatchery for the endangered delta smelt (Fish Conservation and Culture Laboratory, University of California Davis; FCCL). The FCCL population is genetically managed with parentage analysis and the addition of wild fish each year. Molecular monitoring indicates little loss of genetic variation and low differentiation between the wild and conservation populations. Yet, we found an increase in offspring survival to reproductive maturity during the subsequent spawning season (recovery rate) in crosses that included one or both cultured parents. Crosses with higher levels of hatchery ancestry tend to produce a greater number of offspring that are recovered the following year. The recovery rate of a cross decreases when offspring are raised in a tank with fish of high levels of hatchery ancestry. We suggest changes in fish rearing practices at the FCCL to reduce genetic adaptation to captivity, as delta smelt numbers in the wild continue to decline and the use of FCCL fish for reintroduction becomes more likely.

Relationship between effective and demographic population size in continuously distributed populations

Genetic monitoring of wild populations can offer insights into demographic and genetic information simultaneously. However, widespread application of genetic monitoring is hindered by large uncertainty in the estimation and interpretation of target metrics such as contemporary effective population size, Ne . We used four long-term genetic and demographic studies (≥9 years) to evaluate the temporal stability of the relationship between Ne and demographic population size (Nc ). These case studies focused on mammals that are continuously distributed, yet dispersal-limited within the spatial scale of the study. We estimated local, contemporary Ne with single-sample methods (LDNE, Heterozygosity Excess, and Molecular Ancestry) and demographic abundance with either mark-recapture estimates or catch-per-unit effort indices. Estimates of Ne varied widely within each case study suggesting interpretation of estimates is challenging. We found inconsistent correlations and trends both among estimates of Ne and between Ne and Nc suggesting the value of Ne as an indicator of Nc is limited in some cases. In the two case studies with consistent trends between Ne and Nc , FIS was more stable over time and lower, suggesting FIS may be a good indicator that the population was sampled at a spatial scale at which genetic structure is not biasing estimates of Ne . These results suggest that more empirical work on the estimation of Ne in continuous populations is needed to understand the appropriate context to use LDNe as a useful metric in a monitoring programme to detect temporal trends in either Ne or Nc .

Spatio-temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi-arid regions. Fish-parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species' adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio-temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio-temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.

Retention of Ancestral Genetic Variation Across Life-Stages of an Endangered, Long-Lived Iteroparous Fish

As with many endangered, long-lived iteroparous fishes, survival of razorback sucker depends on a management strategy that circumvents recruitment failure that results from predation by non-native fishes. In Lake Mohave, AZ-NV, management of razorback sucker centers on capture of larvae spawned in the lake, rearing them in off-channel habitats, and subsequent release ("repatriation") to the lake when adults are sufficiently large to resist predation. The effects of this strategy on genetic diversity, however, remained uncertain. After correction for differences in sample size among groups, metrics of mitochondrial DNA (mtDNA; number of haplotypes, N H , and haplotype diversity, H D ) and microsatellite (number of alleles, N A , and expected heterozygosity, H E ) diversity did not differ significantly between annual samples of repatriated adults and larval year-classes or among pooled samples of repatriated adults, larvae, and wild fish. These findings indicate that the current management program thus far maintained historical genetic variation of razorback sucker in the lake. Because effective population size, N e , is closely tied to the small census population size (N c = ~1500-3000) of razorback sucker in Lake Mohave, this population will remain at risk from genetic, as well as demographic risk of extinction unless N c is increased substantially.

Biogeography of "Cyprinella lutrensis": intensive genetic sampling from the Pecos River 'melting pot' reveals a dynamic history and phylogenetic complexity

Thorough sampling is necessary to delineate lineage diversity for polytypic "species" such as Cyprinella lutrensis. We conducted extensive mtDNA sampling (cytochrome b and ND4) from the Pecos River, Rio Grande, and South Canadian River, New Mexico. Our study emphasized the Pecos River due to its complex geological history and potential to harbor multiple lineages. We used geometric-morphometric, morphometric, and meristic analyses to test for phenotypic divergence and combined nucDNA with mtDNA to test for cytonuclear disequilibrium and combined our sequences with published data to conduct a phylogenetic re-assessment of the entire C. lutrensis clade. We detected five co-occurring mtDNA lineages in the Pecos River, but no evidence for cytonuclear disequilibrium or phenotypic divergence. Recognized species were interspersed amongst divergent lineages of "C. lutrensis". Allopatric divergence among drainages isolated in the Late Miocene and Pliocene apparently produced several recognized species and major divisions within "C. lutrensis". Pleistocene re-expansion and subsequent re-fragmentation of a centralized lineage founded younger, divergent lineages throughout the Rio Grande basin and Edwards Plateau. There is also evidence of recent introductions to the Rio Grande, Pecos and South Canadian Rivers. Nonetheless, deeply divergent lineages have coexisted since the Pleistocene.

Comparative riverscape genetics reveals reservoirs of genetic diversity for conservation and restoration of Great Plains fishes

We used comparative landscape genetics to examine the relative roles of historical events, intrinsic traits and landscape factors in determining the distribution of genetic diversity of river fishes across the North American Great Plains. Spatial patterns of diversity were overlaid on a patch-based graphical model and then compared within and among three species that co-occurred across five Great Plains watersheds. Species differing in reproductive strategy (benthic vs. pelagic-spawning) were hypothesized to have different patterns of genetic diversity, but the overriding factor shaping contemporary patterns of diversity was the signature of past climates and geological history. Allelic diversity was significantly higher at southern latitudes for Cyprinella lutrensis and Hybognathus placitus, consistent with northward expansion from southern Pleistocene refugia. Within the historical context, all species exhibited lowered occupancy and abundance in heavily fragmented and drier upstream reaches, particularly H. placitus; a pelagic-spawning species, suggesting rates of extirpation have outpaced losses of genetic diversity in this species. Within most tributary basins, genetically diverse populations of each species persisted. Hence, reconnecting genetically diverse populations with those characterized by reduced diversity (regardless of their position within the riverine network) would provide populations with greater genetic and demographic resilience. We discuss cases where cross-basin transfer may be appropriate to enhance genetic diversity and mitigate negative effects of climate change. Overall, striking similarities in genetic patterns and in response to fragmentation and dewatering suggest a common strategy for genetic resource management in this unique riverine fish assemblage.

Clock gene evolution: seasonal timing, phylogenetic signal, or functional constraint?

Genetic determinants of seasonal reproduction are not fully understood but may be important predictors of organism responses to climate change. We used a comparative approach to study the evolution of seasonal timing within a fish community in a natural common garden setting. We tested the hypothesis that allelic length variation in the PolyQ domain of a circadian rhythm gene, Clock1a, corresponded to interspecific differences in seasonal reproductive timing across 5 native and 1 introduced cyprinid fishes (n = 425 individuals) that co-occur in the Rio Grande, NM, USA. Most common allele lengths were longer in native species that initiated reproduction earlier (Spearman's r = -0.70, P = 0.23). Clock1a allele length exhibited strong phylogenetic signal and earlier spawners were evolutionarily derived. Aside from length variation in Clock1a, all other amino acids were identical across native species, suggesting functional constraint over evolutionary time. Interestingly, the endangered Rio Grande silvery minnow (Hybognathus amarus) exhibited less allelic variation in Clock1a and observed heterozygosity was 2- to 6-fold lower than the 5 other (nonimperiled) species. Reduced genetic variation in this functionally important gene may impede this species' capacity to respond to ongoing environmental change.

Time-series analysis reveals genetic responses to intensive management of razorback sucker (Xyrauchen texanus)

Time-series analysis is used widely in ecology to study complex phenomena and may have considerable potential to clarify relationships of genetic and demographic processes in natural and exploited populations. We explored the utility of this approach to evaluate population responses to management in razorback sucker, a long-lived and fecund, but declining freshwater fish species. A core population in Lake Mohave (Arizona-Nevada, USA) has experienced no natural recruitment for decades and is maintained by harvesting naturally produced larvae from the lake, rearing them in protective custody, and repatriating them at sizes less vulnerable to predation. Analyses of mtDNA and 15 microsatellites characterized for sequential larval cohorts collected over a 15-year time series revealed no changes in geographic structuring but indicated significant increase in mtDNA diversity for the entire population over time. Likewise, ratios of annual effective breeders to annual census size (N b /N a) increased significantly despite sevenfold reduction of N a. These results indicated that conservation actions diminished near-term extinction risk due to genetic factors and should now focus on increasing numbers of fish in Lake Mohave to ameliorate longer-term risks. More generally, time-series analysis permitted robust testing of trends in genetic diversity, despite low precision of some metrics.

Multidisciplinary population monitoring when demographic data are sparse: a case study of remote trout populations

The potential of genetic, genomic, and phenotypic metrics for monitoring population trends may be especially high in isolated regions, where traditional demographic monitoring is logistically difficult and only sporadic sampling is possible. This potential, however, is relatively underexplored empirically. Over eleven years, we assessed several such metrics along with traditional ecological knowledge and catch data in a socioeconomically important trout species occupying a large, remote lake. The data revealed largely stable characteristics in two populations over 2–3 generations, but possible contemporary changes in a third population. These potential shifts were suggested by reduced catch rates, reduced body size, and changes in selection implied at one gene-associated single nucleotide polymorphism. A demographic decline in this population, however, was ambiguously supported, based on the apparent lack of temporal change in effective population size, and corresponding traditional knowledge suggesting little change in catch. We illustrate how the pluralistic approach employed has practicality for setting future monitoring efforts of these populations, by guiding monitoring priorities according to the relative merits of different metrics and availability of resources. Our study also considers some advantages and disadvantages to adopting a pluralistic approach to population monitoring where demographic data are not easily obtained.

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species

Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30-year-long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo-genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo-genetic approaches should be the norm in conservation practices.

We combined genetic and demographic data from a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale for conservation purposes. Genetic diversity and effective population sizes are very low, probably due to the strong genetic bottlenecks detected in this study. The species spatial distribution and abundance also decreased during the last decades.

Genetic analysis of captive spawning strategies for the endangered Rio Grande Silvery Minnow

Captive breeding and rearing are central elements in conservation, management, and recovery planning for many endangered species including Rio Grande Silvery Minnow, a North American freshwater cyprinid. Traditionally, the sole purpose of hatcheries was to produce as many fish as feasible for stocking and harvest. Production quotas are also an important consideration in hatchery programs for endangered species, but they must also maintain and maximize genetic diversity of fish produced through implementation of best breeding practices. Here, we assessed genetic outcomes and measures of productivity (number of eggs and larval viability) for three replicates of three mating designs that are used for this small, pelagic-spawning fish. These were 1) monogamous mating, 2) hormone-induced communal spawning, and 3) environmentally cued communal spawning. A total of 180 broodstock and 450 progeny were genotyped. Genetic diversity and egg productivity did not differ significantly among spawning designs (H e : F = 0.52, P = 0.67; H o : F = 0.12, P = 0.89; number of eggs: F = 3.59, P = 0.09), and there was evidence for variance in reproductive success among individuals in all three designs. Allelic richness declined from the broodstock to progeny generation in all breeding designs. There was no significant difference in the genetic effective size (regardless of the method used) among designs. Significantly more viable eggs were produced in environmentally cued communal spawn compared to the alternative strategies (F = 5.72, P = 0.04), but this strategy is the most difficult to implement.

Evaluating the performance of captive breeding techniques for conservation hatcheries: a case study of the delta smelt captive breeding program

The delta smelt, an endangered fish species endemic to the San Francisco Bay-Delta, California, United States, was recently brought into captivity for species preservation. This study retrospectively evaluates the implementation of a genetic management plan for the captive delta smelt population. The captive genetic management plan entails tagging fish, molecular data collection, pedigree reconstruction, relatedness estimation, and recommending fish crosses annually in an effort to minimize the average coancestry in the population and limit inbreeding. We employed 12 microsatellite DNA markers to examine temporal genetic diversity in consecutive, discrete generations to determine the effects of intensive genetic management on the population and to quantify the amount of wild genetic diversity present within each captive generation. Wild fish are incorporated into the captive population each generation to minimize genetic drift, and 91% of the original founders are still represented in the F(3) generation. The average mean kinship in the third generation in captivity was 0.0035. There was no evidence of significant genetic divergence of the captive population from the wild population. The results of this study yield management insights into the practical application of genetic management plans for captive populations and conservation hatcheries, in an attempt to preserve the genetic integrity of endangered species.