The Versatile Roles of the tRNA Epitranscriptome during Cellular Responses to Toxic Exposures and Environmental Stress

GNN Codon Adjacency Tunes Protein Translation

The central dogma treats the ribosome as a molecular machine that reads one mRNA codon at a time as it adds each amino acid to its growing peptide chain. However, this and previous studies suggest that ribosomes actually perceive pairs of adjacent codons as they take three-nucleotide steps along the mRNA. We examined GNN codons, which we find are surprisingly overrepresented in eukaryote protein-coding open reading frames (ORFs), especially immediately after NNU codons. Ribosome profiling experiments in yeast revealed that ribosomes with NNU at their aminoacyl (A) site have particularly elevated densities when NNU is immediately followed (3') by a GNN codon, indicating slower mRNA threading of the NNU codon from the ribosome's A to peptidyl (P) sites. Moreover, if the assessment was limited to ribosomes that have only recently arrived at the next codon, by examining 21-nucleotide ribosome footprints (21-nt RFPs), elevated densities were observed for multiple codon classes when followed by GNN. This striking translation slowdown at adjacent 5'-NNN GNN codon pairs is likely mediated, in part, by the ribosome's CAR surface, which acts as an extension of the A-site tRNA anticodon during ribosome translocation and interacts through hydrogen bonding and pi stacking with the GNN codon. The functional consequences of 5'-NNN GNN codon adjacency are expected to influence the evolution of protein coding sequences.

Deciphering the Diversity in Bacterial Transporters That Salvage Queuosine Precursors

Queuosine (Q) is a modification of the wobble base of tRNA harboring GUN anticodons with roles in decoding accuracy and efficiency. Its synthesis is complex with multiple enzymatic steps, and several pathway intermediates can be salvaged. The only two transporter families known to salvage Q precursors are QPTR/COG1738 and QrtT/QueT. Analyses of the distribution of known Q synthesis and salvage genes in human gut and oral microbiota genomes have suggested that more transporter families remain to be found and that Q precursor exchanges must occur within the structured microenvironments of the mammalian host. Using physical clustering and fusion-based association with Q salvage genes, candidate genes for missing transporters were identified and five were tested experimentally by complementation assays in Escherichia coli. Three genes encoding transporters from three different Pfam families, a ureide permease (PF07168) from Acidobacteriota bacterium, a hemolysin III family protein (PF03006) from Bifidobacterium breve, and a Major Facilitator Superfamily protein (PF07690) from Bartonella henselae, were found to allow the transport of both preQ0 and preQ1 in this heterologous system. This work suggests that many transporter families can evolve to transport Q precursors, reinforcing the concept of transporter plasticity.

Transfer RNA modifications and cellular thermotolerance

RNA molecules are modified post-transcriptionally to acquire their diverse functions. Transfer RNA (tRNA) has the widest variety and largest numbers of RNA modifications. tRNA modifications are pivotal for decoding the genetic code and stabilizing the tertiary structure of tRNA molecules. Alternation of tRNA modifications directly modulates the structure and function of tRNAs and regulates gene expression. Notably, thermophilic organisms exhibit characteristic tRNA modifications that are dynamically regulated in response to varying growth temperatures, thereby bolstering fitness in extreme environments. Here, we review the history and latest findings regarding the functions and biogenesis of several tRNA modifications that contribute to the cellular thermotolerance of thermophiles.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA (mRNA) codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codons' decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we identified that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Direct Nanopore Sequencing for the 17 RNA Modification Types in 36 Locations in the E. coli Ribosome Enables Monitoring of Stress-Dependent Changes

The bacterium Escherichia coli possesses 16S and 23S rRNA strands that have 36 chemical modification sites with 17 different structures. Nanopore direct RNA sequencing using a protein nanopore sensor and helicase brake, which is also a sensor, was applied to the rRNAs. Nanopore current levels, base calling profile, and helicase dwell times for the modifications relative to unmodified synthetic rRNA controls found signatures for nearly all modifications. Signatures for clustered modifications were determined by selective sequencing of writer knockout E. coli and sequencing of synthetic RNAs utilizing some custom-synthesized nucleotide triphosphates for their preparation. The knowledge of each modification's signature, apart from 5-methylcytidine, was used to determine how metabolic and cold-shock stress impact rRNA modifications. Metabolic stress resulted in either no change or a decrease, and one site increased in modification occupancy, while cold-shock stress led to either no change or a decrease. The double modification m4Cm1402 resides in 16S rRNA, and it decreased with both stressors. Using the helicase dwell time, it was determined that the N4 methyl group is lost during both stressors, and the 2'-OMe group remained. In the ribosome, this modification stabilizes binding to the mRNA codon at the P-site resulting in increased translational fidelity that is lost during stress. The E. coli genome has seven rRNA operons (rrn), and the earlier studies aligned the nanopore reads to a single operon (rrnA). Here, the reads were aligned to all seven operons to identify operon-specific changes in the 11 pseudouridines. This study demonstrates that direct sequencing for >16 different RNA modifications in a strand is achievable.

GC-MS based untargeted metabolomics reveals the metabolic response of earthworm (Eudrilus eugeniae) after chronic combinatorial exposure to three different pesticides

In this study GC-MS-based untargeted metabolomics was used to identify the metabolic response of earthworm; Eudrilus eugeniae exposed to sub-lethal concentrations of chlorpyrifos-CHL, cypermethrin-CYP, Glyphosate-GLY, and Combined-C (all three pesticides) at the concentrations of 3, 6, and 12 mg/kg. Principal component analysis of the obtained datasets revealed a clear distinction between the control and treatment groups. The mean weight of the worms in the treated groups decreased significantly (p < 0.05). Among the identified metabolites, oleic acid (~ 93.47%), lysine (~ 92.20%), glutamic acid (~ 91.81%), leucine (~ 90.20%), asparagine (~ 94.20%), methionine (~ 92.27%), malic acid (~ 93.37%), turanose (~ 95.04%), maltose (~ 92.36%), cholesta-3,5-diene (~ 86.11%), galactose (~ 93.20%), cholesterol (~ 91.56%), tocopherol (~ 85.09%), decreased significantly (p < 0.05), whereas myoinositol (~ 83%) and isoleucine (78.09%) increased significantly (p < 0.05) upon exposure to the CHL, CYP, GLY, and C. Overall, the findings suggest that earthworms might be a new entry point for the pesticides into the food chain. The present study highlights that metabolomics can be a reliable approach to understand the effect of different xenobiotics including pesticides on the metabolic response of earthworms.

A precision environmental health approach to prevention of human disease

Human health is determined by the interaction of our environment with the genome, epigenome, and microbiome, which shape the transcriptomic, proteomic, and metabolomic landscape of cells and tissues. Precision environmental health is an emerging field leveraging environmental and system-level ('omic) data to understand underlying environmental causes of disease, identify biomarkers of exposure and response, and develop new prevention and intervention strategies. In this article we provide real-life illustrations of the utility of precision environmental health approaches, identify current challenges in the field, and outline new opportunities to promote health through a precision environmental health framework.

Paternal methotrexate exposure affects sperm small RNA content and causes craniofacial defects in the offspring

Folate is an essential vitamin for vertebrate embryo development. Methotrexate (MTX) is a folate antagonist that is widely prescribed for autoimmune diseases, blood and solid organ malignancies, and dermatologic diseases. Although it is highly contraindicated for pregnant women, because it is associated with an increased risk of multiple birth defects, the effect of paternal MTX exposure on their offspring has been largely unexplored. Here, we found MTX treatment of adult medaka male fish (Oryzias latipes) causes cranial cartilage defects in their offspring. Small non-coding RNA (sncRNAs) sequencing in the sperm of MTX treated males identify differential expression of a subset of tRNAs, with higher abundance for specific 5' tRNA halves. Sperm RNA methylation analysis on MTX treated males shows that m5C is the most abundant and differential modification found in RNAs ranging in size from 50 to 90 nucleotides, predominantly tRNAs, and that it correlates with greater testicular Dnmt2 methyltransferase expression. Injection of sperm small RNA fractions from MTX-treated males into normal fertilized eggs generated cranial cartilage defects in the offspring. Overall, our data suggest that paternal MTX exposure alters sperm sncRNAs expression and modifications that may contribute to developmental defects in their offspring.

Hypoxia induces alterations in tRNA modifications involved in translational control

BACKGROUND

Adaptation to high-altitude hypobaric hypoxia has been shown to require a set of physiological traits enabled by an associated set of genetic modifications, as well as transcriptome regulation. These lead to both lifetime adaptation of individuals to hypoxia at high altitudes and generational evolution of populations as seen for instance in those of Tibet. Additionally, RNA modifications, which are sensitive to environmental exposure, have been shown to play pivotal biological roles in maintaining the physiological functions of organs. However, the dynamic RNA modification landscape and related molecular mechanisms in mouse tissues under hypobaric hypoxia exposure remain to be fully understood. Here, we explore the tissue-specific distribution pattern of multiple RNA modifications across mouse tissues.

RESULTS

By applying an LC-MS/MS-dependent RNA modification detection platform, we identified the distribution of multiple RNA modifications in total RNA, tRNA-enriched fragments, and 17-50-nt sncRNAs across mouse tissues; these patterns were associated with the expression levels of RNA modification modifiers in different tissues. Moreover, the tissue-specific abundance of RNA modifications was sensitively altered across different RNA groups in a simulated high-altitude (over 5500 m) hypobaric hypoxia mouse model with the activation of the hypoxia response in mouse peripheral blood and multiple tissues. RNase digestion experiments revealed that the alteration of RNA modification abundance under hypoxia exposure impacted the molecular stability of tissue total tRNA-enriched fragments and isolated individual tRNAs, such as tRNA^Ala, tRNA^val, tRNA^Glu, and tRNA^Leu. In vitro transfection experiments showed that the transfection of testis total tRNA-enriched fragments from the hypoxia group into GC-2spd cells attenuated the cell proliferation rate and led to a reduction in overall nascent protein synthesis in cells.

CONCLUSIONS

Our results reveal that the abundance of RNA modifications for different classes of RNAs under physiological conditions is tissue-specific and responds to hypobaric hypoxia exposure in a tissue-specific manner. Mechanistically, the dysregulation of tRNA modifications under hypobaric hypoxia attenuated the cell proliferation rate, facilitated the sensitivity of tRNA to RNases, and led to a reduction in overall nascent protein synthesis, suggesting an active role of tRNA epitranscriptome alteration in the adaptive response to environmental hypoxia exposure.

Mitochondrial RNA methyltransferase TRMT61B is a new, potential biomarker and therapeutic target for highly aneuploid cancers

Despite being frequently observed in cancer cells, chromosomal instability (CIN) and its immediate consequence, aneuploidy, trigger adverse effects on cellular homeostasis that need to be overcome by anti-stress mechanisms. As such, these safeguard responses represent a tumor-specific Achilles heel, since CIN and aneuploidy are rarely observed in normal cells. Recent data have revealed that epitranscriptomic marks catalyzed by RNA-modifying enzymes change under various stress insults. However, whether aneuploidy is associated with such RNA modifying pathways remains to be determined. Through an in silico search for aneuploidy biomarkers in cancer cells, we found TRMT61B, a mitochondrial RNA methyltransferase enzyme, to be associated with high levels of aneuploidy. Accordingly, TRMT61B protein levels are increased in tumor cell lines with an imbalanced karyotype as well as in different tumor types when compared to control tissues. Interestingly, while TRMT61B depletion induces senescence in melanoma cell lines with low levels of aneuploidy, it leads to apoptosis in cells with high levels. The therapeutic potential of these results was further validated by targeting TRMT61B in transwell and xenografts assays. We show that TRM61B depletion reduces the expression of several mitochondrial encoded proteins and limits mitochondrial function. Taken together, these results identify a new biomarker of aneuploidy in cancer cells that could potentially be used to selectively target highly aneuploid tumors.

DNA methylation patterns associated with konzo in Sub-Saharan Africa

Konzo, a disease characterized by sudden, irreversible spastic paraparesis, affecting up to 10% of the population in some regions of Sub-Saharan Africa during outbreaks, is strongly associated with dietary exposure to cyanogenic bitter cassava. The molecular mechanisms underlying the development of konzo remain largely unknown. Here, through an analysis of 16 individuals with konzo and matched healthy controls from the same outbreak zones, we identified 117 differentially methylated loci involved in numerous biological processes that may identify cyanogenic-sensitive regions of the genome, providing the first study of epigenomic alterations associated with a clinical phenotype of konzo.

The Role of Nuclear-Encoded Mitochondrial tRNA Charging Enzymes in Human Inherited Disease

Aminoacyl-tRNA synthetases (ARSs) are highly conserved essential enzymes that charge tRNA with cognate amino acids-the first step of protein synthesis. Of the 37 nuclear-encoded human ARS genes, 17 encode enzymes are exclusively targeted to the mitochondria (mt-ARSs). Mutations in nuclear mt-ARS genes are associated with rare, recessive human diseases with a broad range of clinical phenotypes. While the hypothesized disease mechanism is a loss-of-function effect, there is significant clinical heterogeneity among patients that have mutations in different mt-ARS genes and also among patients that have mutations in the same mt-ARS gene. This observation suggests that additional factors are involved in disease etiology. In this review, we present our current understanding of diseases caused by mutations in the genes encoding mt-ARSs and propose explanations for the observed clinical heterogeneity.

The host tRNA epitranscriptome: A new player in RNA virus infections

Viruses completely depend on the host translation machineries to express the viral proteins. Recent data reveal an unprecedented interaction of positive strand RNA ((+)RNA) viruses with the host tRNA epitranscriptome to favor viral protein expression via a specific reprogramming of codon optimality that ultimately favors decoding of the viral codons. We propose that this feature is shared by multiple RNA viruses and that the involved tRNA modifying enzymes represent promising novel targets for the development of broad-spectrum antivirals.

Selenoproteins and the senescence-associated epitranscriptome

Selenium is a naturally found trace element, which provides multiple benefits including antioxidant, anticancer, and antiaging, as well as boosting immunity. One unique feature of selenium is its incorporation as selenocysteine, a rare 21st amino acid, into selenoproteins. Twenty-five human selenoproteins have been discovered, and a majority of these serve as crucial antioxidant enzymes for redox homeostasis. Unlike other amino acids, incorporation of selenocysteine requires a distinctive UGA stop codon recoding mechanism. Although many studies correlating selenium, selenoproteins, aging, and senescence have been performed, it has not yet been explored if the upstream events regulating selenoprotein synthesis play a role in senescence-associated pathologies. The epitranscriptomic writer alkylation repair homolog 8 (ALKBH8) is critical for selenoprotein production, and its deficiency can significantly decrease levels of selenoproteins that are essential for reactive oxygen species (ROS) detoxification, and increase oxidative stress, one of the major drivers of cellular senescence. Here, we review the potential role of epitranscriptomic marks that govern selenocysteine utilization in regulating the senescence program.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

Weight Change Alters the Small RNA Profile of Urinary Extracellular Vesicles in Obesity

INTRODUCTION

Various kidney diseases reportedly show different urinary extracellular vesicle (EV) RNA profiles. Although obesity is one of the main causes of chronic kidney disease, the expression pattern of urinary EV RNA in obesity is uncertain. Our aim was to sequence the small RNA profiles of urinary EVs in obese patients before and after weight reduction and compare them to those of healthy volunteers (HVs).

METHODS

We recruited age-sex-matched obese patients and HVs. The small RNA profiles of urinary EVs were analyzed using RNA sequencing. To evaluate the effect of weight reduction, small RNA profiles of urinary EVs 6 months after bariatric surgery were also analyzed.

RESULTS

The proportion of urinary EVs transfer RNA and microRNA of obese patients differed from that of HVs. Obese patients showed differential expression of 1,343 small RNAs in urinary EVs compared to HVs (fold change ≥2 and p value <0.05). Among those, 61 small RNAs were upregulated in obese patients and downregulated after weight reduction, whereas 167 small RNAs were downregulated in obese patients and upregulated after weight reduction. RNA sequencing revealed the correlation between the specific urinary EV small RNAs and clinical parameters including body weight, low-density lipoprotein cholesterol, triglyceride, high-density lipoprotein cholesterol, serum glucose, estimated glomerular filtration rate, and albuminuria.

CONCLUSION

Obese patients showed distinct urinary EV small RNA profiles compared to HVs. Weight reduction altered urinary EV small-RNA profiles in obese patients.

Dynamic queuosine changes in tRNA couple nutrient levels to codon choice in Trypanosoma brucei

Every type of nucleic acid in cells undergoes programmed chemical post-transcriptional modification. Generally, modification enzymes use substrates derived from intracellular metabolism, one exception is queuine (q)/queuosine (Q), which eukaryotes obtain from their environment; made by bacteria and ultimately taken into eukaryotic cells via currently unknown transport systems. Here, we use a combination of molecular, cell biology and biophysical approaches to show that in Trypanosoma brucei tRNA Q levels change dynamically in response to concentration variations of a sub-set of amino acids in the growth media. Most significant were variations in tyrosine, which at low levels lead to increased Q content for all the natural tRNAs substrates of tRNA-guanine transglycosylase (TGT). Such increase results from longer nuclear dwell time aided by retrograde transport following cytoplasmic splicing. In turn high tyrosine levels lead to rapid decrease in Q content. Importantly, the dynamic changes in Q content of tRNAs have negligible effects on global translation or growth rate but, at least, in the case of tRNA^Tyr it affected codon choice. These observations have implications for the occurrence of other tunable modifications important for ‘normal’ growth, while connecting the intracellular localization of modification enzymes, metabolites and tRNAs to codon selection and implicitly translational output.

Variations in transfer and ribosomal RNA epitranscriptomic status can adapt eukaryote translation to changing physiological and environmental conditions

The timely reprogramming of gene expression in response to internal and external cues is essential to eukaryote development and acclimation to changing environments. Chemically modifying molecular receptors and transducers of these signals is one way to efficiently induce proper physiological responses. Post-translation modifications, regulating protein biological activities, are central to many well-known signal-responding pathways. Recently, messenger RNA (mRNA) chemical (i.e. epitranscriptomic) modifications were also shown to play a key role in these processes. In contrast, transfer RNA (tRNA) and ribosomal RNA (rRNA) chemical modifications, although critical for optimal function of the translation apparatus, and much more diverse and quantitatively important compared to mRNA modifications, were until recently considered as mainly static chemical decorations. We present here recent observations that are challenging this view and supporting the hypothesis that tRNA and rRNA modifications dynamically respond to various cell and environmental conditions and contribute to adapt translation to these conditions.

The occurrence order and cross-talk of different tRNA modifications

Chemical modifications expand the composition of RNA molecules from four standard nucleosides to over 160 modified nucleosides, which greatly increase the complexity and utility of RNAs. Transfer RNAs (tRNAs) are the most heavily modified cellular RNA molecules and contain the largest variety of modifications. Modification of tRNAs is pivotal for protein synthesis and also precisely regulates the noncanonical functions of tRNAs. Defects in tRNA modifications lead to numerous human diseases. Up to now, more than 100 types of modifications have been found in tRNAs. Intriguingly, some modifications occur widely on all tRNAs, while others only occur on a subgroup of tRNAs or even only a specific tRNA. The modification frequency of each tRNA is approximately 7% to 25%, with 5-20 modification sites present on each tRNA. The occurrence and modulation of tRNA modifications are specifically noticeable as plenty of interplays among different sites and modifications have been discovered. In particular, tRNA modifications are responsive to environmental changes, indicating their dynamic and highly organized nature. In this review, we summarized the known occurrence order, cross-talk, and cooperativity of tRNA modifications.

Adaptor Molecules Epitranscriptome Reprograms Bacterial Pathogenicity

The strong decoration of tRNAs with post-transcriptional modifications provides an unprecedented adaptability of this class of non-coding RNAs leading to the regulation of bacterial growth and pathogenicity. Accumulating data indicate that tRNA post-transcriptional modifications possess a central role in both the formation of bacterial cell wall and the modulation of transcription and translation fidelity, but also in the expression of virulence factors. Evolutionary conserved modifications in tRNA nucleosides ensure the proper folding and stability redounding to a totally functional molecule. However, environmental factors including stress conditions can cause various alterations in tRNA modifications, disturbing the pathogen homeostasis. Post-transcriptional modifications adjacent to the anticodon stem-loop, for instance, have been tightly linked to bacterial infectivity. Currently, advances in high throughput methodologies have facilitated the identification and functional investigation of such tRNA modifications offering a broader pool of putative alternative molecular targets and therapeutic avenues against bacterial infections. Herein, we focus on tRNA epitranscriptome shaping regarding modifications with a key role in bacterial infectivity including opportunistic pathogens of the human microbiome.

Lemon exosome-like nanoparticles enhance stress survival of gut bacteria by RNase P-mediated specific tRNA decay

Diet and bile play critical roles in shaping gut microbiota, but the molecular mechanism underlying interplay with intestinal microbiota is unclear. Here, we showed that lemon-derived exosome-like nanoparticles (LELNs) enhance lactobacilli toleration to bile. To decipher the mechanism, we used Lactobacillus rhamnosus GG (LGG) as proof of concept to show that LELNs enhance LGG bile resistance via limiting production of Msp1 and Msp3, resulting in decrease of bile accessibility to cell membrane. Furthermore, we found that decline of Msps protein levels was regulated through specific tRNAser UCC and tRNAser UCG decay. We identified RNase P, an essential housekeeping endonuclease, being responsible for LELNs-induced tRNAser UCC and tRNAser UCG decay. We further identified galacturonic acid-enriched pectin-type polysaccharide as the active factor in LELNs to increase bile resistance and downregulate tRNAser UCC and tRNAser UCG level in the LGG. Our study demonstrates a tRNA-based gene expression regulation mechanism among lactobacilli to increase bile resistance.

Epitranscriptomic profile of Lactobacillus agilis and its adaptation to growth on inulin

Objective

Ribonucleic acids (RNA) are involved in many cellular functions. In general, RNA is made up by only four different ribonucleotides. The modifications of RNA (epitranscriptome) can greatly enhance the structural diversity of RNA, which in turn support some of the RNA functions. To determine whether the epitranscriptome of a specific probiotic is associated with its adaptation to the source of energy, Lactobacillus agilis (YZ050) was selected as a model and its epitranscriptome was profiled and compared by using mass spectrometry.

Results

The L. agilis epitranscriptome (minus rRNA modifications) consists of 17 different RNA modifications. By capturing the L. agilis cells during exponential growth, reproducible profiling was achieved. In a comparative study, the standard source of energy (glucose) in the medium was substituted by a prebiotic inulin, and a downward trend in the L. agilis epitranscriptome was detected. This marks the first report on a system-wide variation of a bacterial epitranscriptome that resulted from adapting to an alternative energy source. No correlation was found between the down-regulated RNA modifications and the expression level of corresponding writer genes. Whereas, the expression level of a specific exonuclease gene, RNase J1, was detected to be higher in cells grown on inulin.

m5U54 tRNA Hypomodification by Lack of TRMT2A Drives the Generation of tRNA-Derived Small RNAs

Transfer RNA (tRNA) molecules contain various post-transcriptional modifications that are crucial for tRNA stability, translation efficiency, and fidelity. Besides their canonical roles in translation, tRNAs also originate tRNA-derived small RNAs (tsRNAs), a class of small non-coding RNAs with regulatory functions ranging from translation regulation to gene expression control and cellular stress response. Recent evidence indicates that tsRNAs are also modified, however, the impact of tRNA epitranscriptome deregulation on tsRNAs generation is only now beginning to be uncovered. The 5-methyluridine (m⁵U) modification at position 54 of cytosolic tRNAs is one of the most common and conserved tRNA modifications among species. The tRNA methyltransferase TRMT2A catalyzes this modification, but its biological role remains mostly unexplored. Here, we show that TRMT2A knockdown in human cells induces m⁵U54 tRNA hypomodification and tsRNA formation. More specifically, m⁵U54 hypomodification is followed by overexpression of the ribonuclease angiogenin (ANG) that cleaves tRNAs near the anticodon, resulting in accumulation of 5′tRNA-derived stress-induced RNAs (5′tiRNAs), namely 5′tiRNA-Gly^GCC and 5′tiRNA-Glu^CTC, among others. Additionally, transcriptomic analysis confirms that down-regulation of TRMT2A and consequently m⁵U54 hypomodification impacts the cellular stress response and RNA stability, which is often correlated with tiRNA generation. Accordingly, exposure to oxidative stress conditions induces TRMT2A down-regulation and tiRNA formation in mammalian cells. These results establish a link between tRNA hypomethylation and ANG-dependent tsRNAs formation and unravel m⁵U54 as a tRNA cleavage protective mark.

Arginine Methylation Regulates Ribosome CAR Function

The ribosome CAR interaction surface is hypothesized to provide a layer of translation regulation through hydrogen-bonding to the +1 mRNA codon that is next to enter the ribosome A site during translocation. The CAR surface consists of three residues, 16S/18S rRNA C1054, A1196 (E. coli 16S numbering), and R146 of yeast ribosomal protein Rps3. R146 can be methylated by the Sfm1 methyltransferase which is downregulated in stressed cells. Through molecular dynamics analysis, we show here that methylation of R146 compromises the integrity of CAR by reducing the cation-pi stacking of the R146 guanidinium group with A1196, leading to reduced CAR hydrogen-bonding with the +1 codon. We propose that ribosomes assembled under stressed conditions have unmethylated R146, resulting in elevated CAR/+1 codon interactions, which tunes translation levels in response to the altered cellular context.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

A Method to Monitor the Introduction of Posttranscriptional Modifications in tRNAs with NMR Spectroscopy

During their biosynthesis, transfer RNAs (tRNAs) are decorated with a large number of posttranscriptional chemical modifications. Methods to directly detect the introduction of posttranscriptional modifications during tRNA maturation are rare and do not provide information on the temporality of modification events. Here, we report a methodology, using NMR as a tool to monitor tRNA maturation in a nondisruptive and continuous fashion in cellular extracts. This method requires the production of substrate tRNA transcripts devoid of modifications and active cell extracts containing the necessary cellular enzymatic activities to modify RNA. The present protocol describes these different aspects of our method and reports the time-resolved NMR monitoring of the yeast tRNA^Phe maturation as an example. The NMR-based methodology presented here could be adapted to investigate diverse features in tRNA maturation.

Broad-range RNA modification analysis of complex biological samples using rapid C18-UPLC-MS

Post-transcriptional RNA modifications play an important role in cellular metabolism with homoeostatic disturbances manifesting as a wide repertoire of phenotypes, reduced stress tolerance and translational perturbation, developmental defects, and diseases, such as type II diabetes, leukaemia, and carcinomas. Hence, there has been an intense effort to develop various methods for investigating RNA modifications and their roles in various organisms, including sequencing-based approaches and, more frequently, liquid chromatography–mass spectrometry (LC-MS)-based methods. Although LC-MS offers numerous advantages, such as being highly sensitive and quantitative over a broad detection range, some stationary phase chemistries struggle to resolve positional isomers. Furthermore, the demand for detailed analyses of complex biological samples often necessitates long separation times, hampering sample-to-sample turnover and making multisample analyses time consuming. To overcome this limitation, we have developed an ultra-performance LC-MS (UPLC-MS) method that uses an octadecyl carbon chain (C18)-bonded silica matrix for the efficient separation of 50 modified ribonucleosides, including positional isomers, in a single 9-min sample-to-sample run. To validate the performance and versatility of our method, we analysed tRNA modification patterns of representative microorganisms from each domain of life, namely Archaea (Methanosarcina acetivorans), Bacteria (Pseudomonas syringae), and Eukarya (Saccharomyces cerevisiae). Additionally, our method is flexible and readily applicable for detection and relative quantification using stable isotope labelling and targeted approaches like multiple reaction monitoring (MRM). In conclusion, this method represents a fast and robust tool for broad-range exploration and quantification of ribonucleosides, facilitating future homoeostasis studies of RNA modification in complex biological samples.

Bacterial tRNA 2'-O-methylation is dynamically regulated under stress conditions and modulates innate immune response

RNA modifications are a well-recognized way of gene expression regulation at the post-transcriptional level. Despite the importance of this level of regulation, current knowledge on modulation of tRNA modification status in response to stress conditions is far from being complete. While it is widely accepted that tRNA modifications are rather dynamic, such variations are mostly assessed in terms of total tRNA, with only a few instances where changes could be traced to single isoacceptor species. Using Escherichia coli as a model system, we explored stress-induced modulation of 2'-O-methylations in tRNAs by RiboMethSeq. This analysis and orthogonal analytical measurements by LC-MS show substantial, but not uniform, increase of the Gm18 level in selected tRNAs under mild bacteriostatic antibiotic stress, while other Nm modifications remain relatively constant. The absence of Gm18 modification in tRNAs leads to moderate alterations in E. coli mRNA transcriptome, but does not affect polysomal association of mRNAs. Interestingly, the subset of motility/chemiotaxis genes is significantly overexpressed in ΔTrmH mutant, this corroborates with increased swarming motility of the mutant strain. The stress-induced increase of tRNA Gm18 level, in turn, reduced immunostimulation properties of bacterial tRNAs, which is concordant with the previous observation that Gm18 is a suppressor of Toll-like receptor 7 (TLR7)-mediated interferon release. This documents an effect of stress induced modulation of tRNA modification that acts outside protein translation.

Enhanced differentiation of isomeric RNA modifications by reducing the size of ions in ion mobility mass spectrometric measurements

With the ability to differentiate different molecular sizes, ion mobility spectrometry (IMS) has great potentials in the analysis of isomeric compounds. However, due to the lack of sensitivity and resolution, IMS has not been commonly used. To address the issue on resolution, the goals of this study are to explore a more effective way to perform IMS by reducing the size of ions prior to the IM measurements, and apply the new approach to the differentiation of isomeric RNA modifications. The size reduction of ribonucleoside ions was effectively accomplished by using the collision-induced dissociation process, in which the N-glycosidic bond in ribonucleoside was cleaved and split the ions into two parts—a smaller nucleobase ion and a neutral molecule of ribose sugar. Since the chemical group that corresponds to most of the RNA modifications makes up a relatively small part of the molecular structure of nucleobases, the differentiation of the dissociated nucleobase ions is expected to require a lower ion mobility resolution than the differentiation of bigger isomeric ribonucleoside ions. By using RNA methylation as a model in this study, the proposed method lowered the required resolution by 16% for the differentiation of 1-methyladenosine and N6-methyladenosine. Similar results were also obtained from the differentiation of methylated cytidine isomers. In comparison to the results obtained from using the conventional tandem mass spectrometric method, there was no significant loss of signals when the proposed method was used. The proposed method is expected to be applicable to other types of isomeric compounds. Also, the same approach is applicable on other IMS platforms.

Characterization of UVA-Induced Alterations to Transfer RNA Sequences

Ultraviolet radiation (UVR) adversely affects the integrity of DNA, RNA, and their nucleoside modifications. By employing liquid chromatography–tandem mass spectrometry (LC–MS/MS)-based RNA modification mapping approaches, we identified the transfer RNA (tRNA) regions most vulnerable to photooxidation. Photooxidative damage to the anticodon and variable loop regions was consistently observed in both modified and unmodified sequences of tRNA upon UVA (λ 370 nm) exposure. The extent of oxidative damage measured in terms of oxidized guanosine, however, was higher in unmodified RNA compared to its modified version, suggesting an auxiliary role for nucleoside modifications. The type of oxidation product formed in the anticodon stem–loop region varied with the modification type, status, and whether the tRNA was inside or outside the cell during exposure. Oligonucleotide-based characterization of tRNA following UVA exposure also revealed the presence of novel photoproducts and stable intermediates not observed by nucleoside analysis alone. This approach provides sequence-specific information revealing potential hotspots for UVA-induced damage in tRNAs.

Regulation and Function of RNA Pseudouridylation in Human Cells

Recent advances in pseudouridine detection reveal a complex pseudouridine landscape that includes messenger RNA and diverse classes of noncoding RNA in human cells. The known molecular functions of pseudouridine, which include stabilizing RNA conformations and destabilizing interactions with varied RNA-binding proteins, suggest that RNA pseudouridylation could have widespread effects on RNA metabolism and gene expression. Here, we emphasize how much remains to be learned about the RNA targets of human pseudouridine synthases, their basis for recognizing distinct RNA sequences, and the mechanisms responsible for regulated RNA pseudouridylation. We also examine the roles of noncoding RNA pseudouridylation in splicing and translation and point out the potential effects of mRNA pseudouridylation on protein production, including in the context of therapeutic mRNAs.

The Genetic Architecture of Carbon Tetrachloride-Induced Liver Fibrosis in Mice

BACKGROUND & AIMS

Liver fibrosis is a multifactorial trait that develops in response to chronic liver injury. Our aim was to characterize the genetic architecture of carbon tetrachloride (CCl₄)-induced liver fibrosis using the Hybrid Mouse Diversity Panel, a panel of more than 100 genetically distinct mouse strains optimized for genome-wide association studies and systems genetics.

METHODS

Chronic liver injury was induced by CCl₄ injections twice weekly for 6 weeks. Four hundred thirty-seven mice received CCl₄ and 256 received vehicle, after which animals were euthanized for liver histology and gene expression. Using automated digital image analysis, we quantified fibrosis as the collagen proportionate area of the whole section, excluding normal collagen.

RESULTS

We discovered broad variation in fibrosis among the Hybrid Mouse Diversity Panel strains, demonstrating a significant genetic influence. Genome-wide association analyses revealed significant and suggestive loci underlying susceptibility to fibrosis, some of which overlapped with loci identified in mouse crosses and human population studies. Liver global gene expression was assessed by RNA sequencing across the strains, and candidate genes were identified using differential expression and expression quantitative trait locus analyses. Gene set enrichment analyses identified the underlying pathways, of which stellate cell involvement was prominent, and coexpression network modeling identified modules associated with fibrosis.

CONCLUSIONS

Our results provide a rich resource for the design of experiments to understand mechanisms underlying fibrosis and for rational strain selection when testing antifibrotic drugs.

Hypomodified tRNA in evolutionarily distant yeasts can trigger rapid tRNA decay to activate the general amino acid control response, but with different consequences

All tRNAs are extensively modified, and modification deficiency often results in growth defects in the budding yeast Saccharomyces cerevisiae and neurological or other disorders in humans. In S. cerevisiae, lack of any of several tRNA body modifications results in rapid tRNA decay (RTD) of certain mature tRNAs by the 5'-3' exonucleases Rat1 and Xrn1. As tRNA quality control decay mechanisms are not extensively studied in other eukaryotes, we studied trm8Δ mutants in the evolutionarily distant fission yeast Schizosaccharomyces pombe, which lack 7-methylguanosine at G46 (m7G46) of their tRNAs. We report here that S. pombe trm8Δ mutants are temperature sensitive primarily due to decay of tRNATyr(GUA) and that spontaneous mutations in the RAT1 ortholog dhp1+ restored temperature resistance and prevented tRNA decay, demonstrating conservation of the RTD pathway. We also report for the first time evidence linking the RTD and the general amino acid control (GAAC) pathways, which we show in both S. pombe and S. cerevisiae. In S. pombe trm8Δ mutants, spontaneous GAAC mutations restored temperature resistance and tRNA levels, and the trm8Δ temperature sensitivity was precisely linked to GAAC activation due to tRNATyr(GUA) decay. Similarly, in the well-studied S. cerevisiae trm8Δ trm4Δ RTD mutant, temperature sensitivity was closely linked to GAAC activation due to tRNAVal(AAC) decay; however, in S. cerevisiae, GAAC mutations increased tRNA loss and exacerbated temperature sensitivity. A similar exacerbated growth defect occurred upon GAAC mutation in S. cerevisiae trm8Δ and other single modification mutants that triggered RTD. Thus, these results demonstrate a conserved GAAC activation coincident with RTD in S. pombe and S. cerevisiae, but an opposite impact of the GAAC response in the two organisms. We speculate that the RTD pathway and its regulation of the GAAC pathway is widely conserved in eukaryotes, extending to other mutants affecting tRNA body modifications.

Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168

Extensive knowledge of both the nature and position of tRNA modifications in all cellular tRNAs has been limited to two bacteria, Escherichia coli and Mycoplasma capricolum. Bacillus subtilis sp subtilis strain 168 is the model Gram-positive bacteria and the list of the genes involved in tRNA modifications in this organism is far from complete. Mass spectrometry analysis of bulk tRNA extracted from B. subtilis, combined with next generation sequencing technologies and comparative genomic analyses, led to the identification of 41 tRNA modification genes with associated confidence scores. Many differences were found in this model Gram-positive bacteria when compared to E. coli. In general, B. subtilis tRNAs are less modified than those in E. coli, even if some modifications, such as m¹A22 or ms²t⁶A, are only found in the model Gram-positive bacteria. Many examples of non-orthologous displacements and of variations in the most complex pathways are described. Paralog issues make uncertain direct annotation transfer from E. coli to B. subtilis based on homology only without further experimental validation. This difficulty was shown with the identification of the B. subtilis enzyme that introduces ψ at positions 31/32 of the tRNAs. This work presents the most up to date list of tRNA modification genes in B. subtilis, identifies the gaps in knowledge, and lays the foundation for further work to decipher the physiological role of tRNA modifications in this important model organism and other bacteria.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Emerging Roles of tRNAs in RNA Virus Infections

Viruses rely on the host cell translation machinery for efficient synthesis of their own proteins. Emerging evidence highlights different roles for host transfer RNAs (tRNAs) in the process of virus replication. For instance, different RNA viruses manipulate host tRNA pools to favor viral protein translation. Interestingly, specific host tRNAs are used as reverse transcription primers and are packaged into retroviral virions. Recent data also demonstrate the formation of tRNA-derived fragments (tRFs) upon infection to facilitate viral replication. Here, we comprehensively discuss how RNA viruses exploit distinct aspects of the host tRNA biology for their benefit. In light of the recent advances in the field, we propose that host tRNA-related pathways and mechanisms represent promising cellular targets for the development of novel antiviral strategies.

Biomarkers of Environmental Toxicants: Exposure and Biological Effects

Biomarkers of environmental toxicants are measures of exposures and effects, some of which can serve to assess disease risk and interindividual susceptibilities.[...].

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the “fifth” ribonucleotide in 1951. Since then, the ever‐increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal–Hreidarsson syndrome, Bowen–Conradi syndrome, or Williams–Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications.

This article is categorized under:

RNA Processing > RNA Editing and Modification

Loss of epitranscriptomic control of selenocysteine utilization engages senescence and mitochondrial reprogramming☆

Critically important to the maintenance of the glutathione (GSH) redox cycle are the activities of many selenocysteine-containing GSH metabolizing enzymes whose translation is controlled by the epitranscriptomic writer alkylation repair homolog 8 (ALKBH8). ALKBH8 is a tRNA methyltransferase that methylates the wobble uridine of specific tRNAs to regulate the synthesis of selenoproteins. Here we demonstrate that a deficiency in the writer ALKBH8 (Alkbh8^def), alters selenoprotein levels and engages senescence, regulates stress response genes and promotes mitochondrial reprogramming. Alkbh8^def mouse embryonic fibroblasts (MEFs) increase many hallmarks of senescence, including senescence associated β-galactosidase, heterocromatic foci, the cyclin dependent kinase inhibitor p16^Ink4a, markers of mitochondrial dynamics as well as the senescence associated secretory phenotype (SASP). Alkbh8^def cells also acquire a stress resistance phenotype that is accompanied by an increase in a number redox-modifying transcripts. In addition, Alkbh8^def MEFs undergo a metabolic shift that is highlighted by a striking increase in the level of uncoupling protein 2 (UCP2) which enhances oxygen consumption and promotes a reliance on glycolytic metabolism. Finally, we have shown that the Alkbh8 deficiency can be exploited and corresponding MEFs are killed by glycolytic inhibition. Our work demonstrates that defects in an epitransciptomic writer promote senescence and mitochondrial reprogramming and unveils a novel adaptive mechanism for coping with defects in selenocysteine utilization.

•

Deficiencies in selenocysteine utilization engages cellular senescence and the senescence associated secretory phenotype.

•

Alkbh8 deficiency promotes mitochondrial elongation, increased oxygen consumption and a reliance on glycolytic metabolism.

•

Cellular adaptions to Alkbh8 deficiency confer stress resistance.

Epitranscriptomic systems regulate the translation of reactive oxygen species detoxifying and disease linked selenoproteins

Here we highlight the role of epitranscriptomic systems in post-transcriptional regulation, with a specific focus on RNA modifying writers required for the incorporation of the 21st amino acid selenocysteine during translation, and the pathologies linked to epitranscriptomic and selenoprotein defects. Epitranscriptomic marks in the form of enzyme-catalyzed modifications to RNA have been shown to be important signals regulating translation, with defects linked to altered development, intellectual impairment, and cancer. Modifications to rRNA, mRNA and tRNA can affect their structure and function, while the levels of these dynamic tRNA-specific epitranscriptomic marks are stress-regulated to control translation. The tRNA for selenocysteine contains five distinct epitranscriptomic marks and the ALKBH8 writer for the wobble uridine (U) has been shown to be vital for the translation of the glutathione peroxidase (GPX) and thioredoxin reductase (TRXR) family of selenoproteins. The reactive oxygen species (ROS) detoxifying selenocysteine containing proteins are a prime examples of how specialized translation can be regulated by specific tRNA modifications working in conjunction with distinct codon usage patterns, RNA binding proteins and specific 3' untranslated region (UTR) signals. We highlight the important role of selenoproteins in detoxifying ROS and provide details on how epitranscriptomic marks and selenoproteins can play key roles in and maintaining mitochondrial function and preventing disease.

Structural Elucidation of Bisulfite Adducts to Pseudouridine That Result in Deletion Signatures during Reverse Transcription of RNA

The recent report of RBS-Seq to map simultaneously the epitranscriptomic modifications N¹-methyladenosine, 5-methylcytosine, and pseudouridine (Ψ) via bisulfite treatment of RNA provides a key advance to locate these important modifications. The locations of Ψ were found by a deletion signature generated during cDNA synthesis after bisulfite treatment for which the chemical details of the reaction are poorly understood. In the present work, the bisulfite reaction with Ψ was explored to identify six isomers of bisulfite adducted to Ψ. We found four of these adducts involved the heterocyclic ring, similar to the reaction with other pyrimidines. The remaining two adducts were bonded to the 1' carbon, which resulted in opening of the ribose ring. The utilization of complementary 1D- and 2D-NMR, Raman, and electronic circular dichroism spectroscopies led to the assignment of the two ribose adducts being the constitutional isomers of an S- and an O-adduct of bisulfite to the ribose, and these are the final products after heating. A mechanistic proposal is provided to rationalize chemically the formation and stereochemistries of all six isomeric bisulfite adducts to Ψ; conversion of intermediate adducts to the two final products is proposed to involve E2, S_N2', and [2,3]-sigmatropic shift reactions. Lastly, a synthetic RNA template with Ψ at a known location was treated with bisulfite, leading to a deletion signature after reverse transcription, supporting the RBS-Seq report. This classical bisulfite reaction used for epigenomic and epitranscriptomic sequencing diverges from the C nucleoside Ψ to form stable bisulfite end products that yield signatures for next-generation sequencing.