Transfer RNA genes experience exceptionally elevated mutation rates

How much is Transcription-associated Mutagenesis Driving tRNA Microevolution?

Transfer RNAs (tRNAs) are among the most highly conserved and frequently transcribed genes. Recent studies have identified transcription-associated mutagenesis (TAM) as a significant contributor to sequence variation around tRNA loci. However, the extent to which TAM drives allelic variation in tRNAs remains unclear, largely due to the confounding effects of strong selection pressures to maintain their structural integrity. This complexity arises because TAM-induced mutations primarily involve nucleotide transitions, which tend to preserve base-pairing stability. To address this dichotomy at the population level, we analyzed tRNA allelic variation in contemporary Caenorhabditis elegans strains. We propose a model of tRNA microevolution driven by TAM and demonstrate that the observed secondary structure characteristics align with our predicted TAM-biased patterns. Furthermore, we developed a continuous Markov substitution model that incorporates TAM-specific mutational biases. This TAM-biased model fits the C. elegans tRNA data more effectively than standard models, such as the general time-reversible (GTR) model. Based on these results, we conclude that TAM plays a significant role in shaping tRNA allelic variation within populations. This finding is consistent with recent experimental studies on tRNA fitness in yeast but challenges prior theoretical and computational analyses that emphasize RNA base-pairing as a primary determinant in genotype-phenotype systems.

Glutamine missense suppressor transfer RNAs inhibit polyglutamine aggregation

Huntington's disease (HD) is caused by polyglutamine (polyQ) repeat expansions in the huntingtin gene. HD-causative polyQ alleles lead to protein aggregation, which is a prerequisite for disease. Translation fidelity modifies protein aggregation, and several studies suggest that mutating one or two glutamine (Gln) residues in polyQ reduces aggregation. Thus, we hypothesized that missense suppression of Gln codons with other amino acids will reduce polyQ aggregate formation in cells. In neuroblastoma cells, we assessed tRNA variants that misread Gln codons with serine (tRNASer C/UUG) or alanine (tRNAAla C/UUG). The tRNAs with the CUG anticodon were more effective at suppressing the CAG repeats in polyQ, and serine and alanine mis-incorporation had differential impacts on polyQ. The expression of tRNASer CUG reduced polyQ protein production as well as both soluble and insoluble aggregate formation. In contrast, cells expressing tRNAAla CUG selectively decreased insoluble polyQ aggregate formation by 2-fold. Mass spectrometry confirmed Ala mis-incorporation at an average level of ∼20% per Gln codon. Cells expressing the missense suppressor tRNAs showed no cytotoxic effects and no defects in growth or global protein synthesis levels. Our findings demonstrate that tRNA-dependent missense suppression of Gln codons is well tolerated in mammalian cells and significantly reduces polyQ levels and aggregates that cause HD.

Isolation and characterization of Salmonella enterica- and Escherichia coli-specific bacteriophages of the genus Epseptimavirus from wastewater in Minnesota

Five lytic bacteriophages specific for Salmonella enterica and Escherichia coli were isolated from wastewater in Minnesota. These phages, designated vB_Sal_EH1, vB_Sal_EH2, vB_Sal_EH3, vB_Sal_EH4, and vB_Sal_EH7, were characterized, and their genomes were sequenced. Phylogenetic analysis showed that they grouped within the genus Epseptimavirus, with genome sizes ranging from 108,554 to 115,218 bp. All five phages exhibited lytic activity against both S. enterica and Shiga-toxin-producing E. coli O157:H7. Transposon mutagenesis of the host genome identified the outer membrane protein BtuB as essential for phage infection, suggesting that it is a putative receptor. Genome sequence comparisons revealed genetic loci that are variable among the isolated phages and potentially influence their host specificity and virulence.

Towards the genomic sequence code of DNA fragility for machine learning

Genomic DNA breakages and the subsequent insertion and deletion mutations are important contributors to genome instability and linked diseases. Unlike the research in point mutations, the relationship between DNA sequence context and the propensity for strand breaks remains elusive. Here, by analyzing the differences and commonalities across myriads of genomic breakage datasets, we extract the sequence-linked rules and patterns behind DNA fragility. We show the overall deconvolution of the sequence influence into short-, mid- and long-range effects, and the stressor-dependent differences in defining the range and compositional effects on DNA fragility. We summarize and release our feature compendium as a library that can be seamlessly incorporated into genomic machine learning procedures, where DNA fragility is of concern, and train a generalized DNA fragility model on cancer-associated breakages. Structural variants (SVs) tend to stabilize regions in which they emerge, with the effect most pronounced for pathogenic SVs. In contrast, the effects of chromothripsis are seen across regions less prone to breakages. We find that viral integration may bring genome fragility, particularly for cancer-associated viruses. Overall, this work offers novel insights into the genomic sequence basis of DNA fragility and presents a powerful machine learning resource to further enhance our understanding of genome (in)stability and evolution.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

De novo variants in the RNU4-2 snRNA cause a frequent neurodevelopmental syndrome

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here we identify the non-coding RNA RNU4-2 as a syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 base pair region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 115 individuals with NDD. Most individuals (77.4%) have the same highly recurrent single base insertion (n.64_65insT). In 54 individuals in whom it could be determined, the de novo variants were all on the maternal allele. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to RNU4-1 and other U4 homologues. Using RNA sequencing, we show how 5' splice-site use is systematically disrupted in individuals with RNU4-2 variants, consistent with the known role of this region during spliceosome activation. Finally, we estimate that variants in this 18 base pair region explain 0.4% of individuals with NDD. This work underscores the importance of non-coding genes in rare disorders and will provide a diagnosis to thousands of individuals with NDD worldwide.

The 'Not-So-Famous Five' in tumorigenesis: tRNAs, tRNA fragments, and tRNA epitranscriptome in concert with AARSs and AIMPs

RNA profiling studies have revealed that ∼75% of the human genome is transcribed to RNA but only a meagre fraction of it is translated to proteins. Majority of transcribed RNA constitute a specialized pool of non-coding RNAs. Human genome contains approximately 506 genes encoding a set of 51 different tRNAs, constituting a unique class of non-coding RNAs that not only have essential housekeeping functions as translator molecules during protein synthesis, but have numerous uncharted regulatory functions. Intriguing findings regarding a variety of non-canonical functions of tRNAs, tRNA derived fragments (tRFs), esoteric epitranscriptomic modifications of tRNAs, along with aminoacyl-tRNA synthetases (AARSs) and ARS-interacting multifunctional proteins (AIMPs), envision a 'peripheral dogma' controlling the flow of genetic information in the backdrop of qualitative information wrung out of the long-live central dogma of molecular biology, to drive cells towards either proliferation or differentiation programs. Our review will substantiate intriguing peculiarities of tRNA gene clusters, atypical tRNA-transcription from internal promoters catalysed by another distinct RNA polymerase enzyme, dynamically diverse tRNA epitranscriptome, intricate mechanism of tRNA-charging by AARSs governing translation fidelity, epigenetic regulation of gene expression by tRNA fragments, and the role of tRNAs and tRNA derived/associated molecules as quantitative determinants of the functional proteome, covertly orchestrating the process of tumorigenesis, through a deregulated tRNA-ome mediating selective codon-biased translation of cancer related gene transcripts.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

A mutation rate model at the basepair resolution identifies the mutagenic effect of polymerase III transcription

De novo mutations occur at substantially different rates depending on genomic location, sequence context and DNA strand. The success of methods to estimate selection intensity, infer demographic history and map rare disease genes, depends strongly on assumptions about the local mutation rate. Here we present Roulette, a genome-wide mutation rate model at basepair resolution that incorporates known determinants of local mutation rate. Roulette is shown to be more accurate than existing models. We use Roulette to refine the estimates of population growth within Europe by incorporating the full range of human mutation rates. The analysis of significant deviations from the model predictions revealed a tenfold increase in mutation rate in nearly all genes transcribed by polymerase III (Pol III), suggesting a new mutagenic mechanism. We also detected an elevated mutation rate within transcription factor binding sites restricted to sites actively used in testis and residing in promoters.

Design, construction, and functional characterization of a tRNA neochromosome in yeast

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Transfer RNAs as dynamic and critical regulators of cancer progression

Transfer RNAs (tRNAs) have been historically viewed as non-dynamic adaptors that decode the genetic code into proteins. Recent work has uncovered dynamic regulatory roles for these fascinating molecules. Advances in tRNA detection methods have revealed that specific tRNAs can become modulated upon DNA copy number and chromatin alterations and can also be perturbed by oncogenic signalling and transcriptional regulators in cancer cells or the tumour microenvironment. Such alterations in the levels of specific tRNAs have been shown to causally impact cancer progression, including metastasis. Moreover, sequencing methods have identified tRNA-derived small RNAs that influence various aspects of cancer progression, such as cell proliferation and invasion, and could serve as diagnostic and prognostic biomarkers or putative therapeutic targets in various cancers. Finally, there is accumulating evidence, including from genetic models, that specific tRNA synthetases - the enzymes responsible for charging tRNAs with amino acids - can either promote or suppress tumour formation. In this Review, we provide an overview of how deregulation of tRNAs influences cancer formation and progression.

The divergence of mutation rates and spectra across the Tree of Life

Owing to advances in genome sequencing, genome stability has become one of the most scrutinized cellular traits across the Tree of Life. Despite its centrality to all things biological, the mutation rate (per nucleotide site per generation) ranges over three orders of magnitude among species and several-fold within individual phylogenetic lineages. Within all major organismal groups, mutation rates scale negatively with the effective population size of a species and with the amount of functional DNA in the genome. This relationship is most parsimoniously explained by the drift-barrier hypothesis, which postulates that natural selection typically operates to reduce mutation rates until further improvement is thwarted by the power of random genetic drift. Despite this constraint, the molecular mechanisms underlying DNA replication fidelity and repair are free to wander, provided the performance of the entire system is maintained at the prevailing level. The evolutionary flexibility of the mutation rate bears on the resolution of several prior conundrums in phylogenetic and population-genetic analysis and raises challenges for future applications in these areas.

Copy number variation in tRNA isodecoder genes impairs mammalian development and balanced translation

The number of tRNA isodecoders has increased dramatically in mammals, but the specific molecular and physiological reasons for this expansion remain elusive. To address this fundamental question we used CRISPR editing to knockout the seven-membered phenylalanine tRNA gene family in mice, both individually and combinatorially. Using ATAC-Seq, RNA-seq, ribo-profiling and proteomics we observed distinct molecular consequences of single tRNA deletions. We show that tRNA-Phe-1-1 is required for neuronal function and its loss is partially compensated by increased expression of other tRNAs but results in mistranslation. In contrast, the other tRNA-Phe isodecoder genes buffer the loss of each of the remaining six tRNA-Phe genes. In the tRNA-Phe gene family, the expression of at least six tRNA-Phe alleles is required for embryonic viability and tRNA-Phe-1-1 is most important for development and survival. Our results reveal that the multi-copy configuration of tRNA genes is required to buffer translation and viability in mammals.

tRNA dysregulation and disease

tRNAs are key adaptor molecules that decipher the genetic code during translation of mRNAs in protein synthesis. In contrast to the traditional view of tRNAs as ubiquitously expressed housekeeping molecules, awareness is now growing that tRNA-encoding genes display tissue-specific and cell type-specific patterns of expression, and that tRNA gene expression and function are both dynamically regulated by post-transcriptional RNA modifications. Moreover, dysregulation of tRNAs, mediated by alterations in either their abundance or function, can have deleterious consequences that contribute to several distinct human diseases, including neurological disorders and cancer. Accumulating evidence shows that reprogramming of mRNA translation through altered tRNA activity can drive pathological processes in a codon-dependent manner. This Review considers the emerging evidence in support of the precise control of functional tRNA levels as an important regulatory mechanism that coordinates mRNA translation and protein expression in physiological cell homeostasis, and highlights key examples of human diseases that are linked directly to tRNA dysregulation.

The Evolutionary Dynamics of the Mitochondrial tRNA in the Cichlid Fish Family

Simple Summary

Cichlids are a unique example of fish diversity and species richness which have been explained by sympatric speciation at different freshwater sources in Africa. The mitochondria contribute to cell vitality by providing energy. It contains a circular genome with an established translation system that is spatially independent of the cytosolic counterpart. The current study aimed to investigate the evolutionary dynamics of the mitochondrial tRNA and its role in cichlids’ diversity. The available cichlid mitogenomes in the public database were filtered, in addition to newly sequenced accessions from a specific cichlid group known as the haplotilapiine lineage that is widely distributed in the Egyptian sector of the Nile River. Based on the comparative analysis of mitogenomic data, we identified 22 tRNA genes, in which a single gene was D-armless, while the cloverleaf secondary structure subdivided into stem-loop formations was predicted and used to define the levels of genetic divergence for the remained tRNAs. Peculiarly, in cichlids, the formation known as “T-arm” showed the lowest polymorphism levels among other structures in contrast to other organisms (e.g., scorpions). Comparing the whole family to the specific haplotilapiine lineage showed that the tryptophan tRNA was the most conserved tRNA, with signatures of possible purifying selection.

Abstract

The mitochondrial transfer RNA genes (tRNAs) attract more attention due to their highly dynamic and rapidly evolving nature. The current study aimed to detect and evaluate the dynamics, characteristic patterns, and variations of mitochondrial tRNAs. The study was conducted in two main parts: first, the published mitogenomic sequences of cichlids mt tRNAs have been filtered. Second, the filtered mitochondrial tRNA and additional new mitogenomes representing the most prevalent Egyptian tilapiine were compared and analyzed. Our results revealed that all 22 tRNAs of cichlids folded into a classical cloverleaf secondary structure with four domains, except for trnS^GCU, missing the D domain in all cichlids. When consensus tRNAs were compared, most of the mutations were observed in the trnP at nucleotide levels (substitutions and indels), in contrast to trnL^UAA. From a structural perspective, the anticodon loop and T-loop formations were the most conserved structures among all parts of the tRNA in contrast to the A-stem and D-loop formations. The trnW was the lowest polymorphic unneutral tRNA among all cichlids (both the family and the haplotilapiine lineage), in contrast with the neutral trnD that was extremely polymorphic among and within the haplotilapiine lineage species compared to other cichlids species. From a phylogenetic perspective, the trnC was extremely hypervariable and neutral tRNA in both haplotilapiine lineage and cichlids but was unable to report correct phylogenetic signal for the cichlids. In contrast to trnI and trnY, less variable neutral tRNAs that were able to cluster the haplotilapiine lineage and cichlids species as previously reported. By observing the DNA polymorphism in the coding DNA sequences (CDS), the highest affected amino acid by non-synonymous mutations was isoleucine and was equally mutated to valine and vice versa; no correlation between mutations in CDS and tRNAs was statistically found. The current study provides an insight into the mitochondrial tRNA evolution and its effect on the cichlid diversity and speciation model at the maternal level.

tRNA derived small RNAs—Small players with big roles

In the past 2 decades, small non-coding RNAs derived from tRNA (tsRNAs or tRNA derived fragments; tRFs) have emerged as new powerful players in the field of small RNA mediated regulation of gene expression, translation, and epigenetic control. tRFs have been identified from evolutionarily divergent organisms from Archaea, the higher plants, to humans. Recent studies have confirmed their roles in cancers and other metabolic disorders in humans and experimental models. They have been implicated in biotic and abiotic stress responses in plants as well. In this review, we summarize the current knowledge on tRFs including types of tRFs, their biogenesis, and mechanisms of action. The review also highlights recent studies involving differential expression profiling of tRFs and elucidation of specific functions of individual tRFs from various species. We also discuss potential considerations while designing experiments involving tRFs identification and characterization and list the available bioinformatics tools for this purpose.

Fitness landscape analysis of a tRNA gene reveals that the wild type allele is sub-optimal, yet mutationally robust

Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 8–10% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust (“flat”), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Epigenetic regulation of human non-coding RNA gene transcription

Recent investigations on the non-protein-coding transcriptome of human cells have revealed previously hidden layers of gene regulation relying on regulatory non-protein-coding (nc) RNAs, including the widespread ncRNA-dependent regulation of epigenetic chromatin states and of mRNA translation and stability. However, despite its centrality, the epigenetic regulation of ncRNA genes has received relatively little attention. In this mini-review, we attempt to provide a synthetic account of recent literature suggesting an unexpected complexity in chromatin-dependent regulation of ncRNA gene transcription by the three human nuclear RNA polymerases. Emerging common features, like the heterogeneity of chromatin states within ncRNA multigene families and their influence on 3D genome organization, point to unexplored issues whose investigation could lead to a better understanding of the whole human epigenomic network.

Comprehensive determination of transcription start sites derived from all RNA polymerases using ReCappable-seq

Determination of eukaryotic transcription start sites (TSSs) has been based on methods that require the cap structure at the 5' end of transcripts derived from Pol II RNA polymerase. Consequently, these methods do not reveal TSSs derived from the other RNA polymerases that also play critical roles in various cell functions. To address this limitation, we developed ReCappable-seq, which comprehensively identifies TSS for both Pol II and non-Pol II transcripts at single-nucleotide resolution. The method relies on specific enzymatic exchange of 5' m⁷G caps and 5' triphosphates with a selectable tag. When applied to human transcriptomes, ReCappable-seq identifies Pol II TSSs that are in agreement with orthogonal methods such as CAGE. Additionally, ReCappable-seq reveals a rich landscape of TSSs associated with Pol III transcripts that have not previously been amenable to study at genome-wide scale. Novel TSS from non-Pol II transcription can be located in the nuclear and mitochondrial genomes. ReCappable-seq interrogates the regulatory landscape of coding and noncoding RNA concurrently and enables the classification of epigenetic profiles associated with Pol II and non-Pol II TSS.

Intragenic tRNA-promoted R-loops orchestrate transcription interference for plant oxidative stress responses

Eukaryotic genomes are transcribed by at least three RNA polymerases, RNAPI, II, and III. Co-transcriptional R-loops play diverse roles in genome regulation and maintenance. However, little is known about how R-loops regulate transcription interference, the transcriptional event that is caused by different RNA polymerases transcribing the same genomic templates. Here, we established that the intragenic transfer RNA (tRNA) genes can promote sense R-loop enrichment (named intra-tR-loops) in Arabidopsis thaliana, and found that intra-tR-loops are decreased in an RNAPIII mutant, NUCLEAR RNA POLYMERASE C, SUBUNIT 7(nrpc7-1). NRPC7 is co-localized with RNAPIIS2P at intragenic tRNA genes and interferes with RNAPIIS2P elongation. Conversely, the binding of NRPC7 at intragenic tRNA genes is increased following inhibition of RNAPII elongation. The transcription of specific tRNA host genes is inhibited by RNAPIII, and the inhibition of tRNA host genes is intra-tR-loop dependent. Moreover, alleviating the inhibition of tRNAPro-induced intra-tR-loops on its host gene AtNUDX1 promotes oxidative stress tolerance in A. thaliana. Our work suggests intra-tR-loops regulate host gene expression by modulating RNA polymerases interference.

Comparative Genomics of Potato Common Scab-Causing Streptomyces spp. Displaying Varying Virulence

Common scab of potato causes important economic losses worldwide following the development of necrotic lesions on tubers. In this study, the genomes of 14 prevalent scab-causing Streptomyces spp. isolated from Prince Edward Island, one of the most important Canadian potato production areas, were sequenced and annotated. Their phylogenomic affiliation was determined, their pan-genome was characterized, and pathogenic determinants involved in their virulence, ranging from weak to aggressive, were compared. 13 out of 14 strains clustered with Streptomyces scabiei, while the last strain clustered with Streptomyces acidiscabies. The toxicogenic and colonization genomic regions were compared, and while some atypical gene organizations were observed, no clear correlation with virulence was observed. The production of the phytotoxin thaxtomin A was also quantified and again, contrary to previous reports in the literature, no clear correlation was found between the amount of thaxtomin A secreted, and the virulence observed. Although no significant differences were observed when comparing the presence/absence of the main virulence factors among the strains of S. scabiei, a distinct profile was observed for S. acidiscabies. Several mutations predicted to affect the functionality of some virulence factors were identified, including one in the bldA gene that correlates with the absence of thaxtomin A production despite the presence of the corresponding biosynthetic gene cluster in S. scabiei LBUM 1485. These novel findings obtained using a large number of scab-causing Streptomyces strains are challenging some assumptions made so far on Streptomyces’ virulence and suggest that other factors, yet to be characterized, are also key contributors.

The genomic loci of specific human tRNA genes exhibit ageing-related DNA hypermethylation

The epigenome has been shown to deteriorate with age, potentially impacting on ageing-related disease. tRNA, while arising from only ˜46 kb (<0.002% genome), is the second most abundant cellular transcript. tRNAs also control metabolic processes known to affect ageing, through core translational and additional regulatory roles. Here, we interrogate the DNA methylation state of the genomic loci of human tRNA. We identify a genomic enrichment for age-related DNA hypermethylation at tRNA loci. Analysis in 4,350 MeDIP-seq peripheral-blood DNA methylomes (16-82 years), identifies 44 and 21 hypermethylating specific tRNAs at study-and genome-wide significance, respectively, contrasting with none hypomethylating. Validation and replication (450k array and independent targeted Bisuphite-sequencing) supported the hypermethylation of this functional unit. Tissue-specificity is a significant driver, although the strongest consistent signals, also independent of major cell-type change, occur in tRNA-iMet-CAT-1-4 and tRNA-Ser-AGA-2-6. This study presents a comprehensive evaluation of the genomic DNA methylation state of human tRNA genes and reveals a discreet hypermethylation with advancing age.

Mfd regulates RNA polymerase association with hard-to-transcribe regions in vivo, especially those with structured RNAs

RNA polymerase (RNAP) encounters various roadblocks during transcription. These obstacles can impede RNAP movement and influence transcription, ultimately necessitating the activity of RNAP-associated factors. One such factor is the bacterial protein Mfd, a highly conserved DNA translocase and evolvability factor that interacts with RNAP. Although Mfd is thought to function primarily in the repair of DNA lesions that stall RNAP, increasing evidence suggests that it may also be important for transcription regulation. However, this is yet to be fully characterized. To shed light on Mfd's in vivo functions, we identified the chromosomal regions where it associates. We analyzed Mfd's impact on RNAP association and transcription regulation genome-wide. We found that Mfd represses RNAP association at many chromosomal regions. We found that these regions show increased RNAP pausing, suggesting that they are hard to transcribe. Interestingly, we noticed that the majority of the regions where Mfd regulates transcription contain highly structured regulatory RNAs. The RNAs identified regulate a myriad of biological processes, ranging from metabolism to transfer RNA regulation to toxin-antitoxin (TA) functions. We found that cells lacking Mfd are highly sensitive to toxin overexpression. Finally, we found that Mfd promotes mutagenesis in at least one toxin gene, suggesting that its function in regulating transcription may promote evolution of certain TA systems and other regions containing strong RNA secondary structures. We conclude that Mfd is an RNAP cofactor that is important, and at times critical, for transcription regulation at hard-to-transcribe regions, especially those that express structured regulatory RNAs.

Whole-Genome Transformation Promotes tRNA Anticodon Suppressor Mutations under Stress

tRNAs are encoded by a large gene family, usually with several isogenic tRNAs interacting with the same codon. Mutations in the anticodon region of other tRNAs can overcome specific tRNA deficiencies. Phylogenetic analysis suggests that such mutations have occurred in evolution, but the driving force is unclear. We show that in yeast suppressor mutations in other tRNAs are able to overcome deficiency of the essential TRT2-encoded tRNAThrCGU at high temperature (40°C). Surprisingly, these tRNA suppressor mutations were obtained after whole-genome transformation with DNA from thermotolerant Kluyveromyces marxianus or Ogataea polymorpha strains but from which the mutations did apparently not originate. We suggest that transient presence of donor DNA in the host facilitates proliferation at high temperature and thus increases the chances for occurrence of spontaneous mutations suppressing defective growth at high temperature. Whole-genome sequence analysis of three transformants revealed only four to five nonsynonymous mutations of which one causing TRT2 anticodon stem stabilization and two anticodon mutations in non-threonyl-tRNAs, tRNALysCUU and tRNAeMetCAU, were causative. Both anticodon mutations suppressed lethality of TRT2 deletion and apparently caused the respective tRNAs to become novel substrates for threonyl-tRNA synthetase. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) data could not detect any significant mistranslation, and reverse transcription-quantitative PCR results contradicted induction of the unfolded protein response. We suggest that stress conditions have been a driving force in evolution for the selection of anticodon-switching mutations in tRNAs as revealed by phylogenetic analysis.IMPORTANCE In this work, we have identified for the first time the causative elements in a eukaryotic organism introduced by applying whole-genome transformation and responsible for the selectable trait of interest, i.e., high temperature tolerance. Surprisingly, the whole-genome transformants contained just a few single nucleotide polymorphisms (SNPs), which were unrelated to the sequence of the donor DNA. In each of three independent transformants, we have identified a SNP in a tRNA, either stabilizing the essential tRNAThrCGU at high temperature or switching the anticodon of tRNALysCUU or tRNAeMetCAU into CGU, which is apparently enough for in vivo recognition by threonyl-tRNA synthetase. LC-MS/MS analysis indeed indicated absence of significant mistranslation. Phylogenetic analysis showed that similar mutations have occurred throughout evolution and we suggest that stress conditions may have been a driving force for their selection. The low number of SNPs introduced by whole-genome transformation may favor its application for improvement of industrial yeast strains.

Manipulation of the human tRNA pool reveals distinct tRNA sets that act in cellular proliferation or cell cycle arrest

Different subsets of the tRNA pool in human cells are expressed in different cellular conditions. The 'proliferation-tRNAs' are induced upon normal and cancerous cell division, while the 'differentiation-tRNAs' are active in non-dividing, differentiated cells. Here we examine the essentiality of the various tRNAs upon cellular growth and arrest. We established a CRISPR-based editing procedure with sgRNAs that each target a tRNA family. We measured tRNA essentiality for cellular growth and found that most proliferation-tRNAs are essential compared to differentiation- tRNAs in rapidly growing cell lines. Yet in more slowly dividing lines, the differentiation-tRNAs were more essential. In addition, we measured the essentiality of each tRNA family upon response to cell cycle arresting signals. Here we detected a more complex behavior with both proliferation-tRNAs and differentiation tRNAs showing various levels of essentiality. These results provide the so-far most comprehensive functional characterization of human tRNAs with intricate roles in various cellular states.

Noncanonical Roles of tRNAs: tRNA Fragments and Beyond

As one of the most abundant and conserved RNA species, transfer RNAs (tRNAs) are well known for their role in reading the codons on messenger RNAs and translating them into proteins. In this review, we discuss the noncanonical functions of tRNAs. These include tRNAs as precursors to novel small RNA molecules derived from tRNAs, also called tRNA-derived fragments, that are abundant across species and have diverse functions in different biological processes, including regulating protein translation, Argonaute-dependent gene silencing, and more. Furthermore, the role of tRNAs in biosynthesis and other regulatory pathways, including nutrient sensing, splicing, transcription, retroelement regulation, immune response, and apoptosis, is reviewed. Genome organization and sequence variation of tRNA genes are also discussed in light of their noncanonical functions. Lastly, we discuss the recent applications of tRNAs in genome editing and microbiome sequencing.

Genetic hitchhiking, mitonuclear coadaptation, and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

Expression of the Neuronal tRNA n-Tr20 Regulates Synaptic Transmission and Seizure Susceptibility

The mammalian genome has hundreds of nuclear-encoded tRNAs, but the contribution of individual tRNA genes to cellular and organismal function remains unknown. Here, we demonstrate that mutations in a neuronally enriched arginine tRNA, n-Tr20, increased seizure threshold and altered synaptic transmission. n-Tr20 expression also modulated seizures caused by an epilepsy-linked mutation in Gabrg2, a gene encoding a GABA_A receptor subunit. Loss of n-Tr20 altered translation initiation by activating the integrated stress response and suppressing mTOR signaling, the latter of which may contribute to altered neurotransmission in mutant mice. Deletion of a highly expressed isoleucine tRNA similarly altered these signaling pathways in the brain, suggesting that regulation of translation initiation is a conserved response to tRNA loss. Our data indicate that loss of a single member of a tRNA family results in multiple cellular phenotypes, highlighting the disease-causing potential of tRNA mutations.

Predicting transfer RNA gene activity from sequence and genome context

Transfer RNA (tRNA) genes are among the most highly transcribed genes in the genome owing to their central role in protein synthesis. However, there is evidence for a broad range of gene expression across tRNA loci. This complexity, combined with difficulty in measuring transcript abundance and high sequence identity across transcripts, has severely limited our collective understanding of tRNA gene expression regulation and evolution. We establish sequence-based correlates to tRNA gene expression and develop a tRNA gene classification method that does not require, but benefits from, comparative genomic information and achieves accuracy comparable to molecular assays. We observe that guanine + cytosine (G + C) content and CpG density surrounding tRNA loci is exceptionally well correlated with tRNA gene activity, supporting a prominent regulatory role of the local genomic context in combination with internal sequence features. We use our tRNA gene activity predictions in conjunction with a comprehensive tRNA gene ortholog set spanning 29 placental mammals to estimate the evolutionary rate of functional changes among orthologs. Our method adds a new dimension to large-scale tRNA functional prediction and will help prioritize characterization of functional tRNA variants. Its simplicity and robustness should enable development of similar approaches for other clades, as well as exploration of functional diversification of members of large gene families.

Enjoy the Silence: Nearly Half of Human tRNA Genes Are Silent

Transfer RNAs (tRNAs) are key components of the translation machinery. They read codons on messenger RNAs (mRNAs) and deliver the appropriate amino acid to the ribosome for protein synthesis. The human genome encodes more than 500 tRNA genes but their individual contribution to the cellular tRNA pool is unclear. In recent years, novel methods were developed to improve the quantification of tRNA gene expression, most of which rely on next-generation sequencing such as small RNA-Seq applied to tRNAs (tRNA-Seq). In a previous study, we presented a bioinformatics strategy to analyse tRNA-Seq datasets that we named 'isodecoder-specific tRNA gene contribution profiling' (Iso-tRNA-CP). Using Iso-tRNA-CP, we showed that tRNA gene expression is cell type- and tissue-specific and that this process can regulate tRNA-derived fragments abundance. An additional observation that stems from that work is that approximately half of human tRNA genes appeared silent or poorly expressed. In this commentary, I discuss this finding in light of the current literature and speculate on potential functions that transcriptionally silent tRNA genes may play. Studying silent tRNA genes may offer a unique opportunity to unravel novel mechanisms of cell regulation associated to tRNA biology.

Pathways to disease from natural variations in human cytoplasmic tRNAs

Perfectly accurate translation of mRNA into protein is not a prerequisite for life. Resulting from errors in protein synthesis, mistranslation occurs in all cells, including human cells. The human genome encodes >600 tRNA genes, providing both the raw material for genetic variation and a buffer to ensure that resulting translation errors occur at tolerable levels. On the basis of data from the 1000 Genomes Project, we highlight the unanticipated prevalence of mistranslating tRNA variants in the human population and review studies on synthetic and natural tRNA mutations that cause mistranslation or de-regulate protein synthesis. Although mitochondrial tRNA variants are well known to drive human diseases, including developmental disorders, few studies have revealed a role for human cytoplasmic tRNA mutants in disease. In the context of the unexpectedly large number of tRNA variants in the human population, the emerging literature suggests that human diseases may be affected by natural tRNA variants that cause mistranslation or de-regulate tRNA expression and nucleotide modification. This review highlights examples relevant to genetic disorders, cancer, and neurodegeneration in which cytoplasmic tRNA variants directly cause or exacerbate disease and disease-linked phenotypes in cells, animal models, and humans. In the near future, tRNAs may be recognized as useful genetic markers to predict the onset or severity of human disease.