11Background: Different tissues express genes with particular codon usage and anticodon 12 tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely 13 understood, as well as its effect on tissue-specific protein levels. 14
Results:We first validated the accuracy of small RNA-seq for tRNA quantification across 15 five human cell lines. We then analyzed tRNA expression in more than 8000 tumor samples 16 from TCGA, together with their paired mRNA-seq and proteomics data, to determine the 17 Relative Translation Efficiency. We thereby elucidate that the dynamic adaptation of the 18 tRNA pool is largely related to the proliferative state across tissues, which determines tissue-19 specific translation efficiency. Furthermore, the aberrant translational efficiency of ProCCA 20 and GlyGGT in cancer, among other codons, which is partly regulated by the tRNA gene 21 copy numbers and their promoter DNA methylation, is associated with poor patient survival. 22Conclusions: The distribution of tissue-specific tRNA pools over the whole cellular 23 translatome affects the subsequent translational efficiency, which functionally determines a 24 condition-specific expression program in tissues both in healthy and tumor states. 25 (14)(15)(16)(17)(18). In this context, The Cancer Genome Atlas (TCGA) has been recently used to 52 investigate the alteration of tRNA gene expression and translational machinery in cancer, 53 which may play a role in driving aberrant translation (19,20). 54To validate the use of small RNA-seq for tRNA quantification, we first compare tRNA levels 55 determined in HEK293 by well-established tRNA sequencing methods (Hydro-tRNAseq and 56 demethylase-tRNA-seq) (12,13,21), with those obtained by small RNA-seq. Then we 57 quantify the tRNA repertoire of five cell lines using Hydro-tRNAseq and perform small RNA-58 seq in parallel. Comparison of the tRNA abundance obtained by both approaches shows that 59 it is possible to accurately estimate relative tRNA abundance of cells and tissues using small 60 RNA-seq. Furthermore, we show that both types of quantification are informative enough to 61 distinguish between the five analyzed human cell lines covering multiple tissue types. In 62 consequence, we apply a tRNA-specific computational pipeline to re-analyze 8,534 small 63 RNA-seq datasets from TCGA (22). We find that the tissue-specificity of tRNA expression is 64 largely proliferation-related, even within healthy tissues. The tRNA quantification of TCGA 65 samples enables their comparison with paired and publicly available mRNA-seq, proteomic, 66 DNA methylation and copy number data, which underscores the role of tRNAs in globally 67 controlling a condition-specific translational program. We discover multiple codons, including 68ProCCA and GlyGGT, whose translational efficiency is compromised and leads to poor 69 prognosis in cancer. Finally, promoter DNA methylation and tRNA gene copy number arise 70 as two regulatory mechanisms controlling tRNA gene expression in cancer. 71
RESULT...