Abstract:Extracting biologically meaningful information from the continuing flood of genomic data is a major challenge in the life sciences. Codon usage bias (CUB) is a general feature of most genomes and is thought to reflect the effects of both natural selection for efficient translation and mutation bias. Here we present a mechanistically interpretable, Bayesian model (ribosome overhead costs Stochastic Evolutionary Model of Protein Production Rate [ROC SEMPPR]) to extract meaningful information from patterns of CUB within a genome. ROC SEMPPR is grounded in population genetics and allows us to separate the contributions of mutational biases and natural selection against translational inefficiency on a gene-by-gene and codon-by-codon basis. Until now, the primary disadvantage of similar approaches was the need for genome scale measurements of gene expression. Here, we demonstrate that it is possible to both extract accurate estimates of codon-specific mutation biases and translational efficiencies while simultaneously generating accurate estimates of gene expression, rather than requiring such information. We demonstrate the utility of ROC SEMPPR using the Saccharomyces cerevisiae S288c genome. When we compare our model fits with previous approaches we observe an exceptionally high agreement between estimates of both codon-specific parameters and gene expression levels ([Formula: see text] in all cases). We also observe strong agreement between our parameter estimates and those derived from alternative data sets. For example, our estimates of mutation bias and those from mutational accumulation experiments are highly correlated ([Formula: see text]). Our estimates of codon-specific translational inefficiencies and tRNA copy number-based estimates of ribosome pausing time ([Formula: see text]), and mRNA and ribosome profiling footprint-based estimates of gene expression ([Formula: see text]) are also highly correlated, thus supporting the hypothesis that selection against translational inefficiency is an important force driving the evolution of CUB. Surprisingly, we find that for particular amino acids, codon usage in highly expressed genes can still be largely driven by mutation bias and that failing to take mutation bias into account can lead to the misidentification of an amino acid's "optimal" codon. In conclusion, our method demonstrates that an enormous amount of biologically important information is encoded within genome scale patterns of codon usage, accessing this information does not require gene expression measurements, but instead carefully formulated biologically interpretable models.

Predicting Gene Expression Level from Codon Usage Bias

Predicting Gene Expression Level from Relative Codon Usage Bias: An Application to Escherichia coli Genome

Codon usage and amino acid usage influence genes expression level

Codon bias as a factor in regulating expression via translation rate in the human genome

Going Beyond Codon Adaptation Index

Estimating Gene Expression and Codon-Specific Translational Efficiencies, Mutation Biases, and Selection Coefficients from Genomic Data Alone

Codon usages in different gene classes of the Escherichia coli genome

Intragenic spatial patterns of codon usage bias in prokaryotic and eukaryotic genomes

Codon usage pattern and predicted gene expression in Arabidopsis thaliana

Transcriptome-wide meta-analysis of codon usage in Escherichia coli

Codon usage bias analysis for the coding sequences of Camellia sinensis and Brassica campestris

CodonO: a New Informatics Method for Measuring Synonymous Codon Usage Bias Within and Across Genomes.

SCUMBLE: a Method for Systematic and Accurate Detection of Codon Usage Bias by Maximum Likelihood Estimation

Relationships among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

Exploring individual codon influence on protein expression: a predictive approach

Selection on codon bias in yeast: a transcriptional hypothesis

Relative Codon Adaptation: A Generic Codon Bias Index for Prediction of Gene Expression

Variations in codon usage are not correlated with heterologous gene expression in Saccharomyces cerevisiae and Escherichia coli

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Gene length and codon usage bias in Drosophila melanogaster, Saccharomyces cerevisiae and Escherichia coli