Disentangling the impact of demographic factors on population differentiation of an endangered freshwater crayfish (Austropotamobius pallipes) using population density and microsatellite data

1. Habitat fragmentation of stream ecosystems often results in decreased connectivity between populations and lower population sizes. Hence, understanding how habitat fragmentation affects genetic erosion is important for the preservation of freshwater biodiversity, in particular, as small populations suffer from loss of genetic diversity through genetic drift and loss of fitness because of inbreeding, increasing the risk of extinction. 2. Here, we assess the impact of demographic factors on population differentiation in the endangered freshwater crayfish Austropotamobius pallipes by analysing population genetic structure, estimating effective population sizes and comparing levels of polymorphism at five microsatellite loci with density estimates of 10 populations within a small French catchment that has become progressively confined to headwaters over the last six decades. 3. Levels of expected heterozygosity and allelic richness per population were relatively low (0.214–0.396 and 1.6–2.6, respectively). We found strong genetic differentiation between these geographically close populations (F_ST = 0.283), with weak statistical evidence for a pattern of isolation by distance. Estimates of effective population size were low (<150) in most populations, but potentially reached several thousands in three populations. 4. Population density and allelic richness were strongly positively correlated. A robust relationship between population density and heterozygosity values was also noted, but only after discarding two populations for which significant genetic signatures of a recent bottleneck were found; these two populations displayed high expected heterozygosity compared with a very low density. Populations with the highest densities of individuals had the highest effective population size estimates and vice versa. 5. Our results clearly show the importance of demographic factors and genetic drift on A. pallipes populations. Furthermore, analysis of genetic and population density data is a pragmatic and efficient approach to corroborate inferences from genetic data and can be particularly useful in the identification of populations experiencing a bottleneck and therefore in conservation genetics studies aiming at identifying priority populations for conservation.     

EFFECTIVE POPULATION SIZE IN A CONTINUOUSLY DISTRIBUTED POPULATION

An individual-based simulation model was created to examine genetic variability, time until fixation and spatial genetic structure in a continuously distributed population. Previous mathematical models for continuously distributed populations have the difficulty that the assumption of independent reproduction and independent dispersal of offspring cause clumped spatial distribution and thus violate an assumption of random spatial distribution. In this study, this problem is avoided by considering the dispersal behavior of offspring. The simulation results showed that the inbreeding effective population size estimated by the rate of decrease of heterozygosity during the first 15 generations corresponds to the neighborhood size calculated by the standard deviation of the dispersal distance (σ_T). This inbreeding effective population size does not greatly change with the area of simulation when the densities and σ_T are the same. However, the inbreeding effective population size estimated by heterozygosity using the first 500 generations is larger than the neighborhood size calculated by the dispersal distance and increases with the area of simulation with the same densities. The variance effective population size, estimated by time until fixation of alleles, increases with dispersal distance (σ_T) and with the area of simulation given the same densities. The inbreeding effective population size and variance effective population size were smaller than the actual population size unless σ_T is sufficiently large (2 σ_T > approximate L/2, where L is a side of the simulation square).     

Short tandem repeat based analysis of genetic variability in Kanarese buffalo of South India

The goal of the present study was assessing genetic diversity within Kanarese buffalo, the dual purpose breed of South India. A total of 48 unrelated animals were genotyped at 23 short tandem repeat (STR markers)loci. The total number of observed alleles was 180 with a mean of 7.83 per locus, which varied from 3 to 12 across different loci. The mean observed and expected heterozygosity in South Kanara buffaloes was estimated to be 0.518 and 0.712 respectively. Within population inbreeding estimate (F _IS) was significantly positive in most of the investigated loci which resulted in significant deviation from Hardy-Weinberg equilibrium at 19 of 23 loci analyzed. Evaluation of South Kanara buffalo population for mutation drift equilibrium revealed no significant heterozygosity excess under three different models of evolution viz. infinite alleles model (IAM), step-wise mutation model (SMM) and two phase model (TPM), thus indicating the absence of any recent genetic bottleneck. The results of the present study will help in formulating rational breeding strategies as well as conservation of this important germplasm.     

Founder effect and bottleneck signatures in an introduced,insular population of elk

The population of elk (Cervus elaphus roosevelti) inhabiting Afognak Island, Alaska, USA arose from an introduction of 8 individuals from an established population in Washington, USA in 1929, and recently peaked at approximately 1,400 individuals. We examined indices of diversity for 15 microsatellite loci in the Afognak population and compared them to levels in the parent population to determine effects of translocation and demography on genetic variation. The Afognak population differed significantly (P < 0.0001) from the source population in both allele and genotype frequencies. Allelic richness, number of private alleles and multilocus heterozygosity, but not percent loci polymorphic, were significantly lower in Afognak elk. Mean inbreeding coefficients within Afognak (f = 0.019) and source (f = −0.006) populations did not differ significantly from zero. Despite the demographic bottleneck, no evidence of a genetic bottleneck in the Afognak population was detected using a test for heterozygosity excess or mode shift of allele frequencies. Simulations indicated that rapid population growth after the translocation resulted in heterozygosity excess for only 8 years. Conversely, a statistic testing for a bottleneck signature in the ratio of allele number to allele size range (M-ratio) was significant for both the Afognak and source populations, suggesting that the Afognak population had effectively undergone serial bottlenecks. Nonetheless, Afognak failed to show a smaller M-ratio than the parent population, suggesting a failure of that statistic to detect the bottleneck associated with introduction. We show that a severe bottleneck followed by rapid population growth may be undetectable using available tests.     

More precisely biased: increasing the number of markers is not a silver bullet in genetic bottleneck testing

In response to our review of the use of genetic bottleneck tests in the conservation literature (Peery et al. 2012, Molecular Ecology, 21 , 3403–3418), Hoban et al. (2013, Molecular Ecology, in press) conducted population genetic simulations to show that the statistical power of genetic bottleneck tests can be increased substantially by sampling large numbers of microsatellite loci, as they suggest is now possible in the age of genomics. While we agree with Hoban and co‐workers in principle, sampling large numbers of microsatellite loci can dramatically increase the probability of committing type 1 errors (i.e. detecting a bottleneck in a stable population) when the mutation model is incorrectly assumed. Using conservative values for mutation model parameters can reduce the probability of committing type 1 errors, but doing so can result in significant losses in statistical power. Moreover, we believe that practical limitations associated with developing large numbers of high‐quality microsatellite loci continue to constrain sample sizes, a belief supported by a literature review of recent studies using next generation sequencing methods to develop microsatellite libraries. conclusion, we maintain that researchers employing genetic bottleneck tests should proceed with caution and carefully assess both statistical power and type 1 error rates associated with their study design.     

Microsatellite analysis of a population crash and bottleneck in the Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense (Asteraceae), and its implications for reintroduction

The Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense, has experienced both a severe population crash associated with an increase in alien ungulate populations on Mauna Kea, and a population bottleneck associated with reintroduction. In this paper, we address the genetic consequences of both demographic events using eight microsatellite loci. The population crash was not accompanied by a significant reduction in number of alleles or heterozygosity. However, the population bottleneck was accompanied by significant reductions in observed number of alleles, effective number of alleles, and expected heterozygosity, though not in observed heterozygosity. The effective size of the population bottleneck was calculated using both observed heterozygosities and allele frequency variances. Both methods corroborated the historical census size of the population bottleneck of at most three individuals. The results suggest that: (i) small populations, even those that result from severe reductions in historical population size and extent, are not necessarily genetically depauperate; and (ii) species reintroduction plans need to be conceived and implemented carefully, with due consideration to the genetic impact of sampling for reintroduction.     

Analysis of molecular genetic diversity and population structure in Amaranthus germplasm using SSR markers

We assessed the molecular genetic diversity and population structure of Amaranthus species accessions using 11 simple sequence repeat markers. A total of 122 alleles were detected, and the number of alleles per marker (N_A) ranged from 6 to 21 with an average of 11.1 alleles. The frequency of major alleles per locus ranged from 0.148 to 0.695, with an average value of 0.496 per marker. The overall polymorphic information content values were 0.436–0.898, with an average value of 0.657. The observed heterozygosity (H_O) and expected heterozygosity (H_E) ranged from 0.056 to 0.876 and from 0.480 to 0.907, with average values of 0.287 and 0.698, respectively. The average H_O (0.240) was lower than the H_E and gene flow (Nm), and showed substantial genetic variability among all populations of amaranth accessions. The sample groupings did not strictly follow the geographic affiliations of the accessions. A similar pattern was obtained using model-based structure analysis without grouping by species type. Knowledge of the genetic diversity and population structure of amaranth can be used to select representative genotypes and manage Amaranthus germplasm breeding programs.     

When are genetic methods useful for estimating contemporary abundance and detecting population trends?

The utility of microsatellite markers for inferring population size and trend has not been rigorously examined, even though these markers are commonly used to monitor the demography of natural populations. We assessed the ability of a linkage disequilibrium estimator of effective population size (N_e) and a simple capture-recapture estimator of abundance (N) to quantify the size and trend of stable or declining populations (true N = 100–10,000), using simulated Wright–Fisher populations. Neither method accurately or precisely estimated abundance at sample sizes of S = 30 individuals, regardless of true N. However, if larger samples of S = 60 or 120 individuals were collected, these methods provided useful insights into abundance and trends for populations of N = 100–500. At small population sizes (N = 100 or 250), precision of the N_e estimates was improved slightly more by a doubling of loci sampled than by a doubling of individuals sampled. In general, monitoring N_e proved a more robust means of identifying stable and declining populations than monitoring N over most of the parameter space we explored, and performance of the N_e estimator is further enhanced if the N_e/N ratio is low. However, at the largest population size (N = 10,000), N estimation outperformed N_e. Both methods generally required ≥ 5 generations to pass between sampling events to correctly identify population trend.     

Novel tetranucleotide microsatellite DNA markers for the wood frog,Rana sylvatica

Fifteen tetranucleotide microsatellite loci were identified and characterized for wood frogs (Rana sylvatica) collected from three vernal pools in the southeastern US. These markers revealed a high degree of genetic diversity (nine to 34 alleles per locus), heterozygosity (30.6–92.3%) and allelic heterogeneity (69% of comparisons were statistically significant). Considerable differentiation among populations was observed as genetic distances (chord) ranged between 0.40 and 0.55 and all F_ST values (0.02–0.05) were statistically significant. Genotypic assignment tests correctly classified 103 of 113 individuals to their respective collection. These markers should prove useful for investigating fine‐scale population structure and metapopulation dynamics.     

Accounting for missing data in the estimation of contemporary genetic effective population size (Ne)

Theoretical models are often applied to population genetic data sets without fully considering the effect of missing data. Researchers can deal with missing data by removing individuals that have failed to yield genotypes and/or by removing loci that have failed to yield allelic determinations, but despite their best efforts, most data sets still contain some missing data. As a consequence, realized sample size differs among loci, and this poses a problem for unbiased methods that must explicitly account for random sampling error. One commonly used solution for the calculation of contemporary effective population size (N_e) is to calculate the effective sample size as an unweighted mean or harmonic mean across loci. This is not ideal because it fails to account for the fact that loci with different numbers of alleles have different information content. Here we consider this problem for genetic estimators of contemporary effective population size (N_e). To evaluate bias and precision of several statistical approaches for dealing with missing data, we simulated populations with known N_e and various degrees of missing data. Across all scenarios, one method of correcting for missing data (fixed‐inverse variance‐weighted harmonic mean) consistently performed the best for both single‐sample and two‐sample (temporal) methods of estimating N_e and outperformed some methods currently in widespread use. The approach adopted here may be a starting point to adjust other population genetics methods that include per‐locus sample size components.     

EVIDENCE FOR AN EXTREME BOTTLENECK IN A RARE MEXICAN PINYON: GENETIC DIVERSITY,DISEQUILIBRIUM, AND THE MATING SYSTEM IN PINUS MAXIMARTINEZII

Maxipiñon (Pinus maximartinezii Rzedowski), which is confined to a single population of approximately 2000 to 2500 mature trees, covers about 400 ha in southern Zacatecas, Mexico. Genetic diversity measured by expected heterozygosity was 0.122, which is moderate for pines. However, percentage polymorphic loci was low, 30.3%. The fixation index (F) of 0.081 indicated only slight heterozygote deficiency. Mating system analysis indicated a significant but low level of selling; the multilocus outcrossing rate, t_m, was 0.816. The mean of single locus estimates, t_s, was smaller (0.761), perhaps suggesting mating among relatives, although the difference between t_m and t_s was not statistically significant. The most striking features of maxipiñon's genetic structure were that no polymorphic locus had more than two alleles and most alleles at polymorphic loci were at intermediate frequencies. This is in contrast to other pines, which often have three to five or more alleles at some loci and in which the distribution of allele frequencies is U-shaped, most alleles being present at frequencies less than 10% or greater than 90%. A population with only two alleles per locus and at intermediate frequencies could occur if the population had been reduced to an extreme bottleneck and then expanded rapidly before random drift modified allele frequencies. A novel origin from a hybridization event would also explain the results. Significant gametic disequilibrium was detected at several pairs of loci in both maternal and paternal gametes. The presence of disequilibrium is in agreement with an origin from an extreme bottleneck, perhaps even a single seed. Furthermore, it demands that the event be relatively recent. The number of generations, as calculated from the observed mean disequilibrium, suggested that maxipiñon derived from an extreme bottleneck four to five generations ago, which is less than 1000 years in this species.     

Bottlenecks and rescue effects in a fluctuating population of golden-mantled ground squirrels (Spermophilus lateralis)

Mammal species characterized by highly fluctuating populations often maintain genetic diversity in response to frequent demographic bottlenecks, suggesting the ameliorating influence of life history and behavioral factors. Immigration in particular is expected to promote genetic recovery and is hypothesized to be the most likely process maintaining genetic diversity in fluctuating mammal populations. Most demographic bottlenecks have been inferred retrospectively, and direct analysis of a natural population before, during, and after a bottleneck is rare. Using a continuous 10-year dataset detailing the complete demographic and genetic history of a fluctuating population of golden-mantled ground squirrels (Spermophilus lateralis), we analyzed the genetic consequences of a 4-year demographic bottleneck that reduced the population to seven adult squirrels, and we evaluated the potential “rescue effect” of immigration. Analysis of six microsatellite loci revealed that, while a decline in allelic richness was observed during the bottleneck, there was no observed excess of heterozygosity, a characteristic bottleneck signature, and no evidence for heterozygote deficiency during the recovery phase. In addition, we found no evidence for inbreeding depression during or after the bottleneck. By identifying immigrants and analyzing their demographic and genetic contributions, we found that immigration promoted demographic recovery and countered the genetic effects of the bottleneck, especially the loss of allelic richness. Within 3 years both population size and genetic variation had recovered to pre-bottleneck levels, supporting the role of immigration in maintaining genetic variation during bottleneck events in fluctuating populations. Our analyses revealed considerable variation among analytical techniques in their ability to detect genetic bottlenecks, suggesting that caution is warranted when evaluating bottleneck events based on one technique.     

Isolation and characterization of novel tetranucleotide microsatellite DNA markers for the spotted salamander,Ambystoma maculatum

Fifteen tetranucleotide microsatellite loci were identified and characterized for spotted salamanders (Ambystoma maculatum) collected from three vernal pools in the south‐eastern USA. These markers revealed a high degree of genetic diversity (7–32 alleles per locus), heterozygosity (31.6–86.3%) and allelic heterogeneity (91% of comparisons were statistically significant). Considerable differentiation among populations was observed as genetic distances (chord) ranged between 0.50 and 0.65, and all F_ST values (0.08–0.14) were statistically significant. Moreover, genotypic assignment tests correctly classified all individuals to their respective collection. These markers should prove useful for investigating fine‐scale population structure and mating system.     

Inference of population subdivision from the VNTR distributions of New Zealanders

A population sample from people of diverse ethnic origins living in New Zealand serves as a database to test methods for inference of population subdivision. The initial null hypothesis, that the population sample is homogeneous across ethnic groups, is easily rejected by likelihood ratio tests. Beyond this, methods for quantifying subdivision can be based on the probability of drawing alleles identical by descent (F _ST ), probabilities of matching multiple locus genotypes, and occurrence of unique alleles. Population genetic theory makes quantitative predictions about the relation betweenF _ST , population sizes, and rates of migration and mutation. Some VNTR loci have mutation rates of 10^–2 per generation, but, contrary to theory, we find no consistent association between the degree of population subdivision and mutation rate. Quantification of population substructure also allows us to relate the magnitudes of genetic distances between ethnic groups in New Zealand to the colonization history of the country. The data suggests that the closest relatives to the Maori are Polynesians, and that no severe genetic bottleneck occurred when the Maori colonized New Zealand. One of the central points of contention regarding the application of VNTR loci in forensics is the appropriate means for estimating match probabilities. Simulations were performed to test the merits of the product rule in the face of subpopulation heterogeneity. Population heterogeneity results in large differences in estimates of multilocus genotype frequencies depending on which subpopulation is used for reference allele frequencies, but, of greater importance for forensic purposes, no five locus genotype had an expected frequency greater than 10^–6. Although this implies that a match with an innocent individual is unlikely, in a large urban area such chance matches are going to occur.Editor's commentsA side-benefit of the collection of DNA data from human populations is the light it may shed on human evolution. The authors discuss the colonization history of New Zealand in the light of such data. From a forensic viewpoint, too much should not be made of the differences between the major ethnic groups within New Zealand, as the forensic community in that country maintains separate databases for Caucasian, Maori and Pacific Islander (Buckletonet al., 1987). It will be of interest in the future to examine subdivision within these groups, as opposed to within the country as a whole. The authors' comments on testing for independence will need to read along with the findings of Zaykinet al. and Maiste and Weir in this volume. The authors had not seen the Budowleet al. (1994) rebuttal to the paper of Kraneet al. (1992).     

The distribution of coalescence times and distances between microsatellite alleles with changing effective population size

We investigate the effects of past changes of the effective population size on the present allelic diversity at a microsatellite marker locus. We first derive the analytical expression of the generating function of the joint probabilities of the time to the Most Recent Common Ancestor for a pair of alleles and of their distance (the difference in allele size). We give analytical solutions in the case of constant population size and the geometrical mutation model. Otherwise, numerical inversion allows the distributions to be calculated in general cases. The effects of population expansion or decrease and the possibility to detect an ancient bottleneck are discussed. The method is extended to samples of three and four alleles, which allows investigating the covariance structure of the frequencies f(k) of pairs of alleles with a size difference of k motifs, and suggesting some approaches to the estimation of past demography.     

The use of agent-based modelling of genetics in conservation genetics studies

Animal Landscape and Man Simulation System a genetically explicit agent-based model was used to obtain measures for the genetic and demographic status of simulated populations. This investigation aimed to test the applicability of this approach for assessing the effect of environmental perturbations on populations’ temporal and spatial dynamics. This was achieved by assessing how three simple scenarios with increasing degree of environmental disturbance, simulated by populations bottlenecks repeated at different intervals, affected the genetic and demographic characteristics of the simulated population. Model outputs from a simplified landscape scenario concurred with theoretical expectations validating the model in a qualitative way. Differences in medians, means and coefficient of variation of the observed (H_o) and expected heterozygosity (H_e), population census size (N), effective population size (N_e), inbreeding coefficient (F) and N_e/N ratio were observed for simulated populations. Impacts occurred rapidly after simulated bottleneck events and genetic estimates were less variable, and therefore more reliable, than demographic estimates. Precise genetic consequences of the bottlenecks repeated at different intervals, and resulting population perturbations, are a complex balance between effects on population sub-structure, size and founding events. Agent-based models are appropriate tools to simulate these interactions, being sufficiently flexible to mimic real population processes under a range of environmental conditions. Such models incorporating explicit genetics provide a promising new approach to evaluate the impact of environmental changes on genetic composition of populations.     

HISTORICAL DEMOGRAPHY AND PRESENT DAY POPULATION STRUCTURE OF THE GREENFINCH,CARDUEUS CHLORIS—AN ANALYSIS OF mtDNA CONTROL-REGION SEQUENCES

Genetic variability within and among 10 geographically distinct populations of Greenfinches (Carduelis chloris) was assayed by directly sequencing a 637 BP part of the mtDNA control region from 194 individuals. Thirteen variable positions defined 18 haplotypes with a maximum sequence divergence of 0.8%. Haplotype (h = 0.28–0.77) and nucleotide (π = 0.058–0.17%) diversities within populations were low, and decreased with increasing latitude (h:r_s = –0.81; π: r_s = –0.89). The distribution of pairwise nucleotide differences fit better with expectations of a “sudden expansion” than of an “equilibrium” model, and the estimates of long term effective population sizes were considerably lower than current census estimates, especially in northern European samples. Selection is an unlikely cause of observed patterns because the distribution of variability conformed to expectations of neutral infinite alleles model and haplotype diversity across populations was positively correlated with heterozygosity (H_E) in nuclear genes (r_s = 0.74, P < 0.05). Hence, a recent bottleneck, followed by serial bottlenecking during the process of post-Pleistocene recolonization of northern Europe, together with recent population expansion provide a plausible explanation for the low genetic diversity in the north. Genetic distances among populations showed a clear pattern of isolation-by-distance, and 14% of the haplotypic variation was among populations, the rest being distributed among individuals within populations. In accordance with allozyme and morphological data, a hierarchical analysis of nucleotide diversity recognized southern European populations as distinct from northern European ones. However, the magnitude of divergence in mtDNA, allozymes and morphology were highly dissimilar (morphology > mtDNA > allozymes).     

On the genetic parameter determining the efficiency of purging: an estimate for Drosophila egg‐to‐pupae viability

The consequences of inbreeding on fitness can be crucial in evolutionary and conservation grounds and depend upon the efficiency of purging against deleterious recessive alleles. Recently, analytical expressions have been derived to predict the evolution of mean fitness, taking into account both inbreeding and purging, which depend on an ‘effective purging coefficient (d_e)’. Here, we explore the validity of that predictive approach and assay the strength of purging by estimating d_e for egg‐to‐pupae viability (EPV) after a drastic reduction in population size in a recently captured base population of Drosophila melanogaster. For this purpose, we first obtained estimates of the inbreeding depression rate (δ) for EPV in the base population, and we found that about 40% was due to segregating recessive lethals. Then, two sets of lines were founded from this base population and were maintained with different effective size throughout the rest of the experiment (N = 6; N = 12), their mean EPV being assayed at different generations. Due to purging, the reductions in mean EPV experienced by these lines were considerably smaller than the corresponding neutral predictions. For the 60% of δ attributable to nonlethal deleterious alleles, our results suggest an effective purging coefficient d_e > 0.02. Similarly, we obtain that d_e > 0.09 is required to roughly account for purging against the pooled inbreeding depression from lethal and nonlethal deleterious alleles. This implies that purging should be efficient for population sizes of the order of a few tens and larger, but might be inefficient against nonlethal deleterious alleles in smaller populations.     

EFFECTS OF POPULATION BOTTLENECKS ON GENETIC DIVERSITY AS MEASURED BY ALLOZYME ELECTROPHORESIS

Low levels of allozyme heterozygosity in populations are often attributed to previous population bottlenecks; however, few experiments have examined the relationship between heterozygosity and bottlenecks under natural conditions. The composition and number of founders of 55 experimental populations of the eastern mosquitofish (Gambusia holbrooki), maintained under simulated field conditions, were manipulated to examine the effects of bottlenecks on three components of allozyme diversity. Correlations between observed and expected values of allozyme heterozygosity, proportions of polymorphic loci, and numbers of alleles per locus were 0.423, 0.602, and 0.772, respectively. The numbers of polymorphic loci and of alleles per locus were more sensitive indicators of differences in genetic diversity between the pre-bottleneck and post-bottleneck populations than was multiple-locus heterozygosity. In many populations, single- and multiple-locus heterozygosity actually increased as a result of the founder event. The weak relationship between a population's heterozygosity and the number and composition of its founders resulted from an increase in the variance of heterozygosity due to drift of allele frequencies. There was little evidence that selection influenced the loss of allozyme variation. When it is not possible to estimate heterozygosity at a large number of polymorphic loci, allozyme surveys attempting to detect founder events and other types of bottlenecks should focus on levels of locus polymorphism and allelic diversity rather than on heterozygosity.     

SMALL POPULATION GENETIC VARIABILITY AT LOCI UNDER STABILIZING SELECTION

Genetic variability at a locus under stabilizing selection in a finite population is investigated using analytic methods and computer simulations. Three measures are examined: the number of alleles k, heterozygosity H, and additive genetic variance Vg. A nearly-neutral theory results. The composite parameter S = NV_M/V_s (where N is the population size, V_M the variance of new mutant allelic effects and V_s the weakness of stabilizing selection) figures prominently in the results. The equilibrium heterozygosity is similar to that of strictly neutral theory, H = 4N_μc/ (1 + 4N_μc), except that μ_c = where c is about 0.5. Simulations corroborate except for very low N. Genetic variability attains similar equilibrium values at both a “lone” locus and at an “embedded” locus. This agrees with my earlier work concerning molecular clock rates. These results modify the neutralist interpretation of data concerning genetic variability and genetic distances between populations. Low H values are proportional not to N but to . This may explain the narrow observed range of H among species. Heterozygosities need not be highly correlated to genetic variances. Genetic variances are not highly dependent on population size except in very small populations which are difficult to sample without bias because the smallest populations go extinct the fastest. Nearly neutral evolution will not be easily distinguished from strictly neutral theory under the Hudson-Kreitman-Aguade inter-/intraspecific variation ratio test, since a similar effective mutation rate holds for genetic distances and D = 2μ_ct, where . As with strictly neutral theory, comparisons across loci should show D and H to be positively correlated because of the shared μ_c. But unlike neutral theory, for a given locus, comparisons across species should show D and H to be negatively correlated. There is no obvious threshold of population size below which genetic variability inevitably declines. Extinction depends on both genetic variation and natural selection. Neither theory nor observation presently indicates the measure of genetic variability (k, H, V_G or other) that best indicates vulnerability of a small population to extinction.