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.

Detecting compensatory covariation signals in protein evolution using reconstructed ancestral sequences

Synonymous and nonsynonymous distances help untangle convergent evolution and recombination

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

How pairwise coevolutionary models capture the collective residue variability in proteins

Accurate Detection of Convergent Amino-Acid Evolution with PCOC

Protein structure prediction from sequence variation

PHACE: Phylogeny-Aware Co-Evolution

Identifying Coevolution Between Amino Acid Residues in Protein Families: Advances in the Improvement and Evaluation of Correlated Mutation Algorithms

What Does Evolution Tell Us About The Structure Of A Functional Amyloid Protein?

From principal component to direct coupling analysis of coevolution in proteins: Low-eigenvalue modes are needed for structure prediction

On the Accuracy of Inferring Energetic Coupling Between Distant Sites in Protein Families from Evolutionary Imprints: Illustrations Using Lattice Model.

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

Parsimonious evolutionary scenario for the origin of allostery and coevolution patterns in proteins

Small Coupling Expansion for Multiple Sequence Alignment

A Novel Algorithm for Detecting Co-Evolutionary Domains in Protein and Nucleotide Sequences.

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

The Roles of Compensatory Evolution and Constraint in Aminoacyl tRNA Synthetase Evolution

Identifying the Genetic Basis of Functional Protein Evolution Using Reconstructed Ancestors

Simultaneous identification of specifically interacting paralogs and inter-protein contacts by Direct-Coupling Analysis

A Maximum Likelihood Method for Detecting Functional Divergence at Individual Codon Sites, with Application to Gene Family Evolution

Shifts in amino acid preferences as proteins evolve: A synthesis of experimental and theoretical work