Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis.

Advances in high-throughput sequencing have facilitated remarkable insights into the diversity and functioning of naturally occurring microbes; however, current sequencing strategies are insufficient to reveal physiological states of microbial communities associated with protein translation dynamics. Transfer RNAs (tRNAs) are core components of protein synthesis machinery, present in all living cells, and are phylogenetically tractable, which make them ideal targets to gain physiological insights into environmental microbes. Here we report a direct sequencing approach, tRNA-seq, and a software suite, tRNA-seq-tools, to recover sequences, abundance profiles, and post-transcriptional modifications of microbial tRNA transcripts. Our analysis of cecal samples using tRNA-seq distinguishes high-fat- and low-fat-fed mice in a comparable fashion to 16S ribosomal RNA gene amplicons, and reveals taxon- and diet-dependent variations in tRNA modifications. Our results provide taxon-specific in situ insights into the dynamics of tRNA gene expression and post-transcriptional modifications within complex environmental microbiomes.

RNA modification landscape of the human mitochondrial tRNALys regulates protein synthesis.

Post-transcriptional RNA modifications play a critical role in the pathogenesis of human mitochondrial disorders, but the mechanisms by which specific modifications affect mitochondrial protein synthesis remain poorly understood. Here we used a quantitative RNA sequencing approach to investigate, at nucleotide resolution, the stoichiometry and methyl modifications of the entire mitochondrial tRNA pool, and establish the relevance to human disease. We discovered that a N1-methyladenosine (m1A) modification is missing at position 58 in the mitochondrial tRNALys of patients with the mitochondrial DNA mutation m.8344 A > G associated with MERRF (myoclonus epilepsy, ragged-red fibers). By restoring the modification on the mitochondrial tRNALys, we demonstrated the importance of the m1A58 to translation elongation and the stability of selected nascent chains. Our data indicates regulation of post-transcriptional modifications on mitochondrial tRNAs is finely tuned for the control of mitochondrial gene expression. Collectively, our findings provide novel insight into the regulation of mitochondrial tRNAs and reveal greater complexity to the molecular pathogenesis of MERRF.

Pseudouridines have context-dependent mutation and stop rates in high-throughput sequencing.

The abundant RNA modification pseudouridine (Ψ) has been mapped transcriptome-wide by chemically modifying pseudouridines with carbodiimide and detecting the resulting reverse transcription stops in high-throughput sequencing. However, these methods have limited sensitivity and specificity, in part due to the use of reverse transcription stops. We sought to use mutations rather than just stops in sequencing data to identify pseudouridine sites. Here, we identify reverse transcription conditions that allow read-through of carbodiimide-modified pseudouridine (CMC-Ψ), and we show that pseudouridines in carbodiimide-treated human ribosomal RNA have context-dependent mutation and stop rates in high-throughput sequencing libraries prepared under these conditions. Furthermore, accounting for the context-dependence of mutation and stop rates can enhance the detection of pseudouridine sites. Similar approaches could contribute to the sequencing-based detection of many RNA modifications.

Determination of tRNA aminoacylation levels by high-throughput sequencing.

Transfer RNA (tRNA) decodes mRNA codons when aminoacylated (charged) with an amino acid at its 3' end. Charged tRNAs turn over rapidly in cells, and variations in charged tRNA fractions are known to be a useful parameter in cellular responses to stress. tRNA charging fractions can be measured for individual tRNA species using acid denaturing gels, or comparatively at the genome level using microarrays. These hybridization-based approaches cannot be used for high resolution analysis of mammalian tRNAs due to their large sequence diversity. Here we develop a high-throughput sequencing method that enables accurate determination of charged tRNA fractions at single-base resolution (Charged DM-tRNA-seq). Our method takes advantage of the recently developed DM-tRNA-seq method, but includes additional chemical steps that specifically remove the 3'A residue in uncharged tRNA. Charging fraction is obtained by counting the fraction of A-ending reads versus A+C-ending reads for each tRNA species in the same sequencing reaction. In HEK293T cells, most cytosolic tRNAs are charged at >80% levels, whereas tRNASer and tRNAThr are charged at lower levels. These low charging levels were validated using acid denaturing gels. Our method should be widely applicable for investigations of tRNA charging as a parameter in biological regulation.

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.

tRNA base methylation identification and quantification via high-throughput sequencing.

Eukaryotic transfer RNAs contain on average 14 modifications. Investigations of their biological functions require the determination of the modification sites and the dynamic variations of the modification fraction. Base methylation represents a major class of tRNA modification. Although many approaches have been used to identify tRNA base methylations, including sequencing, they are generally qualitative and do not report the information on the modification fraction. Dynamic mRNA modifications have been shown to play important biological roles; yet, the extent of tRNA modification fractions has not been reported systemically. Here we take advantage of a recently developed high-throughput sequencing method (DM-tRNA-seq) to identify and quantify tRNA base methylations located at the Watson-Crick face in HEK293T cells at single base resolution. We apply information derived from both base mutations and positional stops from sequencing using a combination of demethylase treatment and cDNA synthesis by a thermophilic reverse transcriptase to compile a quantitative "Modification Index" (MI) for six base methylations in human tRNA and rRNA. MI combines the metrics for mutational and stop components from alignment of sequencing data without demethylase treatment, and the modifications are validated in the sequencing data upon demethylase treatment. We identify many new methylation sites in both human nuclear and mitochondrial-encoded tRNAs not present in the RNA modification databases. The potentially quantitative nature of the MI values obtained from sequencing is validated by primer extension of several tRNAs. Our approach should be widely applicable to identify tRNA methylation sites, analyze comparative fractional modifications, and evaluate the modification dynamics between different samples.

The dynamic N(1)-methyladenosine methylome in eukaryotic messenger RNA.

Gene expression can be regulated post-transcriptionally through dynamic and reversible RNA modifications. A recent noteworthy example is N(6)-methyladenosine (m(6)A), which affects messenger RNA (mRNA) localization, stability, translation and splicing. Here we report on a new mRNA modification, N(1)-methyladenosine (m(1)A), that occurs on thousands of different gene transcripts in eukaryotic cells, from yeast to mammals, at an estimated average transcript stoichiometry of 20% in humans. Employing newly developed sequencing approaches, we show that m(1)A is enriched around the start codon upstream of the first splice site: it preferentially decorates more structured regions around canonical and alternative translation initiation sites, is dynamic in response to physiological conditions, and correlates positively with protein production. These unique features are highly conserved in mouse and human cells, strongly indicating a functional role for m(1)A in promoting translation of methylated mRNA.

A dual fluorescent reporter for the investigation of methionine mistranslation in live cells.

In mammalian cells under oxidative stress, the methionyl-tRNA synthetase (MetRS) misacylates noncognate tRNAs at frequencies as high as 10% distributed among up to 28 tRNA species. Instead of being detrimental for the cell, misincorporation of methionine residues in the proteome reduces the risk of oxidative damage to proteins, which aids the oxidative stress response. tRNA microarrays have been essential for the detection of the full pattern of misacylated tRNAs, but have limited capacity to investigate the misacylation and mistranslation mechanisms in live cells. Here we develop a dual-fluorescence reporter to specifically measure methionine misincorporation at glutamic acid codons GAA and GAG via tRNA(Glu) mismethionylation in human cells. Our method relies on mutating a specific Met codon in the active site of the fluorescent protein mCherry to a Glu codon that renders mCherry nonfluorescent when translation follows the genetic code. Mistranslation utilizing mismethionylated tRNA(Glu) restores fluorescence in proportion to the amount of misacylated tRNA(Glu). This cellular approach works well for both transient transfection and established stable HEK293 lines. It is rapid, straightforward, and well suited for high-throughput activity analysis under a wide range of physiological conditions. As a proof of concept, we apply this method to characterize the effect of human tRNA(Glu) isodecoders on mistranslation and discuss the implications of our findings.

Efficient and quantitative high-throughput tRNA sequencing.

Despite its biological importance, tRNA has not been adequately sequenced by standard methods because of its abundant post-transcriptional modifications and stable structure, which interfere with cDNA synthesis. We achieved efficient and quantitative tRNA sequencing in HEK293T cells by using engineered demethylases to remove base methylations and a highly processive thermostable group II intron reverse transcriptase to overcome these obstacles. Our method, DM-tRNA-seq, should be applicable to investigations of tRNA in all organisms.