A nonlinear method for estimating nucleotide polymorphism from restriction maps

Restriction mapping is used to estimate nucleotide sequence polymorphism when the regions to be studied are too long or too numerous to be sequenced. Restriction mapping is less costly than DNA sequencing, but it does not allow direct measurement of underlying nucleotide polymorphism. It is therefore useful to be able to estimate underlying nucleotide polymorphism from observations of polymorphism in restriction maps, as this offers some of the resolution afforded by DNA sequencing at a reduced cost. Previous estimators of underlying nucleotide polymorphism have assumed that each restriction-enzyme-binding site contains, at most, a single polymorphic nucleotide position (the low-polymorphism-frequency assumption), and this assumption has placed an upper limit on the level of polymorphism that can be resolved by these estimators. The present study documents an estimator which allows relaxation of this assumption. The new estimator more accurately estimates underlying nucleotide polymorphism when the polymorphism level is high enough to falsify the low-polymorphism-frequency assumption. The new estimator therefore yields good results for data sets that are too divergent for analysis by present methods.