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.

Natural Selection on Various Sites of Ribosomal Proteins： a Cladistic View

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

Natural Selection of Protein Structural and Functional Properties: a Single Nucleotide Polymorphism Perspective

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

Selection on protein structure, interaction, and sequence

Evolutionary Selection and Functional Divergence of MT Gene in Mammalian

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

Adaptive evolution and functional significance of the PPARGC1A gene across diverse animal species

The effectiveness of selection in a species affects the direction of amino acid frequency evolution

Codon-Substitution Simulation Models for Detecting Molecular Adaptive Evolution

Codon clusters with biased synonymous codon usage represent hidden functional domains in protein-coding DNA sequences

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

Detection of differential selection pressure and functional-specific sites in subunits of vertebrate neuronal nicotinic acetylcholine receptors.

Recent codon preference reversals in the lineage

Selective pressure on synonymous codon usage in mammalian protein‐encoding genes

Synonymous codon usage and selection on proteins

Detecting Recent Positive Selection with High Accuracy and Reliability by Conditional Coalescent Tree.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

A New Method for Detecting Natural Selection at the Level of Nucleotide Sites

Analysis of the Functional Differentiation Sites in Vertebrate Neuronal Nicotinic Acetylcholine Receptor Subunits

Natural Selection and Functional Potentials of Human Noncoding Elements Revealed by Analysis of Next Generation Sequencing Data.