Correcting for Global Synonymous Selection Improves the Accuracy of Episodic Positive Selection Inference
Hannah VerdonkiD, Alyssa M. PivirottoiD, Jody HeyiD, Sergei L. Kosakovsky PondiD
Abstract
The ratio of nonsynonymous to synonymous substitution rates ({omega}) constitutes a fundamental parameter for inferring adaptive protein evolution, predicated upon the assumption that synonymous substitutions are selectively inert. This premise, however, is increasingly untenable given evidence of selection acting on synonymous substitutions, driven by various biological processes such as translational efficiency and mRNA stability. In this study, we demonstrate that unmodelled synonymous selection introduces substantial bias into{omega} estimation, resulting in elevated false positive rates in tests for positive selection. To rectify this, we present BUSTED+S+MSS, a statistical framework incorporating Multiclass Synonymous Substitution (MSS) models into BUSTED, a method for detecting episodic selection. By partitioning synonymous codons into empirically derived rate classes, this approach accounts for global synonymous constraints. Application to five diverse clades--Drosophila, Caenorhabditis, Enterobacteria, Saccharomyces, and Primates--reveals that the inclusion of MSS components consistently improves model fit and reduces the proportion of genes inferred to be under positive selection. In Enterobacteria, genes retaining significance under the corrected model exhibit weaker constraint on synonymous substitutions (dSs), consistent with the hypothesis that unmodelled purifying selection drives spurious signals of adaptation. Furthermore, an information-theoretic analysis indicates that whilst site-specific variation (SRV) provides the primary correction, global synonymous rate variation (MSS) contributes a distinct second-order correction. In highly divergent alignments, these signals act in concert to improve model fit. The BUSTED+S+MSS framework, especially when coupled with an "error-sink" to absorb alignment artifacts, thus offers a computationally feasible means to disentangle adaptive nonsynonymous substitution from the confounding effects of synonymous constraint.
02 /Specification & Classification
Click any tag to browse matching papers- Citations
- 0
- Organism
- Assay
- Research Focus
- Approach
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- A new comparative framework for estimating selection on synonymous substitutions.
Hannah Verdonk, Alyssa M. Pivirotto, … Sergei L. Kosakovsky Pond·Molecular Biology and EvolutionPeer-reviewed·Mar 23, 2025
- Minus the Error: Testing for Positive Selection in the Presence of Residual Alignment Errors
Avery Selberg, Nathan L Clark, … Sergei L. Kosakovsky Pond·bioRxiv·Nov 13, 2024
- A mutation-selection model of protein evolution under persistent positive selection
A. U. Tamuri, M. dos Reis·bioRxiv·May 18, 2021
- Evolutionary shortcuts via multi-nucleotide substitutions and their impact on natural selection analyses.
Alexander G. Lucaci, Jordan D. Zehr, … Sergei L. Kosakovsky Pond·Molecular Biology and EvolutionPeer-reviewed·Jun 27, 2023
- A codon model of nucleotide substitution with selection on synonymous codon usage
Laura Kubatko, Premal Shah, … Michael Gilchrist·bioRxiv·Aug 12, 2014
- Large-Scale Comparative Analysis of Codon Models Accounting for Protein and Nucleotide Selection
Iakov I. Davydov, Nicolas Salamin, Marc Robinson‐Rechavi·Molecular Biology and EvolutionPeer-reviewed·Mar 4, 2019