Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genome-Wide and Experimental Resolution of Relative Translation Elongation Speed at Individual Gene Level in Human Cells

  • Xinlei Lian ,

    Contributed equally to this work with: Xinlei Lian, Jiahui Guo

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Jiahui Guo ,

    Contributed equally to this work with: Xinlei Lian, Jiahui Guo

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Wei Gu,

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Yizhi Cui,

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Jiayong Zhong,

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Jingjie Jin,

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Qing-Yu He ,

    zhanggong@jnu.edu.cn (GZ); tongwang@jnu.edu.cn (TW); tqyhe@jnu.edu.cn (QYH)

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Tong Wang ,

    zhanggong@jnu.edu.cn (GZ); tongwang@jnu.edu.cn (TW); tqyhe@jnu.edu.cn (QYH)

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯
  • Gong Zhang

    zhanggong@jnu.edu.cn (GZ); tongwang@jnu.edu.cn (TW); tqyhe@jnu.edu.cn (QYH)

    Affiliation Institute of Life and Health Engineering, Key Laboratory of Functional Protein Research of Guangdong Higher Education Institutes, Jinan University, Guangzhou, China

    ⨯

Abstract

In the process of translation, ribosomes first assemble on mRNAs (translation initiation) and then translate along the mRNA (elongation) to synthesize proteins. Elongation pausing is deemed highly relevant to co-translational folding of nascent peptides and the functionality of protein products, which positioned the evaluation of elongation speed as one of the central questions in the field of translational control. By integrating three types of RNA-seq methods, we experimentally and computationally resolved elongation speed, with our proposed elongation velocity index (EVI), a relative measure at individual gene level and under physiological condition in human cells. We successfully distinguished slow-translating genes from the background translatome. We demonstrated that low-EVI genes encoded more stable proteins. We further identified cell-specific slow-translating codons, which might serve as a causal factor of elongation deceleration. As an example for the biological relevance, we showed that the relatively slow-translating genes tended to be associated with the maintenance of malignant phenotypes per pathway analyses. In conclusion, EVI opens a new view to understand why human cells tend to avoid simultaneously speeding up translation initiation and decelerating elongation, and the possible cancer relevance of translating low-EVI genes to gain better protein quality.

Author Summary

In protein synthesis, ribosome assembles to mRNA to initiate translation, followed by the process of elongation to read the codons along the mRNA molecule for polypeptide chain production. It is known that slowing down the elongation speed at certain regions of mRNA is critical for the correct folding of numerous proteins—the so-called “pause-to-fold”. However, it has been an open question to evaluate elongation speed under cellular physiological conditions in genome-wide scale. Here, we used three types of next-generation sequencing approaches to experimentally and computationally address this question. With a new relative measure of elongation velocity index (EVI), we successfully distinguished slow-translating genes. Their protein products are more stable than the background genes. We found that different cell types tended to have distinct slow-translating codons, which might be relevant to the cell/tissue specific tRNA composition. Such elongation deceleration is potentially disease-relevant: cancer cells tend to slow down numerous cancer-favorable genes, and vice versa. Furthermore, we justified that translation initiation and elongation are evolutionarily synergistic as no gene with both high initiation efficiency and low elongation speed was observed: that would cause a traffic jam of ribosomes that should be maximally avoided per evolution.

Citation: Lian X, Guo J, Gu W, Cui Y, Zhong J, Jin J, et al. (2016) Genome-Wide and Experimental Resolution of Relative Translation Elongation Speed at Individual Gene Level in Human Cells. PLoS Genet 12(2): e1005901. https://doi.org/10.1371/journal.pgen.1005901

Editor: Ronald B. Gartenhaus, University of Maryland Medical School, UNITED STATES

Received: October 22, 2015; Accepted: February 5, 2016; Published: February 29, 2016

Copyright: © 2016 Lian et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All sequencing data files are available from GEO database (accession: GSE46613)

Funding: This work was supported by the National Basic Research Program “973” of China (2014CBA02000 to TW and 2011CB910700 to QYH), the National Natural and Science Foundation of China (81372135 to TW; 81322028 and 31300649 to GZ), the Institutional Grant of Excellence of Jinan University, China (50625072) to GZ as well as the Guangdong Natural Science Foundation (2014A030313369) and Guangzhou Key Subject (201504291059373) to TW. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Protein synthesis is a collective outcome of different mechanisms of translation control, including: i) translation initiation, namely the assembly process of methionyl-tRNA loaded and elongation-competent 80S ribosomes at mRNA start (AUG) codons (reviewed in ref. [1–4]), which determines the fraction of mRNA molecules that can be translated [5, 6]; ii) translation elongation, which determines the translating speed of ribosome(s) post-initiation on a single mRNA molecule (reviewed in ref. [7, 8]); and iii) termination that allows the reuse of ribosomes.

Indeed, elongation speed with its biological relevance is one of the most well-known challenges in understanding translational control, especially regarding cells under physiological conditions. If not considering pre-mature termination (ribosome drop-off) and frameshifting, elongation speed is in general independent from the abundance of protein products, at least in mouse embryonic stem cells [9]. This is comparable with numerous reports showing that initiation is the rate-limiting step of translation in eukaryotic cells [5, 10, 11] (reviewed in ref. [2, 7]), which we have computationally modeled the length-dependent initiation by considering the circulation time and its impact on maintaining a functional proteome [12]. Regarding elongation, we previously experimentally validated that the major determinants of elongation speed should include codon selection and cognate tRNA abundance [13]. Interestingly, we and others have reported that the abundance of tRNAs decoding different codons may vary more than one order of magnitude, thus the elongation speed should be non-uniform even on a single mRNA molecule [13–17].

The biological significance of elongation speed control has been emphasized to understand human disease such as cancer, hemophilia and hyperactivity disorder (reviewed in ref. [18]). Kimchi-Sarfaty et al have found that a synonymous mutation of the multi-drug resistance 1 gene (MDR1) can cause structural and functional alteration of the gene product that may be relevant to the timing of co-translational folding as hypothesized by the authors [19]. This notion was supported by our previous experimental findings on the translational pausing caused by clusters of slow-translating codons coordinates protein synthesis and co-translational folding that are critical for protein quality [13] (reviewed in ref. [7, 20–22]).

Although several studies from us and others have monitored the elongation speed of single genes [13, 23–25], the genome-wide elongation speed evaluation on eukaryotic cells under physiological conditions is still an open question. Since decades ago, mathematical indices, including codon bias index (CBI), codon adaptation index (CAI), effective number of codons (Nc) and tRNA adaptation index (tAI) have been widely employed to estimate translation elongation efficiency (reviewed in ref. [26]). Most of these indices are calculated using codon usage and tRNA gene copy number of the whole genome, which persist the same in a certain organism. Here, codon usage is a term with strict definition: the occurrence of a certain codon among all codons used in a certain organism [27]. Indeed, the tRNA concentrations for each codon, including the ratio of cognate to near-cognate, is an important determinant of the elongation speed [28, 29]. However, the tRNA concentration correlates poorly to these indices and can vary dramatically in different tissues and cell types, serving as a primary causal factor of tissue-specific translation profiles (reviewed in ref. [20, 26]). With this regard, we previously established an algorithm to use tRNA abundance to predict the elongation speed and translational pausing sites, which was experimentally validated both in single genes and at genome-wide level [13, 14, 30]. Unfortunately, absolute tRNA abundance information at anticodon resolution is currently unavailable for higher eukaryotes due to the technological hindrance to resolve the high complexity and homology of tRNA molecules (reviewed in [20]). Other than tRNA concentration, it has been investigated whether the mRNA secondary structure, steric effect of tRNAs and amino acid charge are major determinants of elongation velocity. In the view of computational biology, the mRNA secondary structure may affect the elongation speed [16], supported by the pseudoknot-directed translational frameshifting as an example [31]. In contrast, other evidence acquired by the pulse-chase labeling and single ribosome monitoring suggested that stable mRNA secondary structures do not cause any translation delay [24, 25]. Despite of diversified structures of tRNAs, it was proposed that all of the aminoacyl-tRNAs are selected uniformly on the ribosome, suggesting that steric effects have minimum influence on the elongation speed [32]. Interestingly, repeats of positively charged amino acids have been found to cause translational pausing [16, 33, 34]; furthermore, singlets of such amino acids do not necessarily cause translational deceleration [35]. In addition, other specific sequences, such as anti-Shine-Dalgarno sequence and stalling nascent peptides, have been found to pause the ribosomes regarding some specific genes [36, 37].

Due to such complex scenarios, to directly assess the translation elongation speed in a high-throughput manner is necessary. Ribosome profiling provides a detailed snapshot of ribosome occupancy [38], making it feasible to study the ribosome density profile of translation ([9] and reviewed in ref. [8]). Ingolia et al monitored the progression of the average profiles of ribosome footprints (RFPs) and revealed an average translation elongation speed of 5.6 codons/sec in mouse embryonic stem cells; however, this measurement of ribosome elongation has a 60-s delay caused by the harringtonine treatment [9].

We previously reported a strategy to combine the full length sequencing on ribosome nascent-chain complex (RNC) bound mRNA (RNC-mRNA) and total mRNA for the global translation initiation investigation [6, 39]; we showed that the translation ratio (TR, abundance ratio of RNC-mRNA/mRNA for a certain gene) can properly reflect cellular phenotypes. In this study, we integrated three types of current RNA-seq strategies, including mRNA sequencing (mRNA-seq), full-length RNC-mRNA sequencing (RNC-seq) and ribosome profiling (Ribo-seq) (Fig 1A). As an outcome, we resolved global elongation speed by an Elongation Velocity Index (EVI) at individual gene level in human normal and cancer cells under physiological conditions. This allowed us to distinguish slow-translating genes and codons in different human cell lines, respectively. Furthermore, our results favored the hypothesis on the cancer relevance of co-translational folding by providing the experimental and computational evidence on a genome-wide scale.

thumbnail
Download:
Fig 1. Measurement of TR and EVI.

(A) Schematic workflow of mRNA-seq, RNC-seq and Ribo-seq of the same batch of cultured cells. (B) Contribution of translation initiation efficiency and elongation velocity to the ribosome density.

https://doi.org/10.1371/journal.pgen.1005901.g001

Results

Estimation of relative translation velocity by the Elongation Velocity Index

Using reads per kilo base per million (rpkM) as unit, the abundance of mRNA (M), RNC-mRNA (C) and RFP (F) are length-independent. Therefore, the RNC-mRNA ribosome density (Density), which is defined here as F/C, and TR that is defined as C/M [6] can be compared between different genes. Here, RNC-mRNA ribosome density is not a constant, it varies drastically on genome-wide scale either in a single cell type (S1A Fig) or across different cell types (S1B Fig). As the translation initiation is the rate-limiting factor of the entire translation process in eukaryotes [5], TR is a relative measure of the translation initiation efficiency in eukaryotes [6, 12, 39]. Here, we defined the translation initiation efficiency as the ratio of the amount of mRNAs that are being translated to the amount of transcribed mRNAs of a certain gene [6]. RNC-mRNA ribosome density is proportional by the ratio of translation initiation efficiency and elongation speed (Fig 1B). Here, we define the “Elongation Velocity Index” (EVI) as a relative measure of elongation speed to make (1)

Thus: (2)

The TR and EVI values of all of the tested cell lines are provided in S1 Table.

Genome-wide resolution of translation elongation speed at individual gene level under physiological condition

Due to the possible association of elongation speed and protein degradation, plus the protein turnover rate of HeLa cells has been investigated intensively [40], we used this cell line for the EVI evaluation and verification as an initial step. We detected a total of 10,837 genes in HeLa cells that could be quantified from all of the three types of RNA-seq (Fig 1 and S2 Table). Although TR and EVI emphasize on different stages of translational control, we found they significantly correlated to each other, with the Spearman R (Rs) = 0.62 (P<10−38; Fig 2A). This suggests that evolutionarily synergistic roles of translational control may exist to prevent certain detrimental effects in translation. To further probe such roles, we computed the 99% Hotelling's T2 confidence ellipse and genes outside the confidence ellipse were considered as outliers, which were then subjected to clustering analysis based on the Euclidean distance and Ward’s linkage (Fig 2A). Interestingly, we observed two polarized clusters of outlier genes, namely the low EVI cluster (red dots) and the high TR cluster (blue dots). Considering the gene distribution and empirical cumulative density function (CDF) of TR and EVI, we defined a grey quadrant locating genes with the top 1% of TR and the smallest 1% of EVI, simultaneously (Fig 2A and S2 Fig). No gene was found to distribute in such a grey quadrant for HeLa cells (Fig 2A and S2 Fig) as well as HBE, A549 and H1299 cells (Fig 2B). In addition, similar polarized distribution of genes were confirmatively observed in the human lung cell lines (Fig 2B), which helped to rule out the cell-specific bias. Excluding the RFP reads of upstream ORFs (uORFs) almost led no change to the EVI values (S3 Fig).

thumbnail
Download:
Fig 2. Two mutually exclusive and polarized translation control modes.

(A) The TR and EVI of genes in HeLa cells. The green ellipse indicates the 95% confidence ellipse. Genes outside of this ellipse were clustered in two major clusters (red and blue dots) based on their Euclidean distances. The grey region denotes the region with high TR (genes with top 1% TR) and low EVI (genes with lowest 1% EVI), simultaneously. Rs = Spearman R. (B) The TR and EVI of genes of HBE, A549 and H1299 lung cell lines, respectively. (C) Plot matrix of mutual correlation of the TR (lower half triangle) and EVI (upper half triangle) among the four analyzed cell lines. (D) The heatmap of the Spearman R of the panel (C). All of the P-values are less than 10−38. (E,F) Cluster analysis of the TR (E) and EVI (F) of the four cell lines. The lung-derived cell lines are indicated by a green bar.

https://doi.org/10.1371/journal.pgen.1005901.g002

The mutual correlation coefficient of the analyzed four cell lines ranges from 0.42 to 0.83 (all P < 10−38; Fig 2C and 2D) based on TR and EVI, respectively. In comparison, such correlation coefficients calculated based on mRNA and RNC-mRNA ranged from 0.77~0.87, while the RFP of HeLa cells are quite different than the lung-derived cells (Rs = 0.45~0.41; all P < 10−38) (S4 Fig). The hierarchical cluster analysis of the TR and EVI of each gene in the four cell lines showed that the three lung-derived cells are clustered together and deviate from HeLa cells (Fig 2E and 2F), suggesting a tissue-specific pattern both on translation initiation efficiency and elongation speed. In addition, we found 49 overlapping low-EVI genes and 4 overlapping high-TR genes among the three lung-derived cell lines (S4A and S4B Table). In contrast, we observed 4 overlapping low-EVI and no overlapping high-TR genes among all of the four analyzed cell lines (S4C Table). This indicates that the low-EVI genes and high-TR genes are tissue-specific, which corresponds to our previous findings that high-TR genes reflect the specific cellular phenotypes and organ origins [6]. Gene ontology (GO) enrichment analyses showed that the low-EVI genes in four cell lines shared some universal and housekeeping biological processes, e.g. cellular component organization (S5 Fig). At the same time, the GO enrichments also exhibited remarkable tissue specificity (S5 Fig).

Notably, D = TR/EVI spans a wide distribution of 4~5 log units in each cell line (S1A Fig). This range is similar to the span of EVI distribution (Fig 2A and 2B). The correlation of D in HeLa and lung cells are at <0.4 level (S1B Fig), suggesting a high variability of RNC-mRNA ribosome densities across the cells. In comparison, higher correlation of densities were detected among the lung-derived cell lines (S1B Fig), which echoed the tissue specificity of translation elongation. Thus, experimental measurement is necessary to test the EVI and TR correlations.

We then calculated CAI [27], CBI [41] and Nc [42, 43] in all of the analyzed cell lines, respectively. We found that EVI has only minor correlations with CAI (|Rs| < 0.13; all P<10−5; S6 Fig) and CBI (|Rs| < 0.11; all P<10−8; S7 Fig), respectively. Interestingly, EVI has significant negative correlation with Nc (Rs = -0.15~-0.46; all Ps < 10−38; S8 Fig).

Elongation deceleration and protein stability

We have previously validated the phenotypic relevance of high TR [6]; therefore, we next focused on examining whether the slow-translating genes have translational pausing sites and whether they are more stable. Here, we adopted the definition that the highest RFP peaks that higher than 10-fold of average count are translational pausing sites [9].

We found that the RFP distribution of low-EVI genes in HeLa cells, such as NES, showed significant translational pausing sites, and their highest RFP peaks predominantly occur more than 240nt after the start codons (Fig 3A). It is known that the pausing sites after the 240nt threshold predominantly coordinate protein synthesis and co-translational folding [13, 14, 44]. In contrast, some high TR genes in HeLa cells, such as IFIT2, has the highest RFP peaks located next to the respective start codons, indicating its high efficiency of translation initiation and fewer translational pausing sites (Fig 3B). More examples can be found in the S9 Fig. We statistically confirmed these phenomena observed on the randomly chosen genes (Fig 3C). Approximately half of the high TR genes have their highest RFP peaks prior to 240nt away from the start codons, significantly different from the low EVI genes (P = 0.0043~1.8×10−8, Fig 3C). Notably, all of the highest RFP peaks in either the low EVI genes or the high TR genes in the four analyzed cell lines have >10-fold more coverage than the average, thus can be considered as true translational pausing sites but not the sequencing bias as explained by Ingolia et al [9].

thumbnail
Download:
Fig 3. The relevance of low EVI genes to transient translational pausing and protein folding.

(A, B) The RFP coverage of a representative gene with low EVI (A) and high TR (B) in HeLa cells, respectively, are shown with color bars. The highest RFP peaks are indicated by stars. The full length mRNAs are shown in thin lines, while their CDS regions are marked with thick lines. For more examples in all four analyzed cell lines, please refer to S6 Fig. (C) The percentage distribution of the highest RFP peak of each gene relative to the start AUG codon. The highest RFP peaks more than 400 nucleotides are counted in the bin of 400. The P-values of Kolmogorov-Smirnov test between the low EVI genes and high TR genes as well as all genes in each cell line are indicated. The grey dashed lines mark the 240nt position. (D) Correlation between CDS length and EVI in the four cell lines. (E, F) The distribution of genes in HeLa cells according to their protein degradation rate constant (kdeg). Fractions of the genes with low EVI (E) and high TR (F) were indicated in red and blue bars, respectively. For comparison, the fraction distributions of total genes are shown in grey bars as background.

https://doi.org/10.1371/journal.pgen.1005901.g003

In addition, we found that EVI negatively correlated to CDS length for the genes of the four cell lines, respectively, suggesting that longer genes are generally translated slower in these human cells (Fig 3D). This coincides with our previous findings in bacteria that longer proteins possess more translational pausing sites to overcome the folding complexity of multiple structural domains [5, 13, 45]. Thus, these results independently favored the validity of using EVI to resolve gene elongation speeds.

We next acquired the protein degradation rate constant data of HeLa cells from the report by Cambridge et al [40].We found that the proteins encoded by the low EVI outlier genes had 1.3-fold less degradation rate in average than the proteome background (single-tailed t-test, P = 0.029, tested for the null hypothesis that low EVI genes have no less degradation rates than the background genes) in HeLa cells (Fig 3E). As a negative control, the proteins encoded by the high TR outlier genes share similar degradation rates to the proteome background (single-tailed t-test, P = 0.10, tested for the null hypothesis that high TR genes have no less degradation rates than the background genes) (Fig 3F). This suggests that the translational pausing sites of the genes with slow elongation speeds may significantly facilitate the protein folding, thus increases their stability.

Codon preference of low EVI genes in human cells

We next tried to assess the potential causal factor of the slow translation of the low EVI genes. Although the role of mRNA structure in affecting the elongation speed was under strong debate [16, 24, 25, 31], we tried to examine such an influence with our data. We calculated the mRNA fold energy per base for low-EVI and high-TR genes, respectively, using the MFOLD algorithm [46]. In each of the four cell lines, the low-EVI genes were not more stably folded than the high-TR genes (one-tailed KS-test, S10 Fig), showing that the mRNA structure stability did not play a significant role in slowing down the elongation. Thus, we deduct that the amino acid content, such as positively charged amino acids, may be a factor. We found that the fraction of positively charged amino acids (lysine, arginine and histidine, PCAAs) in low EVI genes was close to that in high TR genes, with the maximum difference less than 1.5%. In HBE and A549 cells, low EVI genes contained slightly more PCAAs than high TR genes, while vice versa in HeLa and H1299 cells (Table 1). In the 4 tested cell lines, we found that their low-EVI genes all had significantly different fraction of PCAAs than the genome average (Table 1). This indicated that the positive charge of the nascent peptides may have effects, but it is unlikely to be the primary factor to slow down the translation of low EVI genes, at least in the four analyzed cell lines of this study.

thumbnail
Download:
Table 1. Fraction of positively charged amino acids (lysine, arginine and histidine) in low EVI genes and high TR genes.

https://doi.org/10.1371/journal.pgen.1005901.t001

We previously reported in bacteria that translation pausing sites are created by clustered slow-translating codons in the CDS that paired to low-abundance tRNA species [13, 14]. Here, we further tried to identify the slow-translating codons in these human cells that may serve as a causal factor of the translational pausing sites in the low EVI genes, respectively. For comparison across different cell lines, we quantified the codon preference by the Preference Score (PS) of low EVI genes (PSLow-EVI) in each cell line based on the RSCU of the low EVI genes vs. background genome. The higher the PS of a codon is, the more preferentially used in the low EVI genes than background genome.

We found that the low EVI genes in all of the four tested cells showed a remarkable codon preference for most of the amino acids that was reflected by PSLow-EVI (Fig 4A). In a general view, the composition of preferred and non-preferred codons in HeLa cells were different from the lung cells with a tissue-specific manner (Fig 4A). For example, the GGT and GGA codons that encodes glycine are disfavored in low EVI genes in HeLa cells (P = 0.03, Chi-squared test), but are favored in HBE, A549 and H1299 cells (P < 10−38 in all three cases). The GGG codon is favored in low EVI genes in HeLa cells (P = 0.02, Chi-squared test), but is disfavored in HBE, A549 and H1299 cells (P < 10−38 in all three cases) (Fig 4A). We also detected by the PSLow-EVI analysis that certain synonymous codons are not necessarily and preferentially used. In HeLa cells, the three codons encoding isoleucine are very different in PSLow-EVI values: ATA is highly preferred in the low EVI genes (P = 0.0003, Chi-squared test) (Fig 4A). In contrast, the codons encoding tyrosine and lysine showed minimal preference in the low EVI genes (Fig 4A).

thumbnail
Download:
Fig 4. Codon preference of the slow-translating genes.

(A) Codon preference of low EVI genes in the four cell lines based on Preference Score analysis. For details, please refer to Materials and Methods section. Favored and disfavored codons are shown in red and grey bars, respectively. (B) Mutual correlation of PSLow-EVI of codons among all four analyzed cell lines. (C) Percentage of variance explained in the principle component analysis (PCA) on RSCU. (D) Principle component (PC) based clustering analysis. The PCs with the summed explained variance greater than 80% were taken for clustering analyses based on the standardized Euclidean distances. For each cell line, four groups of genes were analyzed. They were low EVI genes, high TR genes, all genes and a random subset of genes that were inside the 95% confidence ellipse (Fig 2A and 2B). The clusters with multiple gene groups are indicated by grey ellipses.

https://doi.org/10.1371/journal.pgen.1005901.g004

To quantify the tissue-specific codon preference, we compared the codon PSLow-EVI values among all of the tested cell lines (Fig 4B). The lung HBE, A549 and H1299 cells showed strong correlation to each other (Rs = 0.97~0.99), showing that they have significantly similar codon preference in their low EVI genes. On the contrary, the codon preference in HeLa cells are almost not correlated to the three lung-derived cells (Rs = 0.11~0.17, P = 0.18~0.41). These results confirmed the tissue-specific codon preference shown above.

To verify whether PSLow-EVI can distinguish the low EVI genes in an unbiased manner, we examined the codon preference of four groups of genes in each cell line. They were low EVI genes, high TR genes, all genes, and a random subset of genes that locate inside the 99% confidence ellipse in the TR vs. EVI plots (Fig 1C and 1D). We performed principle component analysis (PCA) on the RSCU of all 61 sense codons in each cell and identified the principle component(s) (PC) that can collectively explain more than 80% of the variance (Fig 4C). In HeLa cells, the first two PCs explained 95.7% of the variance, while in the three lung-derived cells, the first PCs alone explained 95.2~98.1% of the variance, respectively (Fig 4C). We next performed hierarchical clustering with the distance metric of standardized Euclidean distance. Algorithm for computing distance between clusters was unweighed average distance (UPGMA). We found that low EVI genes were clustered as an independent group in 2-dimensional (for HeLa cells) or 1-dimensional (for HBE, A549 and H1299 cells) spaces, respectively (Fig 4D). Thus, the low EVI genes have distinct codon preference to slow down their translation in all of the cells examined in this study.

To further rule out the experimental and/or computational bias, we performed similar analyses on the codon preference of high TR genes (PSHigh-TR) (Fig 5A). The PSHigh-TR is also cell-specific, indicating that the tRNA concentration is cell-specific either (S11 Fig). Since the high TR genes are translated faster, we expected a different codon preference of the high TR genes than low EVI genes. Indeed, in HeLa and HBE cells, the PSHigh-TR and PSLow-EVI have no significant correlation (P = 0.53~0.71, Fig 5B). We noted that PSHigh-TR correlates to PSLow-EVI in A549 cells (Rs = -0.78, Fig 5B) and H1299 cells (Rs = 0.89, Fig 5B), respectively; however, the distribution of PSHigh-TR is significantly different from that of PSLow-EVI, with P = 6.99×10−6 and 2.20×10−4 in the two cell lines, respectively (KS-test on distribution).

thumbnail
Download:
Fig 5. Codon preference of the high TR genes.

(A) Preference Score of high TR genes (PSHigh-TR) of the four cell lines. (B) Correlation of the PSHigh-TR and PSLow-EVI.

https://doi.org/10.1371/journal.pgen.1005901.g005

Furthermore, we observed in all of tested cell lines that there was no significant correlations between PS (either PSLow-EVI or PSHigh-TR) and codon usage (S12A Fig), tAI (S12B Fig) or mean codon elongation time (S13 Fig), respectively.

Slow-translating genes and malignant phenotypes

We previously compared lung cancer A549 and H1299 cells with normal HBE cells, respectively, in terms of cancer phenotypes, translation initiation efficiency and proteome [6, 47]. To understand the biological relevance of EVI, we further examined with these cell lines to answer whether the slow-translating gene group is biased on regulating cancer-favorable phenotypes.

First, we found that the TR change and EVI change correlate to each other in A549/HBE (Rp = 0.62; P < 10−38; Fig 6A) and H1299/HBE (Rp = 0.41; P < 10−38; Fig 6B) comparisons, respectively. This echoes the comparisons on the absolute TR and EVI changes (Fig 2B). However, no significant correlation between TR and EVI changes could be found for the EVI up-regulated genes, and only weak TR-EVI change correlation could be found for the EVI down-regulated genes (S14 Fig), indicating that the translation elongation speed is largely independently regulated for these genes.

thumbnail
Download:
Fig 6. Relative EVI analysis and cancer relevance.

(A, B) Correlation of the relative changes of TR and EVI comparing cancer with normal cells. A549/HBE (A) and H1299/HBE (B) results were shown, respectively; the genes with relative EVIs of greater than ±10 folds were indicated by blue dots. (C,D) Comparative proteomics analyses on the low-EVI genes comparing A549 cells with HBE cells. Relative abundance distributions of total proteins (C) and the unfolded proteins (D) are shown. (E, F) Top canonical pathway analysis of low relative EVIs by IPA. Low relative EVI genes (<10 folds) in A549/HBE (E) and H1299/HBE (F) comparisons were subjected to IPA analysis. The top canonical pathways with P < 0.001 (Fisher’s exact test provided by IPA indicated by the threshold line in orange) regulated by of these genes were shown respectively. The complete gene lists for IPA are included in the S2 Table with HGNC gene names.

https://doi.org/10.1371/journal.pgen.1005901.g006

Next, we tried to answer whether the elongation deceleration could significantly affect the protein abundance of these EVI down-regulated genes. Using our previously published comparative proteomics data [6], we found that the protein abundance of EVI down-regulated genes (A549 compared to HBE cell line) obeyed a normal-like distribution, with no significant difference from the proteome background (P = 0.82, two-tailed KS-test) (Fig 6C). Such results suggested that the protein quality of EVI down-regulated genes may be increased by slowing down elongation, which led to no significant impact on their protein abundance. In addition, we have recently reported that A549 cells have more abundant detergent-insoluble proteins than HBE cells [48], suggesting that lung cancer cells are exposed to more severe protein folding problems than normal cells in genome-wide scale. With the comparative proteomics data, we found that the relative abundance distribution of the EVI down-regulated gene product in the unfolded protein fraction had no significant difference comparing A549 with HBE cells (P = 0.85, two-tailed KS-test; Fig 6D). This result indicated that the EVI down-regulation may serve as an important factor to facilitate protein folding in cancer cells.

With IPA knowledgebase, we analyzed the genes with more than 10-fold down-regulation in their EVIs comparing lung cancer cells with normal cells (S5 Table). We found that A549 cells were decelerating the elongation of genes with low EVIs that is significantly relevant to maintain the EGF signaling pathway (P = 1.77×10−4; Fig 6E), which plays important roles in the proliferation and differentiation of lung cancer cells; in addition, systemic lupus erythematosus signaling was also significantly associated (P = 9.07×10−4; Fig 6E). The elongation speed regulation of H1299 cells focuses on the 14-3-3 (P = 1.73×10−5; Fig 6F; S15 Fig), IL-12 (P = 5.07×10−4; Fig 6F) and remodeling of epithelial adherens junctions (P = 8.76×10−4; Fig 6F). These pathways are known to be essential for the survival, active transcription and migrations of cancer cells as suggested by IPA. For example, the protein of 14-3-3 is highly active in lung and breast cancers interacting with AKT and MAPKs (S15 Fig). As a negative control, we found not significant pathways were enriched by analyzing the genes whose EVI were up-regulated for more than 10 folds in cancer cells.

On the contrary, we found that cancer cells tended to accelerate the translation elongation speed of tumor suppressor genes (TSGs). We took the identified TSGs in the TCGA lung adenocarcinoma datasets from TSGene 2.0 Database (http://bioinfo.mc.vanderbilt.edu/TSGene/download.cgi). The EVIs of these TSGs were elevated 1.89 and 1.75-fold in average in A549 and H1299 cells, respectively, as compared with HBE cells. Furthermore, these up-regulations were more significant in the more malignant H1299 cells (P = 0.02, single-tailed Mann-Whitney U-test against the genome-wide background) than in the less malignant A549 cells (P = 0.34, as above).

Discussion

In this study, we realized the evaluation of global elongation speed at individual gene level for human cells under physiological conditions. Together with our previously defined TR analysis for translation initiation [6], we proposed a two-dimensional fine tuning of translational control to form a translation initiation preference on genes with high phenotype relevance, as well as to adjust translation elongation speed that is relevant to protein quality control. Mechanistically, we discerned the codon preference, especially the slow-translating codons, to explain the slowly translated genes in different human cell types. Therefore, we provided the genome-wide evidence, for the first time, to favor the importance of elongation speed deceleration for the maintenance of the malignant phenotypes of cancer cells at steady-state.

Indeed, using EVI to assess the translating speed is the basis and one of the central messages of this study. We tried to provide evidence at various levels, considering multiple independent factors to avoid biased conclusions. First, we demonstrated with RFP coverage analysis that the low EVI genes have significant translational pausing sites that locate mostly after 240 nt. Favorably, we previously found that most of folding-relevant translational pausing sites occur after 240 nt (~10kDa nascent peptide) in prokaryotes [13, 14]. A recent study reported similar phenomena in human cells under heat stress [44]. In addition, we generally ruled out the possible biases in translation initiation calculation caused by the 5’UTR regulatory elements, such as internal ribosomal entry site (IRES) and upstream open reading frames (uORF) [49]. These findings suggested that the EVI analysis in this study distinguished slow-translating genes in human cells. Second, we found that the low EVI genes are more stable in the degradation rate constant analysis. This finding is consistent with and reciprocally supports the pause-to-fold theories [13, 19] (reviewed in ref. [7, 20, 21]). Notably, our EVI evaluation provided genome-wide and individual gene resolution that moves the field forward. Third, as a negative control, we found that the most frequently initiated genes tend to have fast elongation speeds (Fig 2A and 2B). For example, the high TR genes HIST1H1E and HIST1H4C (S9A Fig) have been shown to have very stable and robust structures that are heat-resistant up to 60°C [50], implicating that the correct folding of these proteins may not be co-translational folding dependent. Favorably, O'Brien has proposed that faster translation of some proteins may assist their correct folding by speeding up termination [51].

In addition, we noted that no genes that have both high TR and low EVI properties. As an explanation, it is known that the high initiation efficiency and significant translational pausing sites will increase the risk of ribosome traffic jam [52] that frequently causes frame-shifting and premature drop-off [53], leading to aggregations of aberrant proteins. This is detrimental to cells by wasting ribosome resource, and thus should be maximally avoided per evolution, at least under physiological conditions, which has been evidenced by modeling and experiments [9, 54–56]. In addition, Chu et al have found that slow codons lead to postponed recycle of ribosomes that are the resources of next round translation initiation [57]. This serves as a reasonable explanation on the significant EVI-TR correlation demonstrated in our current study. These evolution relevant rationales biologically supported the validity of using EVI to evaluate elongation speed.

As a mechanistic explanation, we demonstrated the codon preference in terms of translation elongation speed in human cells. Accurate identifications of slow-translating codons relies on the quantitative and codon-wise tRNA abundance data; however, previous efforts found that such data are tissue-specific, time-dependent and difficult to obtain with anticodon-resolution in eukaryotic cells [58–60] (reviewed in ref. [20]). From the available tRNA data [58], we observed that the tRNA abundances between different tissues have general weak or insignificant correlation (|R|<0.4 or P>0.01) for most cases (S16 Fig). Interestingly, based on EVI analysis, we detected slow-translating codons that were tissue/cell type specific as well. With PCA, we found that one or two PC(s) would be sufficient to represent the codon preferences that were capable of distinguishing the low EVI genes from the background genomes. These results emphasized the validity of EVI examination and the codon preference in human cells. Notably, we observed no significant correlations between codon preference and codon usage or tAI, which was comparable with Ingolia et al, proposing that translation elongation speed was independent of codon usage for a remarkable proportion of genes [9].

As expected, we observed that EVI has very weak or no correlations with CAI or CBI, respectively; however, it is significantly correlated to Nc. This can be explained by the codon selection and translational elongation. Lower Nc values suggest that each amino acid is encoded by fewer codons, corresponding to higher homozygosity [42, 43]. In the extreme condition of the lowest Nc of 20, cells should not use the slow-translating codons for each amino acids in that consecutive slow-translating codons will remarkably prolong the translation duration and increase the risk of frameshift, ribosome drop-off and ribosome jamming [53, 61]. Indeed, this is maximally avoided in the real genomes [14]. Thus, the low Nc genes tend to use fast-translating codons, which coincides with their higher translation elongation speed measured in our study (higher EVI).

Indeed, the EVI measurement strategy shown in this study can be applied to more cell types and situations, which will lead to a complete resolution of identifying slow-translating codons in pan eukaryotic systems. Nevertheless, we realized that it is of great importance to independently validate this EVI evaluation on elongation speed by incorporating the tRNA-omics data into this computation per further technological breakthrough.

We demonstrated the possible relevance of co-translational folding to cancer phenotypes to justify the biological relevance of EVI as an example. In one way, as we previously reported, lung cancer cells have the loop-back enhancement of global translation initiation [6]. Hence, high TR should be a productivity-focused translational control mode. In the other way, to increase the quality of tumor-favorable proteins should be demanded by cancer cells. This was supported by the IPA and stability analyses on the slow-translating gene shown in this study. Specifically, EGF pathway over-activation, caused by such as EGFR mutation, has been identified as a prevalent mechanism for the onset and progression of lung cancers [62]. Interestingly, the low relative EVI genes in A549/HBE comparisons included the IPA clustered genes of p39, PI3K and SOS, which are all central messengers for the EGF pathway. Promoting the proper folding of these proteins in A549 cells should be favorable for their malignant phenotypes. Supportive to this notion, similar evidence was also found in the H1299/HBE comparisons.

We realized that the above rationale could not lead to definitive conclusion to link co-translational folding to cancer. But, the significant acceleration of TSG translation elongation especially in highly malignant cancer cells is an interesting finding in this study. Due to the suppression of translational pausing, the multi-domain TSGs would probably become less folded, thus disables their functions. This facilitates the malignancy of the cancer cells. In the contrast, some other genes were decelerated meanwhile for the same purpose. For example, 14-3-3 in H1299 cells is one of the most decelerated genes. 14-3-3 interacts with numerous misfolded proteins and intrinsic disordered proteins and thus plays a key role in aggresome formation [63, 64], helping to clear up the misfolded proteins generated in lung cancer cells, including the non-functional and misfolded TSGs. It also regulates functions of hundreds of proteins by conformational modulation (reviewed in [65]) and involves in cancer development, progression and chemoresistance (reviewed in [66, 67]). Therefore, the proper folding of 14-3-3 is crucial for the malignancy, and the cancer cells consolidated the folding by relatively decelerating 14-3-3 translation elongation. To be noted, we demonstrated that the deceleration of elongation did not statistically change the protein abundance of the EVI down-regulated genes, and stabilizes the protein folding in an environment in which protein misfolding occurs more frequently. These implicate that the protein quality of their gene products may be increased. Although further studies are necessary to fully reveal the detailed mechanisms, our EVI evaluation shed lights on the translational elongation regulation and its synergy with translational initiation regulation in cancer.

Methods

Cell culture

HeLa cells (ATCC, Rockville, MD) were maintained in DMEM (Invitrogen) supplemented with 10% fetal bovine serum (PAA), 1% penicillin/streptomycin and 10 μg/ml ciprofloxacin. Normal human bronchial epithelial (HBE) cells as well as human lung cancer A549 and H1299 cells were maintained as described previously [6]. In brief, the cells were maintained in Dulbecco’s modified Eagle’s medium (Invitrogen, Carlsbad, CA), supplemented with 10% fetal bovine serum (PAA Australia, Weike Biochemical Reagent, Shanghai, China), 1% penicillin/streptomycin and 10 μg/mL ciprofloxacin.

Transcriptome and full length translating mRNA sequencing

The total RNA and RNC-RNA extraction from HeLa cells was performed as previously described [6]. A pooled sample from 3 independent experiments was used for subsequent RNA-seq regarding mRNA and RNC-RNA, respectively. For mRNA and RNC-mRNA, the polyA+ mRNA was selected by RNA Purification Beads (Illumina). The library was constructed by using the Illumina TruSeq RNA sample Prep Kit v2 and sequenced by the Illumina HiSeq-2000 for 50 cycles. High quality reads that passed the Illumina quality filters were kept for the sequence analysis. These sequencing datasets are available at Gene Expression Omnibus database (accession number GSE46613).

Ribosome profiling

It has been debated that the application of cycloheximide in ribosome profiling may cause artifact like the excessive accumulation of RFPs near the start codons [9]. However, recent studies confirmed that the cycloheximide gives comparable result as the no-drug protocol when properly processed [68]. Brief treatment of cells by cycloheximide does not distort the ribosome profiling measurements [9]. In addition, elongation inhibitors like cycloheximide are especially useful to preserve the translational state during sample preparation [9, 69]. Therefore, we chose to use cycloheximide in our ribosome profiling experiment.

Cells were pre-treated with 100 mg/ml cycloheximide for 15 min, followed by 3 washes with pre-chilled phosphate buffered saline (PBS) prior to the addition of 2 ml cell lysis buffer [1% Triton X-100 in ribosome buffer (RB buffer, 20 mM HEPES-KOH (pH 7.4), 15 mM MgCl2, 200 mM KCl, 100 μg/ml cycloheximide and 2mM dithiothreitol)] to each of the T-75 flasks. After 30-min ice-bath, cell lysates were split into 2 pre-chilled 1.5 ml tubes. Cell debris was removed by centrifugation at 16,200 ×g for 10 min at 4°C. Supernatants were transferred into new pre-chilled 1.5 ml tubes with addition of 2 μl Ribolock RNase Inhibitor (40U/μl, Fermentas) in each tube. RNase I (10U/μl, Fermentas) was then added at 0.2 μl per tube, followed by incubation at 37°C for 15min and reaction termination with 1% SDS (1/10 volume per tube). The digested samples were pooled and layered on the surface of 15 ml sucrose buffer (30% sucrose in RB buffer). The ribosomes were pelleted by ultracentrifugation at 185,000 ×g for 5 h at 4°C. RNA extraction was then performed by Trizol method and ribosomal RNA (rRNA) was depleted using Ribo-Zero rRNA Removal Kit (Human/Mouse/Rat) (Epicenter) by following the manufacturer’s instructions. The ~28nt RNA fragments were considered as ribosome footprints (RFPs), which can be visualized in the RFP samples, while are non-detectable in the intact RNC-RNA samples.

The sequencing libraries of RFP were constructed, following the NEBNext Multiplex Small RNA Library Prep Set for Illumina Guide (NEB). The library was resolved by a 6% polyacrylamide gel. The fraction with the insertion size ~28nt was excised and purified from the gel. This fraction was sequenced by an Illumina HiSeq-2000 sequencer for 36 or 50 cycles. The sequencing datasets for HeLa, HBE, A549 and H1299 cells are available at Gene Expression Omnibus database (accession number GSE46613, reviewer access link: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=tfadpuykakkowxi&acc=GSE46613).

We analyzed the read density near the start codon according to the method described in [9]. All 4 cell lines in our study showed the peak of normalized average reads after the start codons are less than 2 (S17 Fig), comparable to the no-drug control in [9]. This evidenced that the cycloheximide in our study did not create artifacts in ribosome profiles.

Sequence analysis

For ribosome profiling datasets, the adapter sequences were removed from all reads. Reads were truncated at their first nucleotides, whose Phred quality scores were less than 10. Reads shorter than 18 nt were then discarded. The rest high quality reads were aligned to the RefSeq-RNA reference sequence (downloaded from http://hgdownload.cse.ucsc.edu/downloads, accessed on Jan. 21st, 2013) using FANSe 2 algorithm [70] with the parameters –L60 –E2 –U1 –S10. For mRNA and RNC-mRNA sequencing datasets, the reads were mapped to RefSeq-RNA reference sequence with FANSe2 algorithm with the parameters #x2013;L55 –E4 –U0 –S10. Alternative splice variants were merged [6]. The expression levels were estimated by using the rpkM unit [71]. The mRNA length information was also acquired from RefSeq.

Hotelling's T2 ellipse analysis

All genes were located in the TR against EVI scatter plot. Principle component analysis followed by multivariate generalization of t-test was performed to compute the Hotelling's T2 statistics [72]. The 99% confidence ellipse was calculated according to Li et al [73]. In brief, for bivariate observation values Z with a sample size of n, considering its unbiased estimate of covariance matrix S, a (1-α) confidence ellipse for prediction is given by the equation Where F2,n-2(1-α) is the (1-α) critical value of an F variate with degrees of freedom 2 and n-2.

Codon preference analysis

The synonymous codon selection preference was calculated using the method of relative synonymous codon usage (RSCU) [74]: where xij is the usage of the j-th codon for the i-th amino acid, which is encoded by ni codons [74].

The Preference Score (PS) for each codon of a certain group of genes was defined as the RSCU of this codon in this group of genes vs. that in the all genes, in log2 scale: . A positive value of PS denotes favored codon in this gene group, while a negative value indicates a disfavored codon for such genes.

The CAI of a gene was calculated according to Sharp et al. [27]: , where L is the number of codons, fi is the frequency of the codon i, and max(fj) is the maximum codon frequency for that amino acid.

The CBI was calculated according to Bennetzen et al. [41]: , where Npfr is the total number of occurrences of preferred codons, Nran is the expected number of the preferred codons if all synonymous codons were used equally, and Ntot is the total number of the 17 amino acids encoded by the preferred codons.

The effective number of codons (Nc) was calculated according to Wright et al. [42]: , where Fn is the average homozygosity for the amino acids having a degeneracy of n codons.

The tAI of each codon was calculated according to dos Reis et al. [75].

Ingenuity Pathway Analysis (IPA)

The genes with low relative EVIs and their respective absolute number of relative fold changes were uploaded and analyzed by IPA (http://www.ingenuity.com/). Core analysis of IPA was then performed as we previously reported, [6, 76]. The likelihood of associated between a set of genes with a pathway in Global Functional Analysis (GFA) and Global Canonical Pathways (GCP) was measured by the P-value, calculated by using the right-tailed Fisher’s Exact Test. P<0.001 was considered statistically significant.

Supporting Information

S1 Table. List of the TR and EVI values of all of the tested cells.

https://doi.org/10.1371/journal.pgen.1005901.s001

(XLSX)

S2 Table. List of gene quantifications of mRNA-seq, RNC-seq and Ribo-seq in HeLa, HBE, A549 and H1299 cells.

Splice variants were merged and HGNC gene names were provided.

https://doi.org/10.1371/journal.pgen.1005901.s002

(XLSX)

S3 Table. The threshold of the top 1% of TR and the lowest 1% of EVI.

https://doi.org/10.1371/journal.pgen.1005901.s003

(DOCX)

S4 Table. Overlapping low-EVI and high-TR genes among different cell lines.

https://doi.org/10.1371/journal.pgen.1005901.s004

(DOCX)

S5 Table. List of genes with the EVI change of more than 10 folds, comparing A549 and H1299 with HBE cells, respectively.

https://doi.org/10.1371/journal.pgen.1005901.s005

(XLSX)

S1 Fig. The distribution of RNC-mRNA ribosome densities (D = TR/EVI).

(A) D distribution in each analyzed cell line. (B) D correlations between cell lines. Rs = Spearman R.

https://doi.org/10.1371/journal.pgen.1005901.s006

(PDF)

S2 Fig. The gene distribution based on TR or EVI.

(A) Empirical cumulative density function (CDF). The grey dashed line denotes the top 1% of TR or the lowest 1% of EVI. The threshold values are listed in the S2 Table. (B) Empirical probability density function (PDF).

https://doi.org/10.1371/journal.pgen.1005901.s007

(PDF)

S3 Fig. The EVI calculated excluding 5’-UTR (x-axes) and the EVI calculated including 5’-UTR (y-axes) in all four tested cell lines.

Rp = Pearson R; Rs = Spearman R.

https://doi.org/10.1371/journal.pgen.1005901.s008

(PDF)

S4 Fig. The mutual correlation of gene expression levels of the 4 analyzed cell lines at the mRNA (A), RNC-mRNA (B) and RFP (C) levels.

The Spearman correlation coefficient (Rs), Pearson correlation coefficient (Rp) and their corresponding P-values are indicated the Figs (D-F).

https://doi.org/10.1371/journal.pgen.1005901.s009

(PDF)

S5 Fig. The Gene ontology enrichment analysis (Biological Processes, BP) for the low-EVI genes in the four cell lines.

The analysis was performed using PANTHER website (http://pantherdb.org/). The significance level was set to P<0.01.

https://doi.org/10.1371/journal.pgen.1005901.s010

(PDF)

S6 Fig. Correlation between EVI and codon adaptation index (CAI).

The Rs and its p-value are shown on the top of each panel.

https://doi.org/10.1371/journal.pgen.1005901.s011

(PDF)

S7 Fig. Correlation between EVI and codon bias index (CBI).

The Rs, Rp and their P-values are indicated on the top of each panel.

https://doi.org/10.1371/journal.pgen.1005901.s012

(PDF)

S8 Fig. Correlation between EVI and Effective Number of Codons (Nc).

The Rs, Rp and their P-values are shown on the top of each panel.

https://doi.org/10.1371/journal.pgen.1005901.s013

(PDF)

S9 Fig. The RFP coverage of representative genes with high TR (blue) and low EVI (red) genes in HeLa (A), HBE (B), A549 (C) and H1299 (D) cells, respectively.

The sequence coverage range of each plot is indicated in the parentheses above each mRNA illustrations. The full length mRNAs are shown in thin lines, and the corresponding CDS regions are marked with thick lines.

https://doi.org/10.1371/journal.pgen.1005901.s014

(PDF)

S10 Fig. The mRNA folding energy of the low-EVI and high-TR genes in four cell lines, calculated using MFOLD algorithm.

One-tailed KS-test were performed for each cell line under the null hypothesize that the low-EVI genes are not more stably folded than the high-TR genes.

https://doi.org/10.1371/journal.pgen.1005901.s015

(PDF)

S11 Fig. Plot matrix of the PSHigh-TR.

The Rs, Rp and their P-values are indicated.

https://doi.org/10.1371/journal.pgen.1005901.s016

(PDF)

S12 Fig. The PSLow-EVI and PSHigh-TR versus codon usage (A) and tAI (B) in HeLa, HBE, A549 and H1299 cells, respectively.

The Rp, Rs and their P-values (Pp and Ps) are listed in the tables.

https://doi.org/10.1371/journal.pgen.1005901.s017

(PDF)

S13 Fig. Correlation between PSLow-EVI (calculated in this study) and mean codon elongation time for HeLa cells.

https://doi.org/10.1371/journal.pgen.1005901.s018

(PDF)

S14 Fig. The EVI-TR correlation of the EVI up- and down-regulated genes in A549 and H1299 cells, respectively.

The Spearman R (Rs) and the P-values (Ps) were indicated.

https://doi.org/10.1371/journal.pgen.1005901.s019

(PDF)

S15 Fig. The top canonical pathway, 14-3-3 mediated signaling pathway, of the slow-translating genes, comparing lung cancer cells H1299 against the normal lung cells HBE.

https://doi.org/10.1371/journal.pgen.1005901.s020

(PDF)

S16 Fig. Mutual correlation of tRNA content of different human tissues.

(A) Plot matrix of mutual tRNA correlation of 7 human tissues normalized by the brain tRNA. (B) Rp and the –log10 p-values of the plot matrix. (C) The Rs and the #x2013;log10 p-values of the plot matrix. For (B,C), the numbers on the axes represent the tissues: 1 = liver, 2 = vulva, 3 = testis, 4 = ovary, 5 = thymus, 6 = lymph node, 7 = spleen.

https://doi.org/10.1371/journal.pgen.1005901.s021

(PDF)

S17 Fig. Metagene analysis of translation initiation of the 4 tested cell lines.

Average ribosome read density profiles of all well-expressed genes with at least 200 RFP reads are shown plotted.

https://doi.org/10.1371/journal.pgen.1005901.s022

(PDF)

Acknowledgments

We thank Qing Wang, Wangling Zhang, Lijuan Yang and Zhipeng Chen, Jinan University, for their technical assistance.

Author Contributions

Conceived and designed the experiments: GZ TW XL. Performed the experiments: XL JG WG YC JZ JJ. Analyzed the data: XL WG TW GZ. Contributed reagents/materials/analysis tools: XL QYH GZ. Wrote the paper: XL QYH TW GZ.

References

  1. 1. Sonenberg N, Hinnebusch AG. Regulation of translation initiation in eukaryotes: mechanisms and biological targets. Cell. 2009;136(4):731–45. pmid:19239892; PubMed Central PMCID: PMCPMC3610329.
  2. 2. Aitken CE, Lorsch JR. A mechanistic overview of translation initiation in eukaryotes. Nat Struct Mol Biol. 2012;19(6):568–76. pmid:22664984.
  3. 3. Larsson O, Tian B, Sonenberg N. Toward a genome-wide landscape of translational control. Cold Spring Harb Perspect Biol. 2013;5(1):a012302. pmid:23209130; PubMed Central PMCID: PMCPMC3579401.
  4. 4. Hinnebusch AG. The scanning mechanism of eukaryotic translation initiation. Annu Rev Biochem. 2014;83:779–812. pmid:24499181.
  5. 5. Arava Y, Wang Y, Storey JD, Liu CL, Brown PO, Herschlag D. Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2003;100(7):3889–94. Epub 2003/03/28. pmid:12660367; PubMed Central PMCID: PMCPMC153018.
  6. 6. Wang T, Cui Y, Jin J, Guo J, Wang G, Yin X, et al. Translating mRNAs strongly correlate to proteins in a multivariate manner and their translation ratios are phenotype specific. Nucleic Acids Res. 2013;41(9):4743–54. Epub 2013/03/23. pmid:23519614; PubMed Central PMCID: PMCPMC3643591.
  7. 7. Zhang G, Ignatova Z. Folding at the birth of the nascent chain: coordinating translation with co-translational folding. Curr Opin Struct Biol. 2011;21(1):25–31. Epub 2010/11/30. pmid:21111607.
  8. 8. Ingolia NT. Ribosome profiling: new views of translation, from single codons to genome scale. Nat Rev Genet. 2014;15(3):205–13. pmid:24468696.
  9. 9. Ingolia NT, Lareau LF, Weissman JS. Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes. Cell. 2011;147(4):789–802. Epub 2011/11/08. pmid:22056041; PubMed Central PMCID: PMCPMC3225288.
  10. 10. Duncan R, Milburn SC, Hershey JW. Regulated phosphorylation and low abundance of HeLa cell initiation factor eIF-4F suggest a role in translational control. Heat shock effects on eIF-4F. J Biol Chem. 1987;262(1):380–8. pmid:3793730.
  11. 11. Shah P, Ding Y, Niemczyk M, Kudla G, Plotkin JB. Rate-limiting steps in yeast protein translation. Cell. 2013;153(7):1589–601. pmid:23791185; PubMed Central PMCID: PMCPMC3694300.
  12. 12. Guo J, Lian X, Zhong J, Wang T, Zhang G. Length-dependent translation initiation benefits the functional proteome of human cells. Mol Biosyst. 2015;11(2):370–8. pmid:25353704.
  13. 13. Zhang G, Hubalewska M, Ignatova Z. Transient ribosomal attenuation coordinates protein synthesis and co-translational folding. Nat Struct Mol Biol. 2009;16(3):274–80. Epub 2009/02/10. pmid:19198590.
  14. 14. Zhang G, Ignatova Z. Generic algorithm to predict the speed of translational elongation: implications for protein biogenesis. PLoS One. 2009;4(4):e5036. pmid:19343177; PubMed Central PMCID: PMCPMC2661179.
  15. 15. Tuller T, Carmi A, Vestsigian K, Navon S, Dorfan Y, Zaborske J, et al. An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell. 2010;141(2):344–54. pmid:20403328.
  16. 16. Tuller T, Veksler-Lublinsky I, Gazit N, Kupiec M, Ruppin E, Ziv-Ukelson M. Composite effects of gene determinants on the translation speed and density of ribosomes. Genome Biol. 2011;12(11):R110. pmid:22050731; PubMed Central PMCID: PMCPMC3334596.
  17. 17. Spencer PS, Siller E, Anderson JF, Barral JM. Silent substitutions predictably alter translation elongation rates and protein folding efficiencies. J Mol Biol. 2012;422(3):328–35. Epub 2012/06/19. pmid:22705285; PubMed Central PMCID: PMCPMC3576719.
  18. 18. Sauna ZE, Kimchi-Sarfaty C. Understanding the contribution of synonymous mutations to human disease. Nat Rev Genet. 2011;12(10):683–91. pmid:21878961.
  19. 19. Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, Ambudkar SV, et al. A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science. 2007;315(5811):525–8. pmid:17185560.
  20. 20. Czech A, Fedyunin I, Zhang G, Ignatova Z. Silent mutations in sight: co-variations in tRNA abundance as a key to unravel consequences of silent mutations. Mol Biosyst. 2010;6(10):1767–72. pmid:20617253.
  21. 21. Plotkin JB, Kudla G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2011;12(1):32–42. pmid:21102527; PubMed Central PMCID: PMCPMC3074964.
  22. 22. Komar AA. A pause for thought along the co-translational folding pathway. Trends Biochem Sci. 2009;34(1):16–24. pmid:18996013.
  23. 23. Pedersen S. Escherichia coli ribosomes translate in vivo with variable rate. EMBO J. 1984;3(12):2895–8. pmid:6396082; PubMed Central PMCID: PMCPMC557784.
  24. 24. Sorensen MA, Kurland CG, Pedersen S. Codon usage determines translation rate in Escherichia coli. J Mol Biol. 1989;207(2):365–77. pmid:2474074.
  25. 25. Wen JD, Lancaster L, Hodges C, Zeri AC, Yoshimura SH, Noller HF, et al. Following translation by single ribosomes one codon at a time. Nature. 2008;452(7187):598–603. pmid:18327250; PubMed Central PMCID: PMCPMC2556548.
  26. 26. Gingold H, Pilpel Y. Determinants of translation efficiency and accuracy. Mol Syst Biol. 2011;7:481. pmid:21487400; PubMed Central PMCID: PMCPMC3101949.
  27. 27. Sharp PM, Cowe E. Synonymous codon usage in Saccharomyces cerevisiae. Yeast. 1991;7(7):657–78. pmid:1776357.
  28. 28. Fluitt A, Pienaar E, Viljoen H. Ribosome kinetics and aa-tRNA competition determine rate and fidelity of peptide synthesis. Comput Biol Chem. 2007;31(5–6):335–46. pmid:17897886; PubMed Central PMCID: PMCPMC2727733.
  29. 29. Chu D, Barnes DJ, von der Haar T. The role of tRNA and ribosome competition in coupling the expression of different mRNAs in Saccharomyces cerevisiae. Nucleic Acids Res. 2011;39(15):6705–14. pmid:21558172; PubMed Central PMCID: PMCPMC3159466.
  30. 30. Fedyunin I, Lehnhardt L, Bohmer N, Kaufmann P, Zhang G, Ignatova Z. tRNA concentration fine tunes protein solubility. FEBS Lett. 2012;586(19):3336–40. Epub 2012/07/24. pmid:22819830.
  31. 31. Tu C, Tzeng TH, Bruenn JA. Ribosomal movement impeded at a pseudoknot required for frameshifting. Proc Natl Acad Sci U S A. 1992;89(18):8636–40. pmid:1528874; PubMed Central PMCID: PMCPMC49975.
  32. 32. Ledoux S, Uhlenbeck OC. Different aa-tRNAs are selected uniformly on the ribosome. Mol Cell. 2008;31(1):114–23. pmid:18614050; PubMed Central PMCID: PMCPMC2709977.
  33. 33. Lu J, Deutsch C. Electrostatics in the ribosomal tunnel modulate chain elongation rates. J Mol Biol. 2008;384(1):73–86. pmid:18822297; PubMed Central PMCID: PMCPMC2655213.
  34. 34. Charneski CA, Hurst LD. Positively charged residues are the major determinants of ribosomal velocity. PLoS Biol. 2013;11(3):e1001508. pmid:23554576; PubMed Central PMCID: PMCPMC3595205.
  35. 35. Curran JF, Yarus M. Rates of aminoacyl-tRNA selection at 29 sense codons in vivo. J Mol Biol. 1989;209(1):65–77. 2478714. pmid:2478714
  36. 36. Wilson DN, Beckmann R. The ribosomal tunnel as a functional environment for nascent polypeptide folding and translational stalling. Curr Opin Struct Biol. 2011;21(2):274–82. pmid:21316217.
  37. 37. Li GW, Oh E, Weissman JS. The anti-Shine-Dalgarno sequence drives translational pausing and codon choice in bacteria. Nature. 2012;484(7395):538–41. pmid:22456704; PubMed Central PMCID: PMCPMC3338875.
  38. 38. Ingolia NT, Ghaemmaghami S, Newman JR, Weissman JS. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science. 2009;324(5924):218–23. Epub 2009/02/14. pmid:19213877; PubMed Central PMCID: PMCPMC2746483.
  39. 39. Zhong J, Cui Y, Guo J, Chen Z, Yang L, He QY, et al. Resolving chromosome-centric human proteome with translating mRNA analysis: a strategic demonstration. J Proteome Res. 2014;13(1):50–9. pmid:24200226.
  40. 40. Cambridge SB, Gnad F, Nguyen C, Bermejo JL, Kruger M, Mann M. Systems-wide proteomic analysis in mammalian cells reveals conserved, functional protein turnover. J Proteome Res. 2011;10(12):5275–84. Epub 2011/11/05. pmid:22050367.
  41. 41. Bennetzen JL, Hall BD. Codon selection in yeast. J Biol Chem. 1982;257(6):3026–31. pmid:7037777.
  42. 42. Wright F. The 'effective number of codons' used in a gene. Gene. 1990;87(1):23–9. pmid:2110097.
  43. 43. Fuglsang A. The 'effective number of codons' revisited. Biochem Biophys Res Commun. 2004;317(3):957–64. pmid:15081433.
  44. 44. Shalgi R, Hurt JA, Krykbaeva I, Taipale M, Lindquist S, Burge CB. Widespread regulation of translation by elongation pausing in heat shock. Mol Cell. 2013;49(3):439–52. pmid:23290915; PubMed Central PMCID: PMCPMC3570722.
  45. 45. Wu G, Nie L, Freeland SJ. The effects of differential gene expression on coding sequence features: analysis by one-way ANOVA. Biochem Biophys Res Commun. 2007;358(4):1108–13. Epub 2007/05/23. pmid:17517370.
  46. 46. Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15. pmid:12824337; PubMed Central PMCID: PMCPMC169194.
  47. 47. Li LP, Lu CH, Chen ZP, Ge F, Wang T, Wang W, et al. Subcellular proteomics revealed the epithelial-mesenchymal transition phenotype in lung cancer. Proteomics. 2011;11(3):429–39. Epub 2011/01/27. pmid:21268272.
  48. 48. Chen Y, Li Y, Zhong J, Zhang J, Chen Z, Yang L, et al. Identification of Missing Proteins Defined by Chromosome-Centric Proteome Project in the Cytoplasmic Detergent-Insoluble Proteins. J Proteome Res. 2015;14(9):3693–709. pmid:26108252.
  49. 49. Fritsch C, Herrmann A, Nothnagel M, Szafranski K, Huse K, Schumann F, et al. Genome-wide search for novel human uORFs and N-terminal protein extensions using ribosomal footprinting. Genome Res. 2012;22(11):2208–18. pmid:22879431; PubMed Central PMCID: PMCPMC3483550.
  50. 50. Caterino TL, Hayes JJ. Structure of the H1 C-terminal domain and function in chromatin condensation. Biochem Cell Biol. 2011;89(1):35–44. pmid:21326361; PubMed Central PMCID: PMCPMC3787537.
  51. 51. O'Brien EP, Vendruscolo M, Dobson CM. Kinetic modelling indicates that fast-translating codons can coordinate cotranslational protein folding by avoiding misfolded intermediates. Nat Commun. 2014;5:2988. pmid:24394622.
  52. 52. Zhang S, Goldman E, Zubay G. Clustering of low usage codons and ribosome movement. J Theor Biol. 1994;170(4):339–54. pmid:7996861.
  53. 53. Zhang G, Fedyunin I, Miekley O, Valleriani A, Moura A, Ignatova Z. Global and local depletion of ternary complex limits translational elongation. Nucleic Acids Res. 2010;38(14):4778–87. Epub 2010/04/03. pmid:20360046; PubMed Central PMCID: PMCPMC2919707.
  54. 54. Racle J, Overney J, Hatzimanikatis V. A computational framework for the design of optimal protein synthesis. Biotechnol Bioeng. 2012;109(8):2127–33. pmid:22334333.
  55. 55. Ciandrini L, Stansfield I, Romano MC. Ribosome traffic on mRNAs maps to gene ontology: genome-wide quantification of translation initiation rates and polysome size regulation. PLoS Comput Biol. 2013;9(1):e1002866. pmid:23382661; PubMed Central PMCID: PMCPMC3561044.
  56. 56. Mitarai N, Pedersen S. Control of ribosome traffic by position-dependent choice of synonymous codons. Phys Biol. 2013;10(5):056011. pmid:24104350.
  57. 57. Chu D, Kazana E, Bellanger N, Singh T, Tuite MF, von der Haar T. Translation elongation can control translation initiation on eukaryotic mRNAs. EMBO J. 2014;33(1):21–34. pmid:24357599; PubMed Central PMCID: PMCPMC3990680.
  58. 58. Dittmar KA, Goodenbour JM, Pan T. Tissue-specific differences in human transfer RNA expression. PLoS Genet. 2006;2(12):e221. pmid:17194224; PubMed Central PMCID: PMCPMC1713254.
  59. 59. Frenkel-Morgenstern M, Danon T, Christian T, Igarashi T, Cohen L, Hou YM, et al. Genes adopt non-optimal codon usage to generate cell cycle-dependent oscillations in protein levels. Mol Syst Biol. 2012;8:572. pmid:22373820; PubMed Central PMCID: PMCPMC3293633.
  60. 60. Pavon-Eternod M, Gomes S, Geslain R, Dai Q, Rosner MR, Pan T. tRNA over-expression in breast cancer and functional consequences. Nucleic Acids Res. 2009;37(21):7268–80. pmid:19783824; PubMed Central PMCID: PMCPMC2790902.
  61. 61. Buchan JR, Stansfield I. Halting a cellular production line: responses to ribosomal pausing during translation. Biol Cell. 2007;99(9):475–87. pmid:17696878.
  62. 62. Paez JG, Janne PA, Lee JC, Tracy S, Greulich H, Gabriel S, et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304(5676):1497–500. pmid:15118125.
  63. 63. Jia B, Wu Y, Zhou Y. 14-3-3 and aggresome formation: implications in neurodegenerative diseases. Prion. 2014;8(2). pmid:24549097; PubMed Central PMCID: PMCPMC4189886.
  64. 64. Bustos DM, Iglesias AA. Intrinsic disorder is a key characteristic in partners that bind 14-3-3 proteins. Proteins. 2006;63(1):35–42. pmid:16444738.
  65. 65. Obsilova V, Kopecka M, Kosek D, Kacirova M, Kylarova S, Rezabkova L, et al. Mechanisms of the 14-3-3 protein function: regulation of protein function through conformational modulation. Physiol Res. 2014;63 Suppl 1:S155–64. pmid:24564655.
  66. 66. Matta A, Siu KW, Ralhan R. 14-3-3 zeta as novel molecular target for cancer therapy. Expert Opin Ther Targets. 2012;16(5):515–23. pmid:22512284.
  67. 67. Reinhardt HC, Yaffe MB. Phospho-Ser/Thr-binding domains: navigating the cell cycle and DNA damage response. Nat Rev Mol Cell Biol. 2013;14(9):563–80. pmid:23969844.
  68. 68. Guttman M, Russell P, Ingolia NT, Weissman JS, Lander ES. Ribosome profiling provides evidence that large noncoding RNAs do not encode proteins. Cell. 2013;154(1):240–51. pmid:23810193; PubMed Central PMCID: PMCPMC3756563.
  69. 69. Oh E, Becker AH, Sandikci A, Huber D, Chaba R, Gloge F, et al. Selective ribosome profiling reveals the cotranslational chaperone action of trigger factor in vivo. Cell. 2011;147(6):1295–308. pmid:22153074; PubMed Central PMCID: PMCPMC3277850.
  70. 70. Xiao CL, Mai ZB, Lian XL, Zhong JY, Jin JJ, He QY, et al. FANSe2: a robust and cost-efficient alignment tool for quantitative next-generation sequencing applications. PLoS One. 2014;9(4):e94250. pmid:24743329; PubMed Central PMCID: PMCPMC3990525.
  71. 71. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. Epub 2008/06/03. pmid:18516045.
  72. 72. Hotelling H. The Generalization of Student's Ratio. Ann Math Statist. 1931;2(3):360–78.
  73. 73. Li H, Wang L, Yan X, Liu Q, Yu C, Wei H, et al. A proton nuclear magnetic resonance metabonomics approach for biomarker discovery in nonalcoholic fatty liver disease. J Proteome Res. 2011;10(6):2797–806. pmid:21563774.
  74. 74. Sharp PM, Li WH. Codon usage in regulatory genes in Escherichia coli does not reflect selection for 'rare' codons. Nucleic Acids Res. 1986;14(19):7737–49. pmid:3534792; PubMed Central PMCID: PMCPMC311793.
  75. 75. dos Reis M, Savva R, Wernisch L. Solving the riddle of codon usage preferences: a test for translational selection. Nucleic Acids Res. 2004;32(17):5036–44. pmid:15448185; PubMed Central PMCID: PMCPMC521650.
  76. 76. Shen S, Guo J, Luo Y, Zhang W, Cui Y, Wang Q, et al. Functional proteomics revealed IL-1beta amplifies TNF downstream protein signals in human synoviocytes in a TNF-independent manner. Biochem Biophys Res Commun. 2014;450(1):538–44. pmid:24928389.