Abstract:Abstract Frequencies of synonymous codons are typically non-uniform, despite the fact that such codons correspond to the same amino acid in the genetic code. This phenomenon, known as codon bias, is broadly believed to be due to a combination of factors including genetic drift, mutational biases, and selection for speed and accuracy of codon translation; however, quantitative modeling of codon bias has been elusive. We have developed a biophysical population genetics model which explains genome-wide codon frequencies observed across 20 organisms. We assume that codons evolve independently of each other under the influence of mutation and selection forces, and that the codon population has reached evolutionary steady state. Our model implements codon-level treatment of mutations with transition/transversion biases, and includes two contributions to codon fitness which describe codon translation speed and accuracy. Furthermore, our model includes wobble pairing – the possibility of codon-anticodon base pairing mismatches at the 3’ nucleotide position of the codon. We find that the observed patterns of genome-wide codon usage are consistent with a strong selective penalty for mistranslated amino acids. Thus codons undergo purifying selection and their relative frequencies are affected in part by mutational robustness. We find that the dependence of codon fitness on translation speed is weaker on average compared to the strength of selection against mistranslation. Although no constraints on codon-anticodon pairing are imposed a priori , a reasonable hierarchy of pairing rates, which conforms to the wobble hypothesis and is consistent with available structural evidence, emerges spontaneously as a model prediction. Finally, treating the translation process explicitly in the context of a finite ribosomal pool has allowed us to estimate mutation rates per nucleotide directly from the coding sequences. Reminiscent of Drake’s observation that mutation rates are inversely correlated with the genome size, we predict that mutation rates are inversely proportional to the number of genes. Overall, our approach offers a unified biophysical and population genetics framework for studying codon bias across all domains of life.

CoLaML: Inferring latent evolutionary modes from heterogeneous gene content

Bayesian hidden Markov tree models for clustering genes with shared evolutionary history

Modeling gene content across a phylogeny to determine when genes become associated

Decoding coalescent hidden Markov models in linear time

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

Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes

A Probabilistic Model for Gene Content Evolution with Duplication, Loss, and Horizontal Transfer

Genomic language model predicts protein co-regulation and function

A Maximum Likelihood Method for Detecting Functional Divergence at Individual Codon Sites, with Application to Gene Family Evolution

Inference of genomic landscapes using ordered Hidden Markov Models with emission densities (oHMMed)

ParallelEvolCCM: Quantifying co-evolutionary patterns among genomic features

Lateral Gene Transfer, Rearrangement and Reconciliation

Hidden Markov Model Variants and their Application

Modelling gene content across a phylogeny to determine when genes become associated

MAAPE: A Modular Approach to Evolutionary Analysis of Protein Embeddings

Coupling hidden Markov models for the discovery of Cis-regulatory modules in multiple species

Bayesian Inference of Dependent Population Dynamics in Coalescent Models

Multimodal learning of noncoding variant effects using genome sequence and chromatin structure

CAFE 5 models variation in evolutionary rates among gene families

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

Linguistic mechanism of the evolution of amino acid frequencies and genomic GC content