tRNA Abundance Research Articles

Mettl1-dependent m7G tRNA modification is essential for maintaining spermatogenesis and fertility in Drosophila melanogaster

Modification of guanosine to N7-methylguanosine (m7G) in the variable loop region of tRNA is catalyzed by the METTL1/WDR4 heterodimer and stabilizes target tRNA. Here, we reveal essential functions of Mettl1 in Drosophila fertility. Knockout of Mettl1 (Mettl1-KO) causes no major effect on the development of non-gonadal tissues, but abolishes the production of elongated spermatids and mature sperm, which is fully rescued by expression of a Mettl1-transgene, but not a catalytic-dead Mettl1 transgene. This demonstrates that Mettl1-dependent m7G is required for spermatogenesis. Mettl1-KO results in a loss of m7G modification on a subset of tRNAs and decreased tRNA abundance. Ribosome profiling shows that Mettl1-KO led to ribosomes stalling at codons decoded by tRNAs that were reduced in abundance. Mettl1-KO also significantly reduces the translation efficiency of genes involved in elongated spermatid formation and sperm stability. Germ cell-specific expression of Mettl1 rescues disrupted m7G tRNA modification and tRNA abundance in Mettl1-KO testes but not in non-gonadal tissues. Ribosome stalling is much less detectable in non-gonadal tissues than in Mettl1-KO testes. These findings reveal a developmental role for m7G tRNA modification and indicate that m7G modification-dependent tRNA abundance differs among tissues.

199 Awardee Talk: Insights into translational regulation through tRNA sequencing, RNA sequencing, and ribosome profiling

Abstract Translation acts as an additional layer of regulation that has an important role in gene expression and function. Due to a phenomenon called codon usage bias, highly expressed genes are thought to be codon-biased to support efficient translation. As more than one codon can code for the same amino acid (synonymous codons), organisms may exhibit preferences for specific codons that facilitate increased expression of important genes due to variations in the availability of corresponding tRNAs. Advances in sequencing technologies have recently permitted the study of the translatome, which refers to the entire population of mRNA associated with ribosomes for protein synthesis and can be investigated through ribosome profiling. This cutting-edge technique allows the identification of actively translated regions, but can also reveal translational pausing events that stem from the presence of SNPs. Our previous research has demonstrated variation in tRNA expression across tissues and different states of health, leading us to consider the connection between tRNA abundance and translational stalling due to SNPs. By coupling ribosome profiling with tRNA sequencing and RNA sequencing, we have investigated all elements of translational machinery to predict translational efficiency, estimate proteome composition, and evaluate tRNA abundance as a source of genetic variation. In this work, we utilized ribosome profiling, tRNAseq, and RNAseq and performed an integrative analysis in bovine tissues (kidney, liver and muscle) as well as murine myoblast cell lines (C2C12). By applying these methods to different bovine tissues, we were able to explore the interplay between tRNA availability and translational stalling events. Moreover, we have identified translationally regulated genes underlying tissue-specific biological processes and found that many upregulated and downregulated genes coincided with high and low translational efficiency respectively. We have also successfully defined stalling sites that depict the regulatory information encoded within the coding sequence of transcripts, which could control translation rate and facilitate proper protein folding. Through the implementation of these next generation sequencing (NGS) methods to cell culture, we were able to evaluate the translatome at various time points (0-min, 30-min, 60-min, and 4-h) post-induction of differentiation. Where most studies have focused on molecular changes between differentiating myoblasts and multinucleated myotubes that are present 7 d after differentiation, we focus on the earliest stages of muscle differentiation that initiate muscle development and therefore kickstart the process of meat production. This work offers an atlas of distinctive stalling sites across bovine tissues and various stages of muscle differentiation, which provides an opportunity to predict codon optimality and understand tissue-specific or time-specific mechanisms of regulating protein synthesis.

Metabolic rewiring during bone development underlies tRNA m7G-associated primordial dwarfism.

Translation of mRNA to protein is tightly regulated by transfer RNAs (tRNAs), which are subject to various chemical modifications that maintain structure, stability, and function. Deficiency of tRNA N7-methylguanosine (m7G) modification in patients causes a type of primordial dwarfism, but the underlying mechanism remains unknown. Here we report that the loss of m7G rewires cellular metabolism, leading to the pathogenesis of primordial dwarfism. Conditional deletion of the catalytic enzyme Mettl1 or missense mutation of the scaffold protein Wdr4 severely impaired endochondral bone formation and bone mass accrual. Mechanistically, Mettl1 knockout decreased abundance of m7G-modified tRNAs and inhibited translation of mRNAs relating to cytoskeleton and Rho GTPase signaling. Meanwhile, Mettl1 knockout enhanced cellular energy metabolism despite incompetent proliferation and osteogenic commitment. Further exploration revealed that impairment of Rho GTPase signaling upregulated the level of branched-chain amino acid transaminase 1 (BCAT1) that rewired cell metabolism and restricted intracellular α-ketoglutarate (αKG). Supplementation of αKG ameliorated the skeletal defect of Mettl1-deficient mice. In addition to the selective translation of metabolism-related mRNAs, we further revealed that Mettl1 knockout globally regulated translation via integrated stress response (ISR) and mammalian target of rapamycin complex 1 (mTORC1) signaling. Restoring translation by targeting either ISR or mTORC1 aggravated bone defects of Mettl1-deficient mice. Overall, our study unveils a critical role of m7G tRNA modification in bone development by regulation of cellular metabolism and indicates suspension of translation initiation as a quality control mechanism in response to tRNA dysregulation.

Impaired 2',3'-cyclic phosphate tRNA repair causes thermo-sensitive genic male sterility in rice.

Hybrid rice, widely planted in Asia, is pathogen resistant and has superior yields, making it a major contributor to global food security. The two-line hybrid rice system, which utilizes mutants exhibiting photo-/thermo-sensitive genic male sterility (P/TGMS), is the leading hybrid rice breeding technology. Mutations in THERMO-SENSITIVE GENIC MALE STERILE 5 (TMS5) accounts for over 95% of current TGMS lines. We previously found that tms5 carries a mutation in ribonuclease ZS1. Despite its importance for breeding robust rice lines, the mechanism underlying tms5-mediated TGMS remains elusive. Here, we demonstrate that TMS5 is a tRNA 2',3'-cyclic phosphatase. The tms5 mutation leads to accumulation of 2',3'-cyclic phosphate (cP)-ΔCCA-tRNAs (tRNAs without 3'CCA ended with cP), which is exacerbated by high temperatures, and reduction in the abundance of mature tRNAs, particularly alanine tRNAs (tRNA-Alas). Overexpression of tRNA-Alas in the tms5 mutant restores male fertility to 70%. Remarkably, male fertility of tms5 mutant is completely restored at high temperatures by knocking out OsVms1 which encodes the enzyme for cP-ΔCCA-tRNA generation. Our study reveals the mechanism underlying tms5-mediated TGMS in rice and provides mechanistic insight into the further improvement of TGMS in hybrid crop development.

Transfer RNA levels are tuned to support differentiation during Drosophila neurogenesis.

Neural differentiation requires a multifaceted program to alter gene expression along the proliferation to differentiation axis. While critical changes occur at the level of transcription, post-transcriptional mechanisms allow fine-tuning of protein output. We investigated the role of tRNAs in regulating gene expression during neural differentiation by quantifying tRNA abundance in neural progenitor-biased and neuron-biased Drosophila larval brains. We found that tRNA profiles are largely consistent between progenitor-biased and neuron-biased brains but significant variation occurs for 10 cytoplasmic isodecoders (individual tRNA genes) and this establishes differential tRNA levels for 8 anticodon groups. We used these tRNA data to investigate relationships between tRNA abundance, codon optimality-mediated mRNA decay, and translation efficiency in progenitors and neurons. Our data reveal that tRNA levels strongly correlate with codon optimality-mediated mRNA decay within each cell type but generally do not explain differences in stabilizing versus destabilizing codons between cell types. Regarding translation efficiency, we found that tRNA expression in neural progenitors preferentially supports translation of mRNAs whose products are in high demand in progenitors, such as those associated with protein synthesis. In neurons, tRNA expression shifts to disfavor translation of proliferation-related transcripts and preferentially support translation of transcripts tied to neuron-specific functions like axon pathfinding and synapse formation. Overall, our analyses reveal that changes in tRNA levels along the neural differentiation axis support optimal gene expression in progenitors and neurons.

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species.

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

Transcriptome and proteome changes triggered by overexpression of the transcriptional regulator Maf1 in the human pathogen Leishmania major.

Maf1, originally described as a repressor of RNA polymerase III (RNAP III) transcription in yeast, participates in multiple functions across eukaryotes. However, the knowledge about Maf1 in protozoan parasites is scarce. To initiate the study of Maf1 in Leishmania major, we generated a cell line that overexpresses this protein. Overexpression of Maf1 led to a significant reduction in the abundance of tRNAs, 5S rRNA, and U4 snRNA, demonstrating that Maf1 regulates RNAP III activity in L. major. To further explore the roles played by Maf1 in this microorganism, global transcriptomic and proteomic changes due to Maf1 overexpression were determined using RNA-sequencing and label-free quantitative mass spectrometry. Compared to wild-type cells, differential expression was observed for 1082 transcripts (615 down-regulated and 467 up-regulated) and 205 proteins (132 down-regulated and 73 up-regulated) in the overexpressing cells. A correlation of 44% was found between transcriptomic and proteomic results. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the differentially expressed genes and proteins are mainly involved in transcription, cell cycle regulation, lipid metabolism and transport, ribosomal biogenesis, carbohydrate metabolism, autophagy, and cytoskeleton modification. Thus, our results suggest the involvement of Maf1 in the regulation of all these processes in L. major, as reported in other species, indicating that the functions performed by Maf1 were established early in eukaryotic evolution. Notably, our data also suggest the participation of L. major Maf1 in mRNA post-transcriptional control, a role that, to the best of our knowledge, has not been described in other organisms.

Predicting the translation efficiency of messenger RNA in mammalian cells.

The degree to which translational control is specified by mRNA sequence is poorly understood in mammalian cells. Here, we constructed and leveraged a compendium of 3,819 ribosomal profiling datasets, distilling them into a transcriptome-wide atlas of translation efficiency (TE) measurements encompassing >140 human and mouse cell types. We subsequently developed RiboNN, a multitask deep convolutional neural network, and classic machine learning models to predict TEs in hundreds of cell types from sequence-encoded mRNA features, achieving state-of-the-art performance (r=0.79 in human and r=0.78 in mouse for mean TE across cell types). While the majority of earlier models solely considered 5' UTR sequence, RiboNN integrates contributions from the full-length mRNA sequence, learning that the 5' UTR, CDS, and 3' UTR respectively possess ~67%, 31%, and 2% per-nucleotide information density in the specification of mammalian TEs. Interpretation of RiboNN revealed that the spatial positioning of low-level di- and tri-nucleotide features (i.e., including codons) largely explain model performance, capturing mechanistic principles such as how ribosomal processivity and tRNA abundance control translational output. RiboNN is predictive of the translational behavior of base-modified therapeutic RNA, and can explain evolutionary selection pressures in human 5' UTRs. Finally, it detects a common language governing mRNA regulatory control and highlights the interconnectedness of mRNA translation, stability, and localization in mammalian organisms.

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 codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codon 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 discovered 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.

TRNA-m1A methylation controls the infection of Magnaporthe oryzae by supporting ergosterol biosynthesis.

Ergosterols are essential components of fungal plasma membranes. Inhibitors targeting ergosterol biosynthesis (ERG) genes are critical for controlling fungal pathogens, including Magnaporthe oryzae, the fungus that causes rice blast. However, the translational mechanisms governing ERG gene expression remain largely unexplored. Here, we show that the Trm6/Trm61 complex catalyzes dynamic N1-methyladenosine at position 58 (m1A58) in 51 transfer RNAs (tRNAs) of M.oryzae, significantly influencing translation at both the initiation and elongation stages. Notably, tRNA m1A58 promotes elongation speed at most cognate codons mainly by enhancing eEF1-tRNA binding rather than affecting tRNA abundance or charging. The absence of m1A58 leads to substantial decreases in the translation of ERG genes, ergosterol production, and, consequently, fungal virulence. Simultaneously targeting the Trm6/Trm61 complex and the ergosterol biosynthesis pathway markedly improves rice blast control. Our findings demonstrate an important role of m1A58-mediated translational regulation in ergosterol production and fungal infection, offering a potential strategy for fungicide development.

Chromosomal DNA sequences of the Pacific saury genome: versatile resources for fishery science and comparative biology.

Pacific saury (Cololabis saira) is a commercially important small pelagic fish species in Asia. In this study, we conducted the first-ever whole genome sequencing of this species, with single molecule, real-time (SMRT) sequencing technology. The obtained high-fidelity (HiFi) long-read sequence data, which amount to ~30-folds of its haploid genome size that was measured with quantitative PCR (1.17 Gb), were assembled into contigs. Scaffolding with Hi-C reads yielded a whole genome assembly containing 24 chromosome-scale sequences, with a scaffold N50 length of 47.7Mb. Screening of repetitive elements including telomeric repeats was performed to characterize possible factors that need to be resolved towards 'telomere-to-telomere' sequencing. The larger genome size than in medaka, a close relative in Beloniformes, is at least partly explained by larger repetitive element quantity, which is reflected in more abundant tRNAs, in the Pacific saury genome. Protein-coding regions were predicted using transcriptome data, which resulted in 22,274 components. Retrieval of Pacific saury homologs of aquaporin (AQP) genes known from other teleost fishes validated high completeness and continuity of the genome assembly. These resources are available at https://treethinkers.nig.ac.jp/saira/ and will assist various molecular-level studies in fishery science and comparative biology.

Near-cognate tRNAs increase the efficiency and precision of pseudouridine-mediated readthrough of premature termination codons.

Programmable RNA pseudouridylation has emerged as a new type of RNA base editor to suppress premature termination codons (PTCs) that can lead to truncated and nonfunctional proteins. However, current methods to correct disease-associated PTCs suffer from low efficiency and limited precision. Here we develop RESTART v3, which uses near-cognate tRNAs to improve the readthrough efficiency of pseudouridine-modified PTCs. We show an average of ~5-fold (range: 2.1- to 9.5-fold) higher editing efficiency than RESTART v2 in cultured cells and achieve functional PTC readthrough in disease cell models of cystic fibrosis and Hurler syndrome. Furthermore, RESTART v3 enables accurate incorporation of the original amino acid for nearly half of the PTC sites, considering the naturally occurring frequencies of sense-to-nonsense codons, without affecting normal termination codons. Although off-target sites were detected, we did not observe changes to the coding information or the expression level of transcripts, and the overall natural tRNA abundance remained constant.

Selection on synonymous sites: the unwanted transcript hypothesis.

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Implications of tRNA abundance on translation elongation across bovine tissues.

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene

BackgroundSingle synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed.ResultsSynonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo.ConclusionIn summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

99 Noncoding Rnas Alter Our Interpretations of Genome to Phenome

Abstract Gene transcription and protein translation are core components in the process of gene expression. Previous research on the regulation of gene expression largely focuses on events prior to translation, including epigenetic regulation, transcription, and RNA processing. However, translation acts as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are thought to be codon-biased to support efficient translation, in which the encoded codons correspond to highly abundant tRNAs. Further, synonymous SNPs were once considered to be silent due to the degeneracy of the genetic code. However, synonymous SNPs may disturb protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Our previous study identified variation in tRNA expression within bovine liver and muscle tissue, suggesting tissue-specific modulation of translation. In addition, advances in sequencing technologies have recently permitted the study of the translatome, which refers to the entire population of mRNA associated with ribosomes for protein synthesis and can be investigated through ribosome profiling. Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNA-seq) and ribosome profiling to investigate the relationship between tRNA expression and translational stalling events. Moreover, we have identified translationally regulated genes underlying tissue-specific biological processes and found that many upregulated and downregulated genes coincided with high and low translational efficiency respectively. We have also successfully defined stalling sites that depict the regulatory information encoded within the coding sequence of transcripts, which could control translation rate and facilitate proper protein folding. This work offers an atlas of distinctive stalling sites across bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation.

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Analyzing the Evolution and Host Adaptation of the Rabies Virus from the Perspective of Codon Usage Bias

Rabies virus (RABV) is a highly pathogenic virus that causes a fatal disease in humans and other mammals, but the mechanism of its evolution, spread, and spillover remains unknown. In this study, we analyzed the codon usage pattern of 2,018 RABV full-length genome sequences from 79 countries collected between 1931 and 2021 to provide an insight into its molecular evolution and unravel its unknown host-adapted pattern. We found that RABV exhibited a weak codon usage bias, with a preference for the codons ending in A (28.10 ± 0.01) or U (26.43 ± 0.02). Moreover, natural selection plays a major role in shaping the codon usage bias of the RABV. Notably, nearly half of the 18 codons in the virus were best matched to the hosts’ most abundant isoacceptor tRNAs, which might account for the wide range of RABV hosts. Furthermore, significant differences were observed in the codon usage patterns of RABV for different host species, suggesting that codon usage bias may be influenced by host-specific factors. In conclusion, our study reveals codon usage patterns of RABV that may help in the development of control strategies and effective vaccines and therapies against this deadly virus.

Mitochondrial Metabolism in the Spotlight: Maintaining Balanced RNAP III Activity Ensures Cellular Homeostasis.

RNA polymerase III (RNAP III) holoenzyme activity and the processing of its products have been linked to several metabolic dysfunctions in lower and higher eukaryotes. Alterations in the activity of RNAP III-driven synthesis of non-coding RNA cause extensive changes in glucose metabolism. Increased RNAP III activity in the S. cerevisiae maf1Δ strain is lethal when grown on a non-fermentable carbon source. This lethal phenotype is suppressed by reducing tRNA synthesis. Neither the cause of the lack of growth nor the underlying molecular mechanism have been deciphered, and this area has been awaiting scientific explanation for a decade. Our previous proteomics data suggested mitochondrial dysfunction in the strain. Using model mutant strains maf1Δ (with increased tRNA abundance) and rpc128-1007 (with reduced tRNA abundance), we collected data showing major changes in the TCA cycle metabolism of the mutants that explain the phenotypic observations. Based on 13C flux data and analysis of TCA enzyme activities, the present study identifies the flux constraints in the mitochondrial metabolic network. The lack of growth is associated with a decrease in TCA cycle activity and downregulation of the flux towards glutamate, aspartate and phosphoenolpyruvate (PEP), the metabolic intermediate feeding the gluconeogenic pathway. rpc128-1007, the strain that is unable to increase tRNA synthesis due to a mutation in the C128 subunit, has increased TCA cycle activity under non-fermentable conditions. To summarize, cells with non-optimal activity of RNAP III undergo substantial adaptation to a new metabolic state, which makes them vulnerable under specific growth conditions. Our results strongly suggest that balanced, non-coding RNA synthesis that is coupled to glucose signaling is a fundamental requirement to sustain a cell's intracellular homeostasis and flexibility under changing growth conditions. The presented results provide insight into the possible role of RNAP III in the mitochondrial metabolism of other cell types.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium.

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.