Estimation of allele frequencies from color-multiplexed electropherograms

Parametric model fitting of unprocessed sequencing-gel trace data and a least-squares optimization algorithm provide a method for accurately determiniug allele freqnenties of single nucleotide substitutions in a population. The method uses trace data corn two homozygous individuals and from either a heterozygous individual or a mixed population of templates. A parametric model is fit to each of the traces to estimate the amount of each of the four flnorescent. dyes that is present at each site. The parameters estimated from each trace are then normalized to account, for scalar variations due to differences in the amount of sample loaded. The parameters estimated from the trace of the heterozygous individual or from the mixture are viewed as a weighted sum of the parameters estimated from the traces of the homozygous individuals. The weights, or allele freqnenties, are estimated by mhhizing the sum of squared errors between the linear combination of homozygous traces and the mixed trace. Comparison of allele frequencies estimated by our method to known frequencies at polymorphic sites in three pools of CEPH individuals show that our method is accurate. Our method is automatic and much less laborintensive than previous approaches.