# Ribo-seq enlightens codon usage bias

Damien Paulet, Alexandre David, Eric Rivals

► To cite this version:

Damien Paulet, Alexandre David, Eric Rivals. Ribo-seq enlightens codon usage bias. DNA Research, 2017, 24 (3), pp.303-210. 10.1093/dnares/dsw062. lirmm-01610054

HAL Id: lirmm-01610054

https://hal-lirmm.ccsd.cnrs.fr/lirmm-01610054

Submitted on 4 Oct 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Full Paper

# Ribo-seq enlightens codon usage bias

Damien Paulet¹,², Alexandre David³, and Eric Rivals¹,²*

¹Department of Computer Science, Laboratoire d'Informatique, de Robotique et de Microélectronique de Montpellier (LIRMM), CNRS et Université de Montpellier, 34095 Montpellier Cedex 5, France, ²Institut de Biologie Computationnelle (IBC); Université de Montpellier, France, and ³Cancer Biology Department, Institut de Génomique Fonctionnelle, CNRS, INSERM, Université de Montpellier, F-34094 Montpellier, France

*To whom correspondence should be addressed. Tel. +0033 4 67 41 86 64. Fax. +0033 4 67 41 85 00.
Email: rivals@lirmm.fr

Edited by Prof. Kenta Nakai

Received 20 June 2016; Editorial decision 14 December 2016; Accepted 15 December 2016

## Abstract

Codon usage is biased between lowly and highly expressed genes in a genome-specific manner. This universal bias has been well assessed in some unicellular species, but remains problematic to assess in more complex species. We propose a new method to compute codon usage bias based on genome wide translational data. A new technique based on sequencing of ribosome protected mRNA fragments (Ribo-seq) allowed us to rank genes and compute codon usage bias with high precision for a great variety of species, including mammals. Genes ranking using Ribo-Seq data confirms the influence of the tRNA pool on codon usage bias and shows a decreasing bias in multicellular species. Ribo-Seq analysis also makes possible to detect preferred codons without information on genes function.

**Key words:** codon usage, high throughput sequencing, synonymous codon, evolution, translation

## 1. Introduction

Usage of synonymous codons is biased in genomes, as some codons are favoured in highly expressed genes.¹,² This phenomenon, called codon usage bias (CUB), has been observed in most genomes, although it seems to vary according to the complexity of genomes.¹ The choice of codons seems to depend mainly on tRNA pool: the main hypothesis is that using codons that have the greatest number of accepting tRNAs increases the efficiency or accuracy of translation of highly expressed genes.³,⁴ In addition to its impact on mRNA translation, CUB has an impact on protein folding.⁵,⁶ Despite its importance, precise measures of CUB were, for years, difficult to obtain at genome wide scale. Sharp et al.⁷ developed a clustering approach from a set of sequenced genes: the RSCU for Relative Synonymous Codon Usage. This measure is inferred with a statistical process: two clusters of genes are determined based on expression levels, in order to maximize the difference in RSCU between the two groups, then the RSCU is computed using the codon frequencies within each

group. Consequently, RSCU accuracy relies on the number of sequenced genes. Here, we propose a new method that permits to directly compute CUB from all translated genes. A recent technique named ribosome profiling (or Ribo-seq) has revolutionized the analysis of translation and permitted to refine the picture of global gene expression control. Ribo-seq is based on deep sequencing of ribosome-protected mRNA fragments (RPF or “footprints”) and permits genome-wide analysis of translated sequences at nucleotide resolution. From a bioinformatic standpoint, Ribo-seq provides a new type of data: indeed, RPF reveal codon occupancy of active ribosomes.⁸ This technique yields a precise picture of translation by quantifying the number of ribosomes at every position in a genome. Researchers often analyse jointly Ribo-seq data and RNA-seq data and compare the number of mRNAs and the number of mRNAs involved in translation, which allows them to measure the efficiency of translation⁹,¹⁰,¹¹,¹²,¹³. We decided to use only Ribo-seq data to measure CUB, based on a simple rationale: the more observed
footprints, the more ribosomes at this position. This means that genes with a high coverage are highly translated, and mutations that lead to a benefit in terms of translation are more prone to be fixed in those genes: as a consequence, Ribo-seq could indicate which genes have a strong CUB. Ribo-seq is often used to measure translation efficiency, as the number of ribosomes per copy of an mRNA species. Here highly translated (which differs from efficiently translated) means that all copies of one mRNA species yield, all together, a large number of proteins.

We ranked genes according to their number of footprints in a precise and limited region of their mRNA (from the 20th to the 200th codons from the start codon). We then split genes into two groups according to their number of footprints and measured codon bias within each group: this way, our two groups are not optimized for RSCU, but derived from translation experiments. Moreover, contrary to previous approaches, we used high-throughput data reflecting the number of ribosomes per sequence and use all translated genes. Therefore, our measure of CUB is directly computed from translational data.

We applied our approach on a variety of genomes from a parasite, like *Plasmodium falciparum*, to multicellular eukaryotes like *Caenorhabditis elegans* and *Homo sapiens*. We corroborate Sharp's results on yeast and *C. elegans*,⁷,¹⁴ but also detect preferred codons in every species. A preferred codon of an amino acid is defined as the codon with the highest frequency in highly expressed genes. Our analysis provides novel insights on the evolution of CUB.

# 2. Materials and methods

## 2.1. ORF selection

Sequences and annotations were collected from NCBI for all organisms. Positions of start codons are crucial in this study. Therefore, we first kept all genes with a unique isoform. Genes with multiple isoforms were rejected if they had different start codons. If all isoforms share the same start codon, we kept the longest common coding sequence. We excluded the mRNAs of mitochondrial genes. Finally, we collected all the spliced mRNA sequences of the selected genes.

## 2.2. Ribo-seq data and mapping

All Ribo-seq data come from already published data (Supplementary Table S1). We chose species belonging to opisthokonts (the Fungi/ Metazoa group) for which Ribo-seq data were available and used *P. falciparum* as an outgroup because of its restricted set of tRNA genes compared to all other species considered in this study. When several datasets were available for one species, we selected datasets according to three conditions: 1/availability of a replicate, 2/ribosome profiling that closely follows the protocol of Ingolia et al.,⁸ and 3/a sufficient sequencing depth after quality filtering. For each selected dataset, we plot the footprint position with respect to the start codon for all genes and determined the majority peak near the start. For all experiments, the peak turns out to be at -12 nucleotides upstream from start codon. As explained in the study by Ingolia et al.,⁸ the “shift” applied to the positions of mapped reads was set to 12 nucleotides.

Reads were mapped with CRAC¹⁵ with a k-mer of length 25: this k-mer is longer than the recommended length for RNA-seq to better fit the short length of RPF. After mapping on mRNA, we took the first position of each read and applied a shift of 12 nucleotides to assign the read to the P-site position of ribosomes.¹⁶

## 2.3. Categorizing genes into lowly and highly translated genes

For each gene, we computed the coverage by ribosome footprints from the 20th to the 200th codons from the start codon (and at least 20 codons before the stop codon). We computed the mean coverage on the selected region after removal all genes having less than 10 footprints. Then, we ranked genes according to this mean coverage. From this ranking, we create two groups that were made uniform in terms of total amount of footprints: the lowly and highly translated genes. To accentuate the contrast between those two groups, we excluded the top 5% of lowly translated genes and the bottom 5% of highly translated genes.

## 2.4. RSCU computing

Finally, we computed the RSCU as defined by Sharp et al.⁷ within each group (lowly and highly translated genes). Instead of computing the bias of codon usage over the whole ORF, we took sequences from the 20th to the 200th codons from the start codon, i.e. the same region as the one used for gene sorting. For clarity, when the RSCU is computed based on Ribo-seq data, we named it RSCU<sub>RS</sub>. We developed a computer program in Java language to perform this computation (see the Availability section).

## 2.5. Copy number of tRNA genes and clustering

We use the number of tRNA genes of a species as a proxy for the size of its tRNA pool.¹⁷,¹⁸ We obtained the copy numbers of tRNA genes from the latest version (2.0 from 2016) of the Genomic tRNA Database (GtRNAbd)¹⁹ for all species, but *Histoplasma capsulatum*. For the latter, we used the Genbank annotation from NCBI database. To cluster species according to their tRNA copy number and according to their average RSCU<sub>RS</sub>, we perform hierarchical clustering with the R software using the Canberra distance with UPGMA algorithm (Unweighted Pair Group Method with Arithmetic Mean). Details are provided in Supplementary Material.

# 3. Results and discussion

Only amino acids that have multiple encoding codons are studied: only 59 codons were used as we excluded ATG (methionine), TGG (tryptophan), and stop codons.

## 3.1. CUB computation

Given a set of mRNA sequences, Sharp has proposed to compute the RSCU for each synonymous codon j of each amino acid i using the formula

$$ \frac{n_i \times x_{ij}}{\sum_{j=1}^{n_i} x_{ij}} $$

where $n_i$ denotes the number of synonymous codons for amino acid i, and $x_{ij}$ denotes the number of occurrences of codon j in the set of sequences. Sharp has selected a subset of mRNAs that were known to be highly expressed. Our proposal is to automatically select mRNA sequences that are highly translated according to their coverage by ribosomes, which is determined by mapping Ribo-seq reads. Then, we use the formula above with counts of codons limited to a range comprised between the 20th and the 200th codons. The choice of this range is explained below.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Comparison of the strength of codon usage bias as measured by the RSCU as originally defined by Sharp et al., and by the RSCU<sub>RS</sub> which is introduced in this work. Both measures are plotted for a subset of highly expressed genes. The RSCU is computed for each codon of an amino acid: it multiplies the relative observed frequency of this codon among all possible codons for the corresponding amino acid by the number of possible codons. In all graphs, the codons on x-axis are ordered by increasing RSCU values computed by Sharp et al. (A) Comparison of RSCU values of Sharp and RSCU<sub>RS</sub> values in *C. elegans* (in two replicates). Both replicate curves for RSCU<sub>RS</sub> are extremely close from each other, and they closely follow that from Sharp. (B) Comparison of RSCU values of Sharp and RSCU<sub>RS</sub> values in *C. albicans*. (C) and (D). Comparison of eight RSCU<sub>RS</sub> curves obtained with codons counts computed in different ranges of codons starting with 20th codon and ending between the 50th and the 400th codons, respectively, in *C. elegans* (C) and in *C. albicans* (D). Apart for the smallest range of codons (i.e. [20-50]) all curves are very close to each other, thereby showing the robustness of with respect to the range of codons taken into account. One observes that RSCU<sub>RS</sub> values reach higher values in the yeast species than in the worm species. Refer to the online version for colors.</figcaption></figure>

Sharp and Bradnam¹⁴ have computed RSCU for *C. elegans*, and we compared his results with our method (RSCU<sub>RS</sub>, for RSCU by Ribo-seq) as illustrated in Fig. 1A. The two methods split genes in two categories: highly and lowly expressed genes, but with completely different inputs: Sharp et al. computed codon bias with a statistical model based on gene sequences, while we used only Ribo-seq footprints. Nevertheless, we obtained a high correlation (>0.97) between Sharp's measure and ours. This correlation is even higher within replicates. In 1992, Lloyd and Sharp<sup>20</sup> published RSCU values for highly expressed genes in *Candida albicans*. Here again, the comparison with RSCU values is shown in Fig. 1B, and we observed a very high correlation between RSCU and RSCU<sub>RS</sub> (>0.98). In both experiments, we used a larger set of genes (Supplementary Table S1A).

We question whether our method, RSCU<sub>RS</sub>, is robust with respect to the range of codons used—by default [20, 200] codons. The lower limit of 20 avoids counting the accumulation of RPF due to translation initiation.⁸ The upper limit of 200 was chosen as a minimum length for including sufficient counts and to maximize the number of mRNAs taken into account. Clearly, the upper limit choice is somehow arbitrary. We computed the RSCU<sub>RS</sub> for eight different upper limits ranging from 50 to 400 codons. The RSCU<sub>RS</sub> curves for each range are shown on a single graph in Fig. 1B for *C. elegans* and

Fig. 1C for *C. albicans*. The agreement among all eight curves is striking for both species. Only the curves of two shortest ranges—[20-50] and [20,100]—depart slightly from the other curves, indicating that a minimum number of codons is necessary to capture a stable signal. All other curves are very close to each other, showing that the RSCU<sub>RS</sub> results are not strongly dependent of the upper limit. Hence, RSCU<sub>RS</sub> appears to be robust with respect to this parameter. Since the curve for range [20,200] is highly correlated with Sharp's curves, logically all eight curves are highly correlated with Sharp's RSCU values as well.

## 3.2. CUB exists in all species

As a first step, we computed the Euclidean distance of RSCU<sub>RS</sub> between highly and lowly translated genes (Fig. 2A) as a mean to evaluate the intensity of CUB for each organism: if this distance is close to 0, it means there are very few differences in codon usage between highly and lowly translated genes. Results gave two clearly distinct species groups. The first group corresponds to species with a lower CUB and comprises mammals, Drosophila, *P. falciparum*, and *H. capsulatum*. The second group of species exhibiting a higher CUB contains three yeasts and *C. elegans*.
<figure><img src="image_3.png" /></figure>

**Figure 2.** Differences in codon usage bias between lowly and highly translated genes. (A) Euclidean distance between RSCU<sub>RS</sub> of lowly and highly translated genes was computed for ten species. This distance clearly partitions the species in two groups: the mean euclidean distance is 1.51 for group 1 (left, in orange), and is 6.64 for group 2 (right, in grey). Group 2 comprises *C. albicans*, *S. cerevisiae*, *S. pombe* and *C. elegans*, while group 1 comprises all vertebrates, *D. melanogaster*, *P. falciparum* and *H. capsulatum*. (B) For each species, comparison of the RSCU<sub>RS</sub> of all codons between lowly and highly translated genes (ltg vs htg). Codons from group 1 species (containing the vertebrates) are shown as orange circles, while codons from group 2 species (that of the budding yeast) as gray circles. A point near the diagonal indicates a similar behavior in ltg and htg. The further apart the point from the diagonal, the higher the bias of that codon. Refer to the online version for colors.

For these species, we compare the RSCU<sub>RS</sub> of all codons between lowly and highly translated genes (ltg vs htg; Fig. 2B). One notices the asymmetry of the diagram with respect to the diagonal: The part of the diagram beyond a RSCU<sub>RS</sub> of two is populated above the diagonal and almost empty below it. This confirms the existence of favoured codons in htg, and their absence in ltg. Codons located far above the diagonal indicate strongly favoured codons in htg. This illustrates the impact of selection pressure linked to the gene expression level. There also exist codons far below the diagonal (which are disfavoured in htg), but their RSCU<sub>RS</sub> remains low (<2): those are either disfavoured or slightly favoured in ltg. Some codons even exhibit a nearly null RSCU<sub>RS</sub> in htg, meaning these are almost forbidden in htg. Apart from the global asymmetry, the difference between the codons of the two species groups is striking. All codons far apart from the diagonal belong to group 2 and very few codons of these species lie near the diagonal: in a species like yeast, the difference of usage between htg and ltg is strong for most codons. Note also that in group 2, some codons are disfavoured in htg and favoured in ltg, suggesting there could be some selection pressure also in ltg, albeit weaker than in htg. In species of group 1, most codons lie around the diagonal, but the bias changes with the RSCU<sub>RS</sub> value: most codons that are disfavoured in htg (dots lying below the horizontal line, i.e. with a RSCU<sub>RS</sub> in htg <1) lie below the diagonal, while most codons favoured in htg (i.e. a RSCU<sub>RS</sub> in htg > 1) lie above the diagonal. This suggests that codon usage preference in htg is small, but is detectable with Ribo-seq data, and thus do exist also in these species.

### 3.3. Patterns of codon preferences

Then, we compared the frequency of codons in highly translated genes. Higher eukaryotes have few strongly preferred codons while

yeasts show marked differences (Fig. 3A): some yeasts present a highly preferred codon for Arginine, Glutamic acid and Cysteine, whereas mammals and Drosophila do not. Glutamine follows remarkable patterns: both mammals and Drosophila strongly prefer CAG, yeasts favour CAA, while *H. capsulatum* and *C. elegans* use CAA/CAG in an equivalent fashion.

Our analysis revealed that some preferred codons are well conserved (Table 1): all species, but *P. falciparum*, have the same preferred codon for phenylalanine, histidine, tyrosine (that is, all aromatic amino acids that can be studied), and asparagine. On the opposite, Leucine has five different preferred codons. Indeed, unique preferred codons are rare: *Schizosaccharomyces pombe* has one only for proline (CCT), and *P. falciparum* has the highest number of unique features with seven unique preferred codons.

### 3.4. Preferred codons and tRNA copy number

It is assumed that a codon with the highest copy number of corresponding tRNAs is prone to be the preferred codon. This result is confirmed in all species, but *H. capsulatum* (Table 2). Nevertheless, species do not equally observe this tendency: in *S. pombe* 15 out 18 amino acids prefer the codon with the highest copy number of tRNA, while only 9-12 amino acids behave the same in mammals. Non-majority preferred codons are shared among species, like species preferring GGT for glycine (*S. pombe*, *Saccharomyces cerevisiae* and *C. albicans*) or preferring GCC for Alanine (*H. sapiens*, *Mus musculus*, *Rattus norvegicus*, *Drosophila melanogaster*, *C. elegans* and *H. capsulatum*).

We inferred two species trees based on the species' proximity either in terms of tRNA copy number, or in terms of average RSCU<sub>RS</sub> over highly translated genes (Fig. 3B and C). These trees are similar, but not equal. tRNA copy number groups the mammals with
<figure><img src="image_4.png" /><figcaption>Comparison of codons' relative frequencies</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Species clustering according to the tRNA gene copy numbers.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Species clustering according to the RSCU<sub>RS</sub></figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Codons lacking tRNA in function of tRNA gene copy numbers.</figcaption></figure>

**Figure 3.** (A) Preferred codons in highly translated genes. To ease the comparison of the bias across codons, we plotted the codon frequency such that each codon has a value comprised in [0,1]. The frequency is the RSCU multiplied by the number of codons for the corresponding amino acid. The range of colors goes from blue for a frequency of 0 until red for a frequency of 1. Same species abbreviations as in Fig. 1. (B) Clustering of the ten species based on tRNA copy numbers. (C) Clustering of the 10 species based on RSCU<sub>RS</sub> in highly translated genes. Both clustering were performed using UPGMA algorithm in R (see Methods section). (D) Variation of the numbers of codons lacking tRNA with the size of the species' tRNA pool. Linear regression lines are solid for codons with no tRNA (R<sup>2</sup> = 0.69) and dotted for preferred codons with no tRNA (R<sup>2</sup> = 0.53). *H. capsulatum* was excluded for both regressions. The regression line for all codons decreases sharply (slope = -0.026), while the regression line for preferred codons remains flat (slope = -0.002). Refer to the online version for colors.

Drosophila and *S. cerevisiae*, and makes a clade with *C. albicans*, *S. pombe* and *C. elegans*, while both *H. capsulatum* and *P. falciparum* form an outgroup (which reflects their small tRNA gene copy numbers; see Table 2). With CUB (Fig. 3C), the clade grouping the mammals and Drosophila remains opposed to the branches containing the yeast species and *P. falciparum*. The striking feature is the position of *H. capsulatum*, which, due to its low tRNA gene number, is located at the base the tRNA-tree, but it is grouped with the mammals in the CUB-tree. Despite its basal position in the tree, *H. capsulatum* is closer to the group of yeasts (Euclidian distance = 29.9) than to the group of mammals (Euclidian distance = 86.5).

### 3.5. Preferred codons, codons lacking tRNA and GC content

Wobbling corresponds to a codon that is decoded without following Watson-Crick base pair rules.²¹ Wobbling is complex to study, but codons lacking tRNA are necessarily subject to it. In Fig. 3D, we compared the numbers of codons lacking tRNA, and of preferred codons lacking tRNA for all species ordered by the size of their tRNA pool. As expected, when the tRNA pool increases, the number of codons lacking tRNA decreases steeply (slope equals -0.026642). However, the number of preferred codons lacking tRNA remains stable and appears not to depend on the size of the tRNA pool (slope dashed curve equals -0.002812).
<table><caption>**Table 1.** Preferred codons per organism in highly (A) and lowly (B) translated genes</caption><thead><tr><th></th><th>Hs</th><th>Mm</th><th>Rn</th><th>Dm</th><th>Hc</th><th>Ce</th><th>Sp</th><th>Sc</th><th>Ca</th><th>Pf</th></tr></thead><tbody><tr><td>Phe</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTT*</td></tr><tr><td>His</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAT</td></tr><tr><td>Tyr</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAT*</td></tr><tr><td>Asn</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAT*</td></tr><tr><td>Lys</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAA</td></tr><tr><td>Cys</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGT*</td><td>TGT*</td><td>TGT</td></tr><tr><td>Glu</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAA</td><td>GAA</td><td>GAA</td></tr><tr><td>Ala</td><td>GCC*</td><td>GCC*</td><td>GCC*</td><td>GCC*</td><td>GCC*</td><td>GCC*</td><td>GCT</td><td>GCT</td><td>GCT</td><td>GCT</td></tr><tr><td>Se</td><td>TCC*</td><td>TCC*</td><td>TCC*</td><td>TCC*</td><td>TCC*</td><td>TCC*</td><td>TCT</td><td>TCT</td><td>TCT</td><td>TCA</td></tr><tr><td>Thr</td><td>ACC*</td><td>ACC*</td><td>ACC*</td><td>ACC*</td><td>ACC*</td><td>ACC*</td><td>ACT</td><td>ACT</td><td>ACT</td><td>ACA</td></tr><tr><td>Ile</td><td>ATC*</td><td>ATC*</td><td>ATC*</td><td>ATC*</td><td>ATC*</td><td>ATC*</td><td>ATT</td><td>ATC*</td><td>ATT</td><td>ATT</td></tr><tr><td>Gly</td><td>GGC</td><td>GGC</td><td>GGC</td><td>GGC</td><td>GGC*</td><td>GGA*</td><td>GGT*</td><td>GGT*</td><td>GGT*</td><td>GGA</td></tr><tr><td>Pro</td><td>CCC</td><td>CCC*</td><td>CCC*</td><td>CCC*</td><td>CCC*</td><td>CCA</td><td>CCT</td><td>CCA</td><td>CCA</td><td>CCA</td></tr><tr><td>Gln</td><td>CAG*</td><td>CAG*</td><td>CAG*</td><td>CAG*</td><td>CAG*</td><td>CAA</td><td>CAA</td><td>CAA</td><td>CAA</td><td>CAA</td></tr><tr><td>Leu</td><td>CTG*</td><td>CTG</td><td>CTG</td><td>CTG*</td><td>CTC*</td><td>CTC*</td><td>TTG*</td><td>TTG</td><td>TTG</td><td>TTA</td></tr><tr><td>Val</td><td>GTG</td><td>GTG</td><td>GTG</td><td>GTG</td><td>GTC</td><td>GTC*</td><td>GTT</td><td>GTT</td><td>GTT</td><td>GTT</td></tr><tr><td>Arg</td><td>CGC*</td><td>CGC*</td><td>CGC*</td><td>CGC*</td><td>CGC*</td><td>CGT</td><td>CGT</td><td>AGA</td><td>AGA</td><td>AGA</td></tr><tr><td>Asp</td><td>GAT*</td><td>GAC</td><td>GAC</td><td>GAC</td><td>GAT*</td><td>GAC</td><td>GAT*</td><td>GAC</td><td>GAT*</td><td>GAT*</td></tr></tbody></table>

B.

<table><thead><tr><th></th><th>Hs</th><th>Ms</th><th>Rn</th><th>Dm</th><th>Hc</th><th>Ce</th><th>Sp</th><th>Sc</th><th>Ca</th><th>Pf</th></tr></thead><tbody><tr><td>Phe</td><td>TTT</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTC</td><td>TTT</td><td>TTT</td><td>TTT</td><td>TTT</td></tr><tr><td>His</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAC</td><td>CAT</td><td>CAT</td><td>CAT</td><td>CAT</td><td>CAT</td><td>CAT</td></tr><tr><td>Tyr</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAC</td><td>TAT</td><td>TAT</td><td>TAT</td><td>TAT</td><td>TAT</td></tr><tr><td>Asn</td><td>AAT</td><td>AAC</td><td>AAC</td><td>AAC</td><td>AAT</td><td>AAT</td><td>AAT</td><td>AAT</td><td>AAT</td><td>AAT</td></tr><tr><td>Lys</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAG</td><td>AAA</td><td>AAA</td><td>AAA</td><td>AAA</td><td>AAA</td></tr><tr><td>Cys</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGC</td><td>TGT</td><td>TGT</td><td>TGT</td><td>TGT</td><td>TGT</td></tr><tr><td>Glu</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAG</td><td>GAA</td><td>GAA</td><td>GAA</td><td>GAA</td><td>GAA</td><td>GAA</td></tr><tr><td>Ala</td><td>GCC</td><td>GCC</td><td>GCC</td><td>GCC</td><td>GCC</td><td>GCT</td><td>GCT</td><td>GCT</td><td>GCT</td><td>GCT</td></tr><tr><td>Ser</td><td>AGC</td><td>AGC</td><td>AGC</td><td>AGC</td><td>TCC</td><td>TCA</td><td>TCT</td><td>TCT</td><td>TCA</td><td>AGT</td></tr><tr><td>Thr</td><td>ACC</td><td>ACC</td><td>ACC</td><td>ACC</td><td>ACC</td><td>ACA</td><td>ACT</td><td>ACT</td><td>ACT</td><td>ACA</td></tr><tr><td>Ile</td><td>ATC</td><td>ATC</td><td>ATC</td><td>ATC</td><td>ATC</td><td>ATT</td><td>ATT</td><td>ATT</td><td>ATT</td><td>ATA</td></tr><tr><td>Gly</td><td>GGC</td><td>GGC</td><td>GGC</td><td>GGC</td><td>GGC</td><td>GGA</td><td>GGT</td><td>GGT</td><td>GGT</td><td>GGA</td></tr><tr><td>Pro</td><td>CCC</td><td>CCC</td><td>CCC</td><td>CCC</td><td>CCA</td><td>CCA</td><td>CCT</td><td>CCA</td><td>CCA</td><td>CCA</td></tr><tr><td>Gln</td><td>CAG</td><td>CAG</td><td>CAG</td><td>CAG</td><td>CAG</td><td>CAA</td><td>CAA</td><td>CAA</td><td>CAA</td><td>CAA</td></tr><tr><td>Leu</td><td>CTG</td><td>CTG</td><td>CTG</td><td>CTG</td><td>CTC</td><td>CTT</td><td>TTA</td><td>TTG</td><td>TTA</td><td>TTA</td></tr><tr><td>Val</td><td>GTG</td><td>GTG</td><td>GTG</td><td>GTG</td><td>GTC</td><td>GTT</td><td>GTT</td><td>GTT</td><td>GTT</td><td>GTA</td></tr><tr><td>Arg</td><td>CGG</td><td>CGG</td><td>CGG</td><td>CGG</td><td>CGG</td><td>AGA</td><td>CGT</td><td>AGA</td><td>AGA</td><td>AGA</td></tr><tr><td>Asp</td><td>GAT</td><td>GAC</td><td>GAC</td><td>GAC</td><td>GAT</td><td>GAT</td><td>GAT</td><td>GAT</td><td>GAT</td><td>GAT</td></tr></tbody></table>

Color code: on a blue background codons shared with *M. musculus*, on a red background, codons shared with *C. albicans*. A: Bolded codons with a star are preferred codons for which a codon with more acceptors exists. B: bolded and underlined codons are identical in lowly and highly translated genes. The following abbreviations of species names are used throughout this article: H. sapiens (Hs), *H. sapiens* (Hs), *M. musculus* (Mm), *R. norvegicus* (Rn), *D. melanogaster* (Dm), *C. elegans* (Ce), *H. capsulatum* (Hc), *S. pombe* (Sp), *S. cerevisiae* (Sc), *C. albicans* (Ca), and *P. falciparum* (Pf).

The last base of a codon is supposed to correlate with the GC content of a genome. This result is partially confirmed in Fig. 4. On one extremity, *P. falciparum* has a very low GC content and no preferred codon ending with a G or C. At the other extremity, all species with a high GC content strongly prefer codons ending with a G or C. For intermediate values of GC content, results are more puzzling: *C. elegans* has an intermediate GC content compared to *C. albicans*, *S. pombe* and *S. cerevisiae*, but has a much stronger preference for codons ending with a G or C.

4. Conclusion and discussion

CUB measures the difference in the frequency of translation of each possible codon for an amino acid. CUB is detected in highly

expressed genes because the selection pressure to choose optimal codons is stronger in those genes. First, it is known that numerous mRNAs are finally not translated²² and have thus no impact on fitness. Hence, to measure CUB, Ribo-seq seems more appropriate to select highly expressed genes than RNA-seq data. Second, various forces and processes (such as mutation bias) impact the nucleotidic composition of a gene.¹⁷ Because it captures the presence of ribosomes, Ribo-seq allows to measure the sole translational advantage of a codon rather than the impact of evolutionary forces unrelated to translation.

Our method based on Ribo-seq—rather than RNA-seq—data allows us to compute a precise and direct measure of CUB. Although Ribo-seq was not originally designed for CUB studies, we show that
<table><caption>**Table 2.** Comparison of favourite codons and tRNA copy numbers</caption><thead><tr><th>Species</th><th>Number of tRNAs</th><th>Number of favourite codons having the highest tRNA copy number (out of 18 cases)</th><th>(%)</th><th>Number of favorite codons lacking tRNA</th><th>Number of codons lacking tRNA</th></tr></thead><tbody><tr><td><i>H. sapiens</i></td><td>589</td><td>9</td><td>50.0</td><td>2</td><td>6</td></tr><tr><td><i>C. elegans</i></td><td>563</td><td>12</td><td>66.7</td><td>6</td><td>15</td></tr><tr><td><i>M. musculus</i></td><td>441</td><td>12</td><td>66.7</td><td>2</td><td>8</td></tr><tr><td><i>R. norvegicus</i></td><td>385</td><td>10</td><td>55.6</td><td>5</td><td>13</td></tr><tr><td><i>D. melanogaster</i></td><td>268</td><td>12</td><td>66.7</td><td>6</td><td>18</td></tr><tr><td><i>S. cerevisiae</i></td><td>259</td><td>15</td><td>83.3</td><td>3</td><td>20</td></tr><tr><td><i>S. pombe</i></td><td>158</td><td>15</td><td>83.3</td><td>2</td><td>16</td></tr><tr><td><i>C. albicans</i></td><td>119</td><td>15</td><td>83.3</td><td>3</td><td>20</td></tr><tr><td><i>H. capsulatum</i></td><td>71</td><td>8</td><td>44.4</td><td>5</td><td>16</td></tr><tr><td><i>P. falciparum</i></td><td>33</td><td>14</td><td>77.8</td><td>7</td><td>26</td></tr></tbody></table>

(The codons for Stop, Met and Trp are excluded). tRNA: transfer RNA. Copy numbers of tRNA genes come from GtRNAbd 2.0¹⁹. The percentage of favourite codons (one per amino acid) that have the highest tRNA copy number lies above 50% for all species but *H. capsulatum*.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Genome percentage of GC bases (GC %) versus GC % of last base of preferred codons. Same species abbreviations as in Fig. 1.</figcaption></figure>

our method assesses the CUB with high precision. Using translation data provides other benefits: first, all translated genes can be included even if functional information is lacking. Second, our method avoids any statistical optimization as CUB is directly computed from translation data. Third, Ribo-seq data allows to compute CUB and to detect preferred codons in a wide variety of species. Fourth, our method may help to measure variation of CUB across conditions or tissues. Nevertheless, the number of Ribo-seq footprints per gene depends on the sequencing depth, which may introduce some bias if depth is low. To correctly classify genes, this method also requires a good sequencing depth in order to get at least some footprints on lowly translated genes: for example, we can grossly estimate that it requires at least 3,000 covered genes to study human's CUB.

Our method was validated by comparison to previous results and our results confirm the current knowledge on CUB. First, the CUB is stronger in highly translated genes.⁷ Second, most higher eukaryotes exhibit a weaker CUB when compared to unicellular species, like yeasts.¹⁷ Third, tRNA copy number has a strong influence on preferred codons, as they indicate preferred codons for most of amino acids.¹⁸ Nevertheless, the tRNA copy number is not the only criterion: the GC content also influences preferred codons, but to a lesser extent than the tRNA pool. The clustering of pluricellular species in Fig. 3B and C suggests that translation can adapt to a changing GC content. Surprisingly, an increasing tRNA pool size does not diminish the number of preferred codons lacking tRNA.

Ribo-seq data on a larger set of species would help to get a better understanding of CUB evolution, but to date, only very few species have been subjected to Ribo-seq experiments. Nevertheless, comparison of codon bias between species seems to reveal a common story that is partially adapted by every organism. In yeast, the level of ribosome occupancy measured using an adequate protocol was found to be negatively correlated with tRNA copy numbers, which is used as a proxy of cognate tRNA abundances.¹³ Our results suggest that the tRNA pool is the key factor for every amino acid in a wide range of species, and that some preferred codons are conserved across the tree of life. However, we clearly distinguished two species groups when comparing the intensities of CUB. The low-intensity group contains mammals and drosophila (plus *P. falciparum* that has a very restricted tRNA pool, which may be due to the difficulty of tRNA gene prediction, and *H. capsulatum* which we discussed later) versus a three-time higher CUB intensity group formed by *C. elegans*, *S. pombe*, *C. albicans* and *S. cerevisiae*. Surprisingly, *C. elegans* does not belong to the low intensity group, while mammals and drosophila do. It has common features with both groups: some amino acids with a yeast-like preferences (like valine or proline), some with a mammals-like preferences (like glutamic acid or aspartic acid) and some with a unique profile (glycine and leucine). A possible scenario could be a decreasing intensity of CUB in multicellular species: nonetheless, this evolution would be achieved progressively, amino acid by amino acid (often through a substitution at the third base of a codon²³—see examples below). Highly translated genes are more subjected to CUB than lowly translated genes and, as shown in Table 1, a general pattern appears from yeasts to mammals: an evolution CAT to CAC for Histidine, from GCT to GCC in Alanine, for example. This general pattern admits some exceptions: for instance, *H. capsulatum* and *C. elegans* favour GTC for valine whereas the common pattern favours an evolution from GTT to GTG. Besides the general tendencies, CUB evolution is to some degree dependent of the species, as suggested by some species specific preferred codons. Including more species would help to further unravel the evolution of CUB, or to investigate the generality of its impacts that were detected only in some species²⁴,²⁵; thus we hope that Ribosome profiling will be applied to an increasing number of species.

The most striking results were obtained for the pathogenic yeast, *H. capsulatum*. Its tRNA pool is closer to that of other yeasts, but its CUB much closer to that of mammals. This suggests that its CUB is
quite independent of its tRNA pool. As *H. capsulatum* is a pathogen of mammals, it is tempting to say that it adapts to its host, but we lack information to conclude. Is *H. capsulatum* able to use the translational machinery of its host (ribosomes and/or tRNAs) or do more tRNA genes remain to be annotated in its genome? Including other pathogenic species in the comparison, such as the fungi Aspergillus, would help to determine whether this feature is specific to *H. capsulatum* or widespread in pathogens. Nevertheless, pathogens could represent an exception in which the CUB depends more on the host's tRNA pool than on their own tRNA pool.

## Availability

The program for computing the codon usage bias from Ribo-seq is freely available at: http://www.lirmm.fr/~rivals/rscu

## Acknowledgements

The authors want to thank France Génomique and the ATGC platform for the logistic support they provided. We thank Nicolas Galtier for his insightful comments on the manuscript.

## Conflict of interest

None declared.

## Supplementary data

Supplementary data are available at www.dnaresearch.oxfordjournals.org.

## Funding

This work is supported by the Fondation pour la Recherche Médicale, grant DBI20131228574, ANR grant (ANR-11-BINF-0002).

## References

1. Grantham, R., Gautier, C., Gouy, M., Mercier, R., and Pavé, A. 1980, Codon catalog usage and the genome hypothesis, *Nucleic Acids Res.*, **8**, r49–62.

2. Gouy, M., and Gautier, C. 1982, Codon usage in bacteria: correlation with gene expressivity, *Nucleic Acids Res.*, **10**, 7055–74.

3. Plotkin, J. B., and Kudla, G. 2011, Synonymous but not the same: the causes and consequences of codon bias. *Nat. Rev. Genet.*, **12**, 32–42.

4. Quax, T. E. F., Claassens, N. J., Söll, D., and van der Oost, J. 2015, Codon bias as a means to fine-tune gene expression, *Mol. Cell*, **59**, 149–61.

5. Pechmann, S., and Frydman, J. 2013, Evolutionary conservation of codon optimality reveals hidden signatures of cotranslational folding, *Nat. Struct. Mol. Biol.*, **20**, 237–43.

6. Yu, C.-H., Dang, Y., Zhou, Z., et al. 2015, Codon usage influences the local rate of translation elongation to regulate co-translational protein folding, *Mol. Cell*, **59**, 744–54.

7. Sharp, P. M., Tuohy, T. M., and Mosurski, K. R. 1986, Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes, *Nucleic Acids Res.*, **14**, 5125–43.

8. Ingolia, N. T., Ghaemmaghami, S., Newman, J. R. S., and Weissman, J. S. 2009, Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling, *Science*, **324**, 218–23.

9. Guo, B.-Y., Zeng, T., and Wu, H.-C. 2014, Recent advances of DNA sequencing via nanopore-based technologies, *Sci. Bull.*, **60**, 287–95.

10. Dunn, J. G., Foo, C. K., Belletier, N. G., Gavis, E. R., and Weissman, J. S. 2013, Ribosome profiling reveals pervasive and regulated stop codon readthrough in Drosophila melanogaster, *eLife*, **2**, e01179.

11. Zupanic, A., Meplan, C., Grellscheid, S. N., et al. 2014, Detecting translational regulation by change point analysis of ribosome profiling data sets, *RNA N. Y. N*, **20**, 1507–18.

12. Gao, X., Wan, J., Liu, B., Ma, M., Shen, B., and Qian, S.-B. 2015, Quantitative profiling of initiating ribosomes in vivo, *Nat. Methods*, **12**, 147–53.

13. Weinberg, D. E., Shah, P., Eichhorn, S. W., Hussmann, J. A., Plotkin, J. B., and Bartel, D. P. 2016, Improved ribosome-footprint and mrna measurements provide insights into dynamics and regulation of yeast translation, *Cell Rep.*, **14**, 1787–99.

14. Sharp, P. M., and Bradnam, K. R. 1997, Appendix 3 Codon Usage in *C. elegans*. In: Riddle, D.L., Meyer, T., Blumenthal, B.J. Al, editors. *C. elegans II*. 2nd Ed. Cold Spring Harbor Laboratory Press: Cold Spring Harbor (NY).

15. Philippe, N., Salson, M., Commes, T., and Rivals, E. 2013, CRAC: an integrated approach to the analysis of RNA-seq reads. *Genome Biol.*, **14**, R30.

16. Stadler, M., and Fire, A. 2011, Wobble base-pairing slows in vivo translation elongation in metazoans, *RNA*, **17**, 2063–73.

17. Lynch, M. 2007, *The Origins of Genome Architecture*. Sinauer Associates, Indiana University at Bloomington.

18. Duret, L. 2000, tRNA gene number and codon usage in the *C. elegans* genome are co-adapted for optimal translation of highly expressed genes, *Trends Genet.*, **16**, 287–9.

19. Chan, P. P., and Lowe, T. M. 2016, GtRNAdb 2.0: an expanded database of transfer RNA genes identified in complete and draft genomes. *Nucleic Acids Res.*, **44**, D184–9.

20. Lloyd, A. T., and Sharp, P. M. 1992, Evolution of codon usage patterns: the extent and nature of divergence between *Candida albicans* and *Saccharomyces cerevisiae*, *Nucleic Acids Res.*, **20**, 5289–95.

21. Crick, F. H. 1966, Codon-anticodon pairing: the wobble hypothesis, *J. Mol. Biol.*, **19**, 548–55.

22. Schwanhäusser, B., Busse, D., Li, N., et al. 2011, Global quantification of mammalian gene expression control, *Nature*, **473**, 337–42.

23. Ikemura, T. 1985, Codon usage and tRNA content in unicellular and multicellular organisms, *Mol. Biol. Evol.*, **2**, 13–34.

24. Zhou, Z., Dang, Y., Zhou, M., et al. 2016, Codon usage is an important determinant of gene expression levels largely through its effects on transcription, *Proc. Natl. Acad. Sci.*, **113**(41), E6117–E6125.

25. Gamble, C. E., Brule, C. E., Dean, K. M., Fields, S., and Grayhack, E. J. 2016, Adjacent codons act in concert to modulate translation efficiency in yeast, *Cell*, **166**, 679–90.