Gammaherpesvirus infection alters transfer RNA splicing and triggers tRNA cleavage

Transfer RNAs (tRNAs) are fundamental for both cellular and viral gene expression during viral infection. Moreover, mounting evidence supports a noncanonical role for tRNA cleavage products in the control of gene expression during diverse conditions of stress and infection. We previously reported that infection with the model murine gammaherpesvirus, MHV68, leads to altered tRNA transcription, suggesting that tRNA regulation may play an important role in mediating viral replication or the host response. To better understand how viral infection alters tRNA expression, we combined Ordered Two Template Relay (OTTR) with tRNA-specific bioinformatic software called tRAX to profile full-length tRNAs and fragmented tRNA-derived RNAs (tDRs) during infection with MHV68. We find that OTTR-tRAX is a powerful sequencing strategy for combined tRNA/tDR profiling and reveals that MHV68 infection triggers pre-tRNA and mature tRNA cleavage, resulting in the accumulation of specific tDRs. Fragments of virally-encoded tRNAs (virtRNAs), as well as virtRNA base modification signatures are also detectable during infection. We present evidence that tRNA splicing factors are involved in the biogenesis of MHV68-induced cleavage products from pre-tRNAs and, in the case of CLP1 kinase, impact infectious virus production. Our data offers new insights into the importance of tRNA processing during gammaherpesvirus infection.

Mistranslating the genetic code with leucine in yeast and mammalian cells

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.

Distinct Stress-Dependent Signatures of Cellular and Extracellular tRNA-Derived Small RNAs

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.

Systematic identification of tRNA genes in Drosophila melanogaster

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.

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.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

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.

R2DT is a framework for predicting and visualising RNA secondary structure using templates

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 .

RNAcentral 2021: secondary structure integration, improved sequence search and new member databases

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.

The UCSC SARS-CoV-2 Genome Browser

Background:

Researchers are generating molecular data pertaining to the SARS-CoV-2 RNA genome and its proteins at an unprecedented rate during the COVID-19 pandemic. As a result, there is a critical need for rapid and continuously updated access to the latest molecular data in a format in which all data can be quickly cross-referenced and compared. We adapted our genome browser visualization tool to the viral genome for this purpose. Molecular data, curated from published studies or from database submissions, are mapped to the viral genome and grouped together into “annotation tracks” where they can be visualized along the linear map of the viral genome sequence and programmatically downloaded in standard format for analysis.

Results:

The UCSC Genome Browser for SARS-CoV-2 (https://genome.ucsc.edu/covid19.html ) provides continuously updated access to the mutations in the many thousands of SARS-CoV-2 genomes deposited in GISAID and the international nucleotide sequencing databases, displayed alongside phylogenetic trees. These data are augmented with alignments of bat, pangolin, and other animal and human coronavirus genomes, including per-base evolutionary rate analysis. All available annotations are cross-referenced on the virus genome, including those from major databases (PDB, RFAM, IEDB, UniProt) as well as up-to-date individual results from preprints. Annotated data include predicted and validated immune epitopes, promising antibodies, RT-PCR and sequencing primers, CRISPR guides (from research, diagnostics, vaccines, and therapies), and points of interaction between human and viral genes. As a community resource, any user can add manual annotations which are quality checked and shared publicly on the browser the next day.

Conclusions:

We invite all investigators to contribute additional data and annotations to this resource to accelerate research and development activities globally. Contact us at genome-www@soe.ucsc.edu with data suggestions or requests for support for adding data. Rapid sharing of data will accelerate SARS-CoV-2 research, especially when researchers take time to integrate their data with those from other labs on a widely-used community browser platform with standardized machine-readable data formats, such as the SARS-CoV-2 Genome Browser.

tRNAviz: explore and visualize tRNA sequence features

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/.

A guide to naming human non-coding RNA genes

Research on non-coding RNA (ncRNA) is a rapidly expanding field. Providing an official gene symbol and name to ncRNA genes brings order to otherwise potential chaos as it allows unambiguous communication about each gene. The HUGO Gene Nomenclature Committee (HGNC, www.genenames.org) is the only group with the authority to approve symbols for human genes. The HGNC works with specialist advisors for different classes of ncRNA to ensure that ncRNA nomenclature is accurate and informative, where possible. Here, we review each major class of ncRNA that is currently annotated in the human genome and describe how each class is assigned a standardised nomenclature.

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.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

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.

Matching tRNA modifications in humans to their known and predicted enzymes

tRNA are post-transcriptionally modified by chemical modifications that affect all aspects of tRNA biology. An increasing number of mutations underlying human genetic diseases map to genes encoding for tRNA modification enzymes. However, our knowledge on human tRNA-modification genes remains fragmentary and the most comprehensive RNA modification database currently contains information on approximately 20% of human cytosolic tRNAs, primarily based on biochemical studies. Recent high-throughput methods such as DM-tRNA-seq now allow annotation of a majority of tRNAs for six specific base modifications. Furthermore, we identified large gaps in knowledge when we predicted all cytosolic and mitochondrial human tRNA modification genes. Only 48% of the candidate cytosolic tRNA modification enzymes have been experimentally validated in mammals (either directly or in a heterologous system). Approximately 23% of the modification genes (cytosolic and mitochondrial combined) remain unknown. We discuss these 'unidentified enzymes' cases in detail and propose candidates whenever possible. Finally, tissue-specific expression analysis shows that modification genes are highly expressed in proliferative tissues like testis and transformed cells, but scarcely in differentiated tissues, with the exception of the cerebellum. Our work provides a comprehensive up to date compilation of human tRNA modifications and their enzymes that can be used as a resource for further studies.

tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences

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 RNA genes experience exceptionally elevated mutation rates

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.

RNAcentral: a comprehensive database of non-coding RNA sequences

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/.

High-Throughput Small RNA Sequencing Enhanced by AlkB-Facilitated RNA de-Methylation (ARM-Seq)

N ¹-methyladenosine (m¹A), N ³-methylcytidine (m³C), and N ¹-methylguanosine (m¹G) 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 m¹A, m³C, or m¹G using the Escherichia coli dealkylating enzyme AlkB, along with instructions for subsequent processing with widely used protocols for small RNA sequencing.

Small RNA Modifications: Integral to Function and Disease

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.

tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes

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/.

GtRNAdb 2.0: an expanded database of transfer RNA genes identified in complete and draft genomes

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/.

ARM-seq: AlkB-facilitated RNA methylation sequencing reveals a complex landscape of modified tRNA fragments

High throughput RNA sequencing has accelerated discovery of the complex regulatory roles of small RNAs, but RNAs containing modified nucleosides may escape detection when those modifications interfere with reverse transcription during RNA-seq library preparation. Here we describe AlkB-facilitated RNA Methylation sequencing (ARM-Seq) which uses pre-treatment with Escherichia coli AlkB to demethylate 1-methyladenosine, 3-methylcytidine, and 1-methylguanosine, all commonly found in transfer RNAs. Comparative methylation analysis using ARM-Seq provides the first detailed, transcriptome-scale map of these modifications, and reveals an abundance of previously undetected, methylated small RNAs derived from tRNAs. ARM-Seq demonstrates that tRNA-derived small RNAs accurately recapitulate the m¹A modification state for well-characterized yeast tRNAs, and generates new predictions for a large number of human tRNAs, including tRNA precursors and mitochondrial tRNAs. Thus, ARM-Seq provides broad utility for identifying previously overlooked methyl-modified RNAs, can efficiently monitor methylation state, and may reveal new roles for tRNA-derived RNAs as biomarkers or signaling molecules.

RNAcentral: an international database of ncRNA sequences

The field of non-coding RNA biology has been hampered by the lack of availability of a comprehensive, up-to-date collection of accessioned RNA sequences. Here we present the first release of RNAcentral, a database that collates and integrates information from an international consortium of established RNA sequence databases. The initial release contains over 8.1 million sequences, including representatives of all major functional classes. A web portal (http://rnacentral.org) provides free access to data, search functionality, cross-references, source code and an integrated genome browser for selected species.

IscR is essential for yersinia pseudotuberculosis type III secretion and virulence

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. coliiron-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.

Bacterial pathogens use regulators that sense environmental cues to enhance their fitness. Here, we identify a transcriptional regulator in the human gut pathogen, Yersinia pseudotuberculosis, which controls a specialized secretion system essential for bacterial growth in mammalian tissues. This regulator was shown in other bacterial species to alter its activity in response to changes in iron concentration and oxidative stress, but has never been studied in Yersinia. Importantly, Y. pseudotuberculosis experiences large changes in iron bioavailability upon transit from the gut to deeper tissues and iron is a critical component in Yersinia virulence, as individuals with iron overload disorders have enhanced susceptibility to systemic Yersinia infections. Our work places this iron-modulated transcriptional regulator within the regulatory network that controls virulence gene expression in Y. pseudotuberculosis, identifying it as a potential new target for antimicrobial agents.

Reclassification of Thermoproteus neutrophilus Stetter and Zillig 1989 as Pyrobaculum neutrophilum comb. nov. based on phylogenetic analysis

The hyperthermophilic crenarchaeon Thermoproteus neutrophilus V24StaT was originally classified before sequence-based phylogenetic analysis became standard for bacterial taxonomy. Subsequent phylogenetic analyses by various groups have shown that strain V24StaT groups more closely with strains of the genus Pyrobaculum than with those in the genus Thermoproteus . Based on phylogenetic comparison of rRNA gene sequences and ribosomal proteins, we propose that strain V24StaT be reclassified as Pyrobaculum neutrophilum comb. nov., with the type strain V24StaT ( = DSM 2338T = JCM 9278T). An emended description of the genus Pyrobaculum is also presented.

Modeling the Thermoproteaceae RNase P RNA

The RNA component of the RNase P complex is found throughout most branches of the tree of life and is principally responsible for removing the 5' leader sequence from pre-tRNA transcripts during tRNA maturation. RNase P RNA has a number of universal core features, however variations in sequence and structure found in homologs across the tree of life require multiple Rfam covariance search models to detect accurately. We describe a new Rfam search model to enable efficient detection of the diminutive archaeal Type T RNase P RNAs, which are missed by existing Rfam models. Using the new model, we establish effective score detection thresholds, and detect four new RNase P RNA genes in recently completed genomes from the crenarchaeal family Thermoproteaceae.

Complete genome sequence of Pyrobaculum oguniense

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.

Comparative genomic and transcriptional analyses of CRISPR systems across the genus Pyrobaculum

Within the domain Archaea, the CRISPR immune system appears to be nearly ubiquitous based on computational genome analyses. Initial studies in bacteria demonstrated that the CRISPR system targets invading plasmid and viral DNA. Recent experiments in the model archaeon Pyrococcus furiosus have uncovered a novel RNA-targeting variant of the CRISPR system. Because our understanding of CRISPR system evolution in other archaea is limited, we have taken a comparative genomic and transcriptomic view of the CRISPR arrays across six diverse species within the crenarchaeal genus Pyrobaculum. We present transcriptional data from each of four species in the genus (P. aerophilum, P. islandicum, P. calidifontis, P. arsenaticum), analyzing mature CRISPR-associated small RNA abundance from over 20 arrays. Within the genus, there is remarkable conservation of CRISPR array structure, as well as unique features that are have not been studied in other archaeal systems. These unique features include: a nearly invariant CRISPR promoter, conservation of direct repeat families, the 5′ polarity of CRISPR-associated small RNA abundance, and a novel CRISPR-specific association with homologues of nurA and herA. These analyses provide a genus-level evolutionary perspective on archaeal CRISPR systems, broadening our understanding beyond existing non-comparative model systems.

Diversity of Antisense and Other Non-Coding RNAs in Archaea Revealed by Comparative Small RNA Sequencing in Four Pyrobaculum Species

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.

Activation of archaeal transcription mediated by recruitment of transcription factor B

Archaeal promoters consist of a TATA box and a purine-rich adjacent upstream sequence (transcription factor B (TFB)-responsive element (BRE)), which are bound by the transcription factors TATA box-binding protein (TBP) and TFB. Currently, only a few activators of archaeal transcription have been experimentally characterized. The best studied activator, Ptr2, mediates activation by recruitment of TBP. Here, we present a detailed biochemical analysis of an archaeal transcriptional activator, PF1088, which was identified in Pyrococcus furiosus by a bioinformatic approach. Operon predictions suggested that an upstream gene, pf1089, is polycistronically transcribed with pf1088. We demonstrate that PF1088 stimulates in vitro transcription by up to 7-fold when the pf1089 promoter is used as a template. By DNase I and hydroxyl radical footprinting experiments, we show that the binding site of PF1088 is located directly upstream of the BRE of pf1089. Mutational analysis indicated that activation requires the presence of the binding site for PF1088. Furthermore, we show that activation of transcription by PF1088 is dependent upon the presence of an imperfect BRE and is abolished when the pf1089 BRE is replaced with a BRE from a strong archaeal promoter. Gel shift experiments showed that TFB recruitment to the pf1089 operon is stimulated by PF1088, and TFB seems to stabilize PF1088 operator binding even in the absence of TBP. Taken together, these results represent the first biochemical evidence for a transcriptional activator working as a TFB recruitment factor in Archaea, for which the designation TFB-RF1 is suggested.

Discovery of Pyrobaculum small RNA families with atypical pseudouridine guide RNA features

In the Eukarya and Archaea, small RNA-guided pseudouridine modification is believed to be an essential step in ribosomal RNA maturation. While readily modeled and identified by computational methods in eukaryotic species, these guide RNAs have not been found in most archaeal genomes. Using high-throughput transcriptome sequencing and comparative genomics, we have identified ten novel small RNA families that appear to function as H/ACA pseudouridylation guide sRNAs, yet surprisingly lack several expected canonical features. The new RNA genes are transcribed and highly conserved across at least six species in the archaeal hyperthermophilic genus Pyrobaculum. The sRNAs exhibit a single hairpin structure interrupted by a conserved kink-turn motif, yet only two of ten families contain the complete canonical structure found in all other H/ACA sRNAs. Half of the sRNAs lack the conserved 3'-terminal ACA sequence, and many contain only a single 3' guide region rather than the canonical 5' and 3' bipartite guides. The predicted sRNA structures contain guide sequences that exhibit strong complementarity to ribosomal RNA or transfer RNA. Most of the predicted targets of pseudouridine modification are structurally equivalent to those known in other species. One sRNA appears capable of guiding pseudouridine modification at positions U54 and U55 in most or all Pyrobaculum tRNAs. We experimentally tested seven predicted pseudouridine modifications in ribosomal RNA, and all but one was confirmed. The structural insights provided by this new set of Pyrobaculum sRNAs will augment existing models and may facilitate the identification and characterization of new guide sRNAs in other archaeal species.

The UCSC Archaeal Genome Browser: 2012 update

The UCSC Archaeal Genome Browser (http://archaea.ucsc.edu) offers a graphical web-based resource for exploration and discovery within archaeal and other selected microbial genomes. By bringing together existing gene annotations, gene expression data, multiple-genome alignments, pre-computed sequence comparisons and other specialized analysis tracks, the genome browser is a powerful aggregator of varied genomic information. The genome browser environment maintains the current look-and-feel of the vertebrate UCSC Genome Browser, but also integrates archaeal and bacterial-specific tracks with a few graphic display enhancements. The browser currently contains 115 archaeal genomes, plus 31 genomes of viruses known to infect archaea. Some of the recently developed or enhanced tracks visualize data from published high-throughput RNA-sequencing studies, the NCBI Conserved Domain Database, sequences from pre-genome sequencing studies, predicted gene boundaries from three different protein gene prediction algorithms, tRNAscan-SE gene predictions with RNA secondary structures and CRISPR locus predictions. We have also developed a companion resource, the Archaeal COG Browser, to provide better search and display of arCOG gene function classifications, including their phylogenetic distribution among available archaeal genomes.

RNAcentral: A vision for an international database of RNA sequences

During the last decade there has been a great increase in the number of noncoding RNA genes identified, including new classes such as microRNAs and piRNAs. There is also a large growth in the amount of experimental characterization of these RNA components. Despite this growth in information, it is still difficult for researchers to access RNA data, because key data resources for noncoding RNAs have not yet been created. The most pressing omission is the lack of a comprehensive RNA sequence database, much like UniProt, which provides a comprehensive set of protein knowledge. In this article we propose the creation of a new open public resource that we term RNAcentral, which will contain a comprehensive collection of RNA sequences and fill an important gap in the provision of biomedical databases. We envision RNA researchers from all over the world joining a federated RNAcentral network, contributing specialized knowledge and databases. RNAcentral would centralize key data that are currently held across a variety of databases, allowing researchers instant access to a single, unified resource. This resource would facilitate the next generation of RNA research and help drive further discoveries, including those that improve food production and human and animal health. We encourage additional RNA database resources and research groups to join this effort. We aim to obtain international network funding to further this endeavor.

Discovery of permuted and recently split transfer RNAs in Archaea

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: tRNA^Asp(GUC) in Aeropyrum pernix and Thermosphaera aggregans, and tRNA^Lys(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: tRNA^iMet(CAU) and tRNA^Tyr(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.

FragSeq: transcriptome-wide RNA structure probing using high-throughput sequencing

Previous efforts to determine structures of non-coding RNA (ncRNA) probed only one RNA at a time with enzymes and chemicals, using gel electrophoresis to identify reactive positions. To accelerate RNA structure inference, we have developed FragSeq, a high-throughput RNA structure probing method that uses high-throughput RNA sequencing on fragments generated by nuclease P1, which specifically cleaves single stranded nucleic acids. In experiments probing the entire mouse nuclear transcriptome, we show that we can accurately and simultaneously map single-stranded regions (ssRNA) in multiple ncRNAs with known structure. We carried out probing in two cell types to demonstrate reproducibility. We also identified and experimentally validated structured regions in ncRNAs never previously probed.

The complete genome sequence of Haloferax volcanii DS2, a model archaeon

BACKGROUND

Haloferax 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 FINDINGS

We 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/SIGNIFICANCE

The completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

GtRNAdb: a database of transfer RNA genes detected in genomic sequence

Transfer RNAs (tRNAs) represent the single largest, best-understood class of non-protein coding RNA genes found in all living organisms. By far, the major source of new tRNAs is computational identification of genes within newly sequenced genomes. To organize the rapidly growing collection and enable systematic analyses, we created the Genomic tRNA Database (GtRNAdb), currently including over 74 000 tRNA genes predicted from 740 species. The web resource provides overview statistics of tRNA genes within each analyzed genome, including information by isotype and genetic locus, easily downloadable primary sequences, graphical secondary structures and multiple sequence alignments. Direct links for each gene to UCSC eukaryotic and microbial genome browsers provide graphical display of tRNA genes in the context of all other local genetic information. The database can be searched by primary sequence similarity, tRNA characteristics or phylogenetic group. The database is publicly available at http://gtrnadb.ucsc.edu.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Microarray analysis of the hyperthermophilic archaeon Pyrococcus furiosus exposed to gamma irradiation

The remarkable survival of the hyperthermophilic archaeon Pyrococcus furiosus to ionizing radiation was previously demonstrated. Using a time course study and whole-genome microarray analyses of mRNA transcript levels, the genes and regulatory pathways involved in the repair of lesions produced by ionizing irradiation (oxidative damage and DNA strand breaks) in P. furiosus were investigated. Data analyses showed that radA, encoding the archaeal homolog of the RecA/Rad51 recombinase, was moderately up regulated by irradiation and that a putative DNA-repair gene cluster was specifically induced by exposure to ionizing radiation. This novel repair system appears to be unique to thermophilic archaea and bacteria and is suspected to be involved in translesion synthesis. Genes that encode for a putative Dps-like iron-chelating protein and two membrane-bound oxidoreductases were differentially expressed following gamma irradiation, potentially in response to oxidative stress. Surprisingly, the many systems involved in oxygen detoxification and redox homeostasis appeared to be constitutively expressed. Finally, we identified several transcriptional regulators and protein kinases highly regulated in response to gamma irradiation.

The UCSC Archaeal Genome Browser

As more archaeal genomes are sequenced, effective research and analysis tools are needed to integrate the diverse information available for any given locus. The feature-rich UCSC Genome Browser, created originally to annotate the human genome, can be applied to any sequenced organism. We have created a UCSC Archaeal Genome Browser, available at , currently with 26 archaeal genomes. It displays G/C content, gene and operon annotation from multiple sources, sequence motifs (promoters and Shine-Dalgarno), microarray data, multi-genome alignments and protein conservation across phylogenetic and habitat categories. We encourage submission of new experimental and bioinformatic analysis from contributors. The purpose of this tool is to aid biological discovery and facilitate greater collaboration within the archaeal research community.

A computational screen for mammalian pseudouridylation guide H/ACA RNAs

The box H/ACA RNA gene family is one of the largest non-protein-coding gene families in eukaryotes and archaea. Recently, we developed snoGPS, a computational screening program for H/ACA snoRNAs, and applied it to Saccharomyces cerevisiae. We report here results of extending our method to screen for H/ACA RNAs in multiple large genomes of related species, and apply it to the human, mouse, and rat genomes. Because of the 250-fold larger search space compared to S. cerevisiae, significant enhancements to our algorithms were required. Complementing extensive cloning experiments performed by others, our findings include the detection and experimental verification of seven new mammalian H/ACA RNAs and the prediction of 23 new H/ACA RNA pseudouridine guide assignments. These assignments include four for H/ACA RNAs previously classified as orphan H/ACA RNAs with no known targets. We also determined systematic syntenic conservation among human and mouse H/ACA RNAs. With this work, 82 of 97 ribosomal RNA pseudouridines and 18 of 32 spliceosomal RNA pseudouridines in mammals have been linked to H/ACA guide RNAs.

The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs

Transfer RNAs (tRNAs) and small nucleolar RNAs (snoRNAs) are two of the largest classes of non-protein-coding RNAs. Conventional gene finders that detect protein-coding genes do not find tRNA and snoRNA genes because they lack the codon structure and statistical signatures of protein-coding genes. Previously, we developed tRNAscan-SE, snoscan and snoGPS for the detection of tRNAs, methylation-guide snoRNAs and pseudouridylation-guide snoRNAs, respectively. tRNAscan-SE is routinely applied to completed genomes, resulting in the identification of thousands of tRNA genes. Snoscan has successfully detected methylation-guide snoRNAs in a variety of eukaryotes and archaea, and snoGPS has identified novel pseudouridylation-guide snoRNAs in yeast and mammals. Although these programs have been quite successful at RNA gene detection, their use has been limited by the need to install and configure the software packages on UNIX workstations. Here, we describe online implementations of these RNA detection tools that make these programs accessible to a wider range of research biologists. The tRNAscan-SE, snoscan and snoGPS servers are available at http://lowelab.ucsc.edu/tRNAscan-SE/, http://lowelab.ucsc.edu/snoscan/ and http://lowelab.ucsc.edu/snoGPS/, respectively.

Genome-wide searching for pseudouridylation guide snoRNAs: analysis of the Saccharomyces cerevisiae genome

One of the largest families of small RNAs in eukaryotes is the H/ACA small nucleolar RNAs (snoRNAs), most of which guide RNA pseudouridine formation. So far, an effective computational method specifically for identifying H/ACA snoRNA gene sequences has not been established. We have developed snoGPS, a program for computationally screening genomic sequences for H/ACA guide snoRNAs. The program implements a deterministic screening algorithm combined with a probabilistic model to score gene candidates. We report here the results of testing snoGPS on the budding yeast Saccharomyces cerevisiae. Six candidate snoRNAs were verified as novel RNA transcripts, and five of these were verified as guides for pseudouridine formation at specific sites in ribosomal RNA. We also predicted 14 new base-pairings between snoRNAs and known pseudouridine sites in S.cerevisiae rRNA, 12 of which were verified by gene disruption and loss of the cognate pseudouridine site. Our findings include the first prediction and verification of snoRNAs that guide pseudouridine modification at more than two sites. With this work, 41 of the 44 known pseudouridine modifications in S.cerevisiae rRNA have been linked with a verified snoRNA, providing the most complete accounting of the H/ACA snoRNAs that guide pseudouridylation in any species.

Plant snoRNA database

The Plant snoRNA database (http://www.scri.sari.ac.uk/plant_snoRNA/) provides information on small nucleolar RNAs from Arabidopsis and eighteen other plant species. Information includes sequences, expression data, methylation and pseudouridylation target modification sites, initial gene organization (polycistronic, single gene and intronic) and the number of gene variants. The Arabidopsis information is divided into box C/D and box H/ACA snoRNAs, and within each of these groups, by target sites in rRNA, snRNA or unknown. Alignments of orthologous genes and gene variants from different plant species are available for many snoRNA genes. Plant snoRNA genes have been given a standard nomenclature, designed wherever possible, to provide a consistent identity with yeast and human orthologues.

Archaeal guide RNAs function in rRNA modification in the eukaryotic nucleus

In eukaryotes, many Box C/D small nucleolar RNAs base pair with ribosomal RNA through short complementary guide sequences, thereby marking up to 100 individual nucleotides of ribosomal RNA for 2'-O-methylation. Function of the eukaryotic Box C/D RNAs depends upon interaction with at least six known proteins. Box C/D RNAs are not known to exist in Bacteria but were recently identified in Archaea by biochemical analysis and computational genomic screens and have likely evolved independently in Archaea and Eukarya for more than 2000 million years. We have microinjected Box C/D RNAs from Pyrococcus furiosus, a hyperthermophilic archaeon, into the nuclei of oocytes from the aquatic frog Xenopus laevis. Our results show that Box C/D RNAs derived from this prokaryote are retained in the nucleus, localize to nucleoli, and interact with the X. laevis Box C/D RNA binding proteins fibrillarin, Nop56, and Nop58. Furthermore, we have demonstrated the ability of archaeal Box C/D RNAs to direct site-specific 2'-O-methylation of ribosomal RNA. Our studies have revealed the remarkable ability of archaeal Box C/D RNAs to assemble into functional RNA-protein complexes in the eukaryotic nucleus.

Widespread use of TATA elements in the core promoters for RNA polymerases III, II, and I in fission yeast

In addition to directing transcription initiation, core promoters integrate input from distal regulatory elements. Except for rare exceptions, it has been generally found that eukaryotic tRNA and rRNA genes do not contain TATA promoter elements and instead use protein-protein interactions to bring the TATA-binding protein (TBP), to the core promoter. Genomewide analysis revealed TATA elements in the core promoters of tRNA and 5S rRNA (Pol III), U1 to U5 snRNA (Pol II), and 37S rRNA (Pol I) genes in Schizosaccharomyces pombe. Using tRNA-dependent suppression and other in vivo assays, as well as in vitro transcription, we demonstrated an obligatory requirement for upstream TATA elements for tRNA and 5S rRNA expression in S. pombe. The Pol III initiation factor Brf is found in complexes with TFIIIC and Pol III in S. pombe, while TBP is not, consistent with independent recruitment of TBP by TATA. Template commitment assays are consistent with this and confirm that the mechanisms of transcription complex assembly and initiation by Pol III in S. pombe differ substantially from those in other model organisms. The results were extended to large-rRNA synthesis, as mutation of the TATA element in the Pol I promoter also abolishes rRNA expression in fission yeast. A survey of other organisms' genomes reveals that a substantial number of eukaryotes may use widespread TATAs for transcription. These results indicate the presence of TATA-unified transcription systems in contemporary eukaryotes and provide insight into the residual need for TBP by all three Pols in other eukaryotes despite a lack of TATA elements in their promoters.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.