Understanding and estimating effective population size for practical application in marine species management

Habitat connectivity drives panda recovery

Globally, the majority of habitat loss is irreversible, and most species will never recover their former ranges. We have learned a great deal about what leads to population decline and extinction, but less about recovery. The recently downlisted giant panda provides a unique opportunity to understand the mechanisms of species recovery. In our study, we estimate giant panda suitable habitats, population density, and gene flow across landscapes to fully investigate the direct and indirect ecological mechanisms underlying bold conservation strategies. We found that the Giant Panda National Survey has modestly but systematically underestimated population size. China's effort to mitigate anthropogenic disturbances was associated with increased panda population density through improving habitat quality and reducing habitat fragmentation. Enhanced landscape connectivity reduced inbreeding via gene flow but indirectly increased inbreeding temporarily due to high local panda density. Although the panda's recovery has been geographically uneven, we provide evidence for improving connectivity and gene flow resulting from conservation efforts. If these processes can be sustained and improved, the panda's path to recovery will be less encumbered by loss of genetic diversity, fostering hope that the present rate of recovery will not be stalled. Findings from this study will not only help guide future giant panda conservation management but also provide a model for how a more mechanistic examination of the genetic processes underlying species recovery can foster the development of more effective strategies for endangered species recovery.

Effective population size in field pea

BACKGROUND

Effective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome.

RESULTS

We calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea.

CONCLUSIONS

Our results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

Estimation of effective number of breeders and effective population size in an abundant and heavily exploited marine teleost

Obtaining reliable estimates of the effective number of breeders (N b) and generational effective population size (N e) for fishery-important species is challenging because they are often iteroparous and highly abundant, which can lead to bias and imprecision. However, recent advances in understanding of these parameters, as well as the development of bias correction methods, have improved the capacity to generate reliable estimates. We utilized samples of both single-cohort young of the year and mixed-age adults from two geographically and genetically isolated stocks of the Australasian snapper (Chrysophrys auratus) to investigate the feasibility of generating reliable N b and N e estimates for a fishery species. Snapper is an abundant, iteroparous broadcast spawning teleost that is heavily exploited by recreational and commercial fisheries. Employing neutral genome-wide SNPs and the linkage-disequilibrium method, we determined that the most reliable N b and N e estimates could be derived by genotyping at least 200 individuals from a single cohort. Although our estimates made from the mixed-age adult samples were generally lower and less precise than those based on a single cohort, they still proved useful for understanding relative differences in genetic effective size between stocks. The correction formulas applied to adjust for biases due to physical linkage of loci and age structure resulted in substantial upward modifications of our estimates, demonstrating the importance of applying these bias corrections. Our findings provide important guidelines for estimating N b and N e for iteroparous species with large populations. This work also highlights the utility of samples originally collected for stock structure and stock assessment work for investigating genetic effective size in fishery-important species.

Temporal genomics help in deciphering neutral and adaptive patterns in the contemporary evolution of kelp populations

The impact of climate change on populations will be contingent upon their contemporary adaptive evolution. In this study, we investigated the contemporary evolution of 4 populations of the cold-water kelp Laminaria digitata by analyzing their spatial and temporal genomic variations using ddRAD-sequencing. These populations were sampled from the center to the southern margin of its north-eastern Atlantic distribution at 2 time points, spanning at least 2 generations. Through genome scans for local adaptation at a single time point, we identified candidate loci that showed clinal variation correlated with changes in sea surface temperature (SST) along latitudinal gradients. This finding suggests that SST may drive the adaptive response of these kelp populations, although factors such as species' demographic history should also be considered. Additionally, we performed a simulation approach to distinguish the effect of selection from genetic drift in allele frequency changes over time. This enabled the detection of loci in the southernmost population that exhibited temporal differentiation beyond what would be expected from genetic drift alone: these are candidate loci which could have evolved under selection over time. In contrast, we did not detect any outlier locus based on temporal differentiation in the population from the North Sea, which also displayed low and decreasing levels of genetic diversity. The diverse evolutionary scenarios observed among populations can be attributed to variations in the prevalence of selection relative to genetic drift across different environments. Therefore, our study highlights the potential of temporal genomics to offer valuable insights into the contemporary evolution of marine foundation species facing climate change.

New insights into population structure, demographic history, and effective population size of the critically endangered blue shark Prionace glauca in the Mediterranean Sea

The blue shark, Prionace glauca, is the most abundant pelagic shark in the open ocean but its vulnerability remains poorly understood while being one of the most fecund sharks. In the Mediterranean Sea, the blue shark is listed as Critically Endangered (CR) by the International Union for Conservation of Nature. The species is facing a strong decline due to fishing, and scientific data regarding its genetic structure and vulnerability are still lacking. Here, we investigated the genetic diversity, demographic history, and population structure of the blue shark within the Mediterranean Sea, from samples of the Gulf of Lion and Malta, using sequences of the mtDNA control region and 22 microsatellite markers. We also compared our mitochondrial data to previous studies to examine the Atlantic-Mediterranean population structure. We assessed the blue shark's genetic vulnerability in the Mediterranean basin by modelling its effective population size. Our results showed a genetic differentiation between the Atlantic and the Mediterranean basins, with limited gene flow between the two areas, and distinct demographic histories making the Mediterranean population an independent management unit. Within the Mediterranean Sea, no sign of population structure was detected, suggesting a single population across the Western and Central parts of the sea. The estimated effective population size was low and highlighted the high vulnerability of the Mediterranean blue shark population, as the estimated size we calculated might not be sufficient to ensure the long-term persistence of the population. Our data also provide additional evidence that the Gulf of Lion area acts as a nursery for P. glauca, where protection is essential for the conservation strategy of the species in the Mediterranean.

Consequences of domestication in eastern oyster: Insights from whole genomic analyses

Selective breeding for production traits has yielded relatively rapid successes with high-fecundity aquaculture species. Discovering the genetic changes associated with selection is an important goal for understanding adaptation and can also facilitate better predictions about the likely fitness of selected strains if they escape aquaculture farms. Here, we hypothesize domestication as a genetic change induced by inadvertent selection in culture. Our premise is that standardized culture protocols generate parallel domestication effects across independent strains. Using eastern oyster as a model and a newly developed 600K SNP array, this study tested for parallel domestication effects in multiple independent selection lines compared with their progenitor wild populations. A single contrast was made between pooled selected strains (1-17 generations in culture) and all wild progenitor samples combined. Population structure analysis indicated rank order levels of differentiation as [wild - wild] < [wild - cultured] < [cultured - cultured]. A genome scan for parallel adaptation to the captive environment applied two methodologically distinct outlier tests to the wild versus selected strain contrast and identified a total of 1174 candidate SNPs. Contrasting wild versus selected strains revealed the early evolutionary consequences of domestication in terms of genomic differentiation, standing genetic diversity, effective population size, relatedness, runs of homozygosity profiles, and genome-wide linkage disequilibrium patterns. Random Forest was used to identify 37 outlier SNPs that had the greatest discriminatory power between bulked wild and selected oysters. The outlier SNPs were in genes enriched for cytoskeletal functions, hinting at possible traits under inadvertent selection during larval culture or pediveliger setting at high density. This study documents rapid genomic changes stemming from hatchery-based cultivation of eastern oysters, identifies candidate loci responding to domestication in parallel among independent aquaculture strains, and provides potentially useful genomic resources for monitoring interbreeding between farm and wild oysters.

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.

Determinants of rapid adaptation in species with large variance in offspring production

The speed of population adaptation to changing biotic and abiotic environments is determined by the interaction between genetic drift, positive selection and linkage effects. Many marine species (fish, crustaceans), invertebrates and pathogens of humans and crops, exhibit sweepstakes reproduction characterized by the production of a very large amount of offspring (fecundity phase) from which only a small fraction may survive to the next generation (viability phase). Using stochastic simulations, we investigate whether the occurrence of sweepstakes reproduction affects the efficiency of a positively selected unlinked locus, and thus, the speed of adaptation since fecundity and/or viability have distinguishable consequences on mutation rate, probability and fixation time of advantageous alleles. We observe that the mean number of mutations at the next generation is always the function of the population size, but the variance increases with stronger sweepstakes reproduction when mutations occur in the parents. On the one hand, stronger sweepstakes reproduction magnifies the effect of genetic drift thus increasing the probability of fixation of neutral allele and decreasing that of selected alleles. On the other hand, the time to fixation of advantageous (as well as neutral) alleles is shortened by stronger sweepstakes reproduction. Importantly, fecundity and viability selection exhibit different probabilities and times to fixation of advantageous alleles under intermediate and weak sweepstakes reproduction. Finally, alleles under both strong fecundity and viability selection display a synergistic efficiency of selection. We conclude that measuring and modelling accurately fecundity and/or viability selection are crucial to predict the adaptive potential of species with sweepstakes reproduction.

Integration of genomic and ecological methods inform management of an undescribed, yet highly exploited, sardine species

Assessing genetic diversity within species is key for conservation strategies in the context of human-induced biotic changes. This is important in marine systems, where many species remain undescribed while being overfished, and conflicts between resource-users and conservation agencies are common. Combining niche modelling with population genomics can contribute to resolving those conflicts by identifying management units and understanding how past climatic cycles resulted in current patterns of genetic diversity. We addressed these issues on an undescribed but already overexploited species of sardine of the genus Harengula. We find that the species distribution is determined by salinity and depth, with a continuous distribution along the Brazilian mainland and two disconnected oceanic archipelagos. Genomic data indicate that such biogeographic barriers are associated with two divergent intraspecific lineages. Changes in habitat availability during the last glacial cycle led to different demographic histories among stocks. One coastal population experienced a 3.6-fold expansion, whereas an island-associated population contracted 3-fold, relative to the size of the ancestral population. Our results indicate that the island population should be managed separately from the coastal population, and that a Marine Protected Area covering part of the island population distribution can support the viability of this lineage.

High germline mutation rates, but not extreme population outbreaks, influence genetic diversity in a keystone coral predator

Lewontin's paradox, the observation that levels of genetic diversity (π) do not scale linearly with census population size (Nc) variation, is an evolutionary conundrum. The most extreme mismatches between π and Nc are found for highly abundant marine invertebrates. Yet, the influences of new mutations on π relative to extrinsic processes such as Nc fluctuations are unknown. Here, we provide the first germline mutation rate (μ) estimate for a marine invertebrate in corallivorous crown-of-thorns sea stars (Acanthaster cf. solaris). We use high-coverage whole-genome sequencing of 14 parent-offspring trios alongside empirical estimates of Nc in Australia's Great Barrier Reef to jointly examine the determinants of π in populations undergoing extreme Nc fluctuations. The A. cf. solaris mean μ was 9.13 x 10-09 mutations per-site per-generation (95% CI: 6.51 x 10-09 to 1.18 x 10-08), exceeding estimates for other invertebrates and showing greater concordance with vertebrate mutation rates. Lower-than-expected Ne (~70,000-180,000) and low Ne/Nc values (0.0047-0.048) indicated weak influences of population outbreaks on long-term π. Our findings are consistent with elevated μ evolving in response to reduced Ne and generation time length, with important implications for explaining high mutational loads and the determinants of genetic diversity in marine invertebrate taxa.

Challenges in estimating effective population sizes from metagenome-assembled genomes

Effective population size (Ne) plays a critical role in shaping the relative efficiency between natural selection and genetic drift, thereby serving as a cornerstone for understanding microbial ecological dynamics. Direct Ne estimation relies on neutral genetic diversity within closely related genomes, which is, however, often constrained by the culturing difficulties for the vast majority of prokaryotic lineages. Metagenome-assembled genomes (MAGs) offer a high-throughput alternative for genomic data acquisition, yet their accuracy in Ne estimation has not been fully verified. This study examines the Thermococcus genus, comprising 66 isolated strains and 29 MAGs, to evaluate the reliability of MAGs in Ne estimation. Despite the even distribution across the Thermococcus phylogeny and the comparable internal average nucleotide identity (ANI) between isolate populations and MAG populations, our results reveal consistently lower Ne estimates from MAG populations. This trend of underestimation is also observed in various MAG populations across three other bacterial genera. The underrepresentation of genetic variation in MAGs, including loss of allele frequency data and variable genomic segments, likely contributes to the underestimation of Ne. Our findings underscore the necessity for caution when employing MAGs for evolutionary studies, which often depend on high-quality genome assemblies and nucleotide-level diversity.

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.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Environmental drivers of size-based population structure, sexual maturity and fecundity: A study of the invasive blue crab Callinectes sapidus (Rathbun, 1896) in the Mediterranean Sea

The blue crab Callinectes sapidus (Rathbun, 1896) is one of the most invasive species in the Mediterranean Sea. Understand how the populations are maintained and how the environment is driving the populations in the areas invaded is the key to an effective future management. This current study is presenting a monthly long-term monitoring of the blue crabs' population structure, body size, sexual maturity, reproduction periods and fecundity, and their relationships with environmental factors in a saltmarshes system in Italy. During winter, high densities (15 ± 8 ind m-2) of early juveniles (< 2 cm) were observed, and their numbers decreased due the population growth until summer. The size-spectra showed that across different temperature (18-26°C) and salinity (24-40 psu) gradients, the growth period for males is faster than for females. Based on abdominal analysis, sexual maturity was defined at ∼12 cm for males and females but the population was in favor of males that were more than 66% of the time mature throughout the year. Copulations periods were identified between spring and autumn when more than 50% of females were matures, and ovigerous females' migrations were observed in late summer. Our study expand our understanding of how the environment interacts to effect physiological and biological processes of alien species and improve our ability to make predictions of how environmental change the distribution of the alien species in the future. Based on our results, we also discuss which population control strategy would be most effective based on the data available in the literature.

Effect of geographic isolation on genetic variation and population structure of Euphrasia nankotaizanensis, a threatened endemic alpine herb in Taiwan

Euphrasia nankotaizanensis (Orobanchaceae) is a rare alpine herb that is endemic to Taiwan. Only four small populations remain in Xue, Nanhu, and Cilai Mountains of Taiwan. The distribution of alpine herbs is severely threatened by climate change, which influences genetic variation and population structure. In this study, we investigated the effects of the natural isolation of alpine habitats on the genetic diversity and geographic structure of populations of E. nankotaizanensis using chloroplast (cp) and nuclear DNA (nrDNA) markers. We found lower levels of genetic diversity in E. nankotaizanensis than in other alpine plants and little to no genetic variation within populations, which could be mainly attributed to the small population size and genetic drift. Only one nrDNA haplotype was present in each population. The lack of monophyly of the four populations in cpDNA probably resulted from lineage sorting or occasional long-distance seed dispersal. Phylogeographic analysis suggested that Nanhu Mountain was probably a refugium over the glacial maxima, agreeing with the potential refugia in central Taiwan. The STRUCTURE and AMOVA analyses revealed significant genetic differentiation in nrDNA among the mountains, which resulted from geographical isolation among these mountains. Estimates of the effective population size (Ne) and demography reflected lower Ne values and a recent population decline, probably implying a greater extinction risk for E. nankotaizanensis. We observed genetic depletion and considerable genetic differentiation among mountain populations, which should be considered in future conservation efforts for this species. In addition, this study provides important insights into the long-term potential of alpine herbs in Taiwan, which are useful for a better prediction of their responses to future climate change.

Fragmented habitat compensates for the adverse effects of genetic bottleneck

In the face of the human-caused biodiversity crisis, understanding the theoretical basis of conservation efforts of endangered species and populations has become increasingly important. According to population genetics theory, population subdivision helps organisms retain genetic diversity, crucial for adaptation in a changing environment. Habitat topography is thought to be important for generating and maintaining population subdivision, but empirical cases are needed to test this assumption. We studied Saimaa ringed seals, landlocked in a labyrinthine lake and recovering from a drastic bottleneck, with additional samples from three other ringed seal subspecies. Using whole-genome sequences of 145 seals, we analyzed the distribution of variation and genetic relatedness among the individuals in relation to the habitat shape. Despite a severe history of genetic bottlenecks with prevalent homozygosity in Saimaa ringed seals, we found evidence for the population structure mirroring the subregions of the lake. Our genome-wide analyses showed that the subpopulations had retained unique variation and largely complementary patterns of homozygosity, highlighting the significance of habitat connectivity in conservation biology and the power of genomic tools in understanding its impact. The central role of the population substructure in preserving genetic diversity at the metapopulation level was confirmed by simulations. Integration of genetic analyses in conservation decisions gives hope to Saimaa ringed seals and other endangered species in fragmented habitats.

Population structure and adaptive differentiation in the sea cucumber Apostichopus californicus and implications for spatial resource management

A growing body of evidence suggests that spatial population structure can develop in marine species despite large population sizes and high gene flow. Characterizing population structure is important for the effective management of exploited species, as it can be used to identify appropriate scales of management in fishery and aquaculture contexts. The California sea cucumber, Apostichopus californicus, is one such exploited species whose management could benefit from further characterization of population structure. Using restriction site-associated DNA (RAD) sequencing, we developed 2075 single nucleotide polymorphisms (SNPs) to quantify genetic structure over a broad section of the species' range along the North American west coast and within the Salish Sea, a region supporting the Washington State A. californicus fishery and developing aquaculture production of the species. We found evidence for population structure (global fixation index (FST) = 0.0068) with limited dispersal driving two patterns of differentiation: isolation-by-distance and a latitudinal gradient of differentiation. Notably, we found detectable population differences among collection sites within the Salish Sea (pairwise FST = 0.001-0.006). Using FST outlier detection and gene-environment association, we identified 10.2% of total SNPs as putatively adaptive. Environmental variables (e.g., temperature, salinity) from the sea surface were more correlated with genetic variation than those same variables measured near the benthos, suggesting that selection on pelagic larvae may drive adaptive differentiation to a greater degree than selection on adults. Our results were consistent with previous estimates of and patterns in population structure for this species in other extents of the range. Additionally, we found that patterns of neutral and adaptive differentiation co-varied, suggesting that adaptive barriers may limit dispersal. Our study provides guidance to decision-makers regarding the designation of management units for A. californicus and adds to the growing body of literature identifying genetic population differentiation in marine species despite large, nominally connected populations.

Minor Genetic Consequences of a Major Mass Mortality: Short-Term Effects in Pisaster ochraceus

AbstractMass mortality events are increasing globally in frequency and magnitude, largely as a result of human-induced change. The effects of these mass mortality events, in both the long and short term, are of imminent concern because of their ecosystem impacts. Genomic data can be used to reveal some of the population-level changes associated with mass mortality events. Here, we use reduced-representation sequencing to identify potential short-term genetic impacts of a mass mortality event associated with a sea star wasting outbreak. We tested for changes in the population for genetic differentiation, diversity, and effective population size between pre-sea star wasting and post-sea star wasting populations of Pisaster ochraceus-a species that suffered high sea star wasting-associated mortality (75%-100% at 80% of sites). We detected no significant population-based genetic differentiation over the spatial scale sampled; however, the post-sea star wasting population tended toward more differentiation across sites than the pre-sea star wasting population. Genetic estimates of effective population size did not detectably change, consistent with theoretical expectations; however, rare alleles were lost. While we were unable to detect significant population-based genetic differentiation or changes in effective population size over this short time period, the genetic burden of this mass mortality event may be borne by future generations, unless widespread recruitment mitigates the population decline. Prior results from P. ochraceus indicated that natural selection played a role in altering allele frequencies following this mass mortality event. In addition to the role of selection found in a previous study on the genomic impacts of sea star wasting on P. ochraceus, our current study highlights the potential role the stochastic loss of many individuals plays in altering how genetic variation is structured across the landscape. Future genetic monitoring is needed to determine long-term genetic impacts in this long-lived species. Given the increased frequency of mass mortality events, it is important to implement demographic and genetic monitoring strategies that capture baselines and background dynamics to better contextualize species' responses to large perturbations.

Population dynamics of Baltic herring since the Viking Age revealed by ancient DNA and genomics

The world's oceans are currently facing major stressors in the form of overexploitation and anthropogenic climate change. The Baltic Sea was home to the first "industrial" fishery ∼800 y ago targeting the Baltic herring, a species that is still economically and culturally important today. Yet, the early origins of marine industries and the long-term ecological consequences of historical and contemporary fisheries remain debated. Here, we study long-term population dynamics of Baltic herring to evaluate the past impacts of humans on the marine environment. We combine modern whole-genome data with ancient DNA (aDNA) to identify the earliest-known long-distance herring trade in the region, illustrating that extensive fish trade began during the Viking Age. We further resolve population structure within the Baltic and observe demographic independence for four local herring stocks over at least 200 generations. It has been suggested that overfishing at Øresund in the 16th century resulted in a demographic shift from autumn-spawning to spring-spawning herring dominance in the Baltic. We show that while the Øresund fishery had a negative impact on the western Baltic herring stock, the demographic shift to spring-spawning dominance did not occur until the 20th century. Instead, demographic reconstructions reveal population trajectories consistent with expected impacts of environmental change and historical reports on shifting fishing targets over time. This study illustrates the joint impact of climate change and human exploitation on marine species as well as the role historical ecology can play in conservation and management policies.

A centenary tale: population genetic insights into the introduction history of the oriental fire-bellied toad (Bombina orientalis) in Beijing

Background

The successful establishment of a species population following a single introduction of a few individuals to a non-native area has been limited. Nevertheless, the oriental fire-bellied toad (Bombina orientalis) population in Beijing is purportedly descended from a single introduction of about 200 individuals translocated from Yantai, Shandong Province, China, in 1927.

Results

To resolve the introduction process and to understand the genetic consequences since that introduction approximately 90 years ago, we investigated the population’s genetic diversity and structure using 261 toads from Beijing and two native Shandong populations and inferred the species’ introduction history using simulation-based approaches. Analysis of mitochondrial DNA (mtDNA) sequences showed the two haplotypes found in Beijing nested within Yantai haplotypes, thus corroborating the historical record of the translocation source. The mtDNA and 11 nuclear microsatellite markers revealed both considerably lower genetic diversity in Beijing than in the source population and strong genetic differentiation between them. Although the current census population in Beijing may be in the range of a few thousand, the effective population size was estimated at only 20–57. Simulations also suggest that this population may have descended from 40–60 founders.

Conclusions

The Beijing population’s genetic patterns were consistent with the consequences of a severe bottleneck during introduction followed by genetic drift. The introduction trajectory constructed for this B. orientalis population reveals the genetic footprints of a small population sustained in isolation for nearly a century. Our results provide an intriguing example of establishment success from limited founders and may inform ex situ conservation efforts as well as the management of biological invasions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02072-z.

Genetic and particle modelling approaches to assessing population connectivity in a deep sea lobster

The emergence of high resolution population genetic techniques, such as genotyping-by-sequencing (GBS), in combination with recent advances in particle modelling of larval dispersal in marine organisms, can deliver powerful new insights to support fisheries conservation and management. In this study, we used this combination to investigate the population connectivity of a commercial deep sea lobster species, the New Zealand scampi, Metanephrops challengeri, which ranges across a vast area of seafloor around New Zealand. This species has limited dispersal capabilities, including larvae with weak swimming abilities and short pelagic duration, while the reptant juvenile/adult stages of the lifecycle are obligate burrow dwellers with limited home ranges. Ninety-one individuals, collected from five scampi fishery management areas around New Zealand, were genotyped using GBS. Using 983 haplotypic genomic loci, three genetically distinct groups were identified: eastern, southern and western. These groups showed significant genetic differentiation with clear source-sink dynamics. The direction of gene flow inferred from the genomic data largely reflected the hydrodynamic particle modelling of ocean current flow around New Zealand. The modelled dispersal during pelagic larval phase highlights the strong connectivity among eastern sampling locations and explains the low genetic differentiation detected among these sampled areas. Our results highlight the value of using a transdisciplinary approach in the inference of connectivity among populations for informing conservation and fishery management.

Population structure and demographic history of the gastropod Thaisella chocolata (Duclos, 1832) from the Southeast Pacific inferred from mitochondrial DNA analyses

The present-day population structure of a species reflects the combination of oceanographic currents, life-history traits, and historical events. However, little is known about the mechanisms that have shaped the gene lineage distribution of marine species inhabiting the Southeast Pacific. Here, we provide a comprehensive phylogeographical study of a species distributed along the Southeast Pacific coastal region by analyzing the endemic gastropod Thaisella chocolata (Duclos, 1832). Sequencing of mitochondrial cytochrome c oxidase subunit 1 (CO1) and 16S rRNA revealed strikingly high haplotypic nucleotide and genetic diversity but a lack of significant population differentiation within the survey area. In addition, a star-shaped phylogeny and significantly negative Tajima's D and Fu's Fs tests of neutrality suggested historical occurrence of rapid demographic expansion. Mismatch distributions and Bayesian inference analyses also confirmed T. chocolata to have undergone two ancestral demographic expansions. Calculations suggested that these expansions began in the lower and middle Pleistocene epoch, likely due to continental shelf development and climatic conditions. These findings could help establish a genetic baseline for T. chocolata as the first step toward sustainable spatial management of this species, as well as understand this species' response to future climate change.

Efficacy of Endangered Pinto Abalone (Haliotis kamtschatkana) Stock Restoration in the Southern Salish Sea From a Genomic Perspective

Declines in abalone populations throughout the world have led to conservation measures including fishery closures and captive breeding programs aimed at stock restoration. Restoration of endangered pinto abalone (Haliotis kamtschatkana) in the southern Salish Sea (Washington State, USA) began in the mid-2000s, and since 2009, nearly 40,000 hatchery-produced juvenile abalone have been outplanted at 21 restoration sites. We used genotyping by sequencing to evaluate the efficacy of this restoration program from a genomic standpoint. Over 49,000 SNPs and 8,000 haplotypes were evaluated across both wild and hatchery-produced abalone. Compared to wild abalone, hatchery-bred abalone had similar heterozygosity, lower allelic richness and effective population size, and higher relatedness. However, more recently bred hatchery progeny sampled prior to outplanting showed intermediate allelic richness, lower relatedness and genomic divergence, and higher effective population size compared to older hatchery outplants sampled after outplanting. We attribute these differences to genetic drift among older hatchery outplants due to mortality and emigration as well as larger numbers of broodstock families in more recent hatchery operations. This suggests that current hatchery practices that combine higher output with larger numbers of broodstock families may more effectively overcome the effects of genetic drift. Effective population size estimates among wild Salish Sea abalone had no upper limit, indicating that these abalone have sufficient genetic diversity to support the restoration program, but also highlighting the need for large broodstock sizes to overcome potential Ryman-Laikre effects. The greatest challenge for the future of the restoration program will be finding sufficient numbers of wild broodstock in the Salish Sea, where pinto abalone remain extremely rare.

Fine Scale Microevolutionary and Demographic Processes Shaping a Wild Metapopulation Dynamics of the South American Fruit Fly Anastrepha fraterculus

Anastrepha fraterculus (Wiedmann) is an important American pest species. Knowledge of its population dynamics is of particular interest for ecology, evolutionary biology, and management programs. In the present study, phenotypic, genotypic, and spatial data were combined, within the frame of landscape genetics, to uncover the spatial population genetic structure (SGS) and demographic processes of an Argentinian local population from the Yungas ecoregion. Eight simple sequence repeats (SSR) loci and six morphometric traits were analysed considering the hierarchical levels: tree/fruit/individual. Genetic variability estimates were high (H_E = 0.72, R_A = 4.39). Multivariate analyses of phenotypic data showed that in average 52.81% of variance is explained by the tree level, followed by between individuals 28.37%. Spatial analysis of morphological traits revealed a negative autocorrelation in all cases. SGS analysis and isolation by distance based on SSR showed no significant autocorrelation for molecular coancestry. The comparison between phenotypic (P_ST) and molecular (F_ST) differentiation identified positive selection in different fruits for all traits. Bayesian analysis revealed a cryptic structure within the population, with three clusters spatially separated. The results of this study showed a metapopulation dynamics. The genetic background of the components of this metapopulation is expected to change through time due to seasonality, repopulation activities, and high gene flow, with an estimated dispersal ability of at least 10 km. Effective population size (Ne) of the metapopulation was estimated in around 800 flies, and within subpopulations (clusters) Ne was associated with the levels of genetic drift experienced by the founding lineages.

A greater yellowstone ecosystem grizzly bear case study: genetic reassessment for managers

There are five grizzly bear (Ursus arctos horribilis) populations in the lower 48 states of the United States. My goal in this Commentary was to ascertain whether genetic diversity is being lost from the isolated GYE grizzly bear population and to better understand any viability implications. I reviewed the scientific literature, including two key genetic studies that the US Fish and Wildlife Service (USFWS) relied upon for their 2007 and current 2017 GYE grizzly bear genetics policy. I discovered that some studies reveal a loss of heterozygosity in the GYE bear population, both historically and in recent decades. Some had a statistically significant depletion rate. My review took place periodically between 2010 and 2021 and indicates that the genome of the GYE grizzly bear population is too small for long-term adaptation. The paper includes a discussion about evolutionary adaptation which invokes time frames rarely considered by nature conservation planners. I also examined genetic statements in the USFWS’s 2017 GYE grizzly bear delisting regulations and highlighted those that seem incongruent with current scientific thought. If this paper is read by some scientists, land managers, administrators, environmentalists, and others with some genetics background, they will better understand some USFWS decisions and policy statements. This case study illustrates that land management agencies can provide a one-sided treatment of some science when writing regulations about genetics.

Genetic structure and diversity of the Chilean flat oyster Ostrea chilensis (Bivalvia: Ostreidae) along its natural distribution from natural beds subject to different fishing histories

Ostrea chilensis (Küster, 1844), the flat oyster, is native to Chile and New Zealand. In Chile, it occurs in a few natural beds, from the northern part of Chiloé Island (41 ºS) to the Guaitecas Archipelago (45 ºS). This bivalve is slow growing, broods its young, and has very limited dispersal potential. The Ostrea chilensis fishery has been over-exploited for a number of decades such that in some locations oysters no longer exist. The aim of this study was to study the genetic diversity of the Chilean flat oyster along its natural distribution to quantify the possible impact of the dredge fishery on wild populations. The genetic structure and diversity of Ostrea chilensis from six natural beds with different histories of fishing activity were estimated. Based on mitochondrial (Cytb) and nuclear (ITS1) DNA sequence variation, our results provide evidence that genetic diversity is different among populations with recent history of wild dredge fishery efforts. We discuss the possible causes of these results. Ultimately, such new information may be used to develop and apply new management measures to promote the sustainable use of this valuable marine resource.

Captive breeding and the conservation of the threatened houbara bustards

Translocation of captive-bred individuals to reinforce wild populations may be an important conservation approach for some species, but can be detrimental when employed to boost exploited wild populations, particularly where repeated long-term reinforcement aims to compensate for repeated unregulated offtake. We review evidence that captive breeding alters multiple physiological, life-history and temperamental traits through founder effects, genetic drift and unintended adaption to captivity; degrades learnt behaviours; and compromises biogeography, population structure and viability through introgression. We highlight these risks for the globally threatened African houbara Chlamydotis undulata and Asian houbara C. macqueenii, 2 bustard species hunted throughout much of their ranges and now subject to multiple large-scale captive-breeding programmes and translocations. In eastern Morocco, annual releases of captive-bred African houbara are 2‒3 times higher than original wild numbers, but no investigation of their potentially deleterious effects has, to our knowledge, been published, although most wild populations may now have been replaced by captive-bred domestic stock, which are reportedly not self-sustaining. Despite multiple decades of reinforcement, we are not aware of any analysis of the contribution of captive breeding to African houbara population dynamics, or of the genomic consequences. Asian houbara release programmes may also be promoting rather than preventing declines, and need to contextualise themselves through rigorous analyses of wild population numbers, demographic rates and threats, maintenance of phylogeographic concordance of released with supplemented populations, profiling of traits crucial to survival and the measurement and modelling of the impacts of reinforcement on physiological and behavioural fitness of wild populations.

A brief history and popularity of methods and tools used to estimate micro-evolutionary forces

Population genetics is a field of research that predates the current generations of sequencing technology. Those approaches, that were established before massively parallel sequencing methods, have been adapted to these new marker systems (in some cases involving the development of new methods) that allow genome-wide estimates of the four major micro-evolutionary forces-mutation, gene flow, genetic drift, and selection. Nevertheless, classic population genetic markers are still commonly used and a plethora of analysis methods and programs is available for these and high-throughput sequencing (HTS) data. These methods employ various and diverse theoretical and statistical frameworks, to varying degrees of success, to estimate similar evolutionary parameters making it difficult to get a concise overview across the available approaches. Presently, reviews on this topic generally focus on a particular class of methods to estimate one or two evolutionary parameters. Here, we provide a brief history of methods and a comprehensive list of available programs for estimating micro-evolutionary forces. We furthermore analyzed their usage within the research community based on popularity (citation bias) and discuss the implications of this bias for the software community. We found that a few programs received the majority of citations, with program success being independent of both the parameters estimated and the computing platform. The only deviation from a model of exponential growth in the number of citations was found for the presence of a graphical user interface (GUI). Interestingly, no relationship was found for the impact factor of the journals, when the tools were published, suggesting accessibility might be more important than visibility.

Ghosts of a Structured Past: Impacts of Ancestral Patterns of Isolation-by-Distance on Divergence-Time Estimation

Isolation-by-distance is a widespread pattern in nature that describes the reduction of genetic correlation between subpopulations with increased geographic distance. In the population ancestral to modern sister species, this pattern may hypothetically inflate population divergence time estimation due to allele frequency differences in subpopulations at the ends of the ancestral population. In this study, we analyze the relationship between the time to the most recent common ancestor and the population divergence time when the ancestral population model is a linear stepping-stone. Using coalescent simulations, we compare the coalescent time to the population divergence time for various ratios of the divergence time over the population size. Next, we simulate whole genomes to obtain single nucleotide polymorphisms (SNPs), and use the Bayesian coalescent program SNAPP to estimate divergence times. We find that as the rate of migration between neighboring demes decreases, the coalescent time becomes significantly greater than the population divergence time when sampled from end demes. Divergence-time overestimation in SNAPP becomes severe when the divergence-to-population size ratio < 10 and migration is low. Finally, we demonstrate the impact of ancestral isolation-by-distance on divergence-time estimation using an empirical dataset of squamates (Tropidurus) endemic to Brazil. We conclude that studies estimating divergence times should be cognizant of the potential ancestral population structure in an explicitly spatial context or risk dramatically overestimating the timing of population splits.

Genomic signatures of spatially divergent selection at clownfish range margins

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

The paradox of retained genetic diversity of Hippocampus guttulatus in the face of demographic decline

Genetic diversity is the raw foundation for evolutionary potential. When genetic diversity is significantly reduced, the risk of extinction is heightened considerably. The long-snouted seahorse (Hippocampus guttulatus) is one of two seahorse species occurring in the North-East Atlantic. The population living in the Ria Formosa (South Portugal) declined dramatically between 2001 and 2008, prompting fears of greatly reduced genetic diversity and reduced effective population size, hallmarks of a genetic bottleneck. This study tests these hypotheses using samples from eight microsatellite loci taken from 2001 and 2013, on either side of the 2008 decline. The data suggest that the population has not lost its genetic diversity, and a genetic bottleneck was not detectable. However, overall relatedness increased between 2001 to 2013, leading to questions of future inbreeding. The effective population size has seemingly increased close to the threshold necessary for the population to retain its evolutionary potential, but whether these results have been affected by sample size is not clear. Several explanations are discussed for these unexpected results, such as gene flow, local decline due to dispersal to other areas of the Ria Formosa, and the potential that the duration of the demographic decline too short to record changes in the genetic diversity. Given the results presented here and recent evidence of a second population decline, the precise estimation of both gene flow and effective population size via more extensive genetic screening will be critical to effective population management.

Population Viability of Sea Turtles in the Context of Global Warming

Sea turtles present a model for the potential impacts of climate change on imperiled species, with projected warming generating concern about their persistence. Various sea turtle life-history traits are affected by temperature; most strikingly, warmer egg incubation temperatures cause female-biased sex ratios and higher embryo mortality. Predictions of sea turtle resilience to climate change are often focused on how resulting male limitation or reduced offspring production may affect populations. In the present article, by reviewing research on sea turtles, we provide an overview of how temperature impacts on incubating eggs may cascade through life history to ultimately affect population viability. We explore how sex-specific patterns in survival and breeding periodicity determine the differences among offspring, adult, and operational sex ratios. We then discuss the implications of skewed sex ratios for male-limited reproduction, consider the negative correlation between sex ratio skew and genetic diversity, and examine consequences for adaptive potential. Our synthesis underscores the importance of considering the effects of climate throughout the life history of any species. Lethal effects (e.g., embryo mortality) are relatively direct impacts, but sublethal effects at immature life-history stages may not alter population growth rates until cohorts reach reproductive maturity. This leaves a lag during which some species transition through several stages subject to distinct biological circumstances and climate impacts. These perspectives will help managers conceptualize the drivers of emergent population dynamics and identify existing knowledge gaps under different scenarios of predicted environmental change.

Anthropogenic river fragmentation reduces long-term viability of the migratory fish Salminus brasiliensis (Characiformes: Bryconidae) populations

Abstract Life-history, geographical barriers, and damming can shape the genetic diversity of freshwater migratory fish, which are particularly vulnerable to anthropogenic impacts. We investigated the genetic diversity of Salminus brasiliensis, a long-distance migratory species that is recognized as an important provider of ecosystem services. We implemented microsatellite analyses to assess genetic diversity and simulate future scenarios for evaluating the long-term viability of dammed and non-dammed populations from the Uruguay River. High levels of genetic diversity were detected for all sampled populations. However, effective population sizes were lower in the uppermost river stretches, where the landscape is highly fragmented. Population structure analysis indicated two spatial genetic populations. It is suggested that this genetic structure preserves populations partially isolated by an ancient natural barrier, instead of being a result of the presence of dams. The simulated genetic scenarios indicated that genetic variability of S. brasiliensis populations from upstream dams could collapse over the years, mainly due to the reduction in the number of alleles. Therefore, besides helping to better understand issues related to the influence of dams on the genetic diversity of migratory fish, our results are especially relevant for driving local fishery policies and management actions for the species conservation.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

Combining genetic and demographic monitoring better informs conservation of an endangered urban snake

Conversion and fragmentation of wildlife habitat often leads to smaller and isolated populations and can reduce a species' ability to disperse across the landscape. As a consequence, genetic drift can quickly lower genetic variation and increase vulnerability to extirpation. For species of conservation concern, quantification of population size and connectivity can clarify the influence of genetic drift in local populations and provides important information for conservation management and recovery strategies. Here, we used genome-wide single nucleotide polymorphism (SNP) data and capture-mark-recapture methods to evaluate the genetic diversity and demography within seven focal sites of the endangered San Francisco gartersnake (Thamnophis sirtalis tetrataenia), a species affected by alteration and isolation of wetland habitats throughout its distribution. The primary goals were to determine the population structure and degree of genetic isolation among T. s. tetrataenia populations and estimate effective size and population abundance within sites to better understand the present and future importance of genetic drift. We also used temporally sampled datasets to examine the magnitude of genetic change over time. We found moderate population genetic structure throughout the San Francisco Peninsula that partitions sites into northern and southern regional clusters. Point estimates of both effective size and population abundance were generally small (≤ 100) for a majority of the sites, and estimates were particularly low in the northern populations. Genetic analyses of temporal datasets indicated an increase in genetic differentiation, especially for the most geographically isolated sites, and decreased genetic diversity over time in at least one site (Pacifica). Our results suggest that drift-mediated processes as a function of small population size and reduced connectivity from neighboring populations may decrease diversity and increase differentiation. Improving genetic diversity and connectivity among T. s. tetrataenia populations could promote persistence of this endangered snake.

Time matters: genetic composition and evaluation of effective population size in temperate coastal fish species

Background

Extensive knowledge on the genetic characterization of marine organisms has been assembled, mainly concerning the spatial distribution and structuring of populations. Temporal monitoring assesses not only the stability in genetic composition but also its trajectory over time, providing critical information for the accurate forecast of changes in genetic diversity of marine populations, particularly important for both fisheries and endangered species management. We assessed fluctuations in genetic composition among different sampling periods in the western Portuguese shore in three fish species.

Methods

White seabream Diplodus sargus, sand smelt Atherina presbyter and shanny Lipophrys pholis were chosen, because of their genetic patterns in distinct ecological environments, insight into historical and contemporary factors influencing population effective size (N_e), and degree of commercial exploitation. Samples were obtained near Lisbon between 2003 and 2014 and screened for genetic variation with mitochondrial and nuclear markers. Analyses included genealogies, genetic diversities, temporal structures and contemporary N_e.

Results

For mtDNA no temporal structure was detected, while for nDNA significant differences were recorded between some sampling periods for the shanny and the sand smelt. Haplotype networks revealed deep genealogies, with various levels of diversification. The shanny revealed a smaller N_e/generation when compared to the other species, which, in turn, revealed no evidence of genetic drift for most study periods. These results highlight the fact that temporal variations in genetic pool composition should be considered when evaluating the population structure of fish species with long distance dispersal, which are more vulnerable to recruitment fluctuations.

Contrasting demographic trends and asymmetric migration rates in a spatially structured amphibian population

Natural populations often persist at the landscape scale as metapopulations, with breeding units (subpopulations) experiencing temporal extinction and recolonization events. Important parameters to forecast population viability in these systems include the ratio of the effective number of breeders (N_b ) to the total number of adults (N_a ) and migration rates among subpopulations. Here, we present the results of a 10-year integrative monitoring program of a metapopulation of the Iberian green frog (Pelophylax perezi) in central Spain. We characterized population dynamics at two main breeding ponds (Gravera and Laguna) using capture-mark-recapture data to estimate N_a in each breeding season, and multilocus genotypes to estimate the effective population size (N_e ), N_b , individual breeding success, and migration rates. Both ponds experienced population decline after a dry season, with Gravera subsequently recovering and Laguna suffering a bottleneck associated with genetic impoverishment. In this subpopulation, average allelic richness and private alleles decreased from 2010 (10.87 and 1.67, respectively) to 2018 (8.0 and 0.20). The N_b /N_a ratio in Laguna in 2018 was twice as high (0.95) than in Gravera (0.41) or in pre-bottleneck Laguna (0.50), suggesting plasticity or genetic compensation through increased individual breeding success. Migration rates were asymmetric between ponds, with a stronger contribution from Gravera to Laguna (29.9% vs. 16.2% in the opposite direction) that may result in a rescue effect. This study emphasizes the importance of integrative demographic approaches for the monitoring of natural populations based on a better understanding of their spatio-temporal dynamics, which provides valuable information for conservation actions.

Estimating effective population size using RADseq: Effects of SNP selection and sample size

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.