Differential expression of human tRNA genes drives the abundance of tRNA-derived fragments.

The human genome encodes hundreds of transfer RNA (tRNA) genes but their individual contribution to the tRNA pool is not fully understood. Deep sequencing of tRNA transcripts (tRNA-Seq) can estimate tRNA abundance at single gene resolution, but tRNA structures and posttranscriptional modifications impair these analyses. Here we present a bioinformatics strategy to investigate differential tRNA gene expression and use it to compare tRNA-Seq datasets from cultured human cells and human brain. We find that sequencing caveats affect quantitation of only a subset of human tRNA genes. Unexpectedly, we detect several cases where the differences in tRNA expression among samples do not involve variations at the level of isoacceptor tRNA sets (tRNAs charged with the same amino acid but using different anticodons), but rather among tRNA genes within the same isodecoder set (tRNAs having the same anticodon sequence). Because isodecoder tRNAs are functionally equal in terms of genetic translation, their differential expression may be related to noncanonical tRNA functions. We show that several instances of differential tRNA gene expression result in changes in the abundance of tRNA-derived fragments (tRFs) but not of mature tRNAs. Examples of differentially expressed tRFs include PIWI-associated RNAs, tRFs present in tissue samples but not in cells cultured in vitro, and somatic tissue-specific tRFs. Our data support that differential expression of tRNA genes regulate noncanonical tRNA functions performed by tRFs.

Association between mitochondrial DNA methylation and internal exposure to polycyclic aromatic hydrocarbons (PAHs), nitrated-PAHs (NPAHs) and oxygenated-PAHs (OPAHs) in young adults from Tianjin, China.

Polycyclic aromatic hydrocarbons (PAHs), nitrated-PAHs (NPAHs) and oxygenated-PAHs (OPAHs) are environmental pollutants with adverse effects on human health. The correlation between the concentrations of PAHs, NPAHs and OPAHs in human plasma and the methylation level of mitochondrial DNA (mtDNA) was investigated using data from 110 plasma samples collected in Tianjin, China. The median concentrations of PAHs, NPAHs and OPAHs were 16.0 (IQR: 14.4-20.7) ng/mL, 82.2 (IQR: 63.1-97.6) ng/mL and 49.6 (IQR: 28.6-53.8) ng/mL, and the mean proportions were 13.4%, 56.5% and 30.1%, respectively. Bisulfite-PCR pyrosequencing was used to measure the methylation level of MT-CO1 and tRNA-Leu. The methylation levels of two mitochondrial genes (MT-CO1, tRNA-Leu) including four CpG sites (MT-CO1-P1, MT-CO1-P2, tRNA-Leu-P1 and tRNA-Leu-P2) were 0.67% ± 1.38%, 13.54% ± 2.59%, 7.23% ± 5.35% and 1.64% ± 2.94%, respectively. To the best of our knowledge, this is the first time that significant correlations were found between PAHs and their derivatives exposure and mtDNA methylation levels.

Full-Range Profiling of tRNA Modifications Using LC-MS/MS at Single-Base Resolution through a Site-Specific Cleavage Strategy.

Transfer RNAs (tRNAs) are the most heavily modified RNA species. Liquid chromatography coupled with mass spectrometry (LC-MS/MS) is a powerful tool for characterizing tRNA modifications, which involves pretreating tRNAs with base-specific ribonucleases to produce smaller oligonucleotides amenable to MS. However, the quality and quantity of products from base-specific digestions are severely impacted by the base composition of tRNAs. This often leads to a loss of sequence information. Here, we report a method for the full-range profiling of tRNA modifications at single-base resolution by combining site-specific RNase H digestion with the LC-MS/MS and RNA-seq techniques. The key steps were designed to generate high-quality products of optimal lengths and ionization properties. A linear correlation between collision energies and the m/z of oligonucleotides significantly improved the information content of collision-induced dissociation (CID) spectra. False positives were eliminated by up to 95% using novel inclusion criteria for collecting a census of modifications. This method is illustrated by the mapping of mouse mitochondrial tRNAHis(GUG) and tRNAVal(UAC), which were hitherto not investigated. The identities and locations of the five species of modifications on these tRNAs were fully characterized. This approach is universally applicable to any tRNA species and provides an experimentally realizable pathway to the de novo sequencing of post-transcriptionally modified tRNAs with high sequence coverage.

Cost-efficiency tradeoff is optimized in various cancer types revealed by genome-wide analysis.

The tradeoff between cost and efficiency is omnipresent in organisms. Specifically, how the evolutionary force shapes the tradeoff between biosynthetic cost and translation efficiency remains unclear. In the cancer community, whether the adjustment of cost-efficiency tradeoff acts as a strategy to facilitate tumor proliferation and contributes to oncogenesis is uninvestigated. To address this issue, we retrieved the gene expression profile in various cancer types and the matched normal samples from The Cancer Genome Atlas (TCGA). We found that the highly expressed genes in cancers generally have higher tAI/nitro ratios than those in normal samples. This is possibly caused by the higher tAI/nitro ratios observed in oncogenes than tumor suppressor genes (TSG). Furthermore, in the cancer samples, derived mutations in oncogenes usually lead to higher tAI/nitro ratios, while those mutations in TSG lead to lower tAI/nitro. For a special case of kidney cancer, we investigated several crucial genes in tumor samples versus normal samples, and discovered that the changes in tAI/nitro ratios are correlated with the changes in translation level. Our study for the first time revealed the optimization of cost-efficiency tradeoff in cancers. The cost-efficiency dilemma is optimized by the tumor cells, and is possibly beneficial for the translation and production of oncogenes, and eventually contributes to proliferation and oncogenesis. Our findings could provide novel perspectives in depicting the cancer genomes and might help unravel the cancer evolution.

Translational efficiency across healthy and tumor tissues is proliferation-related.

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

Improved RNA modification mapping of cellular non-coding RNAs using C- and U-specific RNases.

Locating ribonucleoside modifications within an RNA sequence requires digestion of the RNA into oligoribonucleotides of amenable size for subsequent analysis by LC-MS (liquid chromatography-mass spectrometry). This approach, widely referred to as RNA modification mapping, is facilitated through ribonucleases (RNases) such as T1 (guanosine-specific), U2 (purine-selective) and A (pyrimidine-specific) among others. Sequence coverage by these enzymes depends on positioning of the recognized nucleobase (such as guanine or purine or pyrimidine) in the sequence and its ribonucleotide composition. Using E. coli transfer RNA (tRNA) and ribosomal RNA (rRNA) as model samples, we demonstrate the ability of complementary nucleobase-specific ribonucleases cusativin (C-specific) and MC1 (U-specific) to generate digestion products that facilitate confident mapping of modifications in regions such as G-rich and pyrimidine-rich segments of RNA, and to distinguish C to U sequence differences. These enzymes also increase the number of oligonucleotide digestion products that are unique to a specific RNA sequence. Further, with these additional RNases, multiple modifications can be localized with high confidence in a single set of experiments with minimal dependence on the individual tRNA abundance in a mixture. The sequence overlaps observed with these complementary digestion products and that of RNase T1 improved sequence coverage to 75% or above. A similar level of sequence coverage was also observed for the 2904 nt long 23S rRNA indicating their utility has no dependence on RNA size. Wide-scale adoption of these additional modification mapping tools could help expedite the characterization of modified RNA sequences to understand their structural and functional role in various living systems.

The Origin of tRNA Deduced from Pseudomonas aeruginosa 5' Anticodon-Stem Sequence : Anticodon-stem loop hypothesis.

The riddle of the origin of life is unsolved as yet. One of the best ways to solve the riddle would be to find a vestige of the first life from databases of DNA and/or protein of modern organisms. It would be, especially, important to know the origin of tRNA, because it mediates between genetic information and the amino acid sequence of a protein. Here I attempt to find a vestige of the origin and evolution of tRNA from base sequences of Pseudomonas aeruginosa tRNA gene. It was first perceived that 5' anticodon (AntiC) stem sequences of P. aeruginosa tRNA for translation of G-start codon (GNN) are intimately and mutually related. Then, mutual relations among all of the forty-two 5' AntiC stem sequences of P. aeruginosa tRNA were examined. These relationships imply that P. aeruginosa tRNA originated from four anticodon stem-loops (AntiC-SL) translating GNC codons to the corresponding four amino acids, Gly, Ala, Asp and Val (where N is G, C, A, or T). In contrast to the case of AntiC-stem sequence, a mutual relation map could not be drawn with D-, T- and acceptor-stem sequences of P. aeruginosa tRNA. Thus I conclude that the four AntiC-SLs were the first primeval tRNAs.

Selective Enzymatic Demethylation of N2 ,N2 -Dimethylguanosine in RNA and Its Application in High-Throughput tRNA Sequencing.

The abundant Watson-Crick face methylations in biological RNAs such as N1 -methyladenosine (m1 A), N1 -methylguanosine (m1 G), N3 -methylcytosine (m3 C), and N2 ,N2 -dimethylguanosine (m22 G) cause significant obstacles for high-throughput RNA sequencing by impairing cDNA synthesis. One strategy to overcome this obstacle is to remove the methyl group on these modified bases prior to cDNA synthesis using enzymes. The wild-type E. coli AlkB and its D135S mutant can remove most of m1 A, m1 G, m3 C modifications in transfer RNA (tRNA), but they work poorly on m22 G. Here we report the design and evaluation of a series of AlkB mutants against m22 G-containing model RNA substrates that we synthesize using an improved synthetic method. We show that the AlkB D135S/L118V mutant efficiently and selectively converts m22 G modification to N2 -methylguanosine (m2 G). We also show that this new enzyme improves the efficiency of tRNA sequencing.

Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions.

Semen contains relatively ill-defined regulatory components that likely aid fertilization, but which could also interfere with defense against infection. Each ejaculate contains trillions of exosomes, membrane-enclosed subcellular microvesicles, which have immunosuppressive effects on cells important in the genital mucosa. Exosomes in general are believed to mediate inter-cellular communication, possibly by transferring small RNA molecules. We found that seminal exosome (SE) preparations contain a substantial amount of RNA from 20 to 100 nucleotides (nts) in length. We sequenced 20-40 and 40-100 nt fractions of SE RNA separately from six semen donors. We found various classes of small non-coding RNA, including microRNA (21.7% of the RNA in the 20-40 nt fraction) as well as abundant Y RNAs and tRNAs present in both fractions. Specific RNAs were consistently present in all donors. For example, 10 (of ∼2600 known) microRNAs constituted over 40% of mature microRNA in SE. Additionally, tRNA fragments were strongly enriched for 5'-ends of 18-19 or 30-34 nts in length; such tRNA fragments repress translation. Thus, SE could potentially deliver regulatory signals to the recipient mucosa via transfer of small RNA molecules.

Scp160p is required for translational efficiency of codon-optimized mRNAs in yeast.

The budding yeast multi-K homology domain RNA-binding protein Scp160p binds to >1000 messenger RNAs (mRNAs) and polyribosomes, and its mammalian homolog vigilin binds transfer RNAs (tRNAs) and translation elongation factor EF1alpha. Despite its implication in translation, studies on Scp160p's molecular function are lacking to date. We applied translational profiling approaches and demonstrate that the association of a specific subset of mRNAs with ribosomes or heavy polysomes depends on Scp160p. Interaction of Scp160p with these mRNAs requires the conserved K homology domains 13 and 14. Transfer RNA pairing index analysis of Scp160p target mRNAs indicates a high degree of consecutive use of iso-decoding codons. As shown for one target mRNA encoding the glycoprotein Pry3p, Scp160p depletion results in translational downregulation but increased association with polysomes, suggesting that it is required for efficient translation elongation. Depletion of Scp160p also decreased the relative abundance of ribosome-associated tRNAs whose codons show low potential for autocorrelation on mRNAs. Conversely, tRNAs with highly autocorrelated codons in mRNAs are less impaired. Our data indicate that Scp160p might increase the efficiency of tRNA recharge, or prevent diffusion of discharged tRNAs, both of which were also proposed to be the likely basis for the translational fitness effect of tRNA pairing.

Enhanced detection of post-transcriptional modifications using a mass-exclusion list strategy for RNA modification mapping by LC-MS/MS.

There has been a renewed appreciation for the dynamic nature of ribonucleic acid (RNA) modifications and for the impact of modified RNAs on organism health resulting in an increased emphasis on developing analytical methods capable of detecting modifications within specific RNA sequence contexts. Here we demonstrate that a DNA-based exclusion list enhances data dependent liquid chromatography tandem mass spectrometry (LC-MS/MS) detection of post-transcriptionally modified nucleosides within specific RNA sequences. This approach is possible because all post-transcriptional modifications of RNA, except pseudouridine, result in a mass increase in the canonical nucleoside undergoing chemical modification. Thus, DNA-based sequences reflect the state of the RNA prior to or in the absence of modification. The utility of this exclusion list strategy is demonstrated through the RNA modification mapping of total tRNAs from the bacteria Escherichia coli, Lactococcus lactis, and Streptomyces griseus. Creation of a DNA-based exclusion list is shown to consistently enhance the number of detected modified ribonuclease (RNase) digestion products by ∼20%. All modified RNase digestion products that were detected during standard data dependent acquisition (DDA) LC-MS/MS were also detected when the DNA-based exclusion list was used. Consequently, the increase in detected modified RNase digestion products is attributed to new experimental information only obtained when using the exclusion list. This exclusion list strategy should be broadly applicable to any class of RNA and improves the utility of mass spectrometry approaches for discovery-based analyses of RNA modifications, such as are required for studies of the epitranscriptome.

Mass spectrometry analysis of nucleosides and nucleotides.

Mass spectrometry has been widely utilised in the study of nucleobases, nucleosides and nucleotides as components of nucleic acids and as bioactive metabolites in their own right. In this review, the application of mass spectrometry to such analysis is overviewed in relation to various aspects regarding the analytical mass spectrometric and chromatographic techniques applied and also the various applications of such analysis.

Vaccinia and influenza A viruses select rather than adjust tRNAs to optimize translation.

Transfer RNAs (tRNAs) are central to protein synthesis and impact translational speed and fidelity by their abundance. Here we examine the extent to which viruses manipulate tRNA populations to favor translation of their own genes. We study two very different viruses: influenza A virus (IAV), a medium-sized (13 kB genome) RNA virus; and vaccinia virus (VV), a large (200 kB genome) DNA virus. We show that the total cellular tRNA population remains unchanged following viral infection, whereas the polysome-associated tRNA population changes dramatically in a virus-specific manner. The changes in polysome-associated tRNA levels reflect the codon usage of viral genes, suggesting the existence of local tRNA pools optimized for viral translation.

Selection on codon bias in yeast: a transcriptional hypothesis.

Codons that code for the same amino acid are often used with unequal frequencies. This phenomenon is termed codon bias. Here, we report a computational analysis of codon bias in yeast using experimental and theoretical genome-wide data. We show that the most used codons in highly expressed genes can be predicted by mRNA structural data and that the codon choice at each synonymous site within an mRNA is not random with respect to the local secondary structure. Because we also found that the folding stability of intron sequences is strongly correlated with codon bias and mRNA level, our results suggest that codon bias is linked to mRNA folding structure through a mechanism that, at least partially, operates before pre-mRNA splicing. Consistent with this, we report evidence supporting the adaptation of the tRNA pool to the codon profile of the most expressed genes rather than vice versa. We show that the correlation of codon usage with the gene expression level also includes the stop codons that are normally not decoded by aminoacyl-tRNAs. The results reported here are consistent with a role for transcriptional forces in driving codon usage bias via a mechanism that improves gene expression by optimizing mRNA folding structures.

A covalent approach for site-specific RNA labeling in Mammalian cells.

Advances in RNA research and RNA nanotechnology depend on the ability to manipulate and probe RNA with high precision through chemical approaches, both in vitro and in mammalian cells. However, covalent RNA labeling methods with scope and versatility comparable to those of current protein labeling strategies are underdeveloped. A method is reported for the site- and sequence-specific covalent labeling of RNAs in mammalian cells by using tRNA(Ile2) -agmatidine synthetase (Tias) and click chemistry. The crystal structure of Tias in complex with an azide-bearing agmatine analogue was solved to unravel the structural basis for Tias/substrate recognition. The unique RNA sequence specificity and plastic Tias/substrate recognition enable the site-specific transfer of azide/alkyne groups to an RNA molecule of interest in vitro and in mammalian cells. Subsequent click chemistry reactions facilitate the versatile labeling, functionalization, and visualization of target RNA.

Deep sequencing of RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions.

Cells release RNA-carrying vesicles and membrane-free RNA/protein complexes into the extracellular milieu. Horizontal vesicle-mediated transfer of such shuttle RNA between cells allows dissemination of genetically encoded messages, which may modify the function of target cells. Other studies used array analysis to establish the presence of microRNAs and mRNA in cell-derived vesicles from many sources. Here, we used an unbiased approach by deep sequencing of small RNA released by immune cells. We found a large variety of small non-coding RNA species representing pervasive transcripts or RNA cleavage products overlapping with protein coding regions, repeat sequences or structural RNAs. Many of these RNAs were enriched relative to cellular RNA, indicating that cells destine specific RNAs for extracellular release. Among the most abundant small RNAs in shuttle RNA were sequences derived from vault RNA, Y-RNA and specific tRNAs. Many of the highly abundant small non-coding transcripts in shuttle RNA are evolutionary well-conserved and have previously been associated to gene regulatory functions. These findings allude to a wider range of biological effects that could be mediated by shuttle RNA than previously expected. Moreover, the data present leads for unraveling how cells modify the function of other cells via transfer of specific non-coding RNA species.

Direct monitoring of initiation factor dynamics through formation of 30S and 70S translation-initiation complexes on a quartz crystal microbalance.

Translation initiation is a dynamic and complicated process requiring the building a 70S initiation complex (70S-IC) composed of a ribosome, mRNA, and an initiator tRNA. During the formation of the 70S-IC, initiation factors (IFs: IF1, IF2, and IF3) interact with a ribosome to form a 30S initiation complex (30S-IC) and a 70S-IC. Although some spectroscopic analyses have been performed, the mechanism of binding and dissociation of IFs remains unclear. Here, we employed a 27 MHz quartz crystal microbalance (QCM) to evaluate the process of bacterial IC formation in translation initiation by following frequency changes (mass changes). IFs (IF1, IF2, and IF3), N-terminally fused to biotin carboxyl carrier protein (bio-BCCP), were immobilized on a Neutravidin-covered QCM plate. By using bio-BCCP-IF2 immobilized to the QCM, three steps of the formation of ribosomal initiation complex could be sequentially observed as simple mass changes in real time: binding of a 30S complex to the immobilized IF2, a recruitment of 50S to the 30S-IC, and formation of the 70S-IC. The kinetic parameters implied that the release of IF2 from the 70S-IC could be the rate-limiting step in translation initiation. The IF3-immobilized QCM revealed that the affinity of IF3 for the 30S complex decreased upon the addition of mRNA and fMet-tRNA(fMet) but did not lead to complete dissociation from the 30S-IC. These results suggest that IF3 binds and stays bound to ICs, and its interaction mode is altered during the formation of 30S-IC and 70S-IC and is finally induced to dissociate from ICs by 50S binding. This methodology demonstrated here is applicable to investigate the role of IFs in translation initiation driven by other pathways.

The global identification of tRNA isoacceptors by targeted tandem mass spectrometry.

Transfer ribonucleic acids (tRNA) are a biologically significant class of non-coding ribonucleic acids (ncRNAs) that pose unique analytical challenges for complete characterization. Here we present a robust and simple method for the consistent and accurate identification of individual tRNAs from a pool of total tRNA obtained from cell lysate. Through this method individual isoacceptor tRNAs are identified by the detection of unique oligonucleotide sequences which arise from a single enzymatic digestion. These unique sequences can be detected by monitoring specific transitions from precursor to product ions. Thus, for any pool of known tRNA sequences including posttranscriptional modifications, targeted tandem mass spectrometry can be used for monitoring these specific transitions. The proposed method was developed and validated using a set of known tRNAs from Escherichia coli. This approach was found to identify 41 ± 2 of the predicted 47 isoaccepting tRNAs in E. coli from targeted tandem mass spectrometry using only 24 precursor m/z values. This method should be easily adapted to other bacterial systems for both genomic and posttranscriptional analysis of tRNAs, and is likely suitable for future clinical applications.

Quantitative single cell monitoring of protein synthesis at subcellular resolution using fluorescently labeled tRNA.

We have developed a novel technique of using fluorescent tRNA for translation monitoring (FtTM). FtTM enables the identification and monitoring of active protein synthesis sites within live cells at submicron resolution through quantitative microscopy of transfected bulk uncharged tRNA, fluorescently labeled in the D-loop (fl-tRNA). The localization of fl-tRNA to active translation sites was confirmed through its co-localization with cellular factors and its dynamic alterations upon inhibition of protein synthesis. Moreover, fluorescence resonance energy transfer (FRET) signals, generated when fl-tRNAs, separately labeled as a FRET pair occupy adjacent sites on the ribosome, quantitatively reflect levels of protein synthesis in defined cellular regions. In addition, FRET signals enable detection of intra-populational variability in protein synthesis activity. We demonstrate that FtTM allows quantitative comparison of protein synthesis between different cell types, monitoring effects of antibiotics and stress agents, and characterization of changes in spatial compartmentalization of protein synthesis upon viral infection.