Small Effective Population Size Research Articles

Beyond population size: whole-genome data reveal bottleneck legacies in the peninsular Italian wolf.

Preserving genetic diversity and adaptive potential while avoiding inbreeding depression is crucial for the long-term conservation of natural populations. Despite demographic increases, traces of past bottleneck events at the genomic level should be carefully considered for population management. From this perspective, the peninsular Italian wolf is a paradigmatic case. After being on the brink of extinction in the late 1960s, peninsular Italian wolves rebounded and recolonized most of the peninsula aided by conservation measures, including habitat and legal protection. Notwithstanding their demographic recovery, a comprehensive understanding of the genomic consequences of the historical bottleneck in Italian wolves is still lacking. To fill this gap, we sequenced whole genomes of thirteen individuals sampled in the core historical range of the species in Central Italy to conduct population genomic analyses, including a comparison with wolves from two highly-inbred wolf populations (i.e., Scandinavia and Isle Royale). We found that peninsular Italian wolves, despite their recent recovery, still exhibit relatively low genetic diversity, a small effective population size, signatures of inbreeding, and a non-negligible genetic load. Our findings indicate that the peninsular Italian wolf population is still susceptible to bottleneck legacies, which could lead to local inbreeding depression in case of population reduction or fragmentations. This study emphasizes the importance of considering key genetic parameters to design appropriate long-term conservation management plans.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

The genetic consequences of population marginality: A case study in maritime pine

AbstractAimMarginal tree populations, either those located at the edges of the species' range or in suboptimal environments, are often a valuable genetic resource for biological conservation. However, there is a lack of knowledge about the genetic consequences of population marginality, estimated across entire species' ranges. Our study addresses this gap by providing information about several genetic indicators and their variability in marginal and core populations identified using quantitative marginality indices.LocationSouthwestern Europe and North Africa.MethodsUsing 10,185 SNPs across 82 populations of maritime pine (Pinus pinaster Ait.), a widespread conifer characterised by a fragmented range, we modelled the relationship of seven genetic indicators potentially related to population evolutionary resilience, namely genetic diversity (based on both all SNPs and outlier SNPs), inbreeding, genetic differentiation, recessive genetic load and genomic offset, with population geographical, demo‐historical and ecological marginality (as estimated by nine quantitative indices). Models were constructed for both regional (introducing gene pool as a random factor) and range‐wide spatial scales.ResultsWe showed a trend towards decreasing overall genetic diversity and increasing differentiation with geographic marginality, supporting the centre‐periphery hypothesis (CPH). However, we found no correlation between population inbreeding and marginality, while geographically marginal populations had a lower recessive genetic load (only models without the gene pool effect). Ecologically marginal populations had a higher genomic offset, suggesting higher maladaptation to future climate, albeit some of these populations also had high genetic diversity for climate outliers.Main ConclusionsOverall genetic diversity (but not outlier‐based estimates) and differentiation patterns support the CPH. Ecologically marginal populations and those at the southern edge could be more vulnerable to climate change due to higher climate maladaptation, as predicted by genomic offsets, and/or lower potentially adaptive genetic diversity. This risk is exacerbated by typically small effective population sizes and increasing human impact in marginal populations.

Small towns limit dispersal and reduce genetic diversity in populations of Texas horned lizards.

There is a general expectation that urban populations will be fragmented and the movement of individuals will be restricted leading to low effective population sizes, low genetic diversity, higher inbreeding, and higher differentiation than populations living in more continuous habitat. In this study, we compare the genetic diversity and differentiation of Texas horned lizards that are found in four small towns (Kenedy, Karnes City, Rockdale, and Smithville) in Texas and at Tinker Air Force Base, Oklahoma to populations that occur in 16 natural areas and to an introduced population in South Carolina. We also present more detailed spatial genetic data and home range data for several of the towns. Texas horned lizards (Phrynosoma cornutum) living in small towns have lower genetic diversity, higher differentiation, and smaller effective population sizes than populations located in more natural areas. There was evidence for human-mediated movement of lizards into town; however, it has not been enough to counteract the effects of drift. Home range size is smaller in town than in more natural areas. Genetic patterns suggest dispersal occurs over short distances and is inhibited across areas with a high percent of impervious surface and major roads. These data suggest that effective planning to maintain suitable habitat and corridors to facilitate movement is critical to maintaining small terrestrial species like the Texas horned lizard and must be integrated into the early stages of urban development.

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS)

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species’ habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species’ genetic diversity as a resource for future breeding approaches in this medicinally important species.

Evaluation of heritability partitioning approaches in livestock populations

BackgroundHeritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium.ResultsHere, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively).ConclusionsHeritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used.

Drivers of genomic diversity and phenotypic development in early phases of domestication in Hermetia illucens.

The black soldier fly (BSF), Hermetia illucens, has the ability to efficiently bioremediate organic waste into usable bio-compounds. Understanding the impact of domestication and mass rearing on fitness and production traits is therefore important for sustainable production. This study aimed to assess patterns of genomic diversity and its association to phenotypic development across early generations of mass rearing under two selection strategies: selection for greater larval mass (SEL lines) and no direct artificial selection (NS lines). Genome-wide single nucleotide polymorphism (SNP) data were generated using 2bRAD sequencing, while phenotypic traits relating to production and population fitness were measured. Declining patterns of genomic diversity were observed across three generations of captive breeding, with the lowest diversity recorded for the F3 generation of both selection lines, most likely due to founder effects. The SEL cohort displayed statistically significantly greater larval weight com the NS lines with pronounced genetic and phenotypic directional changes across generations. Furthermore, lower genetic and phenotypic diversity, particularly for fitness traits, were evident for SEL lines, illustrating the trade-off between selecting for mass and the resulting decline in population fitness. SNP-based heritability was significant for growth, but was low or non-significant for fitness traits. Genotype-phenotype correlations were observed for traits, but individual locus effect sizes where small and very few of these loci demonstrated a signature for selection. Pronounced genetic drift, due to small effective population sizes, is likely overshadowing the impacts of selection on genomic diversity and consequently phenotypic development. The results hold particular relevance for genetic management and selective breeding for BSF in future.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Single nucleotide polymorphism profile for quantitative trait nucleotide in populations with small effective size and its impact on mapping and genomic predictions.

Increasing SNP density by incorporating sequence information only marginally increases prediction accuracies of breeding values in livestock. To find out why, we used statistical models and simulations to investigate the shape of distribution of estimated SNP effects (a profile) around quantitative trait nucleotides (QTNs) in populations with a small effective population size (Ne). A QTN profile created by averaging SNP effects around each QTN was similar to the shape of expected pairwise linkage disequilibrium (PLD) based on Ne and genetic distance between SNP, with a distinct peak for the QTN. Populations with smaller Ne showed lower but wider QTN profiles. However, adding more genotyped individuals with phenotypes dragged the profile closer to the QTN. The QTN profile was higher and narrower for populations with larger compared to smaller Ne. Assuming the PLD curve for the QTN profile, 80% of the additive genetic variance explained by each QTN was contained in ± 1/Ne Morgan interval around the QTN, corresponding to 2 Mb in cattle and 5 Mb in pigs and chickens. With such large intervals, identifying QTN is difficult even if all of them are in the data and the assumed genetic architecture is simplistic. Additional complexity in QTN detection arises from confounding of QTN profiles with signals due to relationships, overlapping profiles with closely spaced QTN, and spurious signals. However, small Ne allows for accurate predictions with large data even without QTN identification because QTNs are accounted for by QTN profiles if SNP density is sufficient to saturate the segments.

Conservation genetics of the endangered California Freshwater Shrimp (Syncaris pacifica): watershed and stream networks define gene pool boundaries

Understanding genetic structure and diversity among remnant populations of rare species can inform conservation and recovery actions. We used a population genetic framework to spatially delineate gene pools and estimate gene flow and effective population sizes for the endangered California Freshwater Shrimp Syncaris pacifica. Tissues of 101 individuals were collected from 11 sites in 5 watersheds, using non-lethal tissue sampling. Single Nucleotide Polymorphism markers were developed de novo using ddRAD-seq methods, resulting in 433 unlinked loci scored with high confidence and low missing data. We found evidence for strong genetic structure across the species range. Two hierarchical levels of significant differentiation were observed: (i) five clusters (regional gene pools, FST = 0.38–0.75) isolated by low gene flow were associated with watershed limits and (ii) modest local structure among tributaries within a watershed that are not connected through direct downstream flow (local gene pools, FST = 0.06–0.10). Sampling sites connected with direct upstream-to-downstream water flow were not differentiated. Our analyses suggest that regional watersheds are isolated from one another, with very limited (possibly no) gene flow over recent generations. This isolation is paired with small effective population sizes across regional gene pools (Ne = 62.4–147.1). Genetic diversity was variable across sites and watersheds (He = 0.09–0.22). Those with the highest diversity may have been refugia and are now potential sources of genetic diversity for other populations. These findings highlight which portions of the species range may be most vulnerable to future habitat fragmentation and provide management consideration for maintaining local effective population sizes and genetic connectivity.

Hard edges, soft edges, and species range evolution: A genomic analysis of the Cumberland Plateau salamander

AbstractAimGene flow from central to edge populations is thought to limit population growth at range edges by constraining local adaptation. In this study, we explore the thesis that range edges can differ in their dynamics and be either ‘hard’ (e.g. a river) or ‘soft’ (e.g. ecological gradients). We hypothesize that soft edge populations will have smaller effective population sizes than central populations and that gene flow will be greater from the centre to the edge than vice versa. Conversely, we hypothesize that hard edge populations should have similar effective population sizes to central populations and that gene flow will be equal between the two.LocationKentucky, West Virginia, and Virginia, USA.Taxon Plethodon kentucki (Caudata: Plethodontidae).MethodsWe evaluated landscape suitability using an ecological niche model, then we compared gene flow and effective population sizes between edge and central populations and quantified gene flow between populations. Finally, we characterized landscape genetic variation, testing for isolation by distance and isolation by environment.ResultsWe found continuously decreasing habitat quality along soft edges, with hard edges more variable. Additionally, we found that soft edges had lower effective population sizes than central populations and that gene flow was greater from the centre of the range to the soft edges than the reverse. In hard edges, by contrast, we found effective population sizes in edge populations were similar to central populations, with relatively equal gene flow in both directions.Main ConclusionsUnderstanding why species have range limits is central to investigations of the structure of biodiversity, yet the evolutionary dynamics of range edges remain poorly understood. We show that within a single species with a small range, the evolutionary dynamics operating at range boundaries may depend on the nature of the boundary.

Critically Small Contemporaneous Effective Population Sizes Estimated for Stocks of the African Bonytongue in Western Africa

Critically Small Contemporaneous Effective Population Sizes Estimated for Stocks of the African Bonytongue in Western Africa

Genomics-informed captive breeding can reduce inbreeding depression and the genetic load in zoo populations.

Zoo populations of threatened species are a valuable resource for the restoration of wild populations. However, their small effective population size poses a risk to long-term viability, especially in species with high genetic load. Recent bioinformatic developments can identify harmful genetic variants in genome data. Here, we advance this approach, analysing the genetic load in the threatened pink pigeon (Nesoenas mayeri). We lifted the mutation-impact scores that had been calculated for the chicken (Gallus gallus) to estimate the genetic load in six pink pigeons. Additionally, we perform in silico crossings to predict the genetic load and realized load of potential offspring. We thus identify the optimal mate pairs that are theoretically expected to produce offspring with the least inbreeding depression. We use computer simulations to show how genomics-informed conservation can reduce the genetic load whilst reducing the loss of genome-wide diversity. Genomics-informed management is likely to become instrumental in maintaining the long-term viability of zoo populations.

Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana)

The nonrandom distribution of chromosomal characteristics and functional elements—genomic architecture—impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure—ordination, clustering, and phylogenetic methods—we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

Population Genomic History of the Endangered Anatolian and Cyprian Mouflons in Relation to Worldwide Wild, Feral, and Domestic Sheep Lineages.

Once widespread in their homelands, the Anatolian mouflon (Ovis gmelini anatolica) and the Cyprian mouflon (Ovis gmelini ophion) were driven to near extinction during the 20th century and are currently listed as endangered populations by the International Union for Conservation of Nature. While the exact origins of these lineages remain unclear, they have been suggested to be close relatives of domestic sheep or remnants of proto-domestic sheep. Here, we study whole genome sequences of n = 5 Anatolian mouflons and n = 10 Cyprian mouflons in terms of population history and diversity, comparing them with eight other extant sheep lineages. We find reciprocal genetic affinity between Anatolian and Cyprian mouflons and domestic sheep, higher than all other studied wild sheep genomes, including the Iranian mouflon (O. gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne and mutation load estimates are lower in Cyprian compared with Anatolian mouflons, suggesting the purging of recessive deleterious variants in Cyprian sheep under a small long-term Ne, possibly attributable to founder effects, island isolation, introgression from domestic lineages, or differences in their bottleneck dynamics. Expanding our analyses to worldwide wild and feral Ovis genomes, we observe varying viability metrics among different lineages and a limited consistency between viability metrics and International Union for Conservation of Nature conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Founder populations of Stuckenia pectinata in restored ponds originate from genetically diverse propagule pool: A case study of urban ponds in Brussels, Belgium

Eutrophication and cyanobacteria blooms are considered major problems for the biodiversity and water quality of urban ponds. While biomanipulation techniques such as drawdown with fish removal have great potential to restore turbid ponds to a clearwater status, it remains difficult to predict if and how macrophytes will recover naturally. Here, we used individual genotyping and population genetics based on 20 nuclear microsatellite loci to investigate the recruitment and recolonization strategies of the submerged macrophyte Stuckenia pectinata (L.) Börner. More specifically, we compared the founder genetic diversity of recovering populations just after biomanipulation to the genetic diversity of spontaneous, contiguous populations that settled over an extended period of time and were within the same catchment. Our results showed that turbid ponds may contain a persistent propagule bank that allows for an immediate re-establishment of genetically diverse populations of S. pectinata once a desired clearwater state is restored. Therefore, biomanipulation without sediment removal proved to be successful for founding populations to become immediately integrated with their established populations, thus maintaining the overall diversity of this species within local areas. Additionally, our results demonstrated an excess of heterozygotes in established populations that may be caused by substantial drift in albeit small effective population sizes of this predominantly outbreeding species.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that ‘complex’ organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this ‘drift barrier’ theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans.

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Too simple, too complex, or just right? Advantages, challenges, and guidance for indicators of genetic diversity

Abstract Measuring genetic diversity of wild species using DNA-based data remains resource intensive and time consuming for nearly all species. However, genetic assessments are needed for global conservation commitments, including the Convention on Biological Diversity, and for governments and managers to evaluate conservation progress, as well as prioritizing species and populations to preserve and recover genetic diversity (e.g., via genetic rescue). Recently, indicators were developed for tracking and reporting genetic diversity status and trends for hundreds of species. The indicators quantify two simple proxies of within-population and among-population genetic diversity and adaptive potential: small effective population size (Ne &amp;lt; 500) and the loss of genetically distinct populations. The indicators must balance scientific credibility, practicality, and simplicity. In the present article, we summarize the advantages of these pragmatic indicators, address critiques by scientists for simplifying assumptions and by policymakers for complexity, and propose potential solutions and next steps. We aim to support practitioners putting indicators into policy, action, legislation, and reporting.

Patterns of genetic divergence in the Rio Grande cooter (Pseudemys gorzugi), a riverine turtle inhabiting an arid and anthropogenically modified system.

The lower Rio Grande and Pecos River of the southwest United States have been heavily modified by human activities, profoundly impacting the integrity of their aquatic wildlife. In this context, we focused our study on the population genomics of the Rio Grande Cooter (Pseudemys gorzugi), a freshwater turtle of increasing conservation concern, residing in these two rivers and their tributaries. The genetic data revealed two distinct populations: one in the Pecos and Black Rivers of New Mexico and another in the Rio Grande and Devils River of Texas, with admixed individuals identified at the confluence of the Rio Grande and Pecos River. In addition to having a smaller geographic range, we found lower observed heterozygosity, reduced nucleotide diversity, and a smaller effective population size (Ne) in New Mexico population. Our results depict a significant isolation-by-distance pattern across their distribution, with migration being notably infrequent at river confluences. These findings are pivotal for future conservation and restoration strategies, emphasizing the need to recognize the unique needs of each population.