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.

Modeling substitution and indel processes for AFLP marker evolution and phylogenetic inference

Development of Genome-Wide InDel Markers and Their Integration with SSR, DArT and SNP Markers in Single Barley Map.

mInDel: a high-throughput and efficient pipeline for genome-wide InDel marker development

Random-effects substitution models for phylogenetics via scalable gradient approximations

Accurate reconstruction of insertion-deletion histories by statistical phylogenetics

Single-character insertion-deletion model preserves long indels in ancestral sequence reconstruction

Infinite Mixture Models for Improved Modeling of Across-Site Evolutionary Variation

Optimal Phylogenetic Reconstruction of Insertion and Deletion Events

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment including Long Indels

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

IDSSR: an Efficient Pipeline for Identifying Polymorphic Microsatellites from a Single Genome Sequence

A new hierarchy of phylogenetic models consistent with heterogeneous substitution rates

FamLink2 - A comprehensive tool for likelihood computations in pedigrees analyses involving linked DNA markers accounting for genotype uncertainties

Application of Machine Learning for Ancestry Inference Using Multi-Indel Markers.

A New Method for Modeling Coalescent Processes with Recombination

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution

Indels: computational methods, evolutionary dynamics, and biological applications

An Integer Programming Approach for the Selection of Tag SNPs Using Multi-allelic LD.