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.

Estimating the proportion of beneficial mutations that are not adaptive in mammals

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

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

An estimate of fitness reduction from mutation accumulation in a mammal allows assessment of the consequences of relaxed selection

Beneficial mutation-selection balance and the effect of linkage on positive selection

Genomic mutation rates that neutralize adaptive evolution and natural selection

Rates and Fitness Consequences of New Mutations in Humans

An integrated model of the effects on fitness of beneficial, deleterious, and neutral mutations in sexual populations

Quantification of the effect of mutations using a global probability model of natural sequence variation

Evolving a favorable distribution for mutation effects

Mutation bias alters the distribution of fitness effects of mutations

In-vivo mutation rates and fitness landscape of HIV-1

The rate of beneficial mutations surfing on the wave of a range expansion

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

Mutations in yeast are deleterious on average regardless of the degree of adaptation to the testing environment

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

Accumulation of Spontaneous Mutations in the Ciliate Tetrahymena Thermophila

Genome wide signals of pervasive positive selection in human evolution

Growing Probability of Additional Offspring with a Beneficial Reversal Allele in the Decoupled Continuous-time Mutation-selection Model

Influence of Recombination and GC-biased Gene Conversion on the Adaptive and Nonadaptive Substitution Rate in Mammals versus Birds