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

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.

Aging-induced pseudouridine synthase 10 impairs hematopoietic stem cells

Aged hematopoietic stem cells (HSC) exhibit compromised reconstitution capacity and differentiation-bias towards myeloid lineage, however, the molecular mechanism behind it remains not fully understood. In this study, we observed that the expression of pseudouridine (Ψ) synthase 10 is increased in aged hematopoietic stem and progenitor cells (HSPC) and enforced protein of Ψ synthase 10 (PUS10) recapitulates the phenotype of aged HSC, which is not achieved by its Ψ synthase activity. Consistently, we observed no difference of transcribed RNA pseudouridylation profile between young and aged HSPC. No significant alteration of hematopoietic homeostasis and HSC function is observed in young Pus10-/- mice, while aged Pus10-/- mice exhibit mild alteration of hematopoietic homeostasis and HSC function. Moreover, we observed that PUS10 is ubiquitinated by E3 ubiquitin ligase CRL4DCAF1 complex and the increase of PUS10 in aged HSPC is due to aging-declined CRL4DCAF1- mediated ubiquitination degradation signaling. Taken together, this study for the first time evaluated the role of PUS10 in HSC aging and function, and provided a novel insight into HSC rejuvenation and its clinical application.

H/ACA Small Ribonucleoproteins: Structural and Functional Comparison Between Archaea and Eukaryotes

During ribosome synthesis, ribosomal RNA is modified through the formation of many pseudouridines and methylations which contribute to ribosome function across all domains of life. In archaea and eukaryotes, pseudouridylation of rRNA is catalyzed by H/ACA small ribonucleoproteins (sRNPs) utilizing different H/ACA guide RNAs to identify target uridines for modification. H/ACA sRNPs are conserved in archaea and eukaryotes, as they share a common general architecture and function, but there are also several notable differences between archaeal and eukaryotic H/ACA sRNPs. Due to the higher protein stability in archaea, we have more information on the structure of archaeal H/ACA sRNPs compared to eukaryotic counterparts. However, based on the long history of yeast genetic and other cellular studies, the biological role of H/ACA sRNPs during ribosome biogenesis is better understood in eukaryotes than archaea. Therefore, this review provides an overview of the current knowledge on H/ACA sRNPs from archaea, in particular their structure and function, and relates it to our understanding of the roles of eukaryotic H/ACA sRNP during eukaryotic ribosome synthesis and beyond. Based on this comparison of our current insights into archaeal and eukaryotic H/ACA sRNPs, we discuss what role archaeal H/ACA sRNPs may play in the formation of ribosomes.

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.

Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism

Small RNAs have been studied in detail in domains Bacteria and Eukarya but, in the case of the domain Archaea, the knowledge is scarce and the physiological function of these small RNAs (sRNAs) is still uncertain. To extend the knowledge of sRNAs in the domain Archaea and their possible role in the regulation of the nitrogen assimilation metabolism in haloarchaea, Haloferax mediterranei has been used as a model microorganism. The bioinformatic approach has allowed for the prediction of 295 putative sRNAs genes in the genome of H. mediterranei, 88 of which have been verified by means of RNA-Sequencing (RNA-Seq). The secondary structure of these sRNAs and their possible targets have been identified. Curiously, some of them present as possible target genes relating to nitrogen assimilation, such as glutamate dehydrogenase and the nitrogen regulatory PII protein. Analysis of RNA-Seq data has also revealed differences in the expression pattern of 16 sRNAs according to the nitrogen source. Consequently, RNomic and bioinformatic approaches used in this work have allowed for the identification of new sRNAs in H. mediterranei, some of which show different expression patterns depending on the nitrogen source. This suggests that these sRNAs could be involved in the regulation of nitrogen assimilation and can constitute an important gene regulatory network.

Synthesis, Function, and Heterogeneity of snoRNA-Guided Posttranscriptional Nucleoside Modifications in Eukaryotic Ribosomal RNAs

Ribosomal RNAs contain numerous 2'-O-methylated nucleosides and pseudouridines. Methylation of the 2' oxygen of ribose moieties and isomerization of uridines into pseudouridines are catalyzed by C/D and H/ACA small nucleolar ribonucleoprotein particles, respectively. We review the composition, structure, and mode of action of archaeal and eukaryotic C/D and H/ACA particles. Most rRNA modifications cluster in functionally crucial regions of the rRNAs, suggesting they play important roles in translation. Some of these modifications promote global translation efficiency or modulate translation fidelity. Strikingly, recent quantitative nucleoside modification profiling methods have revealed that a subset of modification sites is not always fully modified. The finding of such ribosome heterogeneity is in line with the concept of specialized ribosomes that could preferentially translate specific mRNAs. This emerging concept is supported by findings that some human diseases are caused by defects in the rRNA modification machinery correlated with a significant alteration of IRES-dependent translation.

Revisiting the structure/function relationships of H/ACA(-like) RNAs: a unified model for Euryarchaea and Crenarchaea

A structural and functional classification of H/ACA and H/ACA-like motifs is obtained from the analysis of the H/ACA guide RNAs which have been identified previously in the genomes of Euryarchaea (Pyrococcus) and Crenarchaea (Pyrobaculum). A unified structure/function model is proposed based on the common structural determinants shared by H/ACA and H/ACA-like motifs in both Euryarchaea and Crenarchaea. Using a computational approach, structural and energetic rules for the guide:target RNA-RNA interactions are derived from structural and functional data on the H/ACA RNP particles. H/ACA(-like) motifs found in Pyrococcus are evaluated through the classification and their biological relevance is discussed. Extra-ribosomal targets found in both Pyrococcus and Pyrobaculum might support the hypothesis of a gene regulation mediated by H/ACA(-like) guide RNAs in archaea.

Capture, Unfolding, and Detection of Individual tRNA Molecules Using a Nanopore Device

Transfer RNAs (tRNA) are the most common RNA molecules in cells and have critical roles as both translators of the genetic code and regulators of protein synthesis. As such, numerous methods have focused on studying tRNA abundance and regulation, with the most widely used methods being RNA-seq and microarrays. Though revolutionary to transcriptomics, these assays are limited by an inability to encode tRNA modifications in the requisite cDNA. These modifications are abundant in tRNA and critical to their function. Here, we describe proof-of-concept experiments where individual tRNA molecules are examined as linear strands using a biological nanopore. This method utilizes an enzymatically ligated synthetic DNA adapter to concentrate tRNA at the lipid bilayer of the nanopore device and efficiently denature individual tRNA molecules, as they are pulled through the α-hemolysin (α-HL) nanopore. Additionally, the DNA adapter provides a loading site for ϕ29 DNA polymerase (ϕ29 DNAP), which acts as a brake on the translocating tRNA. This increases the dwell time of adapted tRNA in the nanopore, allowing us to identify the region of the nanopore signal that is produced by the translocating tRNA itself. Using adapter-modified Escherichia coli tRNA^fMet and tRNA^Lys, we show that the nanopore signal during controlled translocation is dependent on the identity of the tRNA. This confirms that adapter-modified tRNA can translocate end-to-end through nanopores and provide the foundation for future work in direct sequencing of individual transfer RNA with a nanopore-based device.

Evolution of small guide RNA genes in hyperthermophilic archaea

Profiling the RNA production in hyperthermophilic archaea revealed an abundance of small RNA-guided processes near the upper temperature limit of life. Archaea utilize the base-pairing ability of RNA guide sequences to target ribosomal RNAs, transfer RNAs, messenger RNAs, and viral genomes. Cellular processes that are guided by small RNAs include the modification of RNA molecules, trans-splicing, gene regulation, and RNA and DNA degradation. Here, a brief overview of our knowledge on small guide RNA genes in archaeal genomes is provided and examples of their putative roles in genome evolution are described.

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.

Results

Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.

Conclusion

The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-684) contains supplementary material, which is available to authorized users.

Small nucleolar RNAs and RNA-guided post-transcriptional modification

snoRNAs (small nucleolar RNAs) constitute one of the largest and best-studied classes of non-coding RNAs that confer enzymatic specificity. With associated proteins, these snoRNAs form ribonucleoprotein complexes that can direct 2'-O-methylation or pseudouridylation of target non-coding RNAs. Aided by computational methods and high-throughput sequencing, new studies have expanded the diversity of known snoRNA functions. Complexes incorporating snoRNAs have dynamic specificity, and include diverse roles in RNA silencing, telomerase maintenance and regulation of alternative splicing. Evidence that dysregulation of snoRNAs can cause human disease, including cancer, indicates that the full scope of snoRNA roles remains an unfinished story. The diversity in structure, genomic origin and function between snoRNAs found in different complexes and among different phyla illustrates the surprising plasticity of snoRNAs in evolution. The ability of snoRNAs to direct highly specific interactions with other RNAs is a consistent thread in their newly discovered functions. Because they are ubiquitous throughout Eukarya and Archaea, it is likely they were a feature of the last common ancestor of these two domains, placing their origin over two billion years ago. In the present chapter, we focus on recent advances in our understanding of these ancient, but functionally dynamic RNA-processing machines.

Small regulatory RNAs in Archaea

Small regulatory RNAs (sRNAs) are universally distributed in all three domains of life, Archaea, Bacteria, and Eukaryotes. In bacteria, sRNAs typically function by binding near the translation start site of their target mRNAs and thereby inhibit or activate translation. In eukaryotes, miRNAs and siRNAs typically bind to the 3′-untranslated region (3′-UTR) of their target mRNAs and influence translation efficiency and/or mRNA stability. In archaea, sRNAs have been identified in all species investigated using bioinformatic approaches, RNomics, and RNA-Seq. Their size can vary significantly between less than 50 to more than 500 nucleotides. Differential expression of sRNA genes has been studied using northern blot analysis, microarrays, and RNA-Seq. In addition, biological functions have been unraveled by genetic approaches, i.e., by characterization of designed mutants. As in bacteria, it was revealed that archaeal sRNAs are involved in many biological processes, including metabolic regulation, adaptation to extreme conditions, stress responses, and even in regulation of morphology and cellular behavior. Recently, the first target mRNAs were identified in archaea, including one sRNA that binds to the 5′-region of two mRNAs in Methanosarcina mazei Gö1 and a few sRNAs that bind to 3′-UTRs in Sulfolobus solfataricus, three Pyrobaculum species, and Haloferax volcanii, indicating that archaeal sRNAs appear to be able to target both the 5′-UTR or the 3′-UTRs of their respective target mRNAs. In addition, archaea contain tRNA-derived fragments (tRFs), and one tRF has been identified as a major ribosome-binding sRNA in H. volcanii, which downregulates translation in response to stress. Besides regulatory sRNAs, archaea contain further classes of sRNAs, e.g., CRISPR RNAs (crRNAs) and snoRNAs.

RNA-guided isomerization of uridine to pseudouridine--pseudouridylation

Box H/ACA ribonucleoproteins (RNPs), each consisting of one unique guide RNA and 4 common core proteins, constitute a family of complex enzymes that catalyze, in an RNA-guided manner, the isomerization of uridines to pseudouridines (Ψs) in RNAs, a reaction known as pseudouridylation. Over the years, box H/ACA RNPs have been extensively studied revealing many important aspects of these RNA modifying machines. In this review, we focus on the composition, structure, and biogenesis of H/ACA RNPs. We explain the mechanism of how this enzyme family recognizes and specifies its target uridine in a substrate RNA. We discuss the substrates of box H/ACA RNPs, focusing on rRNA (rRNA) and spliceosomal small nuclear RNA (snRNA). We describe the modification product Ψ and its contribution to RNA function. Finally, we consider possible mechanisms of the bone marrow failure syndrome dyskeratosis congenita and of prostate and other cancers linked to mutations in H/ACA RNPs.

RNA at 92 °C: the non-coding transcriptome of the hyperthermophilic archaeon Pyrococcus abyssi

The non-coding transcriptome of the hyperthermophilic archaeon Pyrococcus abyssi is investigated using the RNA-seq technology. A dedicated computational pipeline analyzes RNA-seq reads and prior genome annotation to identify small RNAs, untranslated regions of mRNAs, and cis-encoded antisense transcripts. Unlike other archaea, such as Sulfolobus and Halobacteriales, P. abyssi produces few leaderless mRNA transcripts. Antisense transcription is widespread (215 transcripts) and targets protein-coding genes that appear to evolve more rapidly than average genes. We identify at least three novel H/ACA-like guide RNAs among the newly characterized non-coding RNAs. Long 5′ UTRs in mRNAs of ribosomal proteins and amino-acid biosynthesis genes strongly suggest the presence of cis-regulatory leaders in these mRNAs. We selected a high-interest subset of non-coding RNAs based on their strong promoters, high GC-content, phylogenetic conservation, or abundance. Some of the novel small RNAs and long 5′ UTRs display high GC contents, suggesting unknown structural RNA functions. However, we were surprised to observe that most of the high-interest RNAs are AU-rich, which suggests an absence of stable secondary structure in the high-temperature environment of P. abyssi. Yet, these transcripts display other hallmarks of functionality, such as high expression or high conservation, which leads us to consider possible RNA functions that do not require extensive secondary structure.

RNA pseudouridylation: new insights into an old modification

Pseudouridine is the most abundant post-transcriptionally modified nucleotide in various stable RNAs of all organisms. Pseudouridine is derived from uridine via base-specific isomerization, resulting in an extra hydrogen-bond donor that distinguishes it from other nucleotides. In eukaryotes, uridine-to-pseudouridine isomerization is catalyzed primarily by box H/ACA RNPs, ribonucleoproteins that act as pseudouridylases. When introduced into RNA, pseudouridine contributes significantly to RNA-mediated cellular processes. It was recently discovered that pseudouridylation can be induced by stress, suggesting a regulatory role for pseudouridine. It has also been reported that pseudouridine can be artificially introduced into mRNA by box H/ACA RNPs and that such introduction can mediate nonsense-to-sense codon conversion, thus demonstrating a new means of generating coding or protein diversity.

Ribonucleoproteins in archaeal pre-rRNA processing and modification

Given that ribosomes are one of the most important cellular macromolecular machines, it is not surprising that there is intensive research in ribosome biogenesis. Ribosome biogenesis is a complex process. The maturation of ribosomal RNAs (rRNAs) requires not only the precise cleaving and folding of the pre-rRNA but also extensive nucleotide modifications. At the heart of the processing and modifications of pre-rRNAs in Archaea and Eukarya are ribonucleoprotein (RNP) machines. They are called small RNPs (sRNPs), in Archaea, and small nucleolar RNPs (snoRNPs), in Eukarya. Studies on ribosome biogenesis originally focused on eukaryotic systems. However, recent studies on archaeal sRNPs have provided important insights into the functions of these RNPs. This paper will introduce archaeal rRNA gene organization and pre-rRNA processing, with a particular focus on the discovery of the archaeal sRNP components, their functions in nucleotide modification, and their structures.

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.

RNA processing in the minimal organism Nanoarchaeum equitans

BACKGROUND

The minimal genome of the tiny, hyperthermophilic archaeon Nanoarchaeum equitans contains several fragmented genes and revealed unusual RNA processing pathways. These include the maturation of tRNA molecules via the trans-splicing of tRNA halves and genomic rearrangements to compensate for the absence of RNase P.

RESULTS

Here, the RNA processing events in the N. equitans cell are analyzed using RNA-Seq deep sequencing methodology. All tRNA half precursor and tRNA termini were determined and support the tRNA trans-splicing model. The processing of CRISPR RNAs from two CRISPR clusters was verified. Twenty-seven C/D box small RNAs (sRNAs) and a H/ACA box sRNA were identified. The C/D box sRNAs were found to flank split genes, to form dicistronic tRNA-sRNA precursors and to be encoded within the tRNAMet intron.

CONCLUSIONS

The presented data provide an overview of the production and usage of small RNAs in a cell that has to survive with a highly reduced genome. N. equitans lost many essential metabolic pathways but maintains highly active CRISPR/Cas and rRNA modification systems that appear to play an important role in genome fragmentation.

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.