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.

Site-specific structure and stability constrained substitution models improve phylogenetic inference

Selection among site-dependent structurally constrained substitution models of protein evolution by approximate bayesian computation

A new formulation of protein evolutionary models that account for structural constraints

Substitution models of protein evolution with selection on enzymatic activity

A general substitution matrix for structural phylogenetics.

Exploring evolution to uncover insights into protein mutational stability

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

Robustness of ancestral sequence reconstruction to among-site evolutionary heterogeneity and epistasis

Seq-InSite: sequence supersedes structure for protein interaction site prediction

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

Accurate prediction of site- and amino-acid substitution rates with a mutation-selection model

Improved prediction of site-rates from structure with averaging across homologs

STRUM: structure-based prediction of protein stability changes upon single-point mutation

The Combined Effects Of Amino Acid Substitutions And Indels On The Evolution Of Structure Within Protein Families

Newly developed structure-based methods do not outperform standard sequence-based methods for large-scale phylogenomics

Improved prediction of site‐rates from structure with averaging across homologs

Enhancing predictions of protein stability changes induced by single mutations using MSA-based language models

Selection on protein structure, interaction, and sequence

Structure-based Prediction of the Effects of a Missense Variant on Protein Stability.

A mechanistic stress model of protein evolution accounts for site-specific evolutionary rates and their relationship with packing density and flexibility

The divergence time of protein structures modelled by Markov matrices and its relation to the divergence of sequences