Search

Translation fidelity relies on accurate aminoacylation of transfer RNAs (tRNAs) by aminoacyl-tRNA synthetases (AARSs). AARSs specific for alanine (Ala), leucine (Leu), serine, and pyrrolysine do not recognize the anticodon bases. Single nucleotide anticodon variants in their cognate tRNAs can lead to mistranslation. Human genomes include both rare and more common mistranslating tRNA variants. We investigated three rare human tRNALeu variants that mis-incorporate Leu at phenylalanine or tryptophan codons. Expression of each tRNALeu anticodon variant in neuroblastoma cells caused defects in fluorescent protein production without significantly increased cytotoxicity under normal conditions or in the context of proteasome inhibition. Using tRNA sequencing and mass spectrometry we confirmed that each tRNALeu variant was expressed and generated mistranslation with Leu. To probe the flexibility of the entire genetic code towards Leu mis-incorporation, we created 64 yeast strains to express all possible tRNALeu anticodon variants in a doxycycline-inducible system. While some variants showed mild or no growth defects, many anticodon variants, enriched with G/C at positions 35 and 36, including those replacing Leu for proline, arginine, alanine, or glycine, caused dramatic reductions in growth. Differential phenotypic defects were observed for tRNALeu mutants with synonymous anticodons and for different tRNALeu isoacceptors with the same anticodon. A comparison to tRNAAla anticodon variants demonstrates that Ala mis-incorporation is more tolerable than Leu at nearly every codon. The data show that the nature of the amino acid substitution, the tRNA gene, and the anticodon are each important factors that influence the ability of cells to tolerate mistranslating tRNAs.

The cellular response to stress is an important determinant of disease pathogenesis. Uncovering the molecular fingerprints of distinct stress responses may identify novel biomarkers and key signaling pathways for different diseases. Emerging evidence shows that transfer RNA-derived small RNAs (tDRs) play pivotal roles in stress responses. However, RNA modifications present on tDRs are barriers to accurately quantifying tDRs using traditional small RNA sequencing. Here, AlkB-facilitated methylation sequencing is used to generate a comprehensive landscape of cellular and extracellular tDR abundances in various cell types during different stress responses. Extracellular tDRs are found to have distinct fragmentation signatures from intracellular tDRs and these tDR signatures are better indicators of different stress responses than miRNAs. These distinct extracellular tDR fragmentation patterns and signatures are also observed in plasma from patients on cardiopulmonary bypass. It is additionally demonstrated that angiogenin and RNASE1 are themselves regulated by stressors and contribute to the stress-modulated abundance of sub-populations of cellular and extracellular tDRs. Finally, a sub-population of extracellular tDRs is identified for which AGO2 appears to be required for their expression. Together, these findings provide a detailed profile of stress-responsive tDRs and provide insight about tDR biogenesis and stability in response to cellular stressors.

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.

Transfer RNAs (tRNAs) are ubiquitous adapter molecules that link specific codons in messenger RNA (mRNA) with their corresponding amino acids during protein synthesis. The tRNA genes of Drosophila have been investigated for over half a century but have lacked systematic identification and nomenclature. Here, we review and integrate data within FlyBase and the Genomic tRNA Database (GtRNAdb) to identify the full complement of tRNA genes in the D. melanogaster nuclear and mitochondrial genomes. We apply a logical and informative nomenclature to all tRNA genes, and provide an overview of their characteristics and genomic features.

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE's ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Non-coding RNAs (ncRNA) are essential for all life, and their functions often depend on their secondary (2D) and tertiary structure. Despite the abundance of software for the visualisation of ncRNAs, few automatically generate consistent and recognisable 2D layouts, which makes it challenging for users to construct, compare and analyse structures. Here, we present R2DT, a method for predicting and visualising a wide range of RNA structures in standardised layouts. R2DT is based on a library of 3,647 templates representing the majority of known structured RNAs. R2DT has been applied to ncRNA sequences from the RNAcentral database and produced >13 million diagrams, creating the world's largest RNA 2D structure dataset. The software is amenable to community expansion, and is freely available at https://github.com/rnacentral/R2DT and a web server is found at https://rnacentral.org/r2dt .

Transfer RNAs (tRNAs) are ubiquitous across the tree of life. Although tRNA structure is highly conserved, there is still significant variation in sequence features between clades, isotypes and even isodecoders. This variation not only impacts translation, but as shown by a variety of recent studies, nontranslation-associated functions are also sensitive to small changes in tRNA sequence. Despite the rapidly growing number of sequenced genomes, there is a lack of tools for both small- and large-scale comparative genomics analysis of tRNA sequence features. Here, we have integrated over 150 000 tRNAs spanning all domains of life into tRNAviz, a web application for exploring and visualizing tRNA sequence features. tRNAviz implements a framework for determining consensus sequence features and can generate sequence feature distributions by isotypes, clades and anticodons, among other tRNA properties such as score. All visualizations are interactive and exportable. The web server is publicly available at http://trna.ucsc.edu/tRNAviz/.

Transfer RNAs are the largest, most complex non-coding RNA family, universal to all living organisms. tRNAscan-SE has been the de facto tool for predicting tRNA genes in whole genomes. The newly developed version 2.0 has incorporated advanced methodologies with improved probabilistic search software and a suite of new gene models, enabling better functional classification of predicted genes. This chapter describes the use of the UNIX command-driven and online web versions, illustrating different search modes and options.

Transfer RNAs (tRNAs) are a central component for the biological synthesis of proteins, and they are among the most highly conserved and frequently transcribed genes in all living things. Despite their clear significance for fundamental cellular processes, the forces governing tRNA evolution are poorly understood. We present evidence that transcription-associated mutagenesis and strong purifying selection are key determinants of patterns of sequence variation within and surrounding tRNA genes in humans and diverse model organisms. Remarkably, the mutation rate at broadly expressed cytosolic tRNA loci is likely between 7 and 10 times greater than the nuclear genome average. Furthermore, evolutionary analyses provide strong evidence that tRNA genes, but not their flanking sequences, experience strong purifying selection acting against this elevated mutation rate. We also find a strong correlation between tRNA expression levels and the mutation rates in their immediate flanking regions, suggesting a simple method for estimating individual tRNA gene activity. Collectively, this study illuminates the extreme competing forces in tRNA gene evolution and indicates that mutations at tRNA loci contribute disproportionately to mutational load and have unexplored fitness consequences in human populations.

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2'-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/.

N 1-methyladenosine (m1A), N 3-methylcytidine (m3C), and N 1-methylguanosine (m1G) are common in transfer RNA (tRNA) and tRNA-derived fragments. These modifications alter Watson-Crick base-pairing, and cause pauses or stops during reverse transcription required for most high-throughput RNA sequencing protocols, resulting in inefficient detection of methyl-modified RNAs. Here, we describe a procedure to demethylate RNAs containing m1A, m3C, or m1G using the Escherichia coli dealkylating enzyme AlkB, along with instructions for subsequent processing with widely used protocols for small RNA sequencing.

ABSTRACT Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′- O -methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum , particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches, and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

Small RNAs have the potential to store a secondary layer of labile biological information in the form of modified nucleotides. Emerging evidence has shown that small RNAs including microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs) and tRNA-derived small RNAs (tsRNAs) harbor a diversity of RNA modifications. These findings highlight the importance of RNA modifications in the modulation of basic properties such as RNA stability and other complex physiological processes involved in stress responses, metabolism, immunity, and epigenetic inheritance of environmentally acquired traits, among others. High-resolution, high-throughput methods for detecting, mapping and screening these small RNA modifications now provide opportunities to uncover their diagnostic potential as sensitive disease markers.

High-throughput genome sequencing continues to grow the need for rapid, accurate genome annotation and tRNA genes constitute the largest family of essential, ever-present non-coding RNA genes. Newly developed tRNAscan-SE 2.0 has advanced the state-of-the-art methodology in tRNA gene detection and functional prediction, captured by rich new content of the companion Genomic tRNA Database. Previously, web-server tRNA detection was isolated from knowledge of existing tRNAs and their annotation. In this update of the tRNAscan-SE On-line resource, we tie together improvements in tRNA classification with greatly enhanced biological context via dynamically generated links between web server search results, the most relevant genes in the GtRNAdb and interactive, rich genome context provided by UCSC genome browsers. The tRNAscan-SE On-line web server can be accessed at http://trna.ucsc.edu/tRNAscan-SE/.

Transfer RNAs represent the largest, most ubiquitous class of non-protein coding RNA genes found in all living organisms. The tRNAscan-SE search tool has become the de facto standard for annotating tRNA genes in genomes, and the Genomic tRNA Database (GtRNAdb) was created as a portal for interactive exploration of these gene predictions. Since its published description in 2009, the GtRNAdb has steadily grown in content, and remains the most commonly cited web-based source of tRNA gene information. In this update, we describe not only a major increase in the number of tRNA predictions (>367000) and genomes analyzed (>4370), but more importantly, the integration of new analytic and functional data to improve the quality and biological context of tRNA gene predictions. New information drawn from other sources includes tRNA modification data, epigenetic data, single nucleotide polymorphisms, gene expression and evolutionary conservation. A richer set of analytic data is also presented, including better tRNA functional prediction, non-canonical features, predicted structural impacts from sequence variants and minimum free energy structural predictions. Views of tRNA genes in genomic context are provided via direct links to the UCSC genome browsers. The database can be searched by sequence or gene features, and is available at http://gtrnadb.ucsc.edu/.

The rat kangaroo (long-nosed potoroo, Potorous tridactylus) is a marsupial native to Australia. Cultured rat kangaroo kidney epithelial cells (PtK) are commonly used to study cell biological processes. These mammalian cells are large, adherent, and flat, and contain large and few chromosomes-and are thus ideal for imaging intra-cellular dynamics such as those of mitosis. Despite this, neither the rat kangaroo genome nor transcriptome have been sequenced, creating a challenge for probing the molecular basis of these cellular dynamics. Here, we present the sequencing, assembly and annotation of the draft rat kangaroo de novo transcriptome. We sequenced 679 million reads that mapped to 347,323 Trinity transcripts and 20,079 Unigenes. We present statistics emerging from transcriptome-wide analyses, and analyses suggesting that the transcriptome covers full-length sequences of most genes, many with multiple isoforms. We also validate our findings with a proof-of-concept gene knockdown experiment. We expect that this high quality transcriptome will make rat kangaroo cells a more tractable system for linking molecular-scale function and cellular-scale dynamics.

Type III secretion systems (T3SS) are essential for virulence in dozens of pathogens, but are not required for growth outside the host. Therefore, the T3SS of many bacterial species are under tight regulatory control. To increase our understanding of the molecular mechanisms behind T3SS regulation, we performed a transposon screen to identify genes important for T3SS function in the food-borne pathogen Yersinia pseudotuberculosis. We identified two unique transposon insertions in YPTB2860, a gene that displays 79% identity with the E. coli iron-sulfur cluster regulator, IscR. A Y. pseudotuberculosis iscR in-frame deletion mutant (ΔiscR) was deficient in secretion of Ysc T3SS effector proteins and in targeting macrophages through the T3SS. To determine the mechanism behind IscR control of the Ysc T3SS, we carried out transcriptome and bioinformatic analysis to identify Y. pseudotuberculosis genes regulated by IscR. We discovered a putative IscR binding motif upstream of the Y. pseudotuberculosis yscW-lcrF operon. As LcrF controls transcription of a number of critical T3SS genes in Yersinia, we hypothesized that Yersinia IscR may control the Ysc T3SS through LcrF. Indeed, purified IscR bound to the identified yscW-lcrF promoter motif and mRNA levels of lcrF and 24 other T3SS genes were reduced in Y. pseudotuberculosis in the absence of IscR. Importantly, mice orally infected with the Y. pseudotuberculosis ΔiscR mutant displayed decreased bacterial burden in Peyer's patches, mesenteric lymph nodes, spleens, and livers, indicating an essential role for IscR in Y. pseudotuberculosis virulence. This study presents the first characterization of Yersinia IscR and provides evidence that IscR is critical for virulence and type III secretion through direct regulation of the T3SS master regulator, LcrF.

Pyrobaculum oguniense TE7 is an aerobic hyperthermophilic crenarchaeon isolated from a hot spring in Japan. Here we describe its main chromosome of 2,436,033 bp, with three large-scale inversions and an extra-chromosomal element of 16,887 bp. We have annotated 2,800 protein-coding genes and 145 RNA genes in this genome, including nine H/ACA-like small RNA, 83 predicted C/D box small RNA, and 47 transfer RNA genes. Comparative analyses with the closest known relative, the anaerobe Pyrobaculum arsenaticum from Italy, reveals unexpectedly high synteny and nucleotide identity between these two geographically distant species. Deep sequencing of a mixture of genomic DNA from multiple cells has illuminated some of the genome dynamics potentially shared with other species in this genus.

A great diversity of small, non-coding RNA (ncRNA) molecules with roles in gene regulation and RNA processing have been intensely studied in eukaryotic and bacterial model organisms, yet our knowledge of possible parallel roles for small RNAs (sRNA) in archaea is limited. We employed RNA-seq to identify novel sRNA across multiple species of the hyperthermophilic genus Pyrobaculum, known for unusual RNA gene characteristics. By comparing transcriptional data collected in parallel among four species, we were able to identify conserved RNA genes fitting into known and novel families. Among our findings, we highlight three novel cis-antisense sRNAs encoded opposite to key regulatory (ferric uptake regulator), metabolic (triose-phosphate isomerase), and core transcriptional apparatus genes (transcription factor B). We also found a large increase in the number of conserved C/D box sRNA genes over what had been previously recognized; many of these genes are encoded antisense to protein coding genes. The conserved opposition to orthologous genes across the Pyrobaculum genus suggests similarities to other cis-antisense regulatory systems. Furthermore, the genus-specific nature of these sRNAs indicates they are relatively recent, stable adaptations.

Abstract Background As in eukaryotes, precursor transfer RNAs in Archaea often contain introns that are removed in tRNA maturation. Two unrelated archaeal species display unique pre-tRNA processing complexity in the form of split tRNA genes, in which two to three segments of tRNAs are transcribed from different loci, then trans-spliced to form a mature tRNA. Another rare type of pre-tRNA, found only in eukaryotic algae, is permuted, where the 3' half is encoded upstream of the 5' half, and must be processed to be functional. Results Using an improved version of the gene-finding program tRNAscan-SE, comparative analyses and experimental verifications, we have now identified four novel trans-spliced tRNA genes, each in a different species of the Desulfurococcales branch of the Archaea: tRNAAsp(GUC) in Aeropyrum pernix and Thermosphaera aggregans, and tRNALys(CUU) in Staphylothermus hellenicus and Staphylothermus marinus. Each of these includes features surprisingly similar to previously studied split tRNAs, yet comparative genomic context analysis and phylogenetic distribution suggest several independent, relatively recent splitting events. Additionally, we identified the first examples of permuted tRNA genes in Archaea: tRNAiMet(CAU) and tRNATyr(GUA) in Thermofilum pendens, which appear to be permuted in the same arrangement seen previously in red alga. Conclusions Our findings illustrate that split tRNAs are sporadically spread across a major branch of the Archaea, and that permuted tRNAs are a new shared characteristic between archaeal and eukaryotic species. The split tRNA discoveries also provide new clues to their evolutionary history, supporting hypotheses for recent acquisition via viral or other mobile elements.

BackgroundHaloferax volcanii is an easily culturable moderate halophile that grows on simple defined media, is readily transformable, and has a relatively stable genome. This, in combination with its biochemical and genetic tractability, has made Hfx. volcanii a key model organism, not only for the study of halophilicity, but also for archaeal biology in general.Methodology/principal findingsWe report here the sequencing and analysis of the genome of Hfx. volcanii DS2, the type strain of this species. The genome contains a main 2.848 Mb chromosome, three smaller chromosomes pHV1, 3, 4 (85, 438, 636 kb, respectively) and the pHV2 plasmid (6.4 kb).Conclusions/significanceThe completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

Hyperthermophilic crenarchaea in the genus Pyrobaculum are notable for respiratory versatility, but relatively little is known about the genetics or regulation of crenarchaeal respiratory pathways. We measured global gene expression in Pyrobaculum aerophilum cultured with oxygen, nitrate, arsenate and ferric iron as terminal electron acceptors to identify transcriptional patterns that differentiate these pathways. We also compared genome sequences for four closely related species with diverse respiratory characteristics (Pyrobaculum arsenaticum, Pyrobaculum calidifontis, Pyrobaculum islandicum, and Thermoproteus neutrophilus) to identify genes associated with different respiratory capabilities. Specific patterns of gene expression in P. aerophilum were associated with aerobic respiration, nitrate respiration, arsenate respiration, and anoxia. Functional predictions based on these patterns include separate cytochrome oxidases for aerobic growth and oxygen scavenging, a nitric oxide-responsive transcriptional regulator, a multicopper oxidase involved in denitrification, and an archaeal arsenate respiratory reductase. We were unable to identify specific genes for iron respiration, but P. aerophilum exhibited repressive transcriptional responses to iron remarkably similar to those controlled by the ferric uptake regulator in bacteria. Together, these analyses present a genome-scale view of crenarchaeal respiratory flexibility and support a large number of functional and regulatory predictions for further investigation. The complete gene expression data set can be viewed in genomic context with the Archaeal Genome Browser at archaea.ucsc.edu.

Author(s): The ENCODE Project Consortium; Analysis Coordination; Ewan Birney (corresponding author) [1]; John A. Stamatoyannopoulos (corresponding author) [2]; Anindya Dutta (corresponding author) [3]; Roderic Guigó (corresponding author) [4, 5]; Thomas [...]

Author(s): International Human Genome Sequencing Consortium; Whitehead Institute for Biomedical Research, Center for Genome Research:; Eric S. Lander [1]; Lauren M. Linton [1]; Bruce Birren [1]; Chad Nusbaum [1]; Michael [...]