Accurate and fast methods to estimate the population mutation rate from error prone sequences.

Area of Science:

• Genetics
• Ecology
• Evolutionary Biology

Background:

• Population mutation rate (theta) is a fundamental parameter in genetics, ecology, and evolutionary biology.
• Estimating theta accurately is challenging with error-prone data like expressed sequence tags and draft sequences.
• Random sequence errors manifest as singletons in population datasets.

Purpose of the Study:

• To develop methods for accurate estimation of the population mutation rate (theta) from error-prone sequence data.
• To mitigate the impact of random sequencing errors on genetic parameter estimation.
• To provide reliable tools for analyzing complex population datasets.

Main Methods:

• Implemented an infinite sites model focusing on internal branches of sample genealogy.
• Derived new Watterson and Tajima estimators for theta, specifically for error-prone sequences.
• Modified a maximum-likelihood model to incorporate experimental error and design factors alongside coalescence and mutation.

Main Results:

• Developed three new, accurate, and fast methods for theta estimation from error-prone sequences.
• The new methods independently derived Watterson and Tajima estimators.
• Validated methods using evolutionary simulations and a real dataset from the California seahare.

Conclusions:

• Recommend the use of the three new methods for determining theta from error-prone sequences.
• Advocate for the new maximum likelihood model as a foundation for future coalescent/mutation models.
• Emphasize the importance of accounting for experimental error in genetic analyses.