Estimating genetic variability with restriction endonucleases.

The estimation of the amount of sequence variation in samples of homologous DNA segments is considered. The data are assumed to have been obtained by restriction endonuclease digestion of the segments, from which the numbers and frequencies of the cleavage sites in the sample are determined. An estimator, p, of the proportion of sites that are polymorphic in the sample is derived without assuming any particular population genetic model for the evolution of the population. The estimator is very close to the EWENS, SPIELMAN and HARRIS (1981) estimator that was derived with the symmetric WRIGHT-FISHER neutral mode. ENGELS (1981) has also recently proposed an estimator of the same quantity, and he arrived at his estimator without assuming a particular population genetic model. The sampling variance of p and ENGELS' estimator are derived. It is found that the sampling variance of p is lower than the sampling variance of ENGELS' estimator. Also, the sampling variance of theta, an estimate of theta (=4Nu) is obtained for the symmetric WRIGHT-FISHER neutral model with free recombination and with no recombination.