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.

A Population Genetics-Phylogenetics Approach to Inferring Natural Selection in Coding Sequences

Bayesian Inference of Natural Selection from Allele Frequency Time Series

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

A novel spectral method for inferring general diploid selection from time series genetic data

Inferring demographic and selective histories from population genomic data using a two-step approach in species with coding-sparse genomes: an application to human data

A Genomic Map of the Effects of Linked Selection in Drosophila

Population Genomics on the Fly: Recent Advances in Drosophila

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

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

Adaptive Genic Evolution in the Drosophila Genomes

Inference of Natural Selection from Interspersed Genomic Elements Based on Polymorphism and Divergence

Inferring biophysical models of evolution from genome-wide patterns of codon usage

A path integral approach for allele frequency dynamics under polygenic selection

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion

The influence of GC -biased gene conversion on nonadaptive sequence evolution in short introns of Drosophila melanogaster

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

Accurate inference of population history in the presence of background selection

Detecting Directional Selection from the Polymorphism Frequency Spectrum

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

Coalescence, genetic diversity in sexual populations under selection

Global Landscape of Recent Inferred Darwinian Selection for Homo Sapiens.