Using Maximum Likelihood to Estimate Population Size From Temporal Changes in Allele Frequencies

Corresponding author: Ellen G. Williamson, University of California, Integrative Biology, 3060 VLSB, Berkeley, CA 94720-3141., ellenw{at}socrates.berkeley.edu (E-mail)

Communicating editor: G. A. CHURCHILL

	ABSTRACT

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

We develop a maximum-likelihood framework for using temporal changes in allele frequencies to estimate the number of breeding individuals in a population. We use simulations to compare the performance of this estimator to an F-statistic estimator of variance effective population size. The maximum-likelihood estimator had a lower variance and smaller bias. Taking advantage of the likelihood framework, we extend the model to include exponential growth and show that temporal allele frequency data from three or more sampling events can be used to test for population growth.

BIOLOGISTS are often interested in the number of reproducing individuals in a population, or in the related quantity, the effective population size. If only small numbers of individuals are reproducing successfully, then loss of genetic diversity through genetic drift and increased homozygosity due to unavoidable inbreeding may harm the long-term potential of the population to survive (SOULE 1986 ).

The increase in homozygosity in small populations can be quantified by the increase in variance in allele frequencies among identical loci. In a Wright-Fisher model of diploid size N, the variance in allele frequencies among identical loci with initial allele frequency p is p(1 - p)/2N after one generation. For many other models of population structure, one can define a variance effective population size, N_e, such that after one generation the variance in allele frequency among initially identical loci is p(1 - p)/2N_e (EWENS 1979 ). This variance effective population size is related to the expected value of F, the standardized variance of change in allele frequency. The mathematical relationship between F and the variance effective population size led to the development of a statistic for estimating variance effective population size from temporal samples of allele frequency data (KRIMBAS and TSAKAS 1971 ; NEI and TAJIMA 1981 ; POLLAK 1983 ; WAPLES 1989A ). These methods have been applied to data taken from a wide variety of taxa in both natural and artificial (i.e., hatchery) populations (WAPLES 1990 ; WAPLES and TEEL 1990 ; HEDGECOCK et al. 1992 ; HUSBAND and BARRETT 1992A , HUSBAND and BARRETT 1992B ; JORDAN et al. 1992 ; LESSIOS et al. 1994 ; HEDRICK et al. 1995 ; BURCZYK 1996 ; JORDE and RYMAN 1996 ; RICHARDS and LEBERG 1996 ; MILLER and KAPUSCINSKI 1997 ; SCRIBNER et al. 1997 ; DALLAS et al. 1998 ; LAIKRE et al. 1998 ; LEHMANN et al. 1998 ).

With the advent of high-speed computers, maximum-likelihood methods for estimating population genetic parameters have become increasingly popular. Although these methods can be computationally intensive, maximum-likelihood estimators have many desirable statistical properties. One advantage of using a maximum-likelihood framework is that it is relatively simple to modify the underlying statistical model. In this study we develop a maximum-likelihood framework for estimating population size from temporal changes in allele frequencies. We then compare estimates of population size using maximum likelihood and the F-statistic approach on both real and simulated data. Finally we extend the model to include population growth and show that it is possible in principle to estimate both population size and growth rate simultaneously from temporal allele frequency data, although large numbers of loci are needed.

	THEORY

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

Our likelihood method is based on the Wright-Fisher model of neutral genetic evolution in an isolated population (EWENS 1979 ). The Wright-Fisher model assumes discrete, nonoverlapping generations and constant population size, N. Gametes are chosen randomly and with replacement from the previous generation. Sampling from the parental generation with replacement is equivalent to assuming that the gamete pool is infinitely large.

The likelihood of a discrete parameter value given some data is, by definition, the probability of observing the data given the parameter value (EDWARDS 1992 ). In the present case, the data consist of the counts of each allelic type at several unlinked neutral genetic markers in a random sample taken with replacement from a population at generations 0, t₁, t₂, ... t_f. The configuration of alleles at each locus can be represented by vectors ₀, _t1, _t2, ..._tf, where the ith element in ₀ is the number of alleles of type i. The likelihood of N, the population size, given the data, is equal to the probability of observing the data given N; L(N| ₀, _t1, _t2, ..._tf) = P(₀, _t1, _t2, ..._tf|N). This N is a parameter in the Wright-Fisher model and will not always be the same as the variance effective size estimated by the F-statistic method. However, these quantities are related; when a population is evolving according to the Wright-Fisher model, the variance effective size and N are the same.

In a Wright-Fisher population, the configuration of alleles in each generation has a multinomial probability that depends only on the allele frequencies of the previous generation. Similarly, the probability of the allelic configuration , a sample of size n taken randomly from the gamete pool of a population at a locus with k alleles that has allele frequencies , is also multinomial and is given by

(1)

The joint probability of observing the configuration ₀ in a sample from generation 0, and _t in a second sample taken t generations later depends on the allele frequencies at this locus in the population at generations 0 and t, ₀ and _t. Because the sampling events are statistically independent,

(2)

where P(₀|₀) and P(_t|_t) are multinomial probability mass functions given by Equation 1.

If sample data are from diallelic dominant markers such as RAPDs, then different equations must be used for P(₀|₀) and P(_t|_t). If the markers are in Hardy-Weinberg equilibrium, then one can use diploid phenotype frequencies from dominant diallelic markers instead of allele frequencies and

(3)

instead of Equation 1. Here p_a is the frequency of the recessive allele, n_T is the total number of individuals sampled, and n_a is the number of individuals in the sample that have the recessive phenotype.

Because the parameter ₀ is unknown and not directly of interest, we treat it as a "nuisance" parameter. Similarly, the unobserved random variable _t is unknown and we are not interested in estimating its value. If the joint probability distribution of ₀ and _t is known, then we can compute the marginal likelihood of N by summing P(₀, _t|₀, _t) over all possible values for (₀, _t) weighted by P(_t, ₀|N). The population size, N, determines the probability of _t given ₀, and by the definition of conditional probability, P(_t, ₀|N) = P(_t|₀, N) · P(₀|N). Thus the marginal likelihood of N, given the data, is

(4)

In general, P(₀|N) will be unknown. We assume that any starting frequency, ₀, is equally likely. We choose a uniform distribution for ₀ because no additional parameters need to be specified. We define C to be the number of possible allelic configurations for the population. In the diallelic case, C = N + 1. C_poly = N - 1 of these configurations are polymorphic. If ₀ is uniformly distributed, then P(₀) = . If we assume only polymorphic markers are sampled, then P(₀|N) = .

We now have equations for the probabilities P(₀|₀), P(_t|_t), and P(₀|N) of Equation 3. The remaining probability, P(_t|₀, N), can be computed using a forward transition matrix, M. The number of rows and columns in the matrix is the number of possible allele configurations, C. The population's allelic configurations can be enumerated from 1 to C. A population configuration can easily be converted into an allele frequency vector by dividing the elements of the configuration by N. Each element of M, m_ij, is the probability of a population going from the ith configuration in the parental generation to the jth configuration in the offspring generation. Because we are using the Wright-Fisher model, transition probabilities are multinomial. The m_ij are given by Equation 1 if the sample configuration is replaced by the jth population configuration and the allele frequency vector is replaced by the ith population configuration divided by N. The probability distribution of a population's allelic configuration can be described by a vector of length C, where the ith element in is the probability that a population has the ith configuration. To compute P(_t|₀, N), the transition matrix, M, must be raised to the power t, where t is the number of generations between sampling events. P(_t|₀, N) is equal to M^t multiplied on the left by a row vector ^T₀, representing a population with 100% probability of having initial allele frequencies ₀, and on the right by a column vector _t, representing a population that has allele frequencies _t in generation t,

(5)

.

This method can be easily extended to include more than two sampling times. For example, if there are three samples, then the joint probability of the samples given the population size is

(6)

where ₀, _t1, _t2 and ₀, _t1, _t2 represent sample allele counts and population allele frequencies, respectively, at generations 0, t₁, and t₂.

For loci that have k > 2 alleles, the number of possible configurations, C, is (N + k - 1)!/(N!(k - 1)!) of which C_poly = (N - 1)!/((N - k)!(k - 1)!) are polymorphic (FELLER 1950 ). Unfortunately, working with large transition matrices can be extremely computationally intensive. Because C and C_poly increase rapidly with k, we restrict the analysis in this article to diallelic loci. Likelihood methods for loci with k > 2 will probably require simulation approaches such as Monte Carlo integration or Markov Chain Monte Carlo to approximate the likelihood function (E. C. ANDERSON, unpublished results).

This likelihood function is slightly unusual in that the parameter to be estimated, N, is a discrete parameter that determines the number of summands in the computation of the likelihood (Equation 4 and Equation 6). Thus, it is not practical to differentiate the likelihood function with respect to N either analytically or numerically. The maximum-likelihood estimate of N must be determined on a case-by-case basis from the likelihood function; there seems to be no general equation for the maximum-likelihood estimate of N using temporal allele frequency data.

	ASSESSING THE MAXIMUM-LIKELIHOOD ESTIMATOR

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

Simulation methods:
We compared the maximum-likelihood estimator to the F-statistic estimator with simulation tests. As noted earlier, these two estimators do not always estimate the same quantity. The F-statistic method estimates the variance effective size and the maximum-likelihood approach described above estimates the parameter N of a Wright-Fisher population. In these simulation tests the two estimators estimated the same quantity because our simulated populations evolved according to the Wright-Fisher model. In each simulated replicate we selected starting frequencies for each locus from a uniform distribution. The simulated populations consisted of 50 alleles (25 diploid individuals). We sampled 100 alleles (50 diploid individuals) from the population, with replacement, according to three different sampling regimes (see Table 1). For each replicate we computed both maximum-likelihood and F-statistic estimates of the population size. The F-statistic was computed according to

(7)

from WAPLES 1989A (Equation 9), where p_0,A, p_t,A, p_0,a, and p_t,a are the allele frequencies of allele types A and a in the samples taken at generations 0 and t, respectively. The variance effective size was then computed according to

(8)

also from WAPLES 1989A (Equation 11), where n₀ and n_t are the total number of alleles sampled at generations 0 and t.

When there were three sampling events we used the harmonic mean of the estimates for the two intervals. This estimator was derived by pooling the two F-statistic estimates and using this pooled value to compute N. POLLAK 1983 developed similar estimators for population sampling without replacement (hypergeometric). In the examples we considered, the three estimators derived by POLLAK 1983 are the same as the harmonic mean estimator except larger by multiplicative constants. Because in our simulation tests the F-statistic estimator was biased upward, the Pollak estimators of N were less accurate than our harmonic mean estimator and are not shown.

In all replicates, the number of sampled individuals in each sampling event was twice as large as the size of the simulated population. When sample sizes are small, sampling error rather than genetic drift may be the primary reason for observed changes in allele frequencies (WAPLES 1989B ). We avoided this problem by assuming large samples could be taken from the population. Fortunately, it is frequently possible to sample large numbers of individuals even when the number of breeding individuals is low. Often, juveniles are sampled from populations with high juvenile mortality, and the number of breeding adults may be orders of magnitude lower than the number of juveniles. Thus, our sampling scheme is representative of those found in the literature (WAPLES 1990 ; WAPLES and TEEL 1990 ; HEDGECOCK et al. 1992 ; HUSBAND and BARRETT 1992A , HUSBAND and BARRETT 1992B ; JORDAN et al. 1992 ; LESSIOS et al. 1994 ; HEDRICK et al. 1995 ; BURCZYK 1996 ; JORDE and RYMAN 1996 ; RICHARDS and LEBERG 1996 ; MILLER and KAPUSCINSKI 1997 ; SCRIBNER et al. 1997 ; DALLAS et al. 1998 ; LAIKRE et al. 1998 ; LEHMANN et al. 1998 ).

In an application to real data rather than a simulation test, the F-statistic or the maximum-likelihood estimate can always be computed, but the estimate may be infinity or very large. In our simulation tests, to reduce running time, the program halted all searches for the maximum of the likelihood curve if the maximum-likelihood estimate of population size was determined to be >500. We chose 500 because it is an order of magnitude larger than the true population size, 50. To be fair to both estimators, we discarded all simulation replicates in which either estimator was >500. For each estimator, we computed the mean and standard deviation for the estimates of N using the replicates that were not discarded. As a separate statistic we recorded the number of times each estimator was >500. The number of discarded replicates was the union of these events. These simulation tests were computationally intensive, so we ran many more replicates for the cases in which there were 5, 10, or 15 loci as these are the most realistic values for applications of this method to real data.

Simulation results:
Both estimators tended to overestimate population size and had a high variance for samples with few loci. This upward bias and large variance of the estimators would have been much greater if we had not discarded replicates. However, in all of our simulation tests, the mean maximum-likelihood estimate of population size among the replicates that were not discarded was closer to the simulated population size than was the mean F-statistic estimate (Table 1). Similarly, the variance in the likelihood estimates among these replicates was smaller than the variance in F-statistic estimates (Table 1). Also, the F-statistic estimate was more than an order of magnitude larger than the true value about twice as often as the maximum-likelihood estimate (Table 1). The estimates were positively correlated and the F-statistic estimate was generally >500 whenever the maximum-likelihood estimator was >500. We had similar results when we ran simulation tests with initial allele frequencies drawn from beta distributions with varying parameters (not shown). Thus the maximum-likelihood estimate seems robust to violations of our assumption that allele frequencies are drawn initially from a uniform distribution.

As expected, increasing the number of loci or the number of sampled time-points reduced both the variance and the bias in both estimators (Table 1). When the number of markers was very large, on the order of 100 to 200 loci, then both estimators performed very well. The difference between the two estimators was more obvious with samples of fewer markers. Increasing the number of sampling times improved the estimates for both methods (Table 1). The dramatic improvement in both estimators with multiple sampling (Table 1) may be explained by the initial rapid increase in accuracy in both estimators, particularly the F-statistic estimator, as the number of sampled loci increases. Having samples covering two time intervals of the same length is roughly equivalent to sampling twice as many loci over one time interval.

	ANALYZING A REAL DATA SET

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

We chose data from MILLER and KAPUSCINSKI 1997 to use as an example of the maximum-likelihood approach to estimating population size. These authors typed seven polymorphic microsatellite loci using a historical collection of fish scales taken from a completely isolated population of Northern pike (Esox lucius). Of the seven loci, five had two alleles and two had three alleles. For our analysis in this article, we combined the two least common allelic classes together for the two loci with three alleles to artificially create a completely diallelic data set. The scales had been removed from individuals in the population in 1961, 1977, and 1993. These dates represent two intervals of approximately four generations each. The samples in 1961 and 1993 were taken with replacement to the population. The 1977 sample was taken without replacement, although probably after reproduction for at least some of the sampled individuals. Miller and Kapuscinski used the F-statistic method to estimate variance effective population size two ways, first assuming the 1977 data were generated with replacement, then without. They concluded that the two estimates are not significantly different because the census size is orders of magnitude larger than the estimated number of breeders. In our analysis, for simplicity, we assumed that the samples in the Miller and Kapuscinski study were taken with replacement.

Our results agreed with those of Miller and Kapuscinski. Our maximum-likelihood estimate of population size based on their data was 46, with a support interval, based on a decrease of two in the log-likelihood, of 21 to 112 (Figure 1). In their article, Miller and Kapuscinski estimated variance effective population size using the same F-statistic method that we used for the simulation studies in this article. They estimated the variance effective population size to be 48 with a confidence interval of 19 to 101. The agreement between the two estimators for this data set suggests that the assumptions of the underlying models on which these statistics are based are probably well met by this population, as is discussed by Miller and Kapuscinski in their article.

View larger version (14K): In this window In a new window Download PPT slide   Figure 1. Likelihood function of population size using data from MILLER and KAPUSCINSKI 1997 .

	ESTIMATING POPULATION GROWTH RATE

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

One advantage of the likelihood framework for parameter estimation is the flexibility to modify the underlying model. Here, we modified the Wright-Fisher model to include exponential growth. In this model the population size at generation t, N_t, is the smallest integer less than or equal to r^tN₀, where r is the growth rate and N₀ is the initial population size. Still treating the initial allele frequency as a nuisance parameter, we computed the marginal likelihood of the growth parameter and the starting population size simultaneously. To compute the transition probabilities, the matrix M^t of Equation 5 was replaced by M₁, M₂, ... M_t, where M_i is a matrix with N_i columns and N_i-1 rows. Computing the likelihood of each combination of N₀ and r requires different matrices, M₁, M₂, ... M_t. As a result, this method is computationally intensive, even in the diallelic case, and so we were not able to do extensive simulation tests to evaluate this method.

As an example, we applied the method to the data from the pike population of MILLER and KAPUSCINSKI 1997 . The maximum-likelihood estimate was a growth rate of 120% per generation with an initial diploid population size of 25 in 1961 (Figure 1B). However, with the data available, the likelihood support interval is so large (entire graph in Figure 2) that it is impossible to determine if the population is growing, shrinking, or stable. The census size of this population, estimated from mark-recapture studies, fluctuated without any clear trend in the study interval between 1961 and 1993 (MILLER and KAPUSCINSKI 1997 , Table 4).

View larger version (50K): In this window In a new window Download PPT slide   Figure 2. Likelihood function of initial population size and exponential growth rate parameter using data from MILLER and KAPUSCINSKI 1997 . The shading scale for the log-likelihood is shown on the right.

To test the potential utility of our method for a much larger data set, we generated samples with 150 loci from a simulated population. Three samples were taken with replacement, four generations apart, from the simulated population. The simulated population grew at a rate of 120% per generation with an initial size of 25. These parameters are the maximum-likelihood estimates from the Miller and Kapuscinski data. The maximum-likelihood estimates from the simulated data were equal to the true parameter values (Figure 3). The support interval was much smaller with this many markers, and it could be unambiguously determined that the population was growing (Figure 3). Thus, maximum likelihood can be used to estimate growth rate from temporal changes in allele frequencies, but only with extensive data.

View larger version (25K): In this window In a new window Download PPT slide   Figure 3. Likelihood function of initial population size and exponential growth rate parameter using a simulated data set with 150 loci, initial population size of 25, and growth rate of 1.2. The shading scale for the log-likelihood is shown on the right.

	DISCUSSION AND CONCLUSIONS

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

The F-statistic method of estimating the variance effective population size is known to be biased upward and have a large variance (POLLAK 1983 ; WAPLES 1989A ). Estimates of the number of breeding individuals in a population may be needed to help guide management and conservation decisions for endangered species and populations. Therefore, any method that offers an increase in accuracy in estimating population size from genetic data should be seriously considered. The likelihood estimator for N appears to have a lower variance and smaller bias than the F-statistic estimator of the variance effective size (Table 1).

A disadvantage of using likelihood to estimate N is that it is more computationally intensive than using F-statistics. As a consequence, it is difficult to assess the performance of the likelihood method against the F-statistic approach over a wide range of parameters. Furthermore, even with dramatic improvements in computing speed, it would still be unrealistic to exactly compute the likelihood function for samples of highly polymorphic loci. One approach might be to bin allelic classes together and reduce the number of alleles observed to two or three. Another more desirable approach would be to estimate the likelihood function for the full multiallelic data set by using Monte Carlo simulation or Markov Chain Monte Carlo methods (E. C. ANDERSON, unpublished results). The appropriateness of either approach would have to be determined on a case-by-case basis and would depend on the number of loci, number of alleles, population size, accuracy desired, and limitations in terms of computational resources.

Both the F-statistic and the maximum-likelihood estimator assume that the study population is isolated and that the markers are selectively neutral and not mutating. Selection, migration, and mutation all have the power to change allele frequencies and would therefore change the estimated value of N using temporal allele frequency data. For example, a population receiving a constant stream of immigrants from a large source population every generation may have very stable allele frequencies in spite of a very small size. In this case, population size would be overestimated using a temporal method. On the other hand, a single rare large pulse of immigration from a population between sampling events could cause the population size to be underestimated. It is therefore important to carefully assess the study population to see if it meets these criteria before using temporal changes in allele frequencies to estimate population size or growth rate.

F-statistic estimates of variance effective size are generally improved more dramatically if more loci are sampled rather than more alleles. This is in part due to the fact that very rare alleles can be lost within a generation or two, while intervals between sampling events are often longer than this (LAIKRE et al. 1998 ; LUIKART et al. 1998 ). Also much of the computational difficulty in estimating N with likelihood derives from the inclusion of multiple alleles. Thus, the best allocation of resources may be to gather samples from many diallelic loci. Diallelic dominant markers such as RAPDs can be easily incorporated into the likelihood model and can also be combined with other types of markers. Studies that use diallelic markers should benefit from the increase in accuracy offered by the maximum-likelihood approach developed here.

A primary advantage of the likelihood framework is its flexibility. For example, if one could generate from preexisting data a prior distribution for N, then one could extend the method developed here to a fully Bayesian approach and then estimate a posterior distribution for N. In this article we demonstrated how temporal allele frequency data could be used to estimate population growth rates and population size simultaneously. Temporal allele frequency data could potentially be used to estimate other parameters with likelihood as well. This would require a model relevant to the study population that can be specified with only a few parameters. For example, if we were interested in estimating variance in family size instead of growth rate, a fecundity distribution that could be parameterized with a single parameter would be ideal. From this fecundity distribution, a transition matrix for population allele configurations would need to be derived. Thus, in addition to a more accurate estimate of the population size, the likelihood framework offers the flexibility to modify the working model and the potential to estimate new parameters.

	ACKNOWLEDGMENTS

This work was funded by National Institutes of Health grant GM40282 to M. Slatkin. We thank E. Anderson, J. Beilawski, G. Bertorelle, C. Muirhead, R. Neilsen, B. Rannala, S. Tavaré, and T. Turner for their comments and suggestions.

Manuscript received October 12, 1998; Accepted for publication February 22, 1999.

	LITERATURE CITED

TOP ABSTRACT THEORY ASSESSING THE MAXIMUM-LIKELIHOOD... ANALYZING A REAL DATA... ESTIMATING POPULATION GROWTH... DISCUSSION AND CONCLUSIONS LITERATURE CITED

BURCZYK, J., 1996  Variance effective population size based on multilocus gamete frequencies in coniferous populations: an example of a Scots pine clonal seed orchard. Heredity 77:74-82.

DALLAS, J. F., F. BONHOMME, P. BOURSOT, J. BRITTON-DAVIDIAN, and V. BAUCHAU, 1998  Population genetic structure in a Robertsonian race of house mice: evidence from microsatellite polymorphism. Heredity 80:70-77.

EDWARDS, A. W. F., 1992 Likelihood, Expanded Edition. Johns Hopkins University Press, Baltimore.

EWENS, W. J., 1979 Mathematical Population Genetics. Springer-Verlag, New York.

FELLER, W., 1950 An Introduction to Probability Theory and Its Applications. John Wiley and Sons, New York.

HEDGECOCK, D., V. CHOW, and R. S. WAPLES, 1992  Effective population numbers of shellfish broodstocks estimated from temporal variance in allelic frequencies. Aquaculture 108:215-232.

HEDRICK, P. W., D. HEDGECOCK, and S. HAMELBERG, 1995  Effective population size in winter-run chinook salmon. Conserv. Biol. 9:615-624.

HUSBAND, B. C. and S. C. H. BARRETT, 1992a  Effective population size and genetic drift in tristylous Eichhornia paniculata (Pontederiaceae). Evolution 46:1875-1890.

HUSBAND, B. C. and S. C. H. BARRETT, 1992b  Genetic drift and the maintenance of the style length polymorphism in tristylous populations of Eichhornia paniculata (Pontederiaceae). Heredity 69:440-449.

JORDAN, W. C., A. F. YOUNGSON, D. W. HAY, and A. FERGUSON, 1992  Genetic protein variation in natural populations of Atlantic salmon (Salmo salar) in Scotland: temporal and spatial variation. Can. J. Fish. Aquat. Sci. 49:1863-1872.

JORDE, P. E. and N. RYMAN, 1996  Demographic genetics of brown trout (Salmo trutta) and estimation of effective population size from temporal change of allele frequencies. Genetics 143:1369-1381[Abstract].

KRIMBAS, C. B. and S. TSAKAS, 1971  The genetics of Dacus oleae. V. Changes of esterase polymorphism in a natural population following insecticide control: selection of drift? Evolution 25:454-460.

LAIKRE, L., P. E. JORDE, and N. RYMAN, 1998  Temporal change of mitochondrial DNA haplotype frequencies and female effective size in a brown trout (Salmo trutta) population. Evolution 52:910-915.

LEHMANN, T., W. A. HAWLEY, H. GREBERT, and F. H. COLLINS, 1998  The effective population size of Anopheles gambiae in Kenya: implications for population structure. Mol. Biol. Evol. 15:264-276[Abstract].

LESSIOS, H. A., J. R. WEINBERG, and V. R. STARCZAK, 1994  Temporal variation in populations of the marine isopod Excirolana: how stable are gene frequencies and morphology? Evolution 48:549-563.

LUIKART, G., W. B. SHERWIN, M. STEELE, and F. W. ALLENDORF, 1998  Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change. Mol. Ecol. 7:963-974[Medline].

MILLER, L. M. and A. R. KAPUSCINSKI, 1997  Historical analysis of genetic variation reveals low effective population size in a northern pike (Esox lucius) population. Genetics 147:1249-1258[Abstract].

NEI, M. and F. TAJIMA, 1981  Genetic drift and estimation of effective population size. Genetics 98:625-640[Abstract/Free Full Text].

POLLAK, E., 1983  A new method for estimating the effective population size from allele frequency changes. Genetics 104:531-548[Abstract/Free Full Text].

RICHARDS, C. and P. L. LEBERG, 1996  Temporal changes in allele frequencies and a population's history of severe bottlenecks. Conserv. Biol. 10:832-839.

SCRIBNER, K. T., J. W. ARNTZEN, and T. BURKE, 1997  Effective number of breeding adults in Bufo bufo estimated from age-specific variation at minisatellite loci. Mol. Ecol. 6:701-712[Medline].

SOULÉ, M. E. (Editor), 1986 Conservation Biology: The Science of Scarcity and Diversity. Sinauer Associates, Sunderland, MA.

WAPLES, R. S., 1989a  A generalized approach for estimating effective population size from temporal changes in allele frequency. Genetics 121:379-392[Abstract/Free Full Text].

WAPLES, R. S., 1989b  Temporal variation in allele frequencies: testing the right hypothesis. Evolution 43:1236-1251.

WAPLES, R. S., 1990  Temporal changes of allele frequency in Pacific salmon: implications for mixed-stock fishery analysis. Can. J. Fish. Aquat. Sci. 47:968-976.

WAPLES, R. S. and D. J. TEEL, 1990  Conservation genetics of Pacific salmon: I. Temporal changes in allele frequency. Conserv. Biol. 4:144-156.

This article has been cited by other articles:

				J. P. Bollback, T. L. York, and R. Nielsen Estimation of 2Nes From Temporal Allele Frequency Data Genetics, May 1, 2008; 179(1): 497 - 502. [Abstract] [Full Text] [PDF]

				P. E. Jorde and N. Ryman Unbiased Estimator for Genetic Drift and Effective Population Size Genetics, October 1, 2007; 177(2): 927 - 935. [Abstract] [Full Text] [PDF]

				G. T. Skalski Joint Estimation of Migration Rate and Effective Population Size Using the Island Model Genetics, October 1, 2007; 177(2): 1043 - 1057. [Abstract] [Full Text] [PDF]

				R. S. Waples and M. Yokota Temporal Estimates of Effective Population Size in Species With Overlapping Generations Genetics, January 1, 2007; 175(1): 219 - 233. [Abstract] [Full Text] [PDF]

				T. Kitakado, S. Kitada, Y. Obata, and H. Kishino Simultaneous Estimation of Mixing Rates and Genetic Drift Under Successive Sampling of Genetic Markers With Application to the Mud Crab (Scylla paramamosain) in Japan. Genetics, August 1, 2006; 173(4): 2063 - 2072. [Abstract] [Full Text] [PDF]

				J. Wang Estimation of effective population sizes from data on genetic markers Phil Trans R Soc B, July 29, 2005; 360(1459): 1395 - 1409. [Abstract] [Full Text] [PDF]

				J. Fernandez, J. Galindo, B. Fernandez, A. Perez-Figueroa, A. Caballero, and E. Rolan-Alvarez Genetic Differentiation and Estimation of Effective Population Size and Migration Rates in Two Sympatric Ecotypes of the Marine Snail Littorina saxatilis J. Hered., July 1, 2005; 96(4): 460 - 464. [Abstract] [Full Text] [PDF]

				E. C. Anderson An Efficient Monte Carlo Method for Estimating Ne From Temporally Spaced Samples Using a Coalescent-Based Likelihood Genetics, June 1, 2005; 170(2): 955 - 967. [Abstract] [Full Text] [PDF]

				L. Excoffier, A. Estoup, and J.-M. Cornuet Bayesian Analysis of an Admixture Model With Mutations and Arbitrarily Linked Markers Genetics, March 1, 2005; 169(3): 1727 - 1738. [Abstract] [Full Text] [PDF]

				I. Goldringer and T. Bataillon On the Distribution of Temporal Variations in Allele Frequency: Consequences for the Estimation of Effective Population Size and the Detection of Loci Undergoing Selection Genetics, September 1, 2004; 168(1): 563 - 568. [Abstract] [Full Text] [PDF]

				D. A. Tallmon, G. Luikart, and M. A. Beaumont Comparative Evaluation of a New Effective Population Size Estimator Based on Approximate Bayesian Computation Genetics, June 1, 2004; 167(2): 977 - 988. [Abstract] [Full Text] [PDF]

				M. A. Beaumont Estimation of Population Growth or Decline in Genetically Monitored Populations Genetics, July 1, 2003; 164(3): 1139 - 1160. [Abstract] [Full Text] [PDF]

				G. Laval, M. SanCristobal, and C. Chevalet Maximum-Likelihood and Markov Chain Monte Carlo Approaches to Estimate Inbreeding and Effective Size From Allele Frequency Changes Genetics, July 1, 2003; 164(3): 1189 - 1204. [Abstract] [Full Text] [PDF]

				T. F. Turner, J. P. Wares, and J. R. Gold Genetic Effective Size Is Three Orders of Magnitude Smaller Than Adult Census Size in an Abundant, Estuarine-Dependent Marine Fish (Sciaenops ocellatus) Genetics, November 1, 2002; 162(3): 1329 - 1339. [Abstract] [Full Text] [PDF]

				P. Berthier, M. A. Beaumont, J.-M. Cornuet, and G. Luikart Likelihood-Based Estimation of the Effective Population Size Using Temporal Changes in Allele Frequencies: A Genealogical Approach Genetics, February 1, 2002; 160(2): 741 - 751. [Abstract] [Full Text] [PDF]

				R. Vitalis and D. Couvet Estimation of Effective Population Size and Migration Rate From One- and Two-Locus Identity Measures Genetics, February 1, 2001; 157(2): 911 - 925. [Abstract] [Full Text]

				S. Kitada, T. Hayashi, and H. Kishino Empirical Bayes Procedure for Estimating Genetic Distance Between Populations and Effective Population Size Genetics, December 1, 2000; 156(4): 2063 - 2079. [Abstract] [Full Text]

				E. C. Anderson, E. G. Williamson, and E. A. Thompson Monte Carlo Evaluation of the Likelihood for Ne From Temporally Spaced Samples Genetics, December 1, 2000; 156(4): 2109 - 2118. [Abstract] [Full Text]