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.

Cost-Gain Analysis of Sequence Selection for Nonlinearity Mitigation

Sequence-Selection-Based Constellation Shaping for Nonlinear Channels

Perturbation-based Sequence Selection for Probabilistic Amplitude Shaping

Probabilistic Amplitude Shaping and Nonlinearity Tolerance: Analysis and Sequence Selection Method

Sequence optimisation for compressed sensing CDMA MIMO radar via mutual coherence minimisation

Two-phase optimization of binary sequences with low peak sidelobe level value

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

New Features or Metric on Sequence Comparison

Assessing Economic Viability: A Comparative Analysis of Total Cost of Ownership for Domain-Adapted Large Language Models versus State-of-the-art Counterparts in Chip Design Coding Assistance

SequenceLab: A Comprehensive Benchmark of Computational Methods for Comparing Genomic Sequences

Quadratic optimization for unimodular sequence synthesis and applications

Dual-Step Optimization for Binary Sequences with High Merit Factors

Fast differentiable DNA and protein sequence optimization for molecular design

Methods to Estimate Cryptic Sequence Complexity

Sequence Selection by Pareto Optimization

Optimal and Adaptive Pairwise DNA Sequence Correlation Analysis in Natural Gradient

Improved Small-Signal L2 Gain Analysis for Nonlinear Systems

Future Trends in the Design of Memetic Algorithms: the Case of the Linear Ordering Problem

Nonlinear optimization for a low‐emittance storage ring

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

On the Sequence Evaluation based on Stochastic Processes