BACKGROUND
mrna degradation is a critical step in gene expression, and the decay rates of individual mrnas can vary over two orders of magnitude.
differences in the decay rates of individual mrnas can be specified by several features of the mrnas.
they include sequence motifs that are recognized by trans-acting factors, such as micrornas and rna-binding proteins.
strikingly, in many of these cases, the trans-acting factors can also decrease translation initiation, which suggests a tight coupling of translation initiation and mrna degradation  <cit> .
perturbations of translation elongation can also affect mrna degradation.
for example, strong blocks to translation elongation trigger endonucleolytic cleavage of the mrna in a process called no-go decay  <cit> .
however, until a recent study by coller and colleagues  <cit> , it was not appreciated that subtle differences in the rates of translation elongation due to specific codons, which can be indicated by “codon optimality,” would contribute in a general manner to defining mrna decay rates.
the general model is that “optimal” codons, which are decoded efficiently, are associated with mrna stability, whereas “nonoptimal” codons, which are decoded slowly, are associated with mrna instability.
the study in saccharomyces cerevisiae by coller and colleagues was followed by multiple studies in diverse organisms arguing that codon-mediated mrna decay is a broadly conserved phenomenon .
moreover, data suggest that the codon-mediated mrna decay is accompanied by a reduction in translation efficiency, defined as protein synthesis rates per mrna  <cit> , which may result from a decrease in translation initiation rate  <cit> .
a more recent study in s. cerevisiae by grayhack, fields, and colleagues has demonstrated that adjacent codon pairs also influence translation in a manner distinct from their individual constituent codons  <cit> .
this suggests that optimality of individual codons does not solely define the relationship between codon composition and translation efficiency.
specifically, via a large-scale flow cytometry analysis using gfp reporter variants, the authors have identified  <dig> adjacent codon pairs that act inhibitory on protein expression.
a series of subsequent analyses using the reporter system suggest several important aspects of the inhibition of protein expression mediated by the codon pairs  <cit> .
first, the inhibition is mediated by the codon pairs themselves rather than by the corresponding hexanucleotide sequences, individual constituent codons, or encoded dipeptides.
second, the inhibition occurs during translation elongation largely depending on wobble decoding of either the 5’ and/or 3’ codon.
third, in some instances, the inhibitory codon pairs are associated with a reduction in mrna abundance consistent with the codon pairs eliciting mrna instability.
the findings of additional analyses of genomic data imply that the inhibitory codon pairs may be relevant to the regulation of natural endogenous genes  <cit> .
first, ribosome occupancy is substantially elevated at most of the inhibitory codon pairs, which suggests that the codon pairs slow translation elongation.
second, the inhibitory codon pairs are enriched in genes whose mrna abundance is low.
third, translation efficiency, as assessed by protein abundance per mrna, of genes containing at least one of the  <dig> inhibitory codon pairs are significantly lower than that of genes lacking them.
this tendency persists even when the analysis is controlled for usage bias in individual codons.
the observation that inhibitory codon pairs can reduce translation elongation rates and affect mrna levels raises the possibility that the effects of codon optimality on gene expression parameters could be explained by the presence of inhibitory codon pairs  <cit> .
alternatively, both overall “codon optimality” and inhibitory codon pairs could act in a similar manner to slow elongation and thereby both contribute to changes in gene expression as a downstream read-out of translation elongation rates.
in this work, we examined relationships between the inhibitory codon pairs and additional aspects of gene expression on a genomic scale in s. cerevisiae.
our analysis revealed consistent associations of the inhibitory codon pairs with fast mrna decay in multiple rna kinetic datasets.
the association persisted after controlling for the content of nucleotides, optimality of codons contained by the mrna, and the content of dipeptides encoded by the inhibitory pairs, suggesting that the link between the inhibitory codon pairs and mrna decay rates is not solely due to effects of the covariates.
additional analyses suggest that the presence of the inhibitory codon pairs cannot fully explain the relationship of codon optimality with mrna stability or translation efficiency.
overall, our study has revealed novel aspects concerning the relationship between adjacent codon pairs and the regulation of gene expression.
RESULTS
the inhibitory codon pairs are associated with inefficient synthesis and instability of mrna
it has been shown that genes with the inhibitory codon pairs are enriched in genes whose mrna abundance is low  <cit> .
since mrna abundance is determined by rates of mrna synthesis and decay, we examined whether the previously identified  <dig> inhibitory codon pairs  are associated with mrna synthesis and/or decay.
for this purpose, we computed spearman’s and kendall’s correlation coefficients between the outcome  and predictor  variables .
as an outcome, we used rna kinetic values from two previous studies by cramer and colleagues and one by gresham and colleagues  .
we selected these datasets because the metabolic labeling method, which was used in the studies, has been suggested to be less intrusive than others  <cit> .
to examine whether conclusions from the analyses are independent from the methods of rna kinetic measurements, we also analyzed a dataset from a previous study by coller and colleagues , which was generated via transcription inhibition  <cit> .
as a predictor, we used two types of variables, the fraction of inhibitory codon pairs contained in mrna and a binary variable to indicate the presence or absence of at least one of the  <dig> inhibitory codon pairs.
in this analysis, we aggregated the inhibitory codon pairs rather than analyzing them individually since some of them occur very infrequently in the coding sequences .
for each dataset, the number of genes containing at least one of the inhibitory codon pairs is shown in additional file 1: table s <dig>  although we computed p values by methods described by kim  <cit>  as well as by permutation tests , we primarily used permutation p values with a significant threshold of p <  <dig>  for hypothesis testing.table  <dig> the  <dig> inhibitory codon pairs
shown are the nucleotide sequences of the inhibitory codon pairs, the number of occurrences of each inhibitory codon pair, and the number of genes that contain at least one of the inhibitory codon pairs
the analysis of the four datasets consistently led to two associations.
first, the fraction of the inhibitory codon pairs as well as the binary indicator of the presence thereof was associated with slow mrna synthesis rates, which was suggested by negative correlations .
second, both variables representing the content of the inhibitory codon pairs were also associated with fast mrna decay rates, which was suggested by positive correlations .fig.
<dig> associations of the inhibitory codon pairs with synthesis and decay of mrna and protein.
a boxplot comparing mrna synthesis rates in the “cramer 1” data  between genes containing at least one of the  <dig> inhibitory codon pairs  and those without them .
shown on the top are the spearman’s and kendall’s correlation coefficients and p values  to assess an association of the presence  and absence  of the inhibitory codon pairs with mrna synthesis rate.
b same as  but for the “cramer 2” data.
c same as  but for the “gresham” data.
d same as  but for the “coller” data.
e same as  but for mrna decay rate.
f same as  but for the “cramer 2” data.
g same as  but for the “gresham” data.
h samea s  but for the “coller” data.
i same as  but for protein abundance per mrna.
j same as  but for ribosome occupancy.
k same as  but for protein decay rates
ρ
τ
ρ
τ
ρ
τ
ρ
τ
spearman’s and kendall’s correlation coefficients to assess an association between the fraction of the inhibitory codon pairs and various gene expression variables.
p values obtained according to kim  <cit>  and those based on permutation tests are shown.
same as  but for the presence/absence of the inhibitory codon pairs.
spearman’s and kendall’s partial correlation coefficients controlled for gc content, tai, dipeptide content, coding length to assess an association between the fraction of the inhibitory codon pairs and various gene expression variables.
same as  but for the presence/absence of the inhibitory codon pairs
in principle, the observed associations could be due to the inhibitory codon pairs or to other transcript features since several other transcript features were correlated with the content of the inhibitory codon pairs as well as with mrna synthesis/decay rates in some of the datasets .
these include guanine-cytosine  content, trna adaptation index , which is a metric of codon optimality , the fraction of dinucleotides that are encoded by the inhibitory pairs, and the lengths of coding sequences.
to evaluate contributions of these transcript features to the observed associations, we computed spearman’s and kendall’s partial correlation coefficients between the content of the inhibitory codon pairs and mrna decay/synthesis rates controlling for these confounding factors.
the analysis led to the following two points.
first, for mrna decay rates, the association remained significant when the analysis was individually controlled for gc content, codon optimality , dinucleotide content, and coding lengths with one exception where we analyzed a relationship between the presence/absence of the inhibitory codon pairs and mrna decay rates in the “coller” dataset controlling for codon optimality  using the spearman’s method  .
the association was significant across the analysis methods and datasets when the analysis was controlled for all the covariates .
second, for mrna synthesis rates, the association was no longer consistent across the analysis methods/datasets when the analysis was controlled for codon optimality  or for all covariates .
overall, the results suggest an association of the inhibitory codon pairs with mrna instability on a genomic scale, which appears to be, at least in part, independent of the content of nucleotides, individual constituent codons, and encoded dipeptides as well as of coding lengths.
although we also observed an association between the inhibitory codon pairs and mrna synthesis rates, this could be attributed to usage bias in individual codons.
the association between inhibitory codon pairs and mrna instability is largely dependent on the correct reading frame
the simplest model is that inhibitory codon pairs correlate with mrna decay rates due to their effects on translation elongation and would thus only correlate with mrna decay rates when present in the proper reading frame.
alternatively, it remains possible that the hexanucleotides making up inhibitory codon pairs could affect mrna decay rates directly.
for example, one possibility is that a subset of the sequences might coincidentally match those recognized by trans-acting factors that promote mrna decay.
to address this issue, we examined an association of the occurrence of the hexanucleotide sequences in shifted reading frames as wells as in 3’ untranslated regions  with mrna decay rates.
in the former analysis, we computationally introduced frameshifts by one or two nucleotides to all orfs and repeated otherwise the same correlation analyses as described above.
in the latter analysis, to all orfs, we assigned binary indicators to represent the presence/absence of at least one of the inhibitory codon pairs within 3’ utr based on annotations from previous studies by snyder and colleagues  <cit>  and by steinmetz and colleagues and  <cit> .
the analyses led to the following points.
first, the inhibitory codon pairs in the + <dig> frame were associated with fast mrna decay rates across the analysis methods/datasets except for the “coller” data .
however, the association became inconsistent when gc content, codon optimality, the content of dipeptides encoded by the in-frame inhibitory codon pairs, and coding lengths were individually or simultaneously controlled for .
second, the hexanucleotide sequences in 3’ utr were not consistently associated with mrna decay rates .table  <dig> test for associations of the out-of-frame inhibitory codon pairs with mrna decay rate, protein per mrna, and ribosome occupancy
ρ
ρ
ρ
τ
τ
τ
ρ
ρ
ρ
τ
τ
τ
spearman’s partial correlation coefficients controlled for gc content, tai, dipeptide content, and coding length to assess an association between the fraction of hexanucleotide sequences corresponding to the inhibitory codon pairs in the  <dig>  + <dig>  and + <dig> frames and various gene expression variables.
p values obtained according to kim  <cit>  and those based on permutation tests are shown.
same as  but for kendall’s partial correlation coefficients.
same as  but for the presence/absence of the hexanucleotide sequences.
same as  but for the presence/absence of the hexanucleotide sequences
ρ
τ
ρ
τ
spearman’s and kendall’s correlation coefficients to assess an association between the presence/absence of hexanucleotide sequences corresponding to the inhibitory codon pairs in 3’ utr regions.
the utr annotations are based on a study by snyder and colleagues  <cit> .
p values obtained according to kim  <cit>  and those based on permutation tests are shown.
same as  but for utr annotations based on a study by steinmetz and colleagues  <cit>
in sum, the results are largely consistent with the inhibitory codon pairs affecting mrna decay primarily via its effects on translation elongation kinetics.
the inhibitory codon pairs cannot fully explain the association of codon optimality with mrna decay
the inhibitory codon pairs consist of ten types of codons , all of which are classified as nonoptimal based on tai   <cit> .
therefore, one possibility is that the inhibitory codon pairs could explain the association of codon optimality with mrna decay rates  <cit> .
two observations argue that this is unlikely to be the case.
first, for genes lacking the inhibitory codon pairs, codon optimality was still significantly associated with mrna decay rates .
second, there was a significant partial correlation of codon optimality with mrna decay rates even when the analysis was individually or simultaneously controlled for gc content, the content of inhibitory pairs, dipeptide content, and coding length .fig.
<dig> associations of codon optimality with mrna decay rate and translation efficiency for genes lacking the inhibitory codon pairs.
a scatterplot comparing tai, a metric of codon optimality, and mrna decay rate in the “cramer 1” data .
shown on the top are the spearman’s and kendall’s correlation coefficients and p values .
b same as  but for the “cramer 2” data.
c same as  but for the “gresham” data.
d same as  but for the “coller” data.
e same as  but for protein abundance per mrna.
f same as  but for ribosome occupancy
ρ
τ
ρ
τ
ρ
τ
spearman’s and kendall’s correlation coefficients to assess an association between codon optimality and various gene expression variables.
p values obtained according to kim  <cit>  and those based on permutation tests are shown.
same as  but for partial correlation coefficients controlled for gc content, the fraction of the inhibitory codon pairs, dipeptide content, and coding length.
same as  but for partial correlation coefficients controlled for gc content, the presence/absence of the inhibitory codon pairs, dipeptide content, and coding length
to examine further whether the inhibitory codon pairs can explain effects of codon optimality on mrna decay rates, we took advantage of synonymous reporter systems used in the previous studies that suggested a mechanistic link between codon optimality and mrna stability .
specifically, we compiled and analyzed the content of the inhibitory codon pairs in sequences of reporter constructs used in the studies .
the analyses led to the following two points, which imply that the inhibitory codon pairs can promote mrna instability but cannot explain the relationship between codon optimality and mrna stability.
first, when comparisons are made within a synonymous group, the higher the content of the inhibitory codon pairs, the faster mrna decay.
second, there are multiple examples where synonymous transcripts that differ in codon optimality but not in the fraction of the inhibitory codon pairs exhibit different mrna decay rates.
the latter point is particularly important because, if the effect of codon optimality were solely due to the inhibitory codon pairs, the synonymous transcripts lacking the inhibitory codon pairs would show similar mrna decay rates.table  <dig> number of the inhibitory codon pairs in reporter systems
pgk1
lsm8
his3
naa
his3
naa
naa
naa
naa
naa
naa
naa
naa
naa
naa
ana, not applicable
bthe effects on mrna decay rate have been suggested based on mrna abundance
overall, these results are consistent with the idea that codon optimality affects mrna stability at least in part independently of the inhibitory codon pairs.
the inhibitory codon pairs do not show position effects
a previous reporter-based study has shown that a stretch of nonoptimal codons exhibits an increasing destabilizing effect on mrna with an increasing distance from the start codon  <cit> .
based on the assumption that the stretch of the nonoptimal codons causes ribosome queuing along the upstream region, the observation was interpreted to suggest that the higher the number of slow ribosomes on a transcript the less stable the transcript  <cit> .
to examine whether the inhibitory codon pairs have a similar property, we computed spearman’s and kendall’s correlation coefficients between distances of the inhibitory codon pairs from the start codon and mrna decay rates in the genome-wide rna kinetic data.
in this analysis, we focused on  <dig> orfs that contain one and only one of the inhibitory codon pairs.
if the inhibitory codon pairs had a stronger effect with an increasing distance from the start codon, the distances would be positively correlated with mrna decay rates.
however, we did not observe consistent correlations across the datasets .
this suggests that the inhibitory codon pairs are unlikely to cause a long-range ribosome queuing in upstream regions of natural endogenous mrnas.fig.
<dig> lack of position effect of the inhibitory codon pairs.
a scatterplot comparing mrna decay rate in the “cramer 1” data  and distances of the inhibitory codon pairs from the start codons contained by the mrnas.
shown on the top are the spearman’s and kendall’s correlation coefficients and p values .
b same as  but for the “cramer 2” data.
c same as  but for the “gresham” data.
d same as  but for the “coller” data
codon optimality, adjacent codon pairs, and translation efficiency
it has been shown that codon optimality and the inhibitory codon pairs are associated with translation efficiency  <cit> .
however, another study did not find a correlation between codon optimality and translation efficiency  <cit> .
to reassess this issue, we selected recently published genome-scale data and examined the association between codon optimality and translation efficiency.
translation efficiency, defined as the rate of protein synthesis per mrna, can be obtained by various methods  <cit> .
one among them is to measure protein abundance and normalize it against mrna abundance.
this is based on the assumption that the majority of proteins are stable and that protein abundance is largely determined by mrna abundance and protein synthesis rates, which is consistent with a recent protein half-life measurement in s. cerevisiae  <cit> .
another is to use ribosome occupancy on mrna as a proxy for protein synthesis rates.
this is based on the assumption that the majority of mrna-bound ribosomes are actively engaged in translation, which is consistent with a recent observation in growing s. cerevisiae cells that ribosome occupancy is highly correlated with the rate of translation initiation  <cit> .
for our analyses, we selected recent quantitative proteomic data by mann and colleagues  <cit> , which was used in the previous study by grayhack, fields, and colleagues  <cit> , as well as mrna-seq and ribosome profiling data by weinberg and colleagues  <cit> .
the ribosome profiling data was chosen for two reasons.
first, the protocol used for mrna quantification in the studies does not involve poly enrichment, which could result in 3’ bias, an overestimation of mrna abundance of short mrnas, and thus an underestimation of ribosome occupancy of short mrnas  <cit> .
second, the protocol used for quantification of ribosome-protected rna fragments does not involve cycloheximide treatment, which could introduce multiple artifacts.
using these datasets, we observed that translation efficiency in both metrics was positively correlated with codon optimality .fig.
<dig> associations of codon optimality with protein synthesis and decay.
a scatterplot comparing tai, a metric of codon optimality, and protein abundance per mrna .
shown on the top are the spearman’s and kendall’s correlation coefficients and p values .
b same as  but for ribosome occupancy.
c same as  but for protein decay rates
we then reexamined the association between the inhibitory codon pairs and translation efficiency in the selected data.
for this purpose, we computed spearman’s and kendall’s correlation and partial correlation coefficients controlling for the same set of covariates as we used for the analysis of mrna synthesis/decay rates.
the analysis led to the following observations, which is consistent with the previous study  <cit> .
first, the fraction of the inhibitory codon pairs as well as the binary indicator of the presence thereof was associated with low protein abundance per mrna as well as with low ribosome occupancy .
second, the association remained significant when the analysis was individually controlled for gc content, codon optimality, dinucleotide content, and coding lengths .
third, the association still remained significant when the analysis was controlled for all the covariates .
to assess contributions from the hexanucleotide sequences corresponding to the inhibitory codon pairs, we also examined associations of the sequences in + <dig> and + <dig> frames and in 3’ utrs with translation efficiency.
the analyses led to the following three points.
first, neither inhibitory codon pairs in the + <dig> frame nor those in the + <dig> frame were consistently associated with protein abundance per mrna .
second, although the content of the inhibitory codon pairs in the + <dig> and + <dig> frames were consistently associated with low ribosome occupancy , the association became inconsistent when the covariates were controlled for .
third, the hexanucleotide sequences in 3’ utrs were not associated with protein abundance per mrna or ribosome occupancy .
overall, these analyses confirmed a positive correlation between codon optimality and translation efficiency as well as an association between the inhibitory codon pairs and low translation efficiency, which is at least in part independent of effects of nucleotide content, codon optimality, the content of encoded dipeptides, coding lengths, and corresponding hexanucleotide sequences.
the results are consistent with the inhibitory codon pairs affecting protein synthesis primarily via its effects on translation elongation kinetics.
the inhibitory codon pairs cannot fully explain the association of codon optimality with translation efficiency
having confirmed the associations of codon optimality and the inhibitory codon pairs with translation efficiency, we next examined the possibility that the inhibitory codon pairs can explain the association of codon optimality with translation efficiency.
as was seen for mrna decay rates, two observations argue that this is unlikely to be the case.
first, for genes lacking the inhibitory codon pairs, codon optimality was significantly associated with the both metrics of translation efficiency, protein abundance per mrna  and ribosome occupancy .
second, there was a significant partial correlation of codon optimality with mrna decay rates and translation efficiency even when the analysis was individually or simultaneously controlled for gc content, the content of inhibitory pairs, dipeptide content, and coding length .
codon optimality, adjacent codon pairs, and protein stability
translation kinetics can affect co-translational protein folding, which can, in turn, affect protein function and stability.
we thus examined associations between codon optimality or adjacent codon pair content and protein decay rates.
for this analysis, we used a genome-wide protein half-life measurement by walther and colleagues  <cit> .
the data was generated via metabolic labeling and thus likely to be more accurate than those generated via methods using translation inhibitors.
the analysis led to the following observations.
first, codon optimality was negatively correlated with protein decay rates .
second, the inhibitory codon pairs were significantly associated with fast protein decay rates .
third, however, the association between the inhibitory codon pairs and protein instability became insignificant when tai was controlled for  and insignificant or borderline significant when all covariates were controlled for .
these results suggest that the inhibitory codon pairs are unlikely to be directly linked to protein stability.
analysis of all possible  <dig> adjacent codon pairs
the above analyses are focused on  <dig> inhibitory codon pairs that have been experimentally characterized previously.
however, a recent computational analysis suggests that other codon pairs can be inhibitory to translation  <cit> .
to examine relationships between the adjacent codon pairs and mrna stability in an unbiased manner, we computed spearman’s and kendall’s correlation coefficients between the fraction of each of the possible  <dig> codon pairs and mrna decay rates in the four datasets .
we then compared the correlation coefficients between the  <dig> inhibitory codon pairs and  <dig> other pairs.
we note that the minimum number of orfs containing a given codon pair with a measured decay rate was as small as eight, which is for the cga di-codon, one of the inhibitory codon pairs, in the “cramer” data.
the analysis led to the following two points.
first, in three of the four datasets, the correlations with mrna decay rates are higher for the inhibitory codon pairs than for other pairs, which is roughly consistent with the association of the inhibitory codon pairs with fast mrna decay  .
second, for a substantial number of codon pairs other than the  <dig> inhibitory codon pairs, their content showed significant positive correlations with mrna decay rates, which may be consistent with those codon pairs acting inhibitory on gene expression as has been suggested in a recent study   <cit> .
overall, the results are roughly consistent with the association of the inhibitory codon pairs with fast mrna decay.
however, they also suggest that correlation analyses of individual codon pairs and measured gene expression variables are highly susceptible to experimental noise and/or other confounding factors.
discussion
we have analyzed relationships between adjacent codon pairs and aspects of gene expression in s. cerevisiae.
our results suggest an association of the inhibitory codon pairs with fast rates of mrna decay.
the association is not solely due to gc content, the content of individual codons, the content of encoded dipeptides, or coding lengths.
moreover, there is no consistent association of the hexanucleotide sequences corresponding to the inhibitory codon pairs in shifted frames or in non-coding regions with mrna decay rates, which is largely consistent with the inhibitory codon pairs affecting mrna stability at least in part via their effects on translation elongation kinetics.
although we also observed associations of the inhibitory codon pairs with mrna synthesis rates and protein instability, they could be attributed to bias in the individual constituent codons.
it has been speculated that adjacent codon pairs rather than individual codons may largely underlie the relationship between codon optimality and translation efficiency  <cit> .
our results suggest that neither the effects of codon optimality on translation efficiency nor those on mrna stability can be explained by the  <dig> inhibitory codon pairs.
this observation is more consistent with the model that a sum of subtle effects from individual codons can exert a large impact on gene expression  <cit> .
however, it is still possible that codon pairs rather than individual codons largely define the relationship of codon composition with translation efficiency and mrna stability since there are likely to be other inhibitory codon pairs than those in the high confidence set obtained by the previous study  <cit> .
the genome-wide associations of codon optimality and the inhibitory codon pairs with mrna stability and translation efficiency observed by this work and other studies suggest two non-mutually exclusive possibilities.
the first possibility is that direct mechanistic links largely underlie the association of the codon composition, which can affect speeds and/or accuracy of translation elongation, with mrna stability and translation efficiency.
several non-mutually exclusive mechanisms can link slow/inaccurate translation elongation to fast mrna decay rates and/or low protein synthesis rates.
first, slow translation elongation can cause ribosome queuing in upstream regions, which could, in turn, interfere with translation initiation  <cit> .
second, slow translation elongation may cause ribosome drop-off.
consistent with this idea, the cga di-codon, one of the inhibitory codon pairs, has been linked to a quality control mechanism that can detect stalled ribosomes and cause abortion of translation  <cit> .
third, slow translation elongation may be sensed by a mechanism that can modulate translation initiation rates.
indeed, a recent study in s. cerevisiae has raised the possibility that the dead-box protein dhh <dig> might play a central role in such a mechanism  <cit> .
fourth, nonoptimal codons and/or the inhibitory codon pairs may result in translation repression and mrna instability by compromising translation fidelity.
for example, it is possible that nonoptimal codons and/or the inhibitory codon pairs increase erroneous translation frameshifts.
although there is no overlap between the  <dig> inhibitory codon pairs and previously identified frameshift-inducing sequences  <cit> , it remains to be determined whether any of the inhibitory codon pairs tend to introduce translation frameshifting.
in most cases, translation frameshifts would result in a premature translation termination event at a stop codon in the incorrect frame, which would in turn cause repression of translation initiation and nonsense-mediated mrna decay .
such a mechanism can partly underlie the association of nonoptimal codons and inhibitory codon pairs with mrna instability.
indeed, a recent study suggests that mrnas with high content of nonoptimal codons tend to undergo nonfunctional translation frameshifts and, subsequently, nmd  <cit> .
the second possibility is that the genome-wide associations between codon composition, translation efficiency, and mrna stability are largely due to co-evolution rather than to a mechanistic link, the latter of which has been suggested by studies using artificial reporter systems.
that is, nonoptimal codons and the inhibitory codon pairs may be simply avoided in natural endogenous genes that are highly expressed and efficiently translated.
then, their primary function may be to modulate local translation elongation speed and thereby regulate other processes, such as co-translational folding, which may be largely restricted to a situation where slow ribosomes do not negatively impact overall translation efficiency.
consistent with this view, some studies in unicellular organisms suggest that under physiological conditions translation initiation but not translation elongation mainly defines the rate of protein synthesis  <cit> .
moreover, another study did not find corresponding changes in translation efficiency upon genetic manipulation of trna and thus codon optimality  <cit> .
clearly, further investigation will be needed to rigorously evaluate these two possibilities concerning the codon-mediated gene regulation.
CONCLUSIONS
this study suggests genome-scale associations of the inhibitory codon pairs with mrna decay and translation efficiency, which, in turn, suggest another layer of complexity in the codon-mediated gene regulation.
an important future goal will be to understand whether and how the inhibitory codon pairs mechanistically inhibit protein synthesis and elicit mrna instability.

 mrna degradation is a critical step in gene expression, and the decay rates of individual mrnas can vary over two orders of magnitude. differences in the decay rates of individual mrnas can be specified by several features of the mrnas. they include sequence motifs that are recognized by trans-acting factors, such as micrornas and rna-binding proteins. strikingly, in many of these cases, the trans-acting factors can also decrease translation initiation, which suggests a tight coupling of translation initiation and mrna degradation. perturbations of translation elongation can also affect mrna degradation. for example, strong blocks to translation elongation trigger endonucleolytic cleavage of the mrna in a process called no-go decay. however, until a recent study by coller and colleagues, it was not appreciated that subtle differences in the rates of translation elongation due to specific codons, which can be indicated by “codon optimality,” would contribute in a general manner to defining mrna decay rates. the general model is that “optimal” codons, which are decoded efficiently, are associated with mrna stability, whereas “nonoptimal” codons, which are decoded slowly, are associated with mrna instability. the study in saccharomyces cerevisiae by coller and colleagues was followed by multiple studies in diverse organisms arguing that codon-mediated mrna decay is a broadly conserved phenomenon. moreover, data suggest that the codon-mediated mrna decay is accompanied by a reduction in translation efficiency, defined as protein synthesis rates per mrna, which may result from a decrease in translation initiation rate. a more recent study in s. cerevisiae by grayhack, fields, and colleagues has demonstrated that adjacent codon pairs also influence translation in a manner distinct from their individual constituent codons. this suggests that optimality of individual codons does not solely define the relationship between codon composition and translation efficiency. specifically, via a large-scale flow cytometry analysis using gfp reporter variants, the authors have identified adjacent codon pairs that act inhibitory on protein expression. a series of subsequent analyses using the reporter system suggest several important aspects of the inhibition of protein expression mediated by the codon pairs. first, the inhibition is mediated by the codon pairs themselves rather than by the corresponding hexanucleotide sequences, individual constituent codons, or encoded dipeptides. second, the inhibition occurs during translation elongation largely depending on wobble decoding of either the 5’ and/or 3’ codon. third, in some instances, the inhibitory codon pairs are associated with a reduction in mrna abundance consistent with the codon pairs eliciting mrna instability. the findings of additional analyses of genomic data imply that the inhibitory codon pairs may be relevant to the regulation of natural endogenous genes. first, ribosome occupancy is substantially elevated at most of the inhibitory codon pairs, which suggests that the codon pairs slow translation elongation. second, the inhibitory codon pairs are enriched in genes whose mrna abundance is low. third, translation efficiency, as assessed by protein abundance per mrna, of genes containing at least one of the inhibitory codon pairs are significantly lower than that of genes lacking them. this tendency persists even when the analysis is controlled for usage bias in individual codons. the observation that inhibitory codon pairs can reduce translation elongation rates and affect mrna levels raises the possibility that the effects of codon optimality on gene expression parameters could be explained by the presence of inhibitory codon pairs. alternatively, both overall “codon optimality” and inhibitory codon pairs could act in a similar manner to slow elongation and thereby both contribute to changes in gene expression as a downstream read-out of translation elongation rates. in this work, we examined relationships between the inhibitory codon pairs and additional aspects of gene expression on a genomic scale in s. cerevisiae. our analysis revealed consistent associations of the inhibitory codon pairs with fast mrna decay in multiple rna kinetic datasets. the association persisted after controlling for the content of nucleotides, optimality of codons contained by the mrna, and the content of dipeptides encoded by the inhibitory pairs, suggesting that the link between the inhibitory codon pairs and mrna decay rates is not solely due to effects of the covariates. additional analyses suggest that the presence of the inhibitory codon pairs cannot fully explain the relationship of codon optimality with mrna stability or translation efficiency. overall, our study has revealed novel aspects concerning the relationship between adjacent codon pairs and the regulation of gene expression.  the inhibitory codon pairs are associated with inefficient synthesis and instability of mrna it has been shown that genes with the inhibitory codon pairs are enriched in genes whose mrna abundance is low. since mrna abundance is determined by rates of mrna synthesis and decay, we examined whether the previously identified inhibitory codon pairs are associated with mrna synthesis and/or decay. for this purpose, we computed spearman’s and kendall’s correlation coefficients between the outcome and predictor variables. as an outcome, we used rna kinetic values from two previous studies by cramer and colleagues and one by gresham and colleagues. we selected these datasets because the metabolic labeling method, which was used in the studies, has been suggested to be less intrusive than others. to examine whether conclusions from the analyses are independent from the methods of rna kinetic measurements, we also analyzed a dataset from a previous study by coller and colleagues, which was generated via transcription inhibition. as a predictor, we used two types of variables, the fraction of inhibitory codon pairs contained in mrna and a binary variable to indicate the presence or absence of at least one of the inhibitory codon pairs. in this analysis, we aggregated the inhibitory codon pairs rather than analyzing them individually since some of them occur very infrequently in the coding sequences. for each dataset, the number of genes containing at least one of the inhibitory codon pairs is shown in additional file 1: table s although we computed p values by methods described by kim as well as by permutation tests, we primarily used permutation p values with a significant threshold of p <  for hypothesis testing.table the inhibitory codon pairs shown are the nucleotide sequences of the inhibitory codon pairs, the number of occurrences of each inhibitory codon pair, and the number of genes that contain at least one of the inhibitory codon pairs the analysis of the four datasets consistently led to two associations. first, the fraction of the inhibitory codon pairs as well as the binary indicator of the presence thereof was associated with slow mrna synthesis rates, which was suggested by negative correlations. second, both variables representing the content of the inhibitory codon pairs were also associated with fast mrna decay rates, which was suggested by positive correlations.fig. associations of the inhibitory codon pairs with synthesis and decay of mrna and protein. a boxplot comparing mrna synthesis rates in the “cramer 1” data between genes containing at least one of the inhibitory codon pairs and those without them. shown on the top are the spearman’s and kendall’s correlation coefficients and p values to assess an association of the presence and absence of the inhibitory codon pairs with mrna synthesis rate. b same as but for the “cramer 2” data. c same as but for the “gresham” data. d same as but for the “coller” data. e same as but for mrna decay rate. f same as but for the “cramer 2” data. g same as but for the “gresham” data. h samea s but for the “coller” data. i same as but for protein abundance per mrna. j same as but for ribosome occupancy. k same as but for protein decay rates ρ τ ρ τ ρ τ ρ τ spearman’s and kendall’s correlation coefficients to assess an association between the fraction of the inhibitory codon pairs and various gene expression variables. p values obtained according to kim and those based on permutation tests are shown. same as but for the presence/absence of the inhibitory codon pairs. spearman’s and kendall’s partial correlation coefficients controlled for gc content, tai, dipeptide content, coding length to assess an association between the fraction of the inhibitory codon pairs and various gene expression variables. same as but for the presence/absence of the inhibitory codon pairs in principle, the observed associations could be due to the inhibitory codon pairs or to other transcript features since several other transcript features were correlated with the content of the inhibitory codon pairs as well as with mrna synthesis/decay rates in some of the datasets. these include guanine-cytosine content, trna adaptation index, which is a metric of codon optimality, the fraction of dinucleotides that are encoded by the inhibitory pairs, and the lengths of coding sequences. to evaluate contributions of these transcript features to the observed associations, we computed spearman’s and kendall’s partial correlation coefficients between the content of the inhibitory codon pairs and mrna decay/synthesis rates controlling for these confounding factors. the analysis led to the following two points. first, for mrna decay rates, the association remained significant when the analysis was individually controlled for gc content, codon optimality, dinucleotide content, and coding lengths with one exception where we analyzed a relationship between the presence/absence of the inhibitory codon pairs and mrna decay rates in the “coller” dataset controlling for codon optimality using the spearman’s method. the association was significant across the analysis methods and datasets when the analysis was controlled for all the covariates. second, for mrna synthesis rates, the association was no longer consistent across the analysis methods/datasets when the analysis was controlled for codon optimality or for all covariates. overall, the results suggest an association of the inhibitory codon pairs with mrna instability on a genomic scale, which appears to be, at least in part, independent of the content of nucleotides, individual constituent codons, and encoded dipeptides as well as of coding lengths. although we also observed an association between the inhibitory codon pairs and mrna synthesis rates, this could be attributed to usage bias in individual codons. the association between inhibitory codon pairs and mrna instability is largely dependent on the correct reading frame the simplest model is that inhibitory codon pairs correlate with mrna decay rates due to their effects on translation elongation and would thus only correlate with mrna decay rates when present in the proper reading frame. alternatively, it remains possible that the hexanucleotides making up inhibitory codon pairs could affect mrna decay rates directly. for example, one possibility is that a subset of the sequences might coincidentally match those recognized by trans-acting factors that promote mrna decay. to address this issue, we examined an association of the occurrence of the hexanucleotide sequences in shifted reading frames as wells as in 3’ untranslated regions with mrna decay rates. in the former analysis, we computationally introduced frameshifts by one or two nucleotides to all orfs and repeated otherwise the same correlation analyses as described above. in the latter analysis, to all orfs, we assigned binary indicators to represent the presence/absence of at least one of the inhibitory codon pairs within 3’ utr based on annotations from previous studies by snyder and colleagues and by steinmetz and colleagues and. the analyses led to the following points. first, the inhibitory codon pairs in the + frame were associated with fast mrna decay rates across the analysis methods/datasets except for the “coller” data. however, the association became inconsistent when gc content, codon optimality, the content of dipeptides encoded by the in-frame inhibitory codon pairs, and coding lengths were individually or simultaneously controlled for. second, the hexanucleotide sequences in 3’ utr were not consistently associated with mrna decay rates.table test for associations of the out-of-frame inhibitory codon pairs with mrna decay rate, protein per mrna, and ribosome occupancy ρ ρ ρ τ τ τ ρ ρ ρ τ τ τ spearman’s partial correlation coefficients controlled for gc content, tai, dipeptide content, and coding length to assess an association between the fraction of hexanucleotide sequences corresponding to the inhibitory codon pairs in the + and + frames and various gene expression variables. p values obtained according to kim and those based on permutation tests are shown. same as but for kendall’s partial correlation coefficients. same as but for the presence/absence of the hexanucleotide sequences. same as but for the presence/absence of the hexanucleotide sequences ρ τ ρ τ spearman’s and kendall’s correlation coefficients to assess an association between the presence/absence of hexanucleotide sequences corresponding to the inhibitory codon pairs in 3’ utr regions. the utr annotations are based on a study by snyder and colleagues. p values obtained according to kim and those based on permutation tests are shown. same as but for utr annotations based on a study by steinmetz and colleagues in sum, the results are largely consistent with the inhibitory codon pairs affecting mrna decay primarily via its effects on translation elongation kinetics. the inhibitory codon pairs cannot fully explain the association of codon optimality with mrna decay the inhibitory codon pairs consist of ten types of codons, all of which are classified as nonoptimal based on tai. therefore, one possibility is that the inhibitory codon pairs could explain the association of codon optimality with mrna decay rates. two observations argue that this is unlikely to be the case. first, for genes lacking the inhibitory codon pairs, codon optimality was still significantly associated with mrna decay rates. second, there was a significant partial correlation of codon optimality with mrna decay rates even when the analysis was individually or simultaneously controlled for gc content, the content of inhibitory pairs, dipeptide content, and coding length.fig. associations of codon optimality with mrna decay rate and translation efficiency for genes lacking the inhibitory codon pairs. a scatterplot comparing tai, a metric of codon optimality, and mrna decay rate in the “cramer 1” data. shown on the top are the spearman’s and kendall’s correlation coefficients and p values. b same as but for the “cramer 2” data. c same as but for the “gresham” data. d same as but for the “coller” data. e same as but for protein abundance per mrna. f same as but for ribosome occupancy ρ τ ρ τ ρ τ spearman’s and kendall’s correlation coefficients to assess an association between codon optimality and various gene expression variables. p values obtained according to kim and those based on permutation tests are shown. same as but for partial correlation coefficients controlled for gc content, the fraction of the inhibitory codon pairs, dipeptide content, and coding length. same as but for partial correlation coefficients controlled for gc content, the presence/absence of the inhibitory codon pairs, dipeptide content, and coding length to examine further whether the inhibitory codon pairs can explain effects of codon optimality on mrna decay rates, we took advantage of synonymous reporter systems used in the previous studies that suggested a mechanistic link between codon optimality and mrna stability. specifically, we compiled and analyzed the content of the inhibitory codon pairs in sequences of reporter constructs used in the studies. the analyses led to the following two points, which imply that the inhibitory codon pairs can promote mrna instability but cannot explain the relationship between codon optimality and mrna stability. first, when comparisons are made within a synonymous group, the higher the content of the inhibitory codon pairs, the faster mrna decay. second, there are multiple examples where synonymous transcripts that differ in codon optimality but not in the fraction of the inhibitory codon pairs exhibit different mrna decay rates. the latter point is particularly important because, if the effect of codon optimality were solely due to the inhibitory codon pairs, the synonymous transcripts lacking the inhibitory codon pairs would show similar mrna decay rates.table number of the inhibitory codon pairs in reporter systems pgk1 lsm8 his3 naa his3 naa naa naa naa naa naa naa naa naa naa ana, not applicable bthe effects on mrna decay rate have been suggested based on mrna abundance overall, these results are consistent with the idea that codon optimality affects mrna stability at least in part independently of the inhibitory codon pairs. the inhibitory codon pairs do not show position effects a previous reporter-based study has shown that a stretch of nonoptimal codons exhibits an increasing destabilizing effect on mrna with an increasing distance from the start codon. based on the assumption that the stretch of the nonoptimal codons causes ribosome queuing along the upstream region, the observation was interpreted to suggest that the higher the number of slow ribosomes on a transcript the less stable the transcript. to examine whether the inhibitory codon pairs have a similar property, we computed spearman’s and kendall’s correlation coefficients between distances of the inhibitory codon pairs from the start codon and mrna decay rates in the genome-wide rna kinetic data. in this analysis, we focused on orfs that contain one and only one of the inhibitory codon pairs. if the inhibitory codon pairs had a stronger effect with an increasing distance from the start codon, the distances would be positively correlated with mrna decay rates. however, we did not observe consistent correlations across the datasets. this suggests that the inhibitory codon pairs are unlikely to cause a long-range ribosome queuing in upstream regions of natural endogenous mrnas.fig. lack of position effect of the inhibitory codon pairs. a scatterplot comparing mrna decay rate in the “cramer 1” data and distances of the inhibitory codon pairs from the start codons contained by the mrnas. shown on the top are the spearman’s and kendall’s correlation coefficients and p values. b same as but for the “cramer 2” data. c same as but for the “gresham” data. d same as but for the “coller” data codon optimality, adjacent codon pairs, and translation efficiency it has been shown that codon optimality and the inhibitory codon pairs are associated with translation efficiency. however, another study did not find a correlation between codon optimality and translation efficiency. to reassess this issue, we selected recently published genome-scale data and examined the association between codon optimality and translation efficiency. translation efficiency, defined as the rate of protein synthesis per mrna, can be obtained by various methods. one among them is to measure protein abundance and normalize it against mrna abundance. this is based on the assumption that the majority of proteins are stable and that protein abundance is largely determined by mrna abundance and protein synthesis rates, which is consistent with a recent protein half-life measurement in s. cerevisiae. another is to use ribosome occupancy on mrna as a proxy for protein synthesis rates. this is based on the assumption that the majority of mrna-bound ribosomes are actively engaged in translation, which is consistent with a recent observation in growing s. cerevisiae cells that ribosome occupancy is highly correlated with the rate of translation initiation. for our analyses, we selected recent quantitative proteomic data by mann and colleagues, which was used in the previous study by grayhack, fields, and colleagues, as well as mrna-seq and ribosome profiling data by weinberg and colleagues. the ribosome profiling data was chosen for two reasons. first, the protocol used for mrna quantification in the studies does not involve poly enrichment, which could result in 3’ bias, an overestimation of mrna abundance of short mrnas, and thus an underestimation of ribosome occupancy of short mrnas. second, the protocol used for quantification of ribosome-protected rna fragments does not involve cycloheximide treatment, which could introduce multiple artifacts. using these datasets, we observed that translation efficiency in both metrics was positively correlated with codon optimality.fig. associations of codon optimality with protein synthesis and decay. a scatterplot comparing tai, a metric of codon optimality, and protein abundance per mrna. shown on the top are the spearman’s and kendall’s correlation coefficients and p values. b same as but for ribosome occupancy. c same as but for protein decay rates we then reexamined the association between the inhibitory codon pairs and translation efficiency in the selected data. for this purpose, we computed spearman’s and kendall’s correlation and partial correlation coefficients controlling for the same set of covariates as we used for the analysis of mrna synthesis/decay rates. the analysis led to the following observations, which is consistent with the previous study. first, the fraction of the inhibitory codon pairs as well as the binary indicator of the presence thereof was associated with low protein abundance per mrna as well as with low ribosome occupancy. second, the association remained significant when the analysis was individually controlled for gc content, codon optimality, dinucleotide content, and coding lengths. third, the association still remained significant when the analysis was controlled for all the covariates. to assess contributions from the hexanucleotide sequences corresponding to the inhibitory codon pairs, we also examined associations of the sequences in + and + frames and in 3’ utrs with translation efficiency. the analyses led to the following three points. first, neither inhibitory codon pairs in the + frame nor those in the + frame were consistently associated with protein abundance per mrna. second, although the content of the inhibitory codon pairs in the + and + frames were consistently associated with low ribosome occupancy, the association became inconsistent when the covariates were controlled for. third, the hexanucleotide sequences in 3’ utrs were not associated with protein abundance per mrna or ribosome occupancy. overall, these analyses confirmed a positive correlation between codon optimality and translation efficiency as well as an association between the inhibitory codon pairs and low translation efficiency, which is at least in part independent of effects of nucleotide content, codon optimality, the content of encoded dipeptides, coding lengths, and corresponding hexanucleotide sequences. the results are consistent with the inhibitory codon pairs affecting protein synthesis primarily via its effects on translation elongation kinetics. the inhibitory codon pairs cannot fully explain the association of codon optimality with translation efficiency having confirmed the associations of codon optimality and the inhibitory codon pairs with translation efficiency, we next examined the possibility that the inhibitory codon pairs can explain the association of codon optimality with translation efficiency. as was seen for mrna decay rates, two observations argue that this is unlikely to be the case. first, for genes lacking the inhibitory codon pairs, codon optimality was significantly associated with the both metrics of translation efficiency, protein abundance per mrna and ribosome occupancy. second, there was a significant partial correlation of codon optimality with mrna decay rates and translation efficiency even when the analysis was individually or simultaneously controlled for gc content, the content of inhibitory pairs, dipeptide content, and coding length. codon optimality, adjacent codon pairs, and protein stability translation kinetics can affect co-translational protein folding, which can, in turn, affect protein function and stability. we thus examined associations between codon optimality or adjacent codon pair content and protein decay rates. for this analysis, we used a genome-wide protein half-life measurement by walther and colleagues. the data was generated via metabolic labeling and thus likely to be more accurate than those generated via methods using translation inhibitors. the analysis led to the following observations. first, codon optimality was negatively correlated with protein decay rates. second, the inhibitory codon pairs were significantly associated with fast protein decay rates. third, however, the association between the inhibitory codon pairs and protein instability became insignificant when tai was controlled for and insignificant or borderline significant when all covariates were controlled for. these results suggest that the inhibitory codon pairs are unlikely to be directly linked to protein stability. analysis of all possible adjacent codon pairs the above analyses are focused on inhibitory codon pairs that have been experimentally characterized previously. however, a recent computational analysis suggests that other codon pairs can be inhibitory to translation. to examine relationships between the adjacent codon pairs and mrna stability in an unbiased manner, we computed spearman’s and kendall’s correlation coefficients between the fraction of each of the possible codon pairs and mrna decay rates in the four datasets. we then compared the correlation coefficients between the inhibitory codon pairs and other pairs. we note that the minimum number of orfs containing a given codon pair with a measured decay rate was as small as eight, which is for the cga di-codon, one of the inhibitory codon pairs, in the “cramer” data. the analysis led to the following two points. first, in three of the four datasets, the correlations with mrna decay rates are higher for the inhibitory codon pairs than for other pairs, which is roughly consistent with the association of the inhibitory codon pairs with fast mrna decay. second, for a substantial number of codon pairs other than the inhibitory codon pairs, their content showed significant positive correlations with mrna decay rates, which may be consistent with those codon pairs acting inhibitory on gene expression as has been suggested in a recent study. overall, the results are roughly consistent with the association of the inhibitory codon pairs with fast mrna decay. however, they also suggest that correlation analyses of individual codon pairs and measured gene expression variables are highly susceptible to experimental noise and/or other confounding factors. discussion we have analyzed relationships between adjacent codon pairs and aspects of gene expression in s. cerevisiae. our results suggest an association of the inhibitory codon pairs with fast rates of mrna decay. the association is not solely due to gc content, the content of individual codons, the content of encoded dipeptides, or coding lengths. moreover, there is no consistent association of the hexanucleotide sequences corresponding to the inhibitory codon pairs in shifted frames or in non-coding regions with mrna decay rates, which is largely consistent with the inhibitory codon pairs affecting mrna stability at least in part via their effects on translation elongation kinetics. although we also observed associations of the inhibitory codon pairs with mrna synthesis rates and protein instability, they could be attributed to bias in the individual constituent codons. it has been speculated that adjacent codon pairs rather than individual codons may largely underlie the relationship between codon optimality and translation efficiency. our results suggest that neither the effects of codon optimality on translation efficiency nor those on mrna stability can be explained by the inhibitory codon pairs. this observation is more consistent with the model that a sum of subtle effects from individual codons can exert a large impact on gene expression. however, it is still possible that codon pairs rather than individual codons largely define the relationship of codon composition with translation efficiency and mrna stability since there are likely to be other inhibitory codon pairs than those in the high confidence set obtained by the previous study. the genome-wide associations of codon optimality and the inhibitory codon pairs with mrna stability and translation efficiency observed by this work and other studies suggest two non-mutually exclusive possibilities. the first possibility is that direct mechanistic links largely underlie the association of the codon composition, which can affect speeds and/or accuracy of translation elongation, with mrna stability and translation efficiency. several non-mutually exclusive mechanisms can link slow/inaccurate translation elongation to fast mrna decay rates and/or low protein synthesis rates. first, slow translation elongation can cause ribosome queuing in upstream regions, which could, in turn, interfere with translation initiation. second, slow translation elongation may cause ribosome drop-off. consistent with this idea, the cga di-codon, one of the inhibitory codon pairs, has been linked to a quality control mechanism that can detect stalled ribosomes and cause abortion of translation. third, slow translation elongation may be sensed by a mechanism that can modulate translation initiation rates. indeed, a recent study in s. cerevisiae has raised the possibility that the dead-box protein dhh might play a central role in such a mechanism. fourth, nonoptimal codons and/or the inhibitory codon pairs may result in translation repression and mrna instability by compromising translation fidelity. for example, it is possible that nonoptimal codons and/or the inhibitory codon pairs increase erroneous translation frameshifts. although there is no overlap between the inhibitory codon pairs and previously identified frameshift-inducing sequences, it remains to be determined whether any of the inhibitory codon pairs tend to introduce translation frameshifting. in most cases, translation frameshifts would result in a premature translation termination event at a stop codon in the incorrect frame, which would in turn cause repression of translation initiation and nonsense-mediated mrna decay. such a mechanism can partly underlie the association of nonoptimal codons and inhibitory codon pairs with mrna instability. indeed, a recent study suggests that mrnas with high content of nonoptimal codons tend to undergo nonfunctional translation frameshifts and, subsequently, nmd. the second possibility is that the genome-wide associations between codon composition, translation efficiency, and mrna stability are largely due to co-evolution rather than to a mechanistic link, the latter of which has been suggested by studies using artificial reporter systems. that is, nonoptimal codons and the inhibitory codon pairs may be simply avoided in natural endogenous genes that are highly expressed and efficiently translated. then, their primary function may be to modulate local translation elongation speed and thereby regulate other processes, such as co-translational folding, which may be largely restricted to a situation where slow ribosomes do not negatively impact overall translation efficiency. consistent with this view, some studies in unicellular organisms suggest that under physiological conditions translation initiation but not translation elongation mainly defines the rate of protein synthesis. moreover, another study did not find corresponding changes in translation efficiency upon genetic manipulation of trna and thus codon optimality. clearly, further investigation will be needed to rigorously evaluate these two possibilities concerning the codon-mediated gene regulation.  this study suggests genome-scale associations of the inhibitory codon pairs with mrna decay and translation efficiency, which, in turn, suggest another layer of complexity in the codon-mediated gene regulation. an important future goal will be to understand whether and how the inhibitory codon pairs mechanistically inhibit protein synthesis and elicit mrna instability. 
in saccharomyces cerevisiae, a recent study has identified adjacent codon pairs that mediate strong inhibition of translation elongation. second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mrna stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. finally, although the inhibitory codon pairs are associated with slower mrna synthesis and protein instability, the associations can be attributed to usage bias in individual codons. evidence in diverse organisms suggests that codon optimality is a major determinant of mrna translation and degradation. codon optimality is thought to act by modulating the efficiency of ribosome elongation. first, the inhibitory codon pairs are significantly associated with faster mrna decay. however, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown. to gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in s. cerevisiae. the association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. this study suggests an association of inhibitory codon pairs with mrna stability and thus another layer of complexity in the codon-mediated gene regulation. codon pairmrna translationmrna stabilityhttp://dx.doi.org/ /100000011howard hughes medical instituteissue-copyright-statement© the author 2017
 evidence in diverse organisms suggests that codon optimality is a major determinant of mrna translation and degradation. codon optimality is thought to act by modulating the efficiency of ribosome elongation. in saccharomyces cerevisiae, a recent study has identified adjacent codon pairs that mediate strong inhibition of translation elongation. however, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown.  to gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in s. cerevisiae. our analysis revealed the following points. first, the inhibitory codon pairs are significantly associated with faster mrna decay. the association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mrna stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. finally, although the inhibitory codon pairs are associated with slower mrna synthesis and protein instability, the associations can be attributed to usage bias in individual codons.  this study suggests an association of inhibitory codon pairs with mrna stability and thus another layer of complexity in the codon-mediated gene regulation. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords codon pairmrna translationmrna stabilityhttp://dx.doi.org/ /100000011howard hughes medical instituteissue-copyright-statement© the author 2017 