Combined in silico and experimental identification of the Pyrococcus abyssi H/ACA sRNAs and their target sites in ribosomal RNAs

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.

Emerging Data on the Diversity of Molecular Mechanisms Involving C/D snoRNAs

Box C/D small nucleolar RNAs (C/D snoRNAs) represent an ancient family of small non-coding RNAs that are classically viewed as housekeeping guides for the 2′-O-methylation of ribosomal RNA in Archaea and Eukaryotes. However, an extensive set of studies now argues that they are involved in mechanisms that go well beyond this function. Here, we present these pieces of evidence in light of the current comprehension of the molecular mechanisms that control C/D snoRNA expression and function. From this inventory emerges that an accurate description of these activities at a molecular level is required to let the snoRNA field enter in a second age of maturity.

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.

The human ortholog of archaeal Pus10 produces pseudouridine 54 in select tRNAs where its recognition sequence contains a modified residue

The nearly conserved U54 of tRNA is mostly converted to a version of ribothymidine (T) in Bacteria and eukaryotes and to a version of pseudouridine (Ψ) in Archaea. Conserved U55 is nearly always modified to Ψ55 in all organisms. Orthologs of TrmA and TruB that produce T54 and Ψ55, respectively, in Bacteria and eukaryotes are absent in Archaea. Pus10 produces both Ψ54 and Ψ55 in Archaea. Pus10 orthologs are found in nearly all sequenced archaeal and most eukaryal genomes, but not in yeast and bacteria. This coincides with the presence of Ψ54 in most archaeal tRNAs and some animal tRNAs, but its absence from yeast and bacteria. Moreover, Ψ54 is found in several tRNAs that function as primers for retroviral DNA synthesis. Previously, no eukaryotic tRNA Ψ54 synthase had been identified. We show here that human Pus10 can produce Ψ54 in select tRNAs, including tRNA^Lys3, the primer for HIV reverse transcriptase. This synthase activity of Pus10 is restricted to the cytoplasm and is distinct from nuclear Pus10, which is known to be involved in apoptosis. The sequence GUUCAm¹AAUC (m¹A is 1-methyladenosine) at position 53-61 of tRNA along with a stable acceptor stem results in maximum Ψ54 synthase activity. This recognition sequence is unique for a Ψ synthase in that it contains another modification. In addition to Ψ54, SF9 cells-derived recombinant human Pus10 can also generate Ψ55, even in tRNAs that do not contain the Ψ54 synthase recognition sequence. This activity may be redundant with that of TruB.

Insights into RNA-processing pathways and associated RNA-degrading enzymes in Archaea

RNA-processing pathways are at the centre of regulation of gene expression. All RNA transcripts undergo multiple maturation steps in addition to covalent chemical modifications to become functional in the cell. This includes destroying unnecessary or defective cellular RNAs. In Archaea, information on mechanisms by which RNA species reach their mature forms and associated RNA-modifying enzymes are still fragmentary. To date, most archaeal actors and pathways have been proposed in light of information gathered from Bacteria and Eukarya. In this context, this review provides a state of the art overview of archaeal endoribonucleases and exoribonucleases that cleave and trim RNA species and also of the key small archaeal proteins that bind RNAs. Furthermore, synthetic up-to-date views of processing and biogenesis pathways of archaeal transfer and ribosomal RNAs as well as of maturation of stable small non-coding RNAs such as CRISPR RNAs, small C/D and H/ACA box guide RNAs, and other emerging classes of small RNAs are described. Finally, prospective post-transcriptional mechanisms to control archaeal messenger RNA quality and quantity are discussed.

Contribution of protein Gar1 to the RNA-guided and RNA-independent rRNA:Ψ-synthase activities of the archaeal Cbf5 protein

Archaeal RNA:pseudouridine-synthase (PUS) Cbf5 in complex with proteins L7Ae, Nop10 and Gar1, and guide box H/ACA sRNAs forms ribonucleoprotein (RNP) catalysts that insure the conversion of uridines into pseudouridines (Ψs) in ribosomal RNAs (rRNAs). Nonetheless, in the absence of guide RNA, Cbf5 catalyzes the in vitro formation of Ψ2603 in Pyrococcus abyssi 23S rRNA and of Ψ55 in tRNAs. Using gene-disrupted strains of the hyperthermophilic archaeon Thermococcus kodakarensis, we studied the in vivo contribution of proteins Nop10 and Gar1 to the dual RNA guide-dependent and RNA-independent activities of Cbf5 on 23S rRNA. The single-null mutants of the cbf5, nop10, and gar1 genes are viable, but display a thermosensitive slow growth phenotype. We also generated a single-null mutant of the gene encoding Pus10, which has redundant activity with Cbf5 for in vitro formation of Ψ55 in tRNA. Analysis of the presence of Ψs within the rRNA peptidyl transferase center (PTC) of the mutants demonstrated that Cbf5 but not Pus10 is required for rRNA modification. Our data reveal that, in contrast to Nop10, Gar1 is crucial for in vivo and in vitro RNA guide-independent formation of Ψ2607 (Ψ2603 in P. abyssi) by Cbf5. Furthermore, our data indicate that pseudouridylation at orphan position 2589 (2585 in P. abyssi), for which no PUS or guide sRNA has been identified so far, relies on RNA- and Gar1-dependent activity of Cbf5.

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.

Efficient RNA pseudouridylation by eukaryotic H/ACA ribonucleoproteins requires high affinity binding and correct positioning of guide RNA

H/ACA ribonucleoproteins (H/ACA RNPs) are responsible for introducing many pseudouridines into RNAs, but are also involved in other cellular functions. Utilizing a purified and reconstituted yeast H/ACA RNP system that is active in pseudouridine formation under physiological conditions, we describe here the quantitative characterization of H/ACA RNP formation and function. This analysis reveals a surprisingly tight interaction of H/ACA guide RNA with the Cbf5p-Nop10p-Gar1p trimeric protein complex whereas Nhp2p binds comparably weakly to H/ACA guide RNA. Substrate RNA is bound to H/ACA RNPs with nanomolar affinity which correlates with the GC content in the guide-substrate RNA base pairing. Both Nhp2p and the conserved Box ACA element in guide RNA are required for efficient pseudouridine formation, but not for guide RNA or substrate RNA binding. These results suggest that Nhp2p and the Box ACA motif indirectly facilitate loading of the substrate RNA in the catalytic site of Cbf5p by correctly positioning the upper and lower parts of the H/ACA guide RNA on the H/ACA proteins. In summary, this study provides detailed insight into the molecular mechanism of H/ACA RNPs.

Archaeal fibrillarin-Nop5 heterodimer 2'-O-methylates RNA independently of the C/D guide RNP particle

Archaeal fibrillarin (aFib) is a well-characterized S-adenosyl methionine (SAM)-dependent RNA 2'-O-methyltransferase that is known to act in a large C/D ribonucleoprotein (RNP) complex together with Nop5 and L7Ae proteins and a box C/D guide RNA. In the reaction, the guide RNA serves to direct the methylation reaction to a specific site in tRNA or rRNA by sequence complementarity. Here we show that a Pyrococcus abyssi aFib-Nop5 heterodimer can alone perform SAM-dependent 2'-O-methylation of 16S and 23S ribosomal RNAs in vitro independently of L7Ae and C/D guide RNAs. Using tritium-labeling, mass spectrometry, and reverse transcription analysis, we identified three in vitro 2'-O-methylated positions in the 16S rRNA of P. abyssi, positions lying outside of previously reported pyrococcal C/D RNP methylation sites. This newly discovered stand-alone activity of aFib-Nop5 may provide an example of an ancestral activity retained in enzymes that were recruited to larger complexes during evolution.

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.

Structure-function relationships of archaeal Cbf5 during in vivo RNA-guided pseudouridylation

In Eukarya and Archaea, in addition to protein-only pseudouridine (Ψ) synthases, complexes containing one guide RNA and four proteins can also produce Ψ. Cbf5 protein is the Ψ synthase in the complex. Previously, we showed that Ψ's at positions 1940, 1942, and 2605 of Haloferax volcanii 23S rRNA are absent in a cbf5-deleted strain, and a plasmid-borne copy of cbf5 can rescue the synthesis of these Ψ's. Based on published reports of the structure of archaeal Cbf5 complexed with other proteins and RNAs, we identified several potential residues and structures in H. volcanii Cbf5, which were expected to play important roles in pseudouridylation. We mutated these structures and determined their effects on Ψ production at the three rRNA positions under in vivo conditions. Mutations of several residues in the catalytic domain and certain residues in the thumb loop either abolished Ψ's or produced partial modification; the latter indicates a slower rate of Ψ formation. The universal catalytic aspartate of Ψ synthases could be replaced by glutamate in Cbf5. A conserved histidine, which is common to Cbf5 and TruB is not needed, but another conserved histidine of Cbf5 is required for the in vivo RNA-guided Ψ formation. We also identified a previously unreported novelty in the pseudouridylation activity of Cbf5 where a single stem-loop of a guide H/ACA RNA is used to produce two closely placed Ψ's and mutations of certain residues of Cbf5 abolished one of these two Ψ's. In summary, this first in vivo study identifies several structures of an archaeal Cbf5 protein that are important for its RNA-guided pseudouridylation activity.

Large-scale identification of small noncoding RNA with strand-specific deep sequencing and characterization of a novel virulence-related sRNA in Brucella melitensis

Brucella is the causative agent of brucellosis, a worldwide epidemic zoonosis. Small noncoding RNAs (sRNAs) are important modulators of gene expression and involved in pathogenesis and stress adaptation of Brucella. In this study, using a strand-specific RNA deep-sequencing approach, we identified a global set of sRNAs expressed by B. melitensis 16M. In total, 1321 sRNAs were identified, ranging from 100 to 600 nucleotides. These sRNAs differ in their expression levels and strand and chromosomal distributions. The role of BSR0441, one of these sRNAs, in the virulence of B. melitensis 16M was further characterized. BSR0441 was highly induced during the infection of macrophages and mice. The deletion mutant of BSR0441 showed significantly reduced spleen colonization in the middle and late phases of infection. The expression of the BSR0441 target mRNA genes was also altered in the BSR0441 mutant strain during macrophage and mice infection, which is consistent with its reduced intracellular survival capacity. In summary, Brucella encodes a large number of sRNAs, which may be involved in the stress adaptation and virulence of Brucella. Further investigation of these regulators will extend our understanding of the Brucella pathogenesis mechanism and the interactions between Brucella and its hosts.

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.

Contribution of two conserved histidines to the dual activity of archaeal RNA guide-dependent and -independent pseudouridine synthase Cbf5

In all organisms, several distinct stand-alone pseudouridine synthase (PUS) family enzymes are expressed to isomerize uridine into pseudouridine (Ψ) by specific recognition of RNAs. In addition, Ψs are generated in Archaea and Eukaryotes by PUS enzymes which are organized as ribonucleoprotein particles (RNP)--the box H/ACA s/snoRNPs. For this modification system, a unique TruB-like catalytic PUS subunit is associated with various RNA guides which specifically target and secure substrate RNAs by base-pairing. The archaeal Cbf5 PUS displays the special feature of exhibiting both RNA guide-dependent and -independent activities. Structures of substrate-bound TruB and H/ACA sRNP revealed the importance of histidines in positioning the target uridine in the active site. To analyze the respective role of H60 and H77, we have generated variants carrying alanine substitutions at these positions. The impact of the mutations was analyzed for unguided modifications U(55) in tRNA and U2603 in 23S rRNA, and for activity of the box H/ACA Pab91 sRNP enzyme. H77 (H43 in TruB), but not H60, appeared to be crucial for the RNA guide-independent activity. In contrast to earlier suggestions, H60 was found to be noncritical for the activity of the H/ACA sRNP, but contributes together with H77 to the full activity of H/ACA sRNPs. The data suggest that a similar catalytic process was conserved in the two divergent pseudouridylation systems.

RNA size is a critical factor for U-containing substrate selectivity and permanent pseudouridylated product release during the RNA:Ψ-synthase reaction catalyzed by box H/ACA sRNP enzyme at high temperature

The box H/ACA small ribonucleoprotein particles (H/ACA sRNPs) are RNP enzymes that isomerize uridines (U) into pseudouridines (Ψ) in archaeal RNAs. The RNA component acts as a guide by forming base-pair interactions with the substrate RNA to specify the target nucleotide of the modification to the catalytic subunit Cbf5. Here, we have analyzed association of an H/ACA sRNP enzyme from the hyperthermophilic archaeon Pyrococcus abyssi with synthetic substrate RNAs of different length and with target nucleotide variants, and estimated their turnover at high temperature. In these conditions, we found that a short substrate, which length is restricted to the interaction with RNA guide sequence, has higher turnover rate. However, the longer substrate with additional 5' and 3' sequences non-complementary to the guide RNA is better discriminated by the U to Ψ conversion allowing the RNP enzyme to distinguish the modified product from the substrate. In addition, we identified that the conserved residue Y179 in the catalytic center of Cbf5 is crucial for substrate selectivity.

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.

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.

Role of forefinger and thumb loops in production of Ψ54 and Ψ55 in tRNAs by archaeal Pus10

Pseudouridines (Ψ) are found in structurally and functionally important regions of RNAs. Six families of Ψ synthases, TruA, TruB, TruD, RsuA, RluA, and Pus10 have been identified. Pus10 is present in Archaea and Eukarya. While most archaeal Pus10 produce both tRNA Ψ54 and Ψ55, some produce only Ψ55. Interestingly, human PUS10 has been implicated in apoptosis and Crohn's and Celiac diseases. Homology models of archaeal Pus10 proteins based on the crystal structure of human PUS10 reveal that there are subtle structural differences in all of these Pus10 proteins. These observations suggest that structural changes in homologous proteins may lead to loss, gain, or change of their functions, warranting the need to study the structure-function relationship of these proteins. Using comparison of structural models and a series of mutations, we identified forefinger loop (reminiscent of that of RluA) and an Arg and a Tyr residue of archaeal Pus10 as critical determinants for its Ψ54, but not for its Ψ55 activity. We also found that a Leu residue, in addition to the catalytic Asp, is essential for both activities. Since forefinger loop is needed for both rRNA and tRNA Ψ synthase activities of RluA, but only for tRNA Ψ54 activity of Pus10, archaeal Pus10 proteins must use a different mechanism of recognition for Ψ55 activity. We propose that archaeal Pus10 uses two distinct mechanisms for substrate uridine recognition and binding. However, since we did not observe any mutation that affected only Ψ55 activity, both mechanisms for archaeal Pus10 activities must share some common features.

Comparative study of two box H/ACA ribonucleoprotein pseudouridine-synthases: relation between conformational dynamics of the guide RNA, enzyme assembly and activity

Multiple RNA-guided pseudouridine synthases, H/ACA ribonucleoprotein particles (RNPs) which contain a guide RNA and four proteins, catalyze site-specific post-transcriptional isomerization of uridines into pseudouridines in substrate RNAs. In archaeal particles, the guide small RNA (sRNA) is anchored by the pseudouridine synthase aCBF5 and the ribosomal protein L7Ae. Protein aNOP10 interacts with both aCBF5 and L7Ae. The fourth protein, aGAR1, interacts with aCBF5 and enhances catalytic efficiency. Here, we compared the features of two H/ACA sRNAs, Pab21 and Pab91, from Pyrococcus abyssi. We found that aCBF5 binds much more weakly to Pab91 than to Pab21. Surprisingly, the Pab91 sRNP exhibits a higher catalytic efficiency than the Pab21 sRNP. We thus investigated the molecular basis of the differential efficiencies observed for the assembly and catalytic activity of the two enzymes. For this, we compared profiles of the extent of lead-induced cleavages in these sRNAs during a stepwise reconstitution of the sRNPs, and analyzed the impact of the absence of the aNOP10–L7Ae interaction. Such probing experiments indicated that the sRNAs undergo a series of conformational changes upon RNP assembly. These changes were also evaluated directly by circular dichroism (CD) spectroscopy, a tool highly adapted to analyzing RNA conformational dynamics. In addition, our results reveal that the conformation of helix P1 formed at the base of the H/ACA sRNAs is optimized in Pab21 for efficient aCBF5 binding and RNP assembly. Moreover, P1 swapping improved the assembly of the Pab91 sRNP. Nonetheless, efficient aCBF5 binding probably also relies on the pseudouridylation pocket which is not optimized for high activity in the case of Pab21.

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.

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.

Comparative analysis of RNA families reveals distinct repertoires for each domain of life

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer (HGT), and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.

In cells, DNA carries recipes for making proteins, and proteins perform chemical reactions, including replication of DNA. This interdependency raises questions for early evolution, since one molecule seemingly cannot exist without the other. A resolution to this problem is the RNA world, where RNA is postulated to have been both genetic material and primary catalyst. While artificially selected catalytic RNAs strengthen the chemical plausibility of an RNA world, a biological prediction is that some RNAs should date back to this period. In this study, we ask to what degree RNAs in extant organisms trace back to the common ancestor of cellular life. Using the Rfam RNA families database, we systematically screened genomes spanning the three domains of life (Archaea, Bacteria, Eukarya) for RNA genes, and examined how far back in evolution known RNA families can be traced. We find that 99% of RNA families are restricted to a single domain. Limited conservation within domains implies ongoing emergence of RNA functions during evolution. Of the remaining 1%, half show evidence of horizontal transfer (movement of genes between organisms), and half show an evolutionary history consistent with an RNA world. The oldest RNAs are primarily associated with protein synthesis and export.

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.

2'-O-methylation of the wobble residue of elongator pre-tRNA(Met) in Haloferax volcanii is guided by a box C/D RNA containing unique features

The wobble residue C34 of Haloferax volcanii elongator tRNA(Met) is 2'-O-methylated. Neither a protein enzyme nor a guide RNA for this modification has been described. In this study, we show that this methylation is guided by a box C/D RNA targeting the intron-containing precursor of the tRNA. This guide RNA is starkly different from its homologs. This unique RNA of approximately 75 bases, named sR-tMet, is encoded in the genomes of H. volcanii and several other haloarchaea. A unique feature of sR-tMet is that the mature RNA in H. volcanii is substantially larger than its predicted size, whereas those in other haloarchaea are as predicted. While the 5'-ends of all tested haloarchaeal sR-tMets are equivalent, H. volcanii sR-tMet possesses an additional 51-base extension at its 3' end. This extension is present in the precursor but not in the mature sR-tMet of Halobacterium sp., suggesting differential 3'-end processing of sR-tMet in these two closely related organisms. Archaeal box C/D RNAs mostly contain a K-loop at the C'/D' motif. Another unique feature of sR-tMet is that its C'/D' motif lacks either a conventional K-turn or a K-loop. Instead, it contains two tandem, sheared G•A base pairs and a pyrimidine-pyrimidine pair in the non-canonical stem; the latter may form an alternative K-turn. Gel shift assays indicate that the L7Ae protein can form a stable complex with this unusual C'/D' motif, suggesting a novel RNA structure for L7Ae interaction.

Parallel evolution of transcriptome architecture during genome reorganization

Assembly of genes into operons is generally viewed as an important process during the continual adaptation of microbes to changing environmental challenges. However, the genome reorganization events that drive this process are also the roots of instability for existing operons. We have determined that there exists a statistically significant trend that correlates the proportion of genes encoded in operons in archaea to their phylogenetic lineage. We have further characterized how microbes deal with operon instability by mapping and comparing transcriptome architectures of four phylogenetically diverse extremophiles that span the range of operon stabilities observed across archaeal lineages: a photoheterotrophic halophile (Halobacterium salinarum NRC-1), a hydrogenotrophic methanogen (Methanococcus maripaludis S2), an acidophilic and aerobic thermophile (Sulfolobus solfataricus P2), and an anaerobic hyperthermophile (Pyrococcus furiosus DSM 3638). We demonstrate how the evolution of transcriptional elements (promoters and terminators) generates new operons, restores the coordinated regulation of translocated, inverted, and newly acquired genes, and introduces completely novel regulation for even some of the most conserved operonic genes such as those encoding subunits of the ribosome. The inverse correlation (r=-0.92) between the proportion of operons with such internally located transcriptional elements and the fraction of conserved operons in each of the four archaea reveals an unprecedented view into varying stages of operon evolution. Importantly, our integrated analysis has revealed that organisms adapted to higher growth temperatures have lower tolerance for genome reorganization events that disrupt operon structures.

Pseudouridine formation in archaeal RNAs: The case of Haloferax volcanii

Pseudouridine (Ψ), the isomer of uridine, is commonly found at various positions of noncoding RNAs of all organisms. Ψ residues are formed by a number of single- or multisite specific Ψ synthases, which generally act as stand-alone proteins. In addition, in Eukarya and Archaea, specific ribonucleoprotein complexes, each containing a distinct box H/ACA guide RNA and four core proteins, can produce Ψ at many sites of different cellular RNAs. Cbf5 is the core Ψ synthase in these complexes. Using Haloferax volcanii as an archaeal model organism, we show that, contrary to eukaryotes, the Cbf5 homolog (HVO_2493) is not essential in this archaeon. The Cbf5-deleted strain of H. volcanii completely lacks Ψ at positions 1940, 1942, 2605, and 2591 (Escherichia coli positions 1915, 1917, 2572, and 2586) of its 23S rRNA, and contains reduced steady-state levels of some box H/ACA RNAs. Archaeal Cbf5 is known to have tRNA Ψ55 synthase activity in vitro but we could not confirm this activity in vivo in H. volcanii. Conversely, the Pus10 (previously PsuX) homolog (HVO_1979), which can produce tRNA Ψ55, as well as Ψ54 in vitro, is shown here to be essential in H. volcanii, whereas the corresponding tRNA Ψ55 synthases, Pus4 and TruB, are not essential in yeast and E. coli, respectively. Finally, we demonstrate that HVO_1852, the TruA/Pus3 homolog, is responsible for the pseudouridylation of position 39 in H. volcanii tRNAs and that the corresponding gene is not essential.

Structure of H/ACA RNP protein Nhp2p reveals cis/trans isomerization of a conserved proline at the RNA and Nop10 binding interface

H/ACA small nucleolar and Cajal body ribonucleoproteins (RNPs) function in site-specific pseudouridylation of eukaryotic rRNA and snRNA, rRNA processing, and vertebrate telomerase biogenesis. Nhp2, one of four essential protein components of eukaryotic H/ACA RNPs, forms a core trimer with the pseudouridylase Cbf5 and Nop10 that binds to H/ACA RNAs specifically. Crystal structures of archaeal H/ACA RNPs have revealed how the protein components interact with each other and with the H/ACA RNA. However, in place of Nhp2p, archaeal H/ACA RNPs contain L7Ae, which binds specifically to an RNA K-loop motif absent from eukaryotic H/ACA RNPs, while Nhp2 binds a broader range of RNA structures. We report solution NMR studies of Saccharomyces cerevisiae Nhp2 (Nhp2p), which reveal that Nhp2p exhibits two major conformations in solution due to cis/trans isomerization of the evolutionarily conserved Pro83. The equivalent proline is in the cis conformation in all reported structures of L7Ae and other homologous proteins. Nhp2p has the expected α-β-α fold, but the solution structures of the major conformation of Nhp2p with trans Pro83 and of Nhp2p-S82W with cis Pro83 reveal that Pro83 cis/trans isomerization affects the positions of numerous residues at the Nop10 and RNA binding interface. An S82W substitution, which stabilizes the cis conformation, also stabilizes the association of Nhp2p with H/ACA snoRNPs expressed in vivo. We propose that Pro83 plays a key role in the assembly of the eukaryotic H/ACA RNP, with the cis conformation locking in a stable Cbf5-Nop10-Nhp2 ternary complex and positioning the protein backbone to interact with the H/ACA RNA.

Identification of CRISPR and riboswitch related RNAs among novel noncoding RNAs of the euryarchaeon Pyrococcus abyssi

Background

Noncoding RNA (ncRNA) has been recognized as an important regulator of gene expression networks in Bacteria and Eucaryota. Little is known about ncRNA in thermococcal archaea except for the eukaryotic-like C/D and H/ACA modification guide RNAs.

Results

Using a combination of in silico and experimental approaches, we identified and characterized novel P. abyssi ncRNAs transcribed from 12 intergenic regions, ten of which are conserved throughout the Thermococcales. Several of them accumulate in the late-exponential phase of growth. Analysis of the genomic context and sequence conservation amongst related thermococcal species revealed two novel P. abyssi ncRNA families. The CRISPR family is comprised of crRNAs expressed from two of the four P. abyssi CRISPR cassettes. The 5'UTR derived family includes four conserved ncRNAs, two of which have features similar to known bacterial riboswitches. Several of the novel ncRNAs have sequence similarities to orphan OrfB transposase elements. Based on RNA secondary structure predictions and experimental results, we show that three of the twelve ncRNAs include Kink-turn RNA motifs, arguing for a biological role of these ncRNAs in the cell. Furthermore, our results show that several of the ncRNAs are subjected to processing events by enzymes that remain to be identified and characterized.

Conclusions

This work proposes a revised annotation of CRISPR loci in P. abyssi and expands our knowledge of ncRNAs in the Thermococcales, thus providing a starting point for studies needed to elucidate their biological function.

Regulatory RNAs in Haloferax volcanii

In organisms of all three domains of life, a plethora of sRNAs (small regulatory RNAs) exists in addition to the well-known RNAs such as rRNAs, tRNAs and mRNAs. Although sRNAs have been well studied in eukaryotes and in bacteria, the sRNA population in archaea has just recently been identified and only in a few archaeal species. In the present paper, we summarize our current knowledge about sRNAs and their function in the halophilic archaeon Haloferax volcanii. Using two different experimental approaches, 111 intergenic and 38 antisense sRNAs were identified, as well as 42 tRFs (tRNA-derived fragments). Observation of differential expression under various conditions suggests that these sRNAs might be active as regulators in gene expression like their bacterial and eukaryotic counterparts. The severe phenotypes observed upon deletion and overexpression of sRNA genes revealed that sRNAs are involved in, and important for, a variety of biological functions in H. volcanii and possibly other archaea. Investigation of the Haloferax Lsm protein suggests that this protein is involved in the archaeal sRNA pathway.

Structural and functional evidence of high specificity of Cbf5 for ACA trinucleotide

Cbf5 is the catalytic subunit of the H/ACA small nucleolar/Cajal body ribonucleoprotein particles (RNPs) responsible for site specific isomerization of uridine in ribosomal and small nuclear RNA. Recent evidence from studies on archaeal Cbf5 suggests its second functional role in modifying tRNA U55 independent of guide RNA. In order to act both as a stand-alone and a RNP pseudouridine synthase, Cbf5 must differentiate features in H/ACA RNA from those in tRNA or rRNA. Most H/ACA RNAs contain a hallmark ACA trinucleotide downstream of the H/ACA motif. Here we challenged an archaeal Cbf5 (in the form of a ternary complex with its accessory proteins Nop10 and Gar1) with T-stem-loop RNAs with or without ACA trinucleotide in the stem. Although these substrates were previously shown to be substrates for the bacterial stand-alone pseudouridine synthase TruB, the Cbf5-Nop10-Gar1 complex was only able to modify those without ACA trinucleotide. A crystal structure of Cbf5-Nop10-Gar1 trimer bound with an ACA-containing T-stem-loop revealed that the ACA trinucleotide detracted Cbf5 from the stand-alone binding mode, thereby suggesting that the H/ACA RNP-associated function of Cbf5 likely supersedes its stand-alone function.

Deficiency of the tRNATyr:Psi 35-synthase aPus7 in Archaea of the Sulfolobales order might be rescued by the H/ACA sRNA-guided machinery

Up to now, Psi formation in tRNAs was found to be catalysed by stand-alone enzymes. By computational analysis of archaeal genomes we detected putative H/ACA sRNAs, in four Sulfolobales species and in Aeropyrum pernix, that might guide Psi 35 formation in pre-tRNA(Tyr)(GUA). This modification is achieved by Pus7p in eukarya. The validity of the computational predictions was verified by in vitro reconstitution of H/ACA sRNPs using the identified Sulfolobus solfataricus H/ACA sRNA. Comparison of Pus7-like enzymes encoded by archaeal genomes revealed amino acid substitutions in motifs IIIa and II in Sulfolobales and A. pernix Pus7-like enzymes. These conserved RNA:Psi-synthase- motifs are essential for catalysis. As expected, the recombinant Pyrococcus abyssi aPus7 was fully active and acted at positions 35 and 13 and other positions in tRNAs, while the recombinant S. solfataricus aPus7 was only found to have a poor activity at position 13. We showed that the presence of an A residue 3' to the target U residue is required for P. abyssi aPus7 activity, and that this is not the case for the reconstituted S. solfataricus H/ACA sRNP. In agreement with the possible formation of Psi 35 in tRNA(Tyr)(GUA) by aPus7 in P. abyssi and by an H/ACA sRNP in S. solfataricus, the A36G mutation in the P. abyssi tRNA(Tyr)(GUA) abolished Psi 35 formation when using P. abyssi extract, whereas the A36G substitution in the S. solfataricus pre-tRNA(Tyr) did not affect Psi 35 formation in this RNA when using an S. solfataricus extract.

Small RNAs of the halophilic archaeon Haloferax volcanii

In recent years, sRNAs (small non-coding RNAs) have been found to be abundant in eukaryotes and bacteria and have been recognized as a novel class of gene expression regulators. In contrast, much less is known about sRNAs in archaea, except for snoRNAs (small nucleolar RNAs) that are involved in the modification of bases in stable RNAs. Therefore bioinformatic and experimental RNomics approaches were undertaken to search for the presence of sRNAs in the model archaeon Haloferax volcanii, resulting in more than 150 putative sRNA genes being identified. Northern blot analyses were used to study (differential) expression of sRNA genes. Several chromosomal deletion mutants of sRNA genes were generated and compared with the wild-type. It turned out that two sRNAs are essential for growth at low salt concentrations and high temperatures respectively, and one is involved in the regulation of carbon metabolism. Taken together, it could be shown that sRNAs are as abundant in H. volcanii as they are in well-studied bacterial species and that they fulfil important biological roles under specific conditions.

Archaeal Pus10 proteins can produce both pseudouridine 54 and 55 in tRNA

Pus10, a recently identified pseudouridine (Psi) synthase, does not belong to any of the five commonly identified families of Psi synthases. Pyrococcus furiosus Pus10 has been shown to produce Psi55 in tRNAs. However, in vitro studies have identified another mechanism for tRNA Psi55 production in Archaea, which uses Cbf5 and other core proteins of the H/ACA ribonucleoprotein complex, in a guide RNA-independent manner. Pus10 homologs have been observed in nearly all sequenced archaeal genomes and in some higher eukaryotes, but not in yeast and bacteria. This coincides with the presence of Psi54 in the tRNAs of Archaea and higher eukaryotes and its absence in yeast and bacteria. No tRNA Psi54 synthase has been reported so far. Here, using recombinant Methanocaldococcus jannaschii and P. furiosus Pus10, we show that these proteins can function as synthase for both tRNA Psi54 and Psi55. The two modifications seem to occur independently. Salt concentration dependent variations in these activities of both proteins are observed. The Psi54 synthase activity of M. jannaschii protein is robust, while the same activity of P. furiosus protein is weak. Probable reasons for these differences are discussed. Furthermore, unlike bacterial TruB and yeast Pus4, archaeal Pus10 does not require a U54 x A58 reverse Hoogstein base pair and pyrimidine at position 56 to convert tRNA U55 to Psi55. The homology of eukaryal Pus10 with archaeal Pus10 suggests that the former may also have a tRNA Psi54 synthase activity.

Psiscan: a computational approach to identify H/ACA-like and AGA-like non-coding RNA in trypanosomatid genomes

Background

Detection of non coding RNA (ncRNA) molecules is a major bioinformatics challenge. This challenge is particularly difficult when attempting to detect H/ACA molecules which are involved in converting uridine to pseudouridine on rRNA in trypanosomes, because these organisms have unique H/ACA molecules (termed H/ACA-like) that lack several of the features that characterize H/ACA molecules in most other organisms.

Results

We present here a computational tool called Psiscan, which was designed to detect H/ACA-like molecules in trypanosomes. We started by analyzing known H/ACA-like molecules and characterized their crucial elements both computationally and experimentally.

Next, we set up constraints based on this analysis and additional phylogenic and functional data to rapidly scan three trypanosome genomes (T. brucei, T. cruzi and L. major) for sequences that observe these constraints and are conserved among the species. In the next step, we used minimal energy calculation to select the molecules that are predicted to fold into a lowest energy structure that is consistent with the constraints. In the final computational step, we used a Support Vector Machine that was trained on known H/ACA-like molecules as positive examples and on negative examples of molecules that were identified by the computational analyses but were shown experimentally not to be H/ACA-like molecules. The leading candidate molecules predicted by the SVM model were then subjected to experimental validation.

Conclusion

The experimental validation showed 11 molecules to be expressed (4 out of 25 in the intermediate stage and 7 out of 19 in the final validation after the machine learning stage). Five of these 11 molecules were further shown to be bona fide H/ACA-like molecules. As snoRNA in trypanosomes are organized in clusters, the new H/ACA-like molecules could be used as starting points to manually search for additional molecules in their neighbourhood. All together this study increased our repertoire by fourteen H/ACA-like and six C/D snoRNAs molecules from T. brucei and L. Major. In addition the experimental analysis revealed that six ncRNA molecules that are expressed are not downregulated in CBF5 silenced cells, suggesting that they have structural features of H/ACA-like molecules but do not have their standard function. We termed this novel class of molecules AGA-like, and we are exploring their function.

This study demonstrates the power of tight collaboration between computational and experimental approaches in a combined effort to reveal the repertoire of ncRNA molecles.

RNomics and Modomics in the halophilic archaea Haloferax volcanii: identification of RNA modification genes

Background

Naturally occurring RNAs contain numerous enzymatically altered nucleosides. Differences in RNA populations (RNomics) and pattern of RNA modifications (Modomics) depends on the organism analyzed and are two of the criteria that distinguish the three kingdoms of life. If the genomic sequences of the RNA molecules can be derived from whole genome sequence information, the modification profile cannot and requires or direct sequencing of the RNAs or predictive methods base on the presence or absence of the modifications genes.

Results

By employing a comparative genomics approach, we predicted almost all of the genes coding for the t+rRNA modification enzymes in the mesophilic moderate halophile Haloferax volcanii. These encode both guide RNAs and enzymes. Some are orthologous to previously identified genes in Archaea, Bacteria or in Saccharomyces cerevisiae, but several are original predictions.

Conclusion

The number of modifications in t+rRNAs in the halophilic archaeon is surprisingly low when compared with other Archaea or Bacteria, particularly the hyperthermophilic organisms. This may result from the specific lifestyle of halophiles that require high intracellular salt concentration for survival. This salt content could allow RNA to maintain its functional structural integrity with fewer modifications. We predict that the few modifications present must be particularly important for decoding, accuracy of translation or are modifications that cannot be functionally replaced by the electrostatic interactions provided by the surrounding salt-ions. This analysis also guides future experimental validation work aiming to complete the understanding of the function of RNA modifications in Archaeal translation.

Fisher: a program for the detection of H/ACA snoRNAs using MFE secondary structure prediction and comparative genomics - assessment and update

Background

The H/ACA family of small nucleolar RNAs (snoRNAs) plays a central role in guiding the pseudouridylation of ribosomal RNA (rRNA). In an effort to systematically identify the complete set of rRNA-modifying H/ACA snoRNAs from the genome sequence of the budding yeast, Saccharomyces cerevisiae, we developed a program – Fisher – and previously presented several candidate snoRNAs based on our analysis [1].

Findings

In this report, we provide a brief update of this work, which was aborted after the publication of experimentally-identified snoRNAs [2] identical to candidates we had identified bioinformatically using Fisher. Our motivation for revisiting this work is to report on the status of the candidate snoRNAs described in [1], and secondly, to report that a modified version of Fisher together with the available multiple yeast genome sequences was able to correctly identify several H/ACA snoRNAs for modification sites not identified by the snoGPS program [3]. While we are no longer developing Fisher, we briefly consider the merits of the Fisher algorithm relative to snoGPS, which may be of use for workers considering pursuing a similar search strategy for the identification of small RNAs. The modified source code for Fisher is made available as supplementary material.

Conclusion

Our results confirm the validity of using minimum free energy (MFE) secondary structure prediction to guide comparative genomic screening for RNA families with few sequence constraints.