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.

The McDonald-Kreitman Test and its Extensions under Frequent Adaptation: Problems and Solutions

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

An efficient and robust ABC approach to infer the rate and strength of adaptation

Conflicting signals of adaptive molecular evolution - Where does the neutral theory stand after 50 years?

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

Skepticism toward adaptive signals in DNA sequence comparisons - Is the neutral theory dead yet?

The Dynamics of Genetic Draft in Rapidly Adapting Populations

Genetic draft, selective interference, and population genetics of rapid adaptation

Two Decades of Suspect Evidence for Adaptive Molecular Evolution-Negative Selection Confounding Positive-Selection Signals

Inferring the model and onset of natural selection under varying population size from the site frequency spectrum and haplotype structure

A novel expectation-maximization approach to infer general diploid selection from time-series genetic data

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

Estimating the Parameters of Background Selection and Selective Sweeps in Drosophila in the Presence of Gene Conversion.

Demography-adjusted tests of neutrality based on genome-wide SNP data

Inferring population decline and expansion from microsatellite data: a simulation-based evaluation of the Msvar method

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

Testing for Ancient Selection Using Cross-population Allele Frequency Differentiation

Temporal challenges in detecting balancing selection from population genomic data

The impact of non-neutral synonymous mutations when inferring selection on non-synonymous mutations

Accurate inference of population history in the presence of background selection

Robust identification of local adaptation from allele frequencies