Amounts of genetic drift and the effective size of populations can be estimated from observed temporal shifts in sample allele frequencies. Bias in this so-called temporal method has been noted in ...cases of small sample sizes and when allele frequencies are highly skewed. We characterize bias in commonly applied estimators under different sampling plans and propose an alternative estimator for genetic drift and effective size that weights alleles differently. Numerical evaluations of exact probability distributions and computer simulations verify that this new estimator yields unbiased estimates also when based on a modest number of alleles and loci. At the cost of a larger standard deviation, it thus eliminates the bias associated with earlier estimators. The new estimator should be particularly useful for microsatellite loci and panels of SNPs, representing a large number of alleles, many of which will occur at low frequencies.
Large-scale exploitation of wild animals and plants through fishing, hunting and logging often depends on augmentation through releases of translocated or captively raised individuals. Such releases ...are performed worldwide in vast numbers. Augmentation can be demographically and economically beneficial but can also cause four types of adverse genetic change to wild populations: (1) loss of genetic variation, (2) loss of adaptations, (3) change of population composition, and (4) change of population structure. While adverse genetic impacts are recognized and documented in fisheries, little effort is devoted to actually monitoring them. In forestry and wildlife management, genetic risks associated with releases are largely neglected. We outline key features of programs to effectively monitor consequences of such releases on natural populations.
Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic ...differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.
Atlantic herring is an excellent species for studying the genetic basis of adaptation in geographically distant populations because of its characteristically large population sizes and low genetic ...drift. In this study we compared whole-genome resequencing data of Atlantic herring populations from both sides of the Atlantic Ocean. An important finding was the very low degree of genetic differentiation among geographically distant populations (fixation index = 0.026), suggesting lack of reproductive isolation across the ocean. This feature of the Atlantic herring facilitates the detection of genetic factors affecting adaptation because of the sharp contrast between loci showing genetic differentiation resulting from natural selection and the low background noise resulting from genetic drift. We show that genetic factors associated with timing of reproduction are shared between genetically distinct and geographically distant populations. The genes for thyroid-stimulating hormone receptor (TSHR), the SOX11 transcription factor (SOX11), calmodulin (CALM), and estrogen receptor 2 (ESR2A), all with a significant role in reproductive biology, were among the loci that showed the most consistent association with spawning time throughout the species range. In fact, the same two SNPs located at the 5′ end of TSHR showed the most significant association with spawning time in both the east and west Atlantic. We also identified unexpected haplotype sharing between spring-spawning oceanic herring and autumn-spawning populations across the Atlantic Ocean and the Baltic Sea. The genomic regions showing this pattern are unlikely to control spawning time but may be involved in adaptation to ecological factor(s) shared among these populations.
There are two primary measures of the amount of genetic variation in a population at a locus: heterozygosity and the number of alleles. Effective population size (Ne) provides both an expectation of ...the amount of heterozygosity in a population at drift‐mutation equilibrium and the rate of loss of heterozygosity because of genetic drift. In contrast, the number of alleles in a population at drift‐mutation equilibrium is a function of both Ne and census size (NC). In addition, populations with the same Ne can lose allelic variation at very different rates. Allelic variation is generally much more sensitive to bottlenecks than heterozygosity. Expressions used to adjust for the effects of violations of the ideal population on Ne do not provide good predictions of the loss of allelic variation. These effects are much greater for loci with many alleles, which are often important for adaptation. We show that there is a linear relationship between the reduction of NC and the corresponding reduction of the expected number of alleles at drift‐mutation equilibrium. This makes it possible to predict the expected effect of a bottleneck on allelic variation. Heterozygosity provides good estimates of the rate of adaptive change in the short‐term, but allelic variation provides important information about long‐term adaptive change. The guideline of long‐term Ne being greater than 500 is often used as a primary genetic metric for evaluating conservation status. We recommend that this guideline be expanded to take into account allelic variation as well as heterozygosity.
Population census size (N C) and effective population sizes (N e) are two crucial parameters that influence population viability, wildlife management decisions, and conservation planning. Genetic ...estimators of both N C and N e are increasingly widely used because molecular markers are increasingly available, statistical methods are improving rapidly, and genetic estimators complement or improve upon traditional demographic estimators. We review the kinds and applications of estimators of both N C and N e, and the often undervalued and misunderstood ratio of effective-to-census size (N e /N C). We focus on recently improved and well evaluated methods that are most likely to facilitate conservation. Finally, we outline areas of future research to improve N e and N C estimation in wild populations.
The Atlantic herring (Clupea harengus), one of the most abundant marine fishes in the world, has historically been a critical food source in Northern Europe. It is one of the few marine species that ...can reproduce throughout the brackish salinity gradient of the Baltic Sea. Previous studies based on few genetic markers have revealed a conspicuous lack of genetic differentiation between geographic regions, consistent with huge population sizes and minute genetic drift. Here, we present a cost-effective genome-wide study in a species that lacks a genome sequence. We first assembled a muscle transcriptome and then aligned genomic reads to the transcripts, creating an “exome assembly,” capturing both exons and flanking sequences. We then resequenced pools of fish from a wide geographic range, including the Northeast Atlantic, as well as different regions in the Baltic Sea, aligned the reads to the exome assembly, and identified 440,817 SNPs. The great majority of SNPs showed no appreciable differences in allele frequency among populations; however, several thousand SNPs showed striking differences, some approaching fixation for different alleles. The contrast between low genetic differentiation at most loci and striking differences at others implies that the latter category primarily reflects natural selection. A simulation study confirmed that the distribution of the fixation index F ST deviated significantly from expectation for selectively neutral loci. This study provides insights concerning the population structure of an important marine fish and establishes the Atlantic herring as a model for population genetic studies of adaptation and natural selection.
Population translocations occur for a variety of reasons, from displacement due to climate change to human‐induced transfers. Such actions have adverse effects on genetic variation and understanding ...their microevolutionary consequences requires monitoring. Here, we return to an experimental release of brown trout (Salmo trutta) in order to monitor the genomic effects of population translocations. In 1979, fish from each of two genetically (FST = 0.16) and ecologically separate populations were simultaneously released, at one point in time, to a lake system previously void of brown trout. Here, whole‐genome sequencing of pooled DNA (Pool‐seq) is used to characterize diversity within and divergence between the introduced populations and fish inhabiting two lakes downstream of the release sites, sampled 30 years later (c. 5 generations). Present results suggest that while extensive hybridization has occurred, the two introduced populations are unequally represented in the lakes downstream of the release sites. One population, which is ecologically resident in its original habitat, mainly contributes to the lake closest to the release site. The other population, migratory in its natal habitat, is genetically more represented in the lake further downstream. Genomic regions putatively under directional selection in the new habitat are identified, where allele frequencies in both established populations are more similar to the introduced population stemming from a resident population than the migratory one. Results suggest that the microevolutionary consequences of population translocations, for example, hybridization and adaptation, can be rapid and that Pool‐seq can be used as an initial tool to monitor genome‐wide effects.
The consequence of human induced translocations on natural populations is monitored over c. 3 decades in wild brown trout (Salmo trutta). Two genetically distinct populations released to the same environment are unequally represented in descendant populations, with one introduced population mainly contributing to the lake closest to the release site and the other dominating waters further downstream.
The variance effective population size (
N
eV
) 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
N
eV
of a spatially structured population that are of interest for conservation of species. Since
N
eV
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
N
eV
of the metapopulation is maximized, when subpopulations are weighted proportionally to their long term reproductive contributions. This maximal
N
eV
is the sought for global effective size, since it approximates the gene diversity effective size
N
eGD
, 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
N
eV
, one of which (the backward version of
N
eV
) is most stable, exists for most populations, and is closer to
N
eGD
than the classical notion of
N
eV
. Expressions for the optimal length of the time interval for measuring genetic change are developed, that make it possible to estimate any version of
N
eV
with maximal accuracy.