Estimating the Effective Size of European Wolf Populations

Molecular methods are routinely used to estimate the effective size of populations (N e). However, underlying model assumptions are frequently violated to an unknown extent. Although simulations can detect sources of bias and help to adjust sampling strategies and analyses methods, additional information from empirical data can also be used to calibrate methods and improve molecular N e estimation methods. Here, we take advantage of long-term genetic and ecological monitoring data of the grey wolf (Canis lupus) in Germany, and detailed population genetic studies in Poland, Spain and Portugal to improve N e estimation strategies in this species, and species with similar life history traits. We first calculated N e from average lifetime reproductive success and detailed census data from the German population, which served as a baseline to compare to molecular estimates based on linkage disequilibrium and sibship frequency. This yielded a robust N e/N c estimation that we used to calibrate molecular estimates of German, Polish and Iberian wolf populations. The linkage disequilibrium method was strongly influenced by spatial genetic structure, much more than the sibship frequency method. When N e was estimated in local neighbourhoods, both methods yielded comparable results. Estimates of the metapopulation effective size seemed to correspond generally well with the sum of the estimates of local neighbourhoods. Overall, we found that the number of packs is a good proxy of the effective population size. Using this as a rule of thumb, we evaluated for all European wolf populations the N e 500 indicator and concluded that half of the European wolf populations do not yet fulfil this criterion.

Unraveling the Complexity of the N e/N c Ratio for Conservation of Large and Widespread Pelagic Fish Species: Current Status and Challenges

Estimating and understanding the ratio between effective population size (N e) and census population size (N c) are pivotal in the conservation of large marine pelagic fish species, including bony fish such as tunas and cartilaginous fish such as sharks, given the challenges associated with obtaining accurate estimates of their abundance. The difficulties inherent in capturing and monitoring these species in vast and dynamic marine environments often make direct estimation of their population size challenging. By focusing on N e, it is conceivable in certain cases to approximate census size once the N e/N c ratio is known, although this ratio can vary and does not always increase linearly, as it is influenced by various ecological and evolutionary factors. Thus, this ratio presents challenges and complexities in the context of pelagic species conservation. To delve deeper into these challenges, firstly, we recall the diverse types of effective population sizes, including contemporary and historical sizes, and their implications in conservation biology. Secondly, we outline current knowledge about the influence of life history traits on the N e/N c ratio in the light of examples drawn from large and abundant pelagic fish species. Despite efforts to document an increasing number of marine species using recent technologies and statistical methods, establishing general rules to predict N e/N c remains elusive, necessitating further research and investment. Finally, we recall statistical challenges in relating N e and N c emphasizing the necessity of aligning temporal and spatial scales. This last part discusses the roles of generation and reproductive cycle effective population sizes to predict genetic erosion and guiding management strategies. Collectively, these sections underscore the multifaceted nature of effective population size estimation, crucial for preserving genetic diversity and ensuring the long-term viability of populations. By navigating statistical and theoretical complexities, and addressing methodological challenges, scientists should be able to advance our understanding of the N e/N c ratio.

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

Dealing With Assumptions and Sampling Bias in the Estimation of Effective Population Size: A Case Study in an Amphibian Population

Accurately estimating effective population size (N e) is essential for understanding evolutionary processes and guiding conservation efforts. This study investigates N e estimation methods in spatially structured populations using a population of moor frog (Rana arvalis) as a case study. We assessed the behaviour of N e estimates derived from the linkage disequilibrium (LD) method as we changed the spatial configuration of samples. Moor frog eggs were sampled from 25 breeding patches (i.e., separate vernal ponds, ditches or parts of larger fens) within a single population, revealing an isolation-by-distance pattern and a local spatial genetic structure. Varying buffer sizes around each patch were used to examine the impact of sampling window size on the estimation of effective number of breeders (N b). Our results indicate a downward bias in LD N b estimates with increasing buffer size, suggesting an underestimation of N b. The observed bias is attributed to LD resulting from including genetically divergent individuals (mixture-LD) confounding LD due to drift. This emphasises the significance of considering even subtle spatial genetic patterns. The implications of these findings are discussed, emphasising the need to account for spatial genetic structure to accurately assess population viability and inform conservation efforts. This study contributes to our understanding of the challenges associated with N e estimation in spatially structured populations and underscores the importance of refining methodologies to address population-specific spatial dynamics for effective conservation planning and management.

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.

Global assessment of effective population sizes: Consistent taxonomic differences in meeting the 50/500 rule

Effective population size (Ne) is a particularly useful metric for conservation as it affects genetic drift, inbreeding and adaptive potential within populations. Current guidelines recommend a minimum Ne of 50 and 500 to avoid short-term inbreeding and to preserve long-term adaptive potential respectively. However, the extent to which wild populations reach these thresholds globally has not been investigated, nor has the relationship between Ne and human activities. Through a quantitative review, we generated a dataset with 4610 georeferenced Ne estimates from 3829 populations, extracted from 723 articles. These data show that certain taxonomic groups are less likely to meet 50/500 thresholds and are disproportionately impacted by human activities; plant, mammal and amphibian populations had a <54% probability of reachingN ̂ e = 50 and a <9% probability of reachingN ̂ e = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller medianN ̂ e than unlisted populations, and this was consistent across all taxonomic groups.N ̂ e was reduced in areas with a greater Global Human Footprint, especially for amphibians, birds and mammals, however relationships varied between taxa. We also highlight several considerations for future works, including the role that gene flow and subpopulation structure plays in the estimation ofN ̂ e in wild populations, and the need for finer-scale taxonomic analyses. Our findings provide guidance for more specific thresholds based on Ne and help prioritise assessment of populations from taxa most at risk of failing to meet conservation thresholds.

Low-coverage sequencing and Wahlund effect severely bias estimates of inbreeding, heterozygosity and effective population size in North American wolves

vonHoldt et al. ((2024), Molecular Ecology, 33, e17231) (vH24) used low-coverage (average ~ 7X read depth) restriction site-associated DNA sequence data to estimate individual inbreeding and heterozygosity, and recent effective population size (Ne), in Great Lakes (GL) and Northern Rocky Mountain (RM) wolves. They concluded that RM heterozygosity rapidly declined between 1991 and 2020, and that Ne declined substantially in GL and RM over the last 50 generations. Here, we evaluate the effects of low sequence coverage and sampling strategy on vH24's findings and provide general recommendations for using sequence data to evaluate inbreeding, heterozygosity and Ne. Low-coverage sequencing resulted in downwardly biased estimates of individual inbreeding and heterozygosity, and an erroneous large temporal decline in RM heterozygosity due to declining read depth through time. Additionally, vH24's sampling strategy-which combined individuals from several genetically differentiated populations and across at least eight wolf generations-is expected to have resulted in severe downward bias in estimates of recent Ne for RM. We recommend using high-coverage sequence data (≥ $$ \ge $$ 15-20X) to estimate inbreeding and heterozygosity. Carefully filtering individuals, loci and genotypes, and using genotype imputation or likelihoods can help to minimise bias when low-coverage sequence data must be used. For estimation of contemporary Ne, the marginal benefits of using more than 103-104 loci are small, so aggressive filtering of loci with low average read depth potentially can retain most individuals without sacrificing much precision. Individuals are relatively more valuable than loci because analyses of contemporary Ne should focus on roughly single-generation samples from local breeding populations.

Planning and implementing genetic rescue of an endangered freshwater fish population in a regulated river, where low flow reduces breeding opportunities and may trigger inbreeding depression

Augmenting depleted genetic diversity can improve the fitness and evolutionary potential of wildlife populations, but developing effective management approaches requires genetically monitored test cases. One such case is the small, isolated and inbred Cotter River population of an endangered Australian freshwater fish, the Macquarie perch Macquaria australasica, which over 3 years (2017-2019) received 71 translocated migrants from a closely related, genetically more diverse population. We used genetic monitoring to test whether immigrants bred, interbred with local fish and augmented population genetic diversity. We also investigated whether levels of river flow affected recruitment, inbreeding depression and juvenile dispersal. Fish length was used to estimate the age, birth year cohort and growth of 524 individuals born between 2016 and 2020 under variable flow conditions. DArT genome-wide genotypes were used to assess individual ancestry, heterozygosity, short-term effective population size and identify parent-offspring and full-sibling families. Of 442 individuals born after translocations commenced, only two (0.45%) were of mixed ancestry; these were half-sibs with one translocated parent in common. Numbers of breeders and genetic diversity for five birth year cohorts of the Cotter River fish were low, especially in low-flow years. Additionally, individuals born in the year of lowest flow evidently suffered from inbreeding depression for juvenile growth. The year of highest flow was associated with the largest number of breeders, lowest inbreeding in the offspring and greatest juvenile dispersal distances. Genetic diversity decreased in the upstream direction, flagging restricted access of breeders to the most upstream breeding sites, exacerbated by low river flow. Our results suggest that the effectiveness of translocations could be increased by focussing on upstream sites and moving more individuals per year; using riverine sources should be considered. Our results indicate that river flow sufficient to facilitate fish movement through the system would increase the number of breeders, promote individuals' growth, reduce inbreeding depression and promote genetic rescue.

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

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

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

Population genetic insights into establishment, adaptation, and dispersal of the invasive quagga mussel across perialpine lakes

Human activities have facilitated the invasion of freshwater ecosystems by various organisms. Especially, invasive bivalves such as the quagga mussels, Dreissena bugensis, have the potential to alter ecosystem function as they heavily affect the food web. Quagga mussels occur in high abundance, have a high filtration rate, quickly spread within and between waterbodies via pelagic larvae, and colonize various substrates. They have invaded various waterbodies across the Northern Hemisphere. In Central Europe, they have invaded multiple large and deep perialpine lakes with first recordings in Lake Geneva in 2015 and 2016 in Lake Constance. In the deep perialpine lakes, quagga mussels quickly colonized the littoral zone but are also abundant deeper (>80 m), where they are often thinner and brighter shelled. We analysed 675 quagga mussels using ddRAD sequencing to gain in-depth insights into the genetic population structure of quagga mussels across Central European lakes and across various sites and depth habitats in Lake Constance. We revealed substantial genetic differentiation amongst quagga mussel populations from three unconnected lakes, and all populations showed high genetic diversity and effective population size. In Lake Constance, we detected no genetic differentiation amongst quagga mussels sampled across different sites and depth habitats. We also did not identify any convincing candidate loci evidential for adaptation along a depth gradient and a transplant experiment showed no indications of local adaptation to living in the deep based on investigating growth and survival. Hence, the shallow-water and the deep-water morphotypes seem to be a result of phenotypic plasticity rather than local adaptation to depth. In conclusion, our ddRAD approach revealed insight into the establishment of genetically distinct quagga mussel populations in three perialpine lakes and suggests that phenotypic plasticity and life history traits (broadcast spawner with high fecundity and dispersing pelagic larvae) facilitate the fast spread and colonization of various depth habitats by the quagga mussel.

Practical application of the linkage disequilibrium method for estimating contemporary effective population size: A review

The method to estimate contemporary effective population size (Ne ) based on patterns of linkage disequilibrium (LD) at unlinked loci has been widely applied to natural and managed populations. The underlying model makes many simplifying assumptions, most of which have been evaluated in numerous studies published over the last two decades. Here, these performance evaluations are reviewed and summarized, with a focus on information that facilitates practical application to real populations in nature. Potential sources of bias that are discussed include calculation of r2 (a measure of LD), adjustments for sampling error, physical linkage, age structure, migration and spatial structure, mutation and selection, mating systems, changes in abundance, rare alleles, missing data, genotyping errors, data filtering choices and methods for combining multiple Ne estimates. Factors that affect precision are reviewed, including pseudoreplication that limits the information gained from large genomics datasets, constraints imposed by small samples of individuals, and the challenges in obtaining robust estimates for large populations. Topics that merit further research include the potential to weight r2 values by allele frequency, lump samples of individuals, use genotypic likelihoods rather than called genotypes, prune large LD values and apply the method to species practising partial monogamy.

Managing genomic diversity in conservation programs of Chinese domestic chickens

BACKGROUND

Effective conservation and utilization of farm animals are fundamental for realizing sustainable increases in food production. In situ and ex situ conservation are the two main strategies that are currently used to protect the genetic integrity of Chinese domestic chicken breeds. However, genomic diversity and population structure have not been compared in these conserved populations.

RESULTS

Three hundred and sixty-one individuals from three Chinese domestic chicken breeds were collected from populations conserved in situ and ex situ and genotyped using genotyping-by-sequencing (GBS). First, we used different parameters based on heterozygosity, genomic inbreeding, and linkage disequilibrium to estimate the genomic diversity of these populations, and applied principal component analysis (PCA), neighbor-joining tree, and ADMIXTURE to analyze population structure. We found that the small ex situ conserved populations, which have been maintained in controlled environments, retained less genetic diversity than the in situ conserved populations. In addition, genetic differentiation was detected between the in situ and ex situ conserved populations of the same breed. Next, we analyzed signatures of selection using three statistical methods (fixation index (FST), nucleotide diversity (Pi), and cross-population extended haplotype homozygosity (XP-EHH) to study the genetic footprints that underlie the differentiation between in situ and ex situ conserved populations. We concluded that, in these small populations, differentiation might be caused by genetic drift or by mutations from the original populations. The differentiation observed in the population of Beijing You chicken probably reflects adaptation to environmental changes in temperature and humidity that the animals faced when they were moved from their place of origin to the new site for ex situ conservation.

CONCLUSIONS

Conservation programs of three Chinese domestic chicken breeds have maintained their genomic diversity to a sustainable degree. The small ex situ conserved populations, which are maintained in controlled environments, retain less genetic diversity than populations conserved in situ. In addition, the transfer of populations from their place of origin to another site for conservation purposes results in genetic differentiation, which may be caused by genetic drift or adaptation. This study provides a basis for further optimization of in situ and ex situ conservation programs for domestic chicken breeds in China.

Impact of population structure in the estimation of recent historical effective population size by the software GONE

BACKGROUND

Effective population size (Ne) is a crucial parameter in conservation genetics and animal breeding. A recent method, implemented by the software GONE, has been shown to be rather accurate in estimating recent historical changes in Ne from a single sample of individuals. However, GONE estimations assume that the population being studied has remained isolated for a period of time, that is, without migration or confluence of other populations. If this occurs, the estimates of Ne can be heavily biased. In this paper, we evaluate the impact of migration and admixture on the estimates of historical Ne provided by GONE through a series of computer simulations considering several scenarios: (a) the mixture of two or more ancestral populations; (b) subpopulations that continuously exchange individuals through migration; (c) populations receiving migrants from a large source; and (d) populations with balanced systems of chromosomal inversions, which also generate genetic structure.

RESULTS

Our results indicate that the estimates of historical Ne provided by GONE may be substantially biased when there has been a recent mixture of populations that were previously separated for a long period of time. Similarly, biases may occur when the rate of continued migration between populations is low, or when chromosomal inversions are present at high frequencies. However, some biases due to population structuring can be eliminated by conducting population structure analyses and restricting the estimation to the differentiated groups. In addition, disregarding the genomic regions that are involved in inversions can also remove biases in the estimates of Ne.

CONCLUSIONS

Different kinds of deviations from isolation and panmixia of the populations can generate biases in the recent historical estimates of Ne. Therefore, estimation of past demography could benefit from performing population structure analyses beforehand, by mitigating the impact of these biases on historical Ne estimates.

Reproductive strategies and their consequences for divergence, gene flow, and genetic diversity in three taxa of Clarkia

Differences in reproductive strategies can have important implications for macro- and micro-evolutionary processes. We used a comparative approach through a population genetics lens to evaluate how three distinct reproductive strategies shape patterns of divergence among as well as gene flow and genetic diversity within three closely related taxa in the genus Clarkia. One taxon is a predominantly autonomous self-fertilizer and the other two taxa are predominantly outcrossing but vary in the primary pollinator they attract. In genotyping populations using genotyping-by-sequencing and comparing loci shared across taxa, our results suggest that differences in reproductive strategies in part promote evolutionary divergence among these closely related taxa. Contrary to expectations, we found that the selfing taxon had the highest levels of heterozygosity but a low rate of polymorphism. The high levels of fixed heterozygosity for a subset of loci suggests this pattern is driven by the presence of structural rearrangements in chromosomes common in other Clarkia taxa. In evaluating patterns within taxa, we found a complex interplay between reproductive strategy and geographic distribution. Differences in the mobility of primary pollinators did not translate to a difference in rates of genetic diversity and gene flow within taxa - a pattern likely due to one taxon having a patchier distribution and a less temporally and spatially reliable pollinator. Taken together, this work advances our understanding of the factors that shape gene flow and the distribution of genetic diversity within and among closely related taxa.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.

Landscape genetics reveals contrasting patterns of connectivity in two newt species (Lissotriton montandoni and L. vulgaris)

Ecologically distinct species may respond to landscape changes in different ways. In addition to basic ecological data, the extent of the geographic range has been successfully used as an indicator of species sensitivity to anthropogenic landscapes, with widespread species usually found to be less sensitive compared to range-restricted species. In this study, we investigate connectivity patterns of two closely related but ecologically distinct newt species - the range-restricted, Lissotriton montandoni and the widespread, L. vulgaris - using genomic data, a highly replicated setting (six geographic regions per species), and tools from landscape genetics. Our results show the importance of forest for connectivity in both species, but at the same time suggest differential use of forested habitat, with L. montandoni and L. vulgaris showing the highest connectivity at forest-core and forest-edges, respectively. Anthropogenic landscapes (i.e., higher crop- or urban-cover) increased resistance in both species, but the effect was one to three orders of magnitude stronger in L. montandoni than in L. vulgaris. This result is consistent with a view of L. vulgaris as an ecological generalist. Even so, currently, the negative impact of anthropogenic landscapes is mainly seen in connectivity among L. vulgaris populations, which show significantly stronger isolation and lower effective sizes relative to L. montandoni. Overall, this study emphasizes how habitat destruction is compromising genetic connectivity not only in endemic, range-restricted species of conservation concern but also in widespread generalist species, despite their comparatively lower sensitivity to anthropogenic landscape changes.

Assessment of the Global Variance Effective Size of Subdivided Populations, and Its Relation to Other Effective Sizes

The variance effective population size ([Formula: see text]) is frequently used to quantify the expected rate at which a population's allele frequencies change over time. The purpose of this paper is to find expressions for the global [Formula: see text] of a spatially structured population that are of interest for conservation of species. Since [Formula: see text] depends on allele frequency change, we start by dividing the cause of allele frequency change into genetic drift within subpopulations (I) and a second component mainly due to migration between subpopulations (II). We investigate in detail how these two components depend on the way in which subpopulations are weighted as well as their dependence on parameters of the model such a migration rates, and local effective and census sizes. It is shown that under certain conditions the impact of II is eliminated, and [Formula: see text] of the metapopulation is maximized, when subpopulations are weighted proportionally to their long term reproductive contributions. This maximal [Formula: see text] is the sought for global effective size, since it approximates the gene diversity effective size [Formula: see text], a quantifier of the rate of loss of genetic diversity that is relevant for conservation of species and populations. We also propose two novel versions of [Formula: see text], one of which (the backward version of [Formula: see text]) is most stable, exists for most populations, and is closer to [Formula: see text] than the classical notion of [Formula: see text]. Expressions for the optimal length of the time interval for measuring genetic change are developed, that make it possible to estimate any version of [Formula: see text] with maximal accuracy.

Variance effective population size is affected by census size in sub-structured populations

Measurement of allele frequency shifts between temporally spaced samples has long been used for assessment of effective population size (N_e ), and this 'temporal method' provides estimates of N_e referred to as variance effective size (N_eV ). We show that N_eV of a local population that belongs to a sub-structured population (a metapopulation) is determined not only by genetic drift and migration rate (m), but also by the census size (N_c ). The realized N_eV of a local population can either increase or decrease with increasing m, depending on the relationship between N_e and N_c in isolation. This is shown by explicit mathematical expressions for the factors affecting N_eV derived for an island model of migration. We verify analytical results using high-resolution computer simulations, and show that the phenomenon is not restricted to the island model migration pattern. The effect of N_c on the realized N_eV of a local subpopulation is most pronounced at high migration rates. We show that N_c only affects local N_eV , whereas N_eV for the metapopulation as a whole, inbreeding (N_eI ), and linkage disequilibrium (N_eLD ) effective size are all independent of N_c . Our results provide a possible explanation to the large variation of N_e /N_c ratios reported in the literature, where N_e is frequently estimated by N_eV . They are also important for the interpretation of empirical N_e estimates in genetic management where local N_eV is often used as a substitute for inbreeding effective size, and we suggest an increased focus on metapopulation N_eV as a proxy for N_eI .

Drivers of strong isolation and small effective population size at a leading range edge of a widespread plant

Climate change has influenced species distributions worldwide with upward elevational shifts observed in many systems. Leading range edge populations, like those at upper elevation limits, are crucial for climate change responses but can exhibit low genetic diversity due to founder effects, isolation, or limited outbreeding. These factors can hamper local adaptation at range limits. Using the widespread herb, Argentina anserina, we measured ecological attributes (population density on the landscape, area of population occupancy, and plant and flower density) spanning a 1000 m elevation gradient, with high elevation populations at the range limit. We measured vegetative clonal potential in the greenhouse for populations spanning the gradient. We combined these data with a ddRAD-seq dataset to test the hypotheses that high elevation populations would exhibit ecological and genomic signatures of leading range edge populations. We found that population density on the landscape declined towards the high elevation limit, as is expected towards range edges. However, plant density was elevated within edge populations. In the greenhouse, high elevation plants exhibited stronger clonal potential than low elevation plants, likely explaining increased plant density in the field. Phylogeographic analysis supported more recent colonization of high elevation populations which were also more genetically isolated, had more extreme heterozygote excess and had smaller effective population size than low. Results support that colonization of high elevations was likely accompanied by increased asexuality, contributing to a decline in effective population size. Despite high plant density in leading edge populations, their small effective size, isolation and clonality could constrain adaptive potential.

Eco-evolutionary factors that influence its demographic oscillations in Prochilodus costatus (Actinopterygii: Characiformes) populations evidenced through a genetic spatial-temporal evaluation

The human activity impact on wild animal populations is indicated by eco-evolutionary and demographic processes, along with their survival and capacity to evolve; consequently, such data can contribute toward enhancing genetic-based conservation programs. In this context, knowledge on the life-history and the eco-evolutionary processes is required to understand extant patterns of population structure in Prochilodus costatus a Neotropical migratory fish that has been threatened due to loss and fragmentation of its natural habitat since 1960s promoted by the expansion of hydroelectric power plant construction programs. This study evaluated the eco-evolutionary parameters that cause oscillations in the demography and structure of P. costatus populations. An integrated approach was used, including temporal and spatial sampling, next-generation sequencing of eight microsatellite loci, multivariate genetic analysis, and demographic life-history reconstruction. The results provided evidence of the complex interplay of ecological-evolutionary and human-interference events on the life history of this species in the upper basin. In particular, spawning wave behavior might have ecological triggers resulting in an overlapping of distinct genetic generations, and arising distinct migratory and nonmigratory genetic patterns living in the same area. An abrupt decrease in the effective population size of the P. costatus populations in the recent past (1960-80) was likely driven by environment fragmentation promoted by the construction of the Três Marias hydropower dam. The low allelic diversity that resulted from this event is still detected today; thus, active stocking programs are not effective at expanding the genetic diversity of this species in the river basin. Finally, this study highlights the importance of using mixed methods to understand spatial and temporal variation in genetic structure for effective mitigation and conservation programs for threatened species that are directly affected by human actions.

Effective population size in a partially clonal plant is not predicted by the number of genetic individuals

Estimating effective population size (N _e) is important for theoretical and practical applications in evolutionary biology and conservation. Nevertheless, estimates of N _e in organisms with complex life-history traits remain scarce because of the challenges associated with estimation methods. Partially clonal plants capable of both vegetative (clonal) growth and sexual reproduction are a common group of organisms for which the discrepancy between the apparent number of individuals (ramets) and the number of genetic individuals (genets) can be striking, and it is unclear how this discrepancy relates to N _e. In this study, we analysed two populations of the orchid Cypripedium calceolus to understand how the rate of clonal versus sexual reproduction affected N _e. We genotyped >1000 ramets at microsatellite and SNP loci, and estimated contemporary N _e with the linkage disequilibrium method, starting from the theoretical expectation that variance in reproductive success among individuals caused by clonal reproduction and by constraints on sexual reproduction would lower N _e. We considered factors potentially affecting our estimates, including different marker types and sampling strategies, and the influence of pseudoreplication in genomic data sets on N _e confidence intervals. The magnitude of N _e/N _ramets and N _e/N _genets ratios we provide may be used as reference points for other species with similar life-history traits. Our findings demonstrate that N _e in partially clonal plants cannot be predicted based on the number of genets generated by sexual reproduction, because demographic changes over time can strongly influence N _e. This is especially relevant in species of conservation concern in which population declines may not be detected by only ascertaining the number of genets.

Genetic diversity, structure, and effective population size of an endangered, endemic hoary bat, 'ōpe'ape'a, across the Hawaiian Islands

Island bat species are disproportionately at risk of extinction, and Hawai'i's only native terrestrial land mammal, the Hawaiian hoary bat (Lasiurus semotus) locally known as 'ōpe'ape'a, is no exception. To effectively manage this bat species with an archipelago-wide distribution, it is important to determine the population size on each island and connectivity between islands. We used 18 nuclear microsatellite loci and one mitochondrial gene from 339 individuals collected from 1988-2020 to evaluate genetic diversity, population structure and estimate effective population size on the Islands of Hawai'i, Maui, O'ahu, and Kaua'i. Genetic differentiation occurred between Hawai'i and Maui, both of which were differentiated from O'ahu and Kaua'i. The population on Maui presents the greatest per-island genetic diversity, consistent with their hypothesized status as the original founding population. A signature of isolation by distance was detected between islands, with contemporary migration analyses indicating limited gene flow in recent generations, and male-biased sex dispersal within Maui. Historical and long-term estimates of genetic effective population sizes were generally larger than contemporary estimates, although estimates of contemporary genetic effective population size lacked upper bounds in confidence intervals for Hawai'i and Kaua'i. Contemporary genetic effective population sizes were smaller on O'ahu and Maui. We also detected evidence of past bottlenecks on all islands with the exception of Hawai'i. Our study provides population-level estimates for the genetic diversity and geographic structure of 'ōpe'ape'a, that could be used by agencies tasked with wildlife conservation in Hawai'i.

Genetic variation of endangered Jankowski’s Bunting (Emberiza jankowskii): High connectivity and a moderate history of demographic decline

Introduction

Continued discovery of “mismatch” patterns between population size and genetic diversity, involving wild species such as insects, amphibians, birds, mammals, and others, has raised issues about how population history, especially recent dynamics under human disturbance, affects currently standing genetic variation. Previous studies have revealed high genetic diversity in endangered Jankowski’s Bunting. However, it is unclear how the demographic history and recent habitat changes shape the genetic variation of Jankowski’s Bunting.

Methods

To explore the formation and maintenance of high genetic diversity in endangered Jankowski’s Bunting, we used a mitochondrial control region (partial mtDNA CR) and 15 nuclear microsatellite markers to explore the recent demographic history of Jankowski’s Bunting, and we compared the historical and contemporary gene flows between populations to reveal the impact of habitat change on population connectivity. Specifically, we aimed to test the following hypotheses: (1) Jankowski’s Bunting has a large historical Ne and a moderate demographic history; and (2) recent habitat change might have no significant impact on the species’ population connectivity.

Results

The results suggested that large historical effective population size, as well as severe but slow population decline, may partially explain the high observable genetic diversity. Comparison of historical (over the past 4Ne generations) and contemporary (1–3 generations) gene flow indicated that the connectivity between five local populations was only marginally affected by landscape changes.

Discussion

Our results suggest that high population connectivity and a moderate history of demographic decline are powerful explanations for the rich genetic variation in Jankowski’s Bunting. Although there is no evidence that the genetic health of Jankowski’s Bunting is threatened, the time-lag effects on the genetic response to recent environmental changes is a reminder to be cautious about the current genetic characteristics of this species. Where possible, factors influencing genetic variation should be integrated into a systematic framework for conducting robust population health assessments. Given the small contemporary population size, inbreeding, and ecological specialization, we recommend that habitat protection be maintained to maximize the genetic diversity and population connectivity of Jankowski’s Bunting.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

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.

Simulation-Based Evaluation of Methods, Data Types, and Temporal Sampling Schemes for Detecting Recent Population Declines

Understanding recent population trends is critical to quantifying species vulnerability and implementing effective management strategies. To evaluate the accuracy of genomic methods for quantifying recent declines (beginning <120 generations ago), we simulated genomic data using forward-time methods (SLiM) coupled with coalescent simulations (msprime) under a number of demographic scenarios. We evaluated both site frequency spectrum (SFS)-based methods (momi2, Stairway Plot) and methods that employ linkage disequilibrium information (NeEstimator, GONE) with a range of sampling schemes (contemporary-only samples, sampling two time points, and serial sampling) and data types (RAD-like data and whole-genome sequencing). GONE and momi2 performed best overall, with >80% power to detect severe declines with large sample sizes. Two-sample and serial sampling schemes could accurately reconstruct changes in population size, and serial sampling was particularly valuable for making accurate inferences when genotyping errors or minor allele frequency cutoffs distort the SFS or under model mis-specification. However, sampling only contemporary individuals provided reliable inferences about contemporary size and size change using either site frequency or linkage-based methods, especially when large sample sizes or whole genomes from contemporary populations were available. These findings provide a guide for researchers designing genomics studies to evaluate recent demographic declines.

The estimates of effective population size based on linkage disequilibrium are virtually unaffected by natural selection

The effective population size (Ne) is a key parameter to quantify the magnitude of genetic drift and inbreeding, with important implications in human evolution. The increasing availability of high-density genetic markers allows the estimation of historical changes in Ne across time using measures of genome diversity or linkage disequilibrium between markers. Directional selection is expected to reduce diversity and Ne, and this reduction is modulated by the heterogeneity of the genome in terms of recombination rate. Here we investigate by computer simulations the consequences of selection (both positive and negative) and recombination rate heterogeneity in the estimation of historical Ne. We also investigate the relationship between diversity parameters and Ne across the different regions of the genome using human marker data. We show that the estimates of historical Ne obtained from linkage disequilibrium between markers (NeLD) are virtually unaffected by selection. In contrast, those estimates obtained by coalescence mutation-recombination-based methods can be strongly affected by it, which could have important consequences for the estimation of human demography. The simulation results are supported by the analysis of human data. The estimates of NeLD obtained for particular genomic regions do not correlate, or they do it very weakly, with recombination rate, nucleotide diversity, proportion of polymorphic sites, background selection statistic, minor allele frequency of SNPs, loss of function and missense variants and gene density. This suggests that NeLD measures mainly reflect demographic changes in population size across generations.

Relative Precision of the Sibship and LD Methods for Estimating Effective Population Size With Genomics-Scale Datasets

Computer simulations were used to compare relative precision of 2 widely used single-sample methods for estimating effective population size (Ne)-the sibship method and the linkage disequilibrium (LD) method. Emphasis is on performance when thousands of gene loci are used, which now can easily be achieved even for nonmodel species. Results show that unless Ne is very small, if at least 500-2000 diallelic loci are used, precision of the LD method is higher than the maximum possible precision for the sibship method, which occurs when all sibling relationships have been correctly identified. Results also show that when precision is high for both methods, their estimates of Ne are highly and positively correlated, which limits additional gains in precision that might be obtained by combining information from the 2 estimators.

Global Commitments to Conserving and Monitoring Genetic Diversity Are Now Necessary and Feasible

Global conservation policy and action have largely neglected protecting and monitoring genetic diversity—one of the three main pillars of biodiversity. Genetic diversity (diversity within species) underlies species’ adaptation and survival, ecosystem resilience, and societal innovation. The low priority given to genetic diversity has largely been due to knowledge gaps in key areas, including the importance of genetic diversity and the trends in genetic diversity change; the perceived high expense and low availability and the scattered nature of genetic data; and complicated concepts and information that are inaccessible to policymakers. However, numerous recent advances in knowledge, technology, databases, practice, and capacity have now set the stage for better integration of genetic diversity in policy instruments and conservation efforts. We review these developments and explore how they can support improved consideration of genetic diversity in global conservation policy commitments and enable countries to monitor, report on, and take action to maintain or restore genetic diversity.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species

Effective population size affects the efficacy of selection, rate of evolution by drift, and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides buttery populations maintain low contemporary (variance) effective population sizes (e.g., ~200 individuals) and thus evolve rapidly by genetic drift. In contrast, populations harbored high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational time scales. Consistent with this hypothesis, we detected low but non-trivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations, and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

Genomic inference of contemporary effective population size in a large island population of collared flycatchers (Ficedula albicollis)

Due to its central importance to many aspects of evolutionary biology and population genetics, the long-term effective population size (N_e ) has been estimated for numerous species and populations. However, estimating contemporary N_e is difficult and in practice this parameter is often unknown. In principle, contemporary N_e can be estimated using either analyses of temporal changes in allele frequencies, or the extent of linkage disequilibrium (LD) between unlinked markers. We applied these approaches to estimate contemporary N_e of a relatively recently founded island population of collared flycatchers (Ficedula albicollis). We sequenced the genomes of 85 birds sampled in 1993 and 2015, and applied several temporal methods to estimate N_e at a few thousand (4000-7000). The approach based on LD provided higher estimates of N_e (20,000-32,000) and was associated with high variance, often resulting in infinite N_e . We conclude that whole-genome sequencing data offers new possibilities to estimate high (>1000) contemporary N_e , but also note that such estimates remain challenging, in particular for LD-based methods for contemporary N_e estimation.

Muskrats as a bellwether of a drying delta

Wetlands worldwide are under threat from anthropogenic impacts. In large protected North American areas such as Yellowstone and Wood Buffalo National Parks, aquatic habitats are disappearing and wetland-dependent fauna are in decline^1-3. Here we investigate population dynamics of an indicator species in Canada's Peace-Athabasca Delta ("the delta"), a World Heritage Site. Based on population surveys, habitat mapping and genetic data from 288 muskrats, we use agent-based modeling and genetic analyses to explain population expansion and decline of the semi-aquatic muskrat (Ondatra zibethicus). Simulations quantify a large population (~500,000 individuals) following flood-induced habitat gains, with decreased size (~10,000 individuals) during drying. Genetic analyses show extremely low long-term effective population size (N_e: 60-127), supporting a legacy of population bottlenecks. Our simulations indicate that the muskrat population in the delta is a metapopulation with individuals migrating preferentially along riparian pathways. Related individuals found over 40 km apart imply dispersal distances far greater than their typical home range (130 m). Rapid metapopulation recovery is achieved via riparian corridor migration and passive flood-transport of individuals. Source-sink dynamics show wetland loss impacts on the muskrat metapopulation's spatial extent. Dramatic landscape change is underway, devastating local fauna, including this generalist species even in a protected ecosystem.

Pedigree-based assessment of recent population connectivity in a threatened rattlesnake

Managing endangered species in fragmented landscapes requires estimating dispersal rates between populations over contemporary timescales. Here, we developed a new method for quantifying recent dispersal using genetic pedigree data for close and distant kin. Specifically, we describe an approach that infers missing shared ancestors between pairs of kin in habitat patches across a fragmented landscape. We then applied a stepping-stone model to assign unsampled individuals in the pedigree to probable locations based on minimizing the number of movements required to produce the observed locations in sampled kin pairs. Finally, we used all pairs of reconstructed parent-offspring sets to estimate dispersal rates between habitat patches under a Bayesian model. Our approach measures connectivity over the timescale represented by the small number of generations contained within the pedigree and so is appropriate for estimating the impacts of recent habitat changes due to human activity. We used our method to estimate recent movement between newly discovered populations of threatened Eastern Massasauga rattlesnakes (Sistrurus catenatus) using data from 2996 RAD-based genetic loci. Our pedigree analyses found no evidence for contemporary connectivity between five genetic groups, but, as validation of our approach, showed high dispersal rates between sample sites within a single genetic cluster. We conclude that these five genetic clusters of Eastern Massasauga rattlesnakes have small numbers of resident snakes and are demographically isolated conservation units. More broadly, our methodology can be widely applied to determine contemporary connectivity rates, independent of bias from shared genetic similarity due to ancestry that impacts other approaches.

Genetic structure of recently fragmented suburban populations of European stag beetle

Habitat loss and fragmentation due to urbanization can negatively affect metapopulation persistence when gene flow among populations is reduced and population sizes decrease. Inference of patterns and processes of population connectivity derived from spatial genetic analysis has proven invaluable for conservation and management. However, a more complete account of population dynamics may be obtained by combining spatial and temporal sampling. We, therefore, performed a genetic study on European stag beetle (Lucanus cervus L.) populations in a suburban context using samples collected in three locations and during the period 2002-2016. The sampling area has seen recent landscape changes which resulted in population declines. Through the use of a suite of F ST, clustering analysis, individual assignment, and relatedness analysis, we assessed fine scale spatiotemporal genetic variation within and among habitat patches using 283 individuals successfully genotyped at 17 microsatellites. Our findings suggested the three locations to hold demographically independent populations, at least over time scales of relevance to conservation, though with higher levels of gene flow in the past. Contrary to expectation from tagging studies, dispersal appeared to be mainly female-biased. Although the life cycle of stag beetle suggests its generations to be discrete, no clear temporal structure was identified, which could be attributed to the varying duration of larval development. Since population bottlenecks were detected and estimates of effective number of breeders were low, conservation actions are eminent which should include the establishment of suitable dead wood for oviposition on both local and regional scales to increase (re)colonization success and connectivity among current populations.

Insights from Population Genomics to Enhance and Sustain Biological Control of Insect Pests

Biological control-the use of organisms (e.g., nematodes, arthropods, bacteria, fungi, viruses) for the suppression of insect pest species-is a well-established, ecologically sound and economically profitable tactic for crop protection. This approach has served as a sustainable solution for many insect pest problems for over a century in North America. However, all pest management tactics have associated risks. Specifically, the ecological non-target effects of biological control have been examined in numerous systems. In contrast, the need to understand the short- and long-term evolutionary consequences of human-mediated manipulation of biological control organisms for importation, augmentation and conservation biological control has only recently been acknowledged. Particularly, population genomics presents exceptional opportunities to study adaptive evolution and invasiveness of pests and biological control organisms. Population genomics also provides insights into (1) long-term biological consequences of releases, (2) the ecological success and sustainability of this pest management tactic and (3) non-target effects on native species, populations and ecosystems. Recent advances in genomic sequencing technology and model-based statistical methods to analyze population-scale genomic data provide a much needed impetus for biological control programs to benefit by incorporating a consideration of evolutionary consequences. Here, we review current technology and methods in population genomics and their applications to biological control and include basic guidelines for biological control researchers for implementing genomic technology and statistical modeling.

Genetic signature of disease epizootic and reintroduction history in an endangered carnivore

Emerging infectious diseases have recently increased in wildlife and can result in population declines and the loss of genetic diversity in susceptible populations. As populations of impacted species decline, genetic diversity can be lost, with ramifications including reduced effective population size and increased population structuring. For species of conservation concern, which may already have low genetic diversity, the loss of genetic diversity can be especially important. To investigate the impacts of a novel pathogen on genetic diversity in a genetically depauperate endangered species, we assessed the ramifications of a sylvatic plague-induced bottleneck in black-footed ferrets (Mustela nigripes). Following a plague epizootic, we genotyped 184 ferrets from Conata Basin and Badlands National Park, South Dakota, at seven microsatellite loci. We compared our results to pre-plague studies in the same population. We observed population substructuring into three genetic clusters. These clusters reflect founder effects from ferret reintroduction events followed by genetic drift. Compared to the pre-plague population, we observed losses of allelic diversity in all clusters, as well as significantly reduced heterozygosity in one cluster. These results indicate that disease epizootics may reduce population size and also genetic diversity. Our results suggest the importance of early and sustained management in mitigating disease epizootics in naïve populations for the maintenance of genetic diversity.

Population viability analysis of small population: a case study for Asian elephant in China

Small populations are at risk of extinction from deterministic and stochastic factors. Less than 250 Asian elephants (Elephas maximus) remain in China, and are distributed in a few isolated areas; yet, population viability analyses of this endangered population have not been conducted. Here, the current genetic status of the Pu'Er-Mengyang Asian elephant populations in China was analyzed, and the risk of extinction was predicted over the next 500 years. Factors affecting the viability of this population were determined through simulations. The genetic diversity of the population was very low (mean allele number: 3.1; expected heterozygosity: 0.463), even though a recent population bottleneck was not detected. The effective population size was approximately 24.1 adult elephants. Enough adult breeding individuals exist to maintain population viability. VORTEX simulation model showed that this population would not go extinct in the next 500 years. However, illegal poaching and harvesting could negatively affect population size. A sensitivity analysis showed that the mean stochastic growth rate of the study population is sensitive to sex ratio, number of breeding females, mortality of females of different age classes, carrying capacity, and lethal equivalents. Based on our results, we suggest that action should be taken to alleviate inbreeding and any further loss of genetic diversity, by connecting fragmented elephant habitat or by translocating individual elephants. In addition, human-elephant conflict should be mitigated using various modern approaches, including crop guarding techniques, and by encouraging farmers to switch to crops and income sources not vulnerable to elephant raids.

Variation in brown rat cranial shape shows directional selection over 120 years in New York City

Urbanization exposes species to novel environments and selection pressures that may change morphological traits within a population. We investigated how the shape and size of crania and mandibles changed over time within a population of brown rats (Rattus norvegicus) living in Manhattan, New York, USA, a highly urbanized environment. We measured 3D landmarks on the cranium and mandible of 62 adult individuals sampled in the 1890s and 2010s. Static allometry explained approximately 22% of shape variation in crania and mandible datasets, while time accounted for approximately 14% of variation. We did not observe significant changes in skull size through time or between the sexes. Estimating the P-matrix revealed that directional selection explained temporal change of the crania but not the mandible. Specifically, rats from the 2010s had longer noses and shorter upper molar tooth rows, traits identified as adaptive to colder environments and higher quality or softer diets, respectively. Our results highlight the continual evolution to selection pressures. We acknowledge that urban selection pressures impacting cranial shape likely began in Europe prior to the introduction of rats to Manhattan. Yet, our study period spanned changes in intensity of artificial lighting, human population density, and human diet, thereby altering various aspects of rat ecology and hence pressures on the skull.

Long-term stocking practices threaten the original genetic diversity of the southernmost European populations of Atlantic salmon Salmo salar

Many Atlantic salmon Salmo salar populations in Europe are threatened by previous stocking with foreign hatchery strains. Temporal patterns of genetic characteristics of salmon from northern Spain, the southernmost European populations, were compared before and after species decline and heavy stocking with specimens from northern Europe. Eleven microsatellite loci were analysed in archival (scales from 1958-1960) and contemporary (2007-2008) samples from the River Sella. Temporal analyses revealed a similar heterozygosity between archival and contemporary samples, despite a drastic decrease in population abundance, while the contemporary sample showed a higher allelic richness due to the occurrence of foreign alleles. Considering only the alleles with at least 4% frequency in the archival sample, 2 alleles exclusive to the River Sella were absent in the contemporary sample, and 14 alleles showed a decrease of at least 4% frequency. Four alleles common in Scotland showed a high occurrence in the contemporary sample, so they are good candidates as markers of introgression of foreign genes. The heavy stocking with non-native Scottish broodstocks between 1970 and 1990 caused the introgression found in the contemporary sample when compared with the pristine population. An abrupt decrease was evident when the estimates of effective number of breeders were adjusted to take into account overlapping generations (NbAdj), effective population size (NeAdj) estimated from NbAdj, and number of breeders estimated using the sibship assignment method (NbSIB). The very low effective size values found in the contemporary sample, together with the detrimental synergy between genetic drift and high rates of introgression, represent a severe risk for the conservation of native salmon.