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.

Two decades of suspect evidence for adaptive DNA-sequence evolution – Less negative selection misconstrued as positive selection

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

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

Genome wide signals of pervasive positive selection in human evolution

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

Variation In The Analysis Of Positively Selected Sites Using Nonsynonymous/Synonymous Rate Ratios: An Example Using Influenza Virus

Adaptive Genic Evolution in the Drosophila Genomes

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

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

On the unfounded enthusiasm for soft selective sweeps II: Examining recent evidence from humans, flies, and viruses

Rapid Evolution of Male Reproductive Genes in the Descent of Man

The Genome-wide Signature of Short-term Temporal Selection

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

A Chronological Atlas of Natural Selection in the Human Genome during the Past Half-million Years

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

Evolutionary Insights from a Large-Scale Survey of Population-Genomic Variation

Weak Selection Revealed by the Whole-Genome Comparison of the X Chromosome and Autosomes of Human and Chimpanzee

No evidence that natural selection has been less effective at removing deleterious mutations in Europeans than in West Africans

The Origin of Additive Genetic Variance Driven by Positive Selection.

Global Landscape of Recent Inferred Darwinian Selection for Homo Sapiens.

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