CodABC: a computational framework to coestimate recombination, substitution, and molecular adaptation rates by...

Area of Science:

• Computational biology
• Molecular evolution
• Bioinformatics

Background:

• Estimating substitution and recombination rates is crucial for understanding molecular evolution of protein-coding sequences.
• Existing methods often rely on simplified models or struggle with joint estimation of multiple evolutionary parameters.

Discussion:

• CodABC is a novel computational framework designed for the joint estimation of recombination, substitution, and synonymous/nonsynonymous rates from coding sequence data.
• It employs approximate Bayesian computation (ABC) with and without regression adjustment, incorporating various codon models, intracodon recombination, and longitudinal sampling.
• The framework accommodates nuisance parameters like codon frequencies, transition matrices, and site heterogeneity, enhancing its applicability.

Key Insights:

• CodABC provides accurate joint parameter estimates for recombining coding sequences.
• It demonstrates superior performance compared to maximum-likelihood methods that utilize more approximate models.
• The flexibility in handling nuisance parameters increases the robustness of the estimates.

Outlook:

• CodABC offers a powerful tool for detailed molecular evolution studies.
• Its availability with a GUI, documentation, and parallel processing capabilities facilitates broader adoption in bioinformatics research.
• Future work may involve further refinement of codon models and integration with other evolutionary analyses.