A generalized mechanistic codon model.

Area of Science:

• Computational Biology
• Evolutionary Biology
• Genomics

Background:

• Codon evolution models are crucial for understanding sequence evolution and detecting selection forces.
• Existing models often require numerous parameters, limiting their efficiency and interpretability.

Purpose of the Study:

• To develop a novel, mechanistic model for codon evolution using Kronecker products.
• To reduce the number of parameters required for estimating codon substitution rates.
• To improve the fit and interpretability of codon evolution models.

Main Methods:

• A new codon substitution rate matrix was constructed using the Kronecker product of nucleotide substitution matrices.
• The model incorporates equilibrium codon frequencies and a selection rate parameter.
• Parameters were optimized using maximum likelihood estimation.
• Model performance was evaluated using computer simulations and empirical data sets, comparing AICc values.

Main Results:

• The proposed model estimates the codon substitution matrix with significantly fewer parameters (19) compared to traditional methods (3,721).
• The model demonstrates a better fit to biological data across various datasets compared to existing codon evolution models.
• The parameters derived from the model are interpretable, similar to exchangeability rates in empirical models.

Conclusions:

• The Kronecker product-based model offers a more parsimonious and accurate approach to modeling codon evolution.
• This novel method enhances the ability to detect selection forces and understand sequence evolution.
• The model's improved fit and interpretability make it a valuable tool for evolutionary and genomic studies.