Landscape of RNA pseudouridylation in archaeon Sulfolobus islandicus.

Pseudouridine, one of the most abundant RNA modifications, is synthesized by stand-alone or RNA-guided pseudouridine synthases. Here, we comprehensively mapped pseudouridines in rRNAs, tRNAs and small RNAs in the archaeon Sulfolobus islandicus and identified Cbf5-associated H/ACA RNAs. Through genetic deletion and in vitro modification assays, we determined the responsible enzymes for these modifications. The pseudouridylation machinery in S. islandicus consists of the stand-alone enzymes aPus7 and aPus10, and six H/ACA RNA-guided enzymes that account for all identified pseudouridines. These H/ACA RNAs guide the modification of all eleven sites in rRNAs, two sites in tRNAs, and two sites in CRISPR RNAs. One H/ACA RNA shows exceptional versatility by targeting eight different sites. aPus7 and aPus10 are responsible for modifying positions 13, 54 and 55 in tRNAs. We identified four atypical H/ACA RNAs that lack the lower stem and the ACA motif and confirmed their function both in vivo and in vitro. Intriguingly, atypical H/ACA RNAs can be modified by Cbf5 in a guide-independent manner. Our data provide the first global view of pseudouridylation in archaea and reveal unexpected structures, substrates, and activities of archaeal H/ACA RNPs.

The archaeal Lsm protein from Pyrococcus furiosus binds co-transcriptionally to poly(U)-rich target RNAs.

Posttranscriptional processes in Bacteria include the association of small regulatory RNAs (sRNA) with a target mRNA. The sRNA/mRNA annealing process is often mediated by an RNA chaperone called Hfq. The functional role of bacterial and eukaryotic Lsm proteins is partially understood, whereas knowledge about archaeal Lsm proteins is scarce. Here, we used the genetically tractable archaeal hyperthermophile Pyrococcus furiosus to identify the protein interaction partners of the archaeal Sm-like proteins (PfuSmAP1) using mass spectrometry and performed a transcriptome-wide binding site analysis of PfuSmAP1. Most of the protein interaction partners we found are part of the RNA homoeostasis network in Archaea including ribosomal proteins, the exosome, RNA-modifying enzymes, but also RNA polymerase subunits, and transcription factors. We show that PfuSmAP1 preferentially binds messenger RNAs and antisense RNAs recognizing a gapped poly(U) sequence with high affinity. Furthermore, we found that SmAP1 co-transcriptionally associates with target RNAs. Our study reveals that in contrast to bacterial Hfq, PfuSmAP1 does not affect the transcriptional activity or the pausing behaviour of archaeal RNA polymerases. We propose that PfuSmAP1 recruits antisense RNAs to target mRNAs and thereby executes its putative regulatory function on the posttranscriptional level.

Reversible RNA phosphorylation stabilizes tRNA for cellular thermotolerance.

Post-transcriptional modifications have critical roles in tRNA stability and function1-4. In thermophiles, tRNAs are heavily modified to maintain their thermal stability under extreme growth temperatures5,6. Here we identified 2'-phosphouridine (Up) at position 47 of tRNAs from thermophilic archaea. Up47 confers thermal stability and nuclease resistance to tRNAs. Atomic structures of native archaeal tRNA showed a unique metastable core structure stabilized by Up47. The 2'-phosphate of Up47 protrudes from the tRNA core and prevents backbone rotation during thermal denaturation. In addition, we identified the arkI gene, which encodes an archaeal RNA kinase responsible for Up47 formation. Structural studies showed that ArkI has a non-canonical kinase motif surrounded by a positively charged patch for tRNA binding. A knockout strain of arkI grew slowly at high temperatures and exhibited a synthetic growth defect when a second tRNA-modifying enzyme was depleted. We also identified an archaeal homologue of KptA as an eraser that efficiently dephosphorylates Up47 in vitro and in vivo. Taken together, our findings show that Up47 is a reversible RNA modification mediated by ArkI and KptA that fine-tunes the structural rigidity of tRNAs under extreme environmental conditions.

Insights into synthesis and function of KsgA/Dim1-dependent rRNA modifications in archaea.

Ribosomes are intricate molecular machines ensuring proper protein synthesis in every cell. Ribosome biogenesis is a complex process which has been intensively analyzed in bacteria and eukaryotes. In contrast, our understanding of the in vivo archaeal ribosome biogenesis pathway remains less characterized. Here, we have analyzed the in vivo role of the almost universally conserved ribosomal RNA dimethyltransferase KsgA/Dim1 homolog in archaea. Our study reveals that KsgA/Dim1-dependent 16S rRNA dimethylation is dispensable for the cellular growth of phylogenetically distant archaea. However, proteomics and functional analyses suggest that archaeal KsgA/Dim1 and its rRNA modification activity (i) influence the expression of a subset of proteins and (ii) contribute to archaeal cellular fitness and adaptation. In addition, our study reveals an unexpected KsgA/Dim1-dependent variability of rRNA modifications within the archaeal phylum. Combining structure-based functional studies across evolutionary divergent organisms, we provide evidence on how rRNA structure sequence variability (re-)shapes the KsgA/Dim1-dependent rRNA modification status. Finally, our results suggest an uncoupling between the KsgA/Dim1-dependent rRNA modification completion and its release from the nascent small ribosomal subunit. Collectively, our study provides additional understandings into principles of molecular functional adaptation, and further evolutionary and mechanistic insights into an almost universally conserved step of ribosome synthesis.

Expansion segments in bacterial and archaeal 5S ribosomal RNAs.

The large ribosomal RNAs of eukaryotes frequently contain expansion sequences that add to the size of the rRNAs but do not affect their overall structural layout and are compatible with major ribosomal function as an mRNA translation machine. The expansion of prokaryotic ribosomal RNAs is much less explored. In order to obtain more insight into the structural variability of these conserved molecules, we herein report the results of a comprehensive search for the expansion sequences in prokaryotic 5S rRNAs. Overall, 89 expanded 5S rRNAs of 15 structural types were identified in 15 archaeal and 36 bacterial genomes. Expansion segments ranging in length from 13 to 109 residues were found to be distributed among 17 insertion sites. The strains harboring the expanded 5S rRNAs belong to the bacterial orders Clostridiales, Halanaerobiales, Thermoanaerobacterales, and Alteromonadales as well as the archael order Halobacterales When several copies of a 5S rRNA gene are present in a genome, the expanded versions may coexist with normal 5S rRNA genes. The insertion sequences are typically capable of forming extended helices, which do not seemingly interfere with folding of the conserved core. The expanded 5S rRNAs have largely been overlooked in 5S rRNA databases.

Structural Studies of HNA Substrate Specificity in Mutants of an Archaeal DNA Polymerase Obtained by Directed Evolution.

Archaeal DNA polymerases from the B-family (polB) have found essential applications in biotechnology. In addition, some of their variants can accept a wide range of modified nucleotides or xenobiotic nucleotides, such as 1,5-anhydrohexitol nucleic acid (HNA), which has the unique ability to selectively cross-pair with DNA and RNA. This capacity is essential to allow the transmission of information between different chemistries of nucleic acid molecules. Variants of the archaeal polymerase from Thermococcus gorgonarius, TgoT, that can either generate HNA from DNA (TgoT_6G12) or DNA from HNA (TgoT_RT521) have been previously identified. To understand how DNA and HNA are recognized and selected by these two laboratory-evolved polymerases, we report six X-ray structures of these variants, as well as an in silico model of a ternary complex with HNA. Structural comparisons of the apo form of TgoT_6G12 together with its binary and ternary complexes with a DNA duplex highlight an ensemble of interactions and conformational changes required to promote DNA or HNA synthesis. MD simulations of the ternary complex suggest that the HNA-DNA hybrid duplex remains stable in the A-DNA helical form and help explain the presence of mutations in regions that would normally not be in contact with the DNA if it were not in the A-helical form. One complex with two incorporated HNA nucleotides is surprisingly found in a one nucleotide-backtracked form, which is new for a DNA polymerase. This information can be used for engineering a new generation of more efficient HNA polymerase variants.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea.

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac4C nucleotides in thermophilic archaeal tRNAs, the presence of m7G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

A newly identified duplex RNA unwinding activity of archaeal RNase J depends on processive exoribonucleolysis coupled steric occlusion by its structural archaeal loops.

RNase J is a prokaryotic 5'-3' exo/endoribonuclease that functions in mRNA decay and rRNA maturation. Here, we report a novel duplex unwinding activity of mpy-RNase J, an archaeal RNase J from Methanolobus psychrophilus, which enables it to degrade duplex RNAs with hairpins up to 40 bp when linking a 5' single-stranded overhangs of ≥ 7 nt, corresponding to the RNA channel length. A 6-nt RNA-mpy-RNase J-S247A structure reveals the RNA-interacting residues and a steric barrier at the RNA channel entrance comprising two archaeal loops and two helices. Mutagenesis of the residues key to either exoribonucleolysis or RNA translocation diminished the duplex unwinding activity. Substitution of the residues in the steric barrier yielded stalled degradation intermediates at the duplex RNA regions. Thus, an exoribonucleolysis-driven and steric occlusion-based duplex unwinding mechanism was identified. The duplex unwinding activity confers mpy-RNase J the capability of degrading highly structured RNAs, including the bacterial REP RNA, and archaeal mRNAs, rRNAs, tRNAs, SRPs, RNase P and CD-box RNAs, providing an indicative of the potential key roles of mpy-RNase J in pleiotropic RNA metabolisms. Hydrolysis-coupled duplex unwinding activity was also detected in a bacterial RNase J, which may use a shared but slightly different unwinding mechanism from archaeal RNase Js, indicating that duplex unwinding is a common property of the prokaryotic RNase Js.

Dynamic RNA acetylation revealed by quantitative cross-evolutionary mapping.

N4-acetylcytidine (ac4C) is an ancient and highly conserved RNA modification that is present on tRNA and rRNA and has recently been investigated in eukaryotic mRNA1-3. However, the distribution, dynamics and functions of cytidine acetylation have yet to be fully elucidated. Here we report ac4C-seq, a chemical genomic method for the transcriptome-wide quantitative mapping of ac4C at single-nucleotide resolution. In human and yeast mRNAs, ac4C sites are not detected but can be induced-at a conserved sequence motif-via the ectopic overexpression of eukaryotic acetyltransferase complexes. By contrast, cross-evolutionary profiling revealed unprecedented levels of ac4C across hundreds of residues in rRNA, tRNA, non-coding RNA and mRNA from hyperthermophilic archaea. Ac4C is markedly induced in response to increases in temperature, and acetyltransferase-deficient archaeal strains exhibit temperature-dependent growth defects. Visualization of wild-type and acetyltransferase-deficient archaeal ribosomes by cryo-electron microscopy provided structural insights into the temperature-dependent distribution of ac4C and its potential thermoadaptive role. Our studies quantitatively define the ac4C landscape, providing a technical and conceptual foundation for elucidating the role of this modification in biology and disease4-6.

Archaeosine Modification of Archaeal tRNA: Role in Structural Stabilization.

Archaeosine (G+) is a structurally complex modified nucleoside found quasi-universally in the tRNA of Archaea and located at position 15 in the dihydrouridine loop, a site not modified in any tRNA outside the Archaea G+ is characterized by an unusual 7-deazaguanosine core structure with a formamidine group at the 7-position. The location of G+ at position 15, coupled with its novel molecular structure, led to a hypothesis that G+ stabilizes tRNA tertiary structure through several distinct mechanisms. To test whether G+ contributes to tRNA stability and define the biological role of G+, we investigated the consequences of introducing targeted mutations that disrupt the biosynthesis of G+ into the genome of the hyperthermophilic archaeon Thermococcus kodakarensis and the mesophilic archaeon Methanosarcina mazei, resulting in modification of the tRNA with the G+ precursor 7-cyano-7-deazaguansine (preQ0) (deletion of arcS) or no modification at position 15 (deletion of tgtA). Assays of tRNA stability from in vitro-prepared and enzymatically modified tRNA transcripts, as well as tRNA isolated from the T. kodakarensis mutant strains, demonstrate that G+ at position 15 imparts stability to tRNAs that varies depending on the overall modification state of the tRNA and the concentration of magnesium chloride and that when absent results in profound deficiencies in the thermophily of T. kodakarensisIMPORTANCE Archaeosine is ubiquitous in archaeal tRNA, where it is located at position 15. Based on its molecular structure, it was proposed to stabilize tRNA, and we show that loss of archaeosine in Thermococcus kodakarensis results in a strong temperature-sensitive phenotype, while there is no detectable phenotype when it is lost in Methanosarcina mazei Measurements of tRNA stability show that archaeosine stabilizes the tRNA structure but that this effect is much greater when it is present in otherwise unmodified tRNA transcripts than in the context of fully modified tRNA, suggesting that it may be especially important during the early stages of tRNA processing and maturation in thermophiles. Our results demonstrate how small changes in the stability of structural RNAs can be manifested in significant biological-fitness changes.

Role of disordered regions in transferring tyrosine to its cognate tRNA.

Aminoacyl tRNA synthetase (AARS) plays an important role in transferring each amino acid to its cognate tRNA. Specifically, tyrosyl tRNA synthetase (TyrRS) is involved in various functions including protection from DNA damage due to oxidative stress, protein synthesis and cell signaling and can be an attractive target for controlling the pathogens by early inhibition of translation. TyrRS has two disordered regions, which lack a stable 3D structure in solution, and are involved in tRNA synthetase catalysis and stability. One of the disordered regions undergoes disorder-to-order transition (DOT) upon complex formation with tRNA whereas the other remains disordered (DR). In this work, we have explored the importance of these disordered regions using molecular dynamics simulations of both free and RNA-complexed states. We observed that the DOT and DR regions of the first subunit acts as a flap and interact with the acceptor arm of the tRNA. The DOT-DR flap closes when tyrosine (TyrRSTyr) is present at the active site of the complex and opens in the presence of tyrosine monophosphate (TyrRSYMP). The DOT and DR regions of the second subunit interact with the anticodon stem as well as D-loop of the tRNA, which might be involved in stabilizing the complex. The anticodon loop of the tRNA binds to the structured region present in the C-terminal of the protein, which is observed to be flexible during simulations. Detailed energy calculations also show that TyrRSTyr complex has stronger binding energy between tRNA and protein compared to TyrRSYMP; on the contrary, the anticodon is strongly bound in TyrRSYMP. The results obtained in the present study provide additional insights for understanding catalysis and the involvement of disordered regions in Tyr transfer to cognate tRNA.

In Vivo RNA Chemical Footprinting Analysis in Archaea.

RNA structural conformation and dynamics govern the functional properties of all RNA/RNP. Accordingly, defining changes of RNA structure and dynamics in various conditions may provide detailed insight into how RNA structural properties regulate the function of RNA/RNP. Traditional chemical footprinting analysis using chemical modifiers allows to sample the dynamics and conformation landscape of diverse RNA/RNP. However, many chemical modifiers are limited in their capacity to provide unbiased information reflecting the in vivo RNA/RNP structural landscape. In the recent years, the development of selective-2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology that uses powerful new chemical modifiers has significantly improved in vitro and in vivo structural probing of secondary and tertiary interactions of diverse RNA species at the single nucleotide level.Although the original discovery of Archaea as an independent domain of life is intimately linked to the technological development of RNA analysis, our understanding of in vivo RNA structural conformation and dynamics in this domain of life remains scarce.This protocol describes the in vivo use of SHAPE chemistry in two evolutionary divergent model Archaea, Sulfolobus acidocaldarius and Haloferax volcanii.

The presence of the ACA box in archaeal H/ACA guide RNAs promotes atypical pseudouridylation.

Archaea and eukaryotes, in addition to protein-only enzymes, also possess ribonucleoproteins containing an H/ACA guide RNA plus four proteins that produce pseudouridine (Ψ). Although typical conditions for these RNA-guided reactions are known, certain variant conditions allow pseudouridylation. We used mutants of the two stem-loops of the Haloferax volcanii sR-h45 RNA that guides three pseudouridylations in 23S rRNA and their target RNAs to characterize modifications under various atypical conditions. The 5' stem-loop produces Ψ2605 and the 3' stem-loop produces Ψ1940 and Ψ1942. The latter two modifications require unpaired "UVUN" (V = A, C, or G) in the target and ACA box in the guide. Ψ1942 modification requires the presence of U1940 (or Ψ1940). Ψ1940 is not produced in the Ψ1942-containing substrate, suggesting a sequential modification of the two residues. The ACA box of a single stem-loop guide is not required when typically unpaired "UN" is up to 17 bases from its position in the guide, but is needed when the distance increases to 19 bases or the N is paired. However, ANA of the H box of the double stem-loop guide is needed even for the 5' typical pseudouridylation. The most 5' unpaired U in a string of U's is converted to Ψ, and in the absence of an unpaired U, a paired U can also be modified. Certain mutants of the Cbf5 protein affect pseudouridylation by the two stem-loops of sR-h45 differently. This study will help elucidate the conditions for production of nonconstitutive Ψ's, determine functions for orphan H/ACA RNAs and in target designing.

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

RNA stabilization in hyperthermophilic archaea.

Analyses of the RNA metabolism of hyperthermophilic archaea highlight the efficiency of regulatory RNAs and RNA-guided processes at extreme temperatures. These organisms must overcome the intrinsic thermolability of RNAs. Elevated levels of RNA modifications and structured GC-rich regions are observed for many universal noncoding RNA families. Guide RNAs are often protected from degradation by their presence within ribonucleoprotein complexes. Modification and ligation of RNA termini can be employed to impair exonucleolytic degradation. Finally, antisense strand transcription promotes the formation of RNA duplexes and can be used to stabilize RNA regions. In our review, we provide examples of these RNA stabilization mechanisms that have been observed in hyperthermophilic archaeal model organisms.

The role of RNA structure in translational regulation by L7Ae protein in archaea.

A recent study has shown that archaeal L7Ae binds to a putative k-turn structure in the 5'-leader of the mRNA of its structural gene to regulate translation. To function as a regulator, the RNA should be unstructured in the absence of protein, but it should adopt a k-turn-containing stem-loop on binding L7Ae. Sequence analysis of UTR sequences indicates that their k-turn elements will be unable to fold in the absence of L7Ae, and we have demonstrated this experimentally in solution using FRET for the Archaeoglobus fulgidus sequence. We have solved the X-ray crystal structure of the complex of the A. fulgidus RNA bound to its cognate L7Ae protein. The RNA adopts a standard k-turn conformation that is specifically recognized by the L7Ae protein, so stabilizing the stem-loop. In-line probing of the natural-sequence UTR shows that the RNA is unstructured in the absence of L7Ae binding, but folds on binding the protein such that the ribosome binding site is occluded. Thus, L7Ae regulates its own translation by switching the conformation of the RNA to alter accessibility.

The multistructural forms of box C/D ribonucleoprotein particles.

Structural biology studies of archaeal and yeast box C/D ribonucleoprotein particles (RNPs) reveal a surprisingly wide range of forms. If form ever follows function, the different structures of box C/D small ribonucleoprotein particles (snoRNPs) may reflect their versatile functional roles beyond what has been recognized. A large majority of box C/D RNPs serve to site-specifically methylate the ribosomal RNA, typically as independent complexes. Select members of the box C/D snoRNPs also are essential components of the megadalton RNP enzyme, the small subunit processome that is responsible for processing ribosomal RNA. Other box C/D RNPs continue to be uncovered with either unexpected or unknown functions. We summarize currently known box C/D RNP structures in this review and identify the Nop56/58 and box C/D RNA subunits as the key elements underlying the observed structural diversity, and likely, the diverse functional roles of box C/D RNPs.

The trimeric Hef-associated nuclease HAN is a 3'→5' exonuclease and is probably involved in DNA repair.

Nucleases play important roles in nucleic acid metabolism. Some archaea encode a conserved protein known as Hef-associated nuclease (HAN). In addition to its C-terminal DHH nuclease domain, HAN also has three N-terminal domains, including a DnaJ-Zinc-finger, ribosomal protein S1-like, and oligonucleotide/oligosaccharide-binding fold. To further understand HAN's function, we biochemically characterized the enzymatic properties of HAN from Pyrococcus furiosus (PfuHAN), solved the crystal structure of its DHH nuclease domain, and examined its role in DNA repair. Our results show that PfuHAN is a Mn2+-dependent 3'-exonuclease specific to ssDNA and ssRNA with no activity on blunt and 3'-recessive double-stranded DNA. Domain truncation confirmed that the intrinsic nuclease activity is dependent on the C-terminal DHH nuclease domain. The crystal structure of the DHH nuclease domain adopts a trimeric topology, with each subunit adopting a classical DHH phosphoesterase fold. Yeast two hybrid assay confirmed that the DHH domain interacts with the IDR peptide of Hef nuclease. Knockout of the han gene or its C-terminal DHH nuclease domain in Haloferax volcanii resulted in increased sensitivity to the DNA damage reagent MMS. Our results imply that HAN nuclease might be involved in repairing stalled replication forks in archaea.

Epitranscriptomics: RNA Modifications in Bacteria and Archaea.

The increasingly complex functionality of RNA is contrasted by its simple chemical composition. RNA is generally built from only four different nucleotides (adenine, guanine, cytosine, and uracil). To date, >160 chemical modifications are known to decorate RNA molecules and thereby alter their function or stability. Many RNA modifications are conserved throughout bacteria, archaea, and eukaryotes, while some are unique to each branch of life. Most known modifications occur at internal positions, while there is limited diversity at the termini. The dynamic nature of RNA modifications and newly discovered regulatory functions of some of these RNA modifications gave birth to a new field, now often referred to as "epitranscriptomics." This review highlights the major developments in this field and summarizes detection principles for internal as well as 5'-terminal mRNA modifications in prokaryotes and archaea to investigate their biological significance.

Group I introns are widespread in archaea.

Group I catalytic introns have been found in bacterial, viral, organellar, and some eukaryotic genomes, but not in archaea. All known archaeal introns are bulge-helix-bulge (BHB) introns, with the exception of a few group II introns. It has been proposed that BHB introns arose from extinct group I intron ancestors, much like eukaryotic spliceosomal introns are thought to have descended from group II introns. However, group I introns have little sequence conservation, making them difficult to detect with standard sequence similarity searches. Taking advantage of recent improvements in a computational homology search method that accounts for both conserved sequence and RNA secondary structure, we have identified 39 group I introns in a wide range of archaeal phyla, including examples of group I introns and BHB introns in the same host gene.