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.

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

Polymorphism-aware models in RevBayes: Species trees, disentangling Balancing Selection and CG-biased gene conversion

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

Balancing selection and the functional effects of shared polymorphism in cryptic species

A nearly-neutral biallelic Moran model with biased mutation and linear and quadratic selection

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

How does plant chemodiversity evolve? Testing five hypotheses in one population genetic model

Automatic Discovery of Optimal Discrete Character Models

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

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

Background selection as null hypothesis in population genomics: insights and challenges from Drosophila studies

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

The best of both worlds: Combining population genetic and quantitative genetic models

Bridging the gap between the evolutionary dynamics and the molecular mechanisms of meiosis : a model based exploration of the intra-genomic Red Queen

A model-based approach for identifying signatures of balancing selection in genetic data

Temporal challenges in detecting balancing selection from population genomic data

Modelling survival and allele complementation in the evolution of genomes with polymorphic loci

A Coalescent Model for a Sweep of a Unique Standing Variant

Genomic evidence of rapid and stable adaptive oscillations over seasonal time scales in Drosophila

Direct benefits of genetic mosaicism and intraorganismal selection: modeling coevolution between a long-lived tree and a short-lived herbivore

Complex structural variation and behavioral interactions underpin a balanced sexual mimicry polymorphism