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.

MODELTEST: testing the model of DNA substitution

Statistical tests of models of DNA substitution

DnaSAM: Software to perform neutrality testing for large datasets with complex null models

A new hierarchy of phylogenetic models consistent with heterogeneous substitution rates

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

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

bModelTest: Bayesian phylogenetic site model averaging and model comparison

Simulation And Analysis Of In Vitro Dna Evolution

Uncovering differential tolerance to deletions versus substitutions with a protein language model

data to detect putative modifications to DNA bases Modeling kinetic rate variation in third generation DNA sequencing

DnaSP, DNA polymorphism analyses by the coalescent and other methods

Markov Model Applied to Gene Evolution

Modeling DNA methylation dynamics with approaches from phylogenetics

Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes

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

SENCA: A Multilayered Codon Model to Study the Origins and Dynamics of Codon Usage

Detecting DNA Modifications from SMRT Sequencing Data by Modeling Sequence Context Dependence of Polymerase Kinetic.

Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

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

Modeling Kinetic Rate Variation in Third Generation DNA Sequencing Data to Detect Putative Modifications to DNA Bases

A Simple, General Result for the Variance of Substitution Number in Molecular Evolution