Back

Correcting for Global Synonymous Selection Improves the Accuracy of Episodic Positive Selection Inference

Verdonk, H. E.; Pivirotto, A.; Hey, J.; Kosakovsky Pond, S. L.

2026-06-06 bioinformatics
10.64898/2026.06.02.729680 bioRxiv
Show 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.

Matching journals

The top 5 journals account for 50% of the predicted probability mass.

1
Molecular Biology and Evolution
542 papers in training set
Top 0.3%
18.6%
2
PLOS Computational Biology
1863 papers in training set
Top 3%
12.0%
3
Systematic Biology
144 papers in training set
Top 0.2%
9.8%
4
BMC Evolutionary Biology
18 papers in training set
Top 0.1%
7.9%
5
PeerJ
308 papers in training set
Top 1%
4.4%
50% of probability mass above
6
Virus Evolution
155 papers in training set
Top 0.5%
3.5%
7
Bioinformatics
1204 papers in training set
Top 5%
3.2%
8
PLOS ONE
5266 papers in training set
Top 37%
3.2%
9
GENETICS
483 papers in training set
Top 2%
2.5%
10
BMC Bioinformatics
457 papers in training set
Top 3%
2.5%
11
Genome Biology and Evolution
338 papers in training set
Top 2%
2.5%
12
Nucleic Acids Research
1281 papers in training set
Top 9%
1.7%
13
Proceedings of the National Academy of Sciences
2444 papers in training set
Top 28%
1.7%
14
Nature Communications
5641 papers in training set
Top 46%
1.7%
15
eLife
5828 papers in training set
Top 54%
1.3%
16
Scientific Reports
3612 papers in training set
Top 61%
1.3%
17
NAR Genomics and Bioinformatics
242 papers in training set
Top 3%
1.1%
18
Journal of Computational Biology
48 papers in training set
Top 0.9%
1.1%
19
Cell Systems
201 papers in training set
Top 4%
1.0%
20
Bioinformatics Advances
203 papers in training set
Top 4%
1.0%
21
Nature Ecology & Evolution
18 papers in training set
Top 0.3%
1.0%
22
BMC Genomics
406 papers in training set
Top 7%
1.0%
23
Communications Biology
993 papers in training set
Top 26%
1.0%
24
PLOS Genetics
862 papers in training set
Top 12%
0.8%
25
Genome Biology
637 papers in training set
Top 8%
0.8%
26
BioData Mining
22 papers in training set
Top 1%
0.6%
27
Biophysical Journal
631 papers in training set
Top 5%
0.6%
28
Molecular Phylogenetics and Evolution
69 papers in training set
Top 0.9%
0.6%
29
Algorithms for Molecular Biology
17 papers in training set
Top 0.2%
0.6%
30
New Phytologist
346 papers in training set
Top 5%
0.6%