Abstract:Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

Conflicting signals of adaptive molecular evolution - Where does the neutral theory stand after 50 years?

Skepticism toward adaptive signals in DNA sequence comparisons - Is the neutral theory dead yet?

Two Decades of Suspect Evidence for Adaptive Molecular Evolution-Negative Selection Confounding Positive-Selection Signals

Two decades of suspect evidence for adaptive DNA-sequence evolution – Less negative selection misconstrued as positive selection

From Drift to Draft: How Much Do Beneficial Mutations Actually Contribute to Predictions of Ohta’s Slightly Deleterious Model of Molecular Evolution?

Genome wide signals of pervasive positive selection in human evolution

Estimating the rate of adaptive molecular evolution in the presence of slightly deleterious mutations and population size change

Adaptive Genic Evolution in the Drosophila Genomes

The McDonald-Kreitman Test and its Extensions under Frequent Adaptation: Problems and Solutions

Adaptive protein evolution at the Adh locus in Drosophila

Detecting Natural Selection at the DNA Level

Towards an improved understanding of molecular evolution: the relative roles of selection, drift, and everything in between

A novel test for selection on cis-regulatory elements reveals positive and negative selection acting on mammalian transcriptional enhancers

Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

L-Shaped Distributions of the Relative Substitution Rates (c/µ) in SARS-COV-2 with or without Molecular Clocks, Challenging Mainstream Evolutionary Theories

Beginner's Guide on the Use of PAML to Detect Positive Selection

Detecting Positively Selected Sites from Amino Acid Sequences: an Implicit Codon Model

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Applying ecological models to communities of genetic elements: the case of neutral theory

Does Positive Selection Drive Transcription Factor Binding Site Turnover? A Test with Drosophila Cis-Regulatory Modules.

Neutrality Tests Using Dna Polymorphism From Multiple Samples