Identification of CRISPR and riboswitch related RNAs among novel noncoding RNAs of the euryarchaeon Pyrococcus abyssi
© Phok et al; licensee BioMed Central Ltd. 2011
Received: 15 March 2011
Accepted: 13 June 2011
Published: 13 June 2011
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.
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.
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.
A plethora of noncoding RNAs (ncRNAs), including small RNAs that bind to proteins or base pair with target RNAs, have been found to operate at all levels of gene regulation ranging from the control of enzymatic activity to the regulation of the initiation of transcription and translation [1–3]. However, whole genome analyses of both prokaryotic and eukaryotic organisms have generally disregarded ncRNA genes. In Bacteria, systematic searches for functional intergenic regions have led to the discovery of more than 200 bacterial trans-acting ncRNAs of about 50 to 500 nucleotides mostly in E. coli[1, 4, 5], but also in other pathogenic species [2, 6–8]. Functional analysis of these ncRNAs identified many of them as regulators of bacterial stress responses. Furthermore, the recent discovery of riboswitches  and RNA-based thermosensors  in bacterial 5' untranslated regions (5'UTRs) has enlarged the range of posttranscriptional control of gene expression. In Archaea, computational and experimental analysis lead to the identification and characterization of the homologues of the eukaryotic box C/D- and H/ACA guide snoRNAs, which are involved in modification and maturation of tRNAs and rRNAs in Crenarchaea and Euryarchaea [11–14].
RNomics studies have revealed the occurrence of stable antisense RNAs and the expression of ncRNAs from intergenic regions in Sulfolobus solfataricus and Archeoglobulus fulgidus[15–17]. In silico searches using the Pyrococcus furiosus and Methanocaldococcus janashii genomes, predicted several novel ncRNAs originating from intergenic regions [18, 19]. More recently, transcriptome analysis has revealed several dozen ncRNAs in the halophilic euryarchaeon Haloferax volcannii[20, 21] and more than two hundred in the methanogenic crenarchaeon Methanosarcina mazeï. Among the ncRNAs in Archeoglobulus fulgidus, Sulfolobus solfataricus, Sulfolobus acidocaldarius[16, 23] and Pyroccocus furiosus, several correspond to ladder-like transcripts issued from repeated genomic sequences or CRISPR loci, [25, 26]. The CRISPR/Cas defense system (for review ), identified in most archaeal genomes as well as many bacterial genomes, provides acquired immunity against viruses and plasmids by targeting nucleic acid in a sequence-specific manner [28, 29]. The anatomy of a CRISPR locus has been defined as an array of short direct repeats of 20 to 50 base pairs, often containing palindromic sequences [30, 31]. Irrespective of the precise mechanism of the defensive action of CRISPR/Cas systems, there is a consensus that transcription of the CRISPR cassette initiates in or near the leader sequence followed by processing by the Cas proteins of the RNA precursor into fragments (crRNAs) corresponding to the interval of the repeats .
Several stable RNAs in the Archaea, including the box C/D and H/ACA guide RNAs, the ribosomal RNAs and RNase P, form ribonucleoprotein complexes (RNPs) with the multifunctional L7Ae protein that binds to the Kink-turn (K-turn) [17, 32, 33] or the related Kink-loop (K-loop) RNA structural motif . These widespread motifs  provide a platform for assembly of RNPs or, in the case of the S-adenosylmethionine and lysine riboswitches, orient strands that base pair to form pseudoknots [36, 37].
The goal of this study was to identify and characterize novel ncRNAs in Pyrococcus abyssi, a thermococcal archaeon (hyperthermophilic and anaerobic member of the euryarchaeal phylum), which is one of the first archaea whose genome was sequenced . Stable noncoding RNAs that have been identified in P. abyssi are ribosomal RNAs (1 16S, 1 23S, and 2 5S genes), RNase P (1 gene), tRNAs (46 genes), 7S RNA (1 gene), H/ACA guide RNAs (7 genes) and C/D guide RNAs (59 genes). Most of these genes have a significantly higher (G+C) content compared to the rest of the genome, which is AT rich. The (G+C) content of the P. abyssi genome is 44% compared to 66% and 70% in rRNA and tRNA genes, respectively. Considering the availability of several related thermococcal genomes and the AT rich character of P. abyssi genome, we performed computational searches for novel ncRNAs. We clustered InterGenic Regions (IGRs) based on primary and secondary structure features, and sequence conservation in other thermococcal genomes. Northern blotting showed that 24 out of the 82 selected IGRs are transcribed. Additional primer extension and C ircular R apid A mplification of c DNA E nds (C-RACE or CR-RT-PCR) experiments and in silico analysis showed that twelve of these transcripts have characteristics of regulatory ncRNA: three are from CRISPR cassettes, three are from mRNA 5'UTRs and six are from intergenic regions. Altogether, this study allowed us to define two novel families of ncRNA in P. abyssi, the CRIPSR and the 5'UTR derived ncRNAs.
Detection of novel ncRNAs
Novel ncRNAs validated in this study
Cons (53 nt)
Cons (53 nt)
Conserved unique locus
All except Tsi
Conserved repetitive locus
Pho, Pfu, Tko, Tsi
TATA (60 nt)
Pho, Pfu, Tko, Tsi
Specific repetitive locus
Annotation and expression of P. abyssi CRISPR loci
No signals were detected by Northern blotting with strand-specific direct or reverse probes against sequence matching spacers 2 and 7 of CRISPR 2 and spacer 2 of CRISPR 3 suggesting that these loci are not transcribed (data not shown). This tentative conclusion was confirmed by more sensitive tests involving primer extension analysis to map 5' ends of the CRISPR precursors and to identify promoters (see below). In contrast, probes against sequences matching spacer 1 of CRISPR 1 and spacers 4 and 12 of CRISPR 4 allowed the detection of transcripts in all three growth phases (Figure 1B; Table 1). The approximately 60 nt length of CRISPR-derived RNAs, named hereafter Cr1-1, Cr4-1 and Cr4-12, are of the size reported for crRNAs corresponding to a spacer sequence and part of the direct repeat [16, 24]. Moreover, primer extensions with total RNA extracted from cells in entry-stationary phase (Additional file 3, Figure S2A) accurately identified their 5' ends (Figure 1C), which correspond to the dicing site reported for the P. furiosus endoribonuclease Cas 6 . Since only single reverse transcription arrests are observed (Additional file 3, Figure S2A) even though multiple Northern blot signals are detected (Figure 1B), we propose that length heterogeneity results from partial or incomplete 3'end processing. In conclusion, the CRISPR 1 and 4 cassettes are transcribed to produce small crRNAs in P. abyssi. No signals corresponding to precursor or other intermediates were detected by Northern hybridization (Figure 1B). Nevertheless, primer extensions with oligonucleotides hybridizing to the 5' junction of their respective first direct repeat resulted in the detection of an RNA precursor and thus permitted the identification of the transcription starts of the CRISPR 1 and 4 RNA precursors (preCr1 and 4, respectively) (Figure 1C & Additional file 3Figure S2A). Similar experiments with specific primers of CRISPR 2 and 3 did not reveal any RNA precursor confirming that CRISPR 2 and 3 are not expressed. Signals were not detected in Northern hybridization with probes against the ultimate spacer 22 of CRISPR 1 and spacer 28 of CRISPR 4, which harbor degenerate repeats. These results suggest that the direct repeats are important for the processing and/or stability of the mature crRNAs. To identify transcription signals, we analyzed the 200 bp region upstream of the first direct repeat of each P. abyssi CRISPR cassette. This analysis recognized conserved AT rich motifs arguing for the presence of promoters. The CRISPR 2 leader sequence has three AT-rich and one A-rich element upstream of the first direct repeat, not corresponding to any typical P. abyssi consensus promoter sequence as defined in . In contrast, the CRISPR 1 and 4 cassettes have consensus promoter sequences 10 nt upstream of the experimentally determined 5' starts of preCr1 and 4 (Figure 1C). Two additional short elements at position +30 and -100 relative to the transcription starts were detected in CRISPR 1 and 4. These sequences are also observed in other thermococcal CRISPR cassettes at the same distance from the first direct repeat (data not shown). It should be noted that the general organization of the CRISPR 3 leader is similar to those of CRISPR 1 and 4 except for the distance between the +25 A-stretch and the first direct repeat (Figure 1C). This feature might account for the absence of transcription of the CRISPR 3 locus.
Expression of unique intergenic regions conserved among thermococcal genomes
Expression of repetitive intergenic regions conserved among thermococcal genomes
Among the three ncRNAs clustered in this category (Table 1), sRk49 shows sequence similarities with two additional regions in P. abyssi (49.2 and 49.3) and three in P. horikoshii(Additional file 5, Figure S4A). Sequence similarities preclude the design of specific probes to distinguish between sRk49 and 49.2. Northern blot signals provided evidence for a growth stage-specific transcription pattern with the most abundant species migrating about 150 nt. These transcripts could arise from either of the loci harboring an upstream TATA-like promoter (sRk49 and 49.2) (Additional file 5, Figure S4B & C). CR-RT-PCR experiments failed to detect primary transcripts, but suggested that antisense transcripts are expressed from the sRk49 operon (data not shown). Specific RNA structures or RNA modifications could have impeded the accurate amplification of the primary transcripts. Moreover, related sequences in the P. horikoshii genome are all clustered with CRISPR cassettes showing an organization comparable to that observed for the 49.2 region. All sRk49 variants are atypical in that they contain multiple short polyA stretches making RNA folding into stable helical structures implausible.
Expression of repetitive intergenic regions specific to P. abyssi
As mentioned earlier, the sRkB and sRkC loci are part of a set of six similar sequences that had no counterparts in other thermococcal genomes and are therefore specific to P. abyssi. Consistent with their lack of expression, the BRE/TATA promoter sequences are missing in the other four other P. abyssi loci (Additional file 4, Figure S3B). Strikingly, one of these loci corresponds to the 3' end of the PAB1452 open reading frame. This ORF has been identified as a 'single OrfB element' of the IS605/IS200 family of transposases .
Discovery of novel ncRNAs by combined bioinformatic approaches
The rationale for carrying out a search of ncRNAs in P. abyssi, the best studied Thermococcale, was that with the exception of the box H/ACA  and the box C/D  noncoding guide RNA families, few thermococcal ncRNA families have been described. Our study provides evidence for the growth-regulated expression of 12 novel ncRNAs in the P. abyssi genome consisting of three crRNAs from CRISPR loci, four ncRNAs from unique loci conserved throughout Thermococcales of which one appears to be a homologue of the P. furiosus SccA RNA , three ncRNAs from repetitive loci conserved throughout Thermococcales, and two ncRNAs from repetitive loci specific to P. abyssi. Since small proteins (sproteins) of less than 50 amino acids have been largely disregarded in genome annotations , we cannot exclude the possibility that the ncRNAs identified here encode sproteins or peptides. We therefore searched for open reading frames within the ncRNA sequences starting with AUG, GUG or UUG, and ending with UGA, UAG or UAA. For several of the ncRNAs, it was possible to find small ORFs of 25 to 40 amino acids. The only significant ORF identified corresponds to the distal portion of a putative transposase CDS (see below). We also analyzed the GC content of the ncRNAs. We found that it ranges from 33% to 66%. This large spectrum is consistent with the GC content for known ncRNAs (low for C/D guide RNAs, high for tRNAs and rRNAs). Since the computational analyses were restricted to intergenic regions, no cis-encoded antisense ncRNA complementary to ORFs, as identified in transcriptome analyses of Sulfolobus solfataricus and Methanosarcina mazei[22, 45] could be predicted. Below, experimental and genomic features of the new P. abyssi ncRNA families are discussed.
Processing of ncRNAs
In this study, transcripts of different length were identified at one locus suggesting that RNA processing occurs in P. abyssi. The sRkB, sRkC, sRk48 and sRk52 primary transcripts with triphoshorylated 5' ends are processed into shorter 5' end monophosphorylated transcripts through maturation reactions (Figure 3C and Figure 5C). For example, the 3' end heterogeneities observed for sRkB, sRkC, sRk48 and sRk52 RNAs could arise from trimming of the primary transcript by the 3' to 5' exonucleolytic activity carried by the Pyrococcus exosome . The 5' monophosphorylated RNAs could result from endonucleolytic or pyrophosphohydrolytic activity. Previous examples of processing in Pyrococcus were limited to tRNA processing by RNase P  and crRNA maturation by the Cas 6 endonuclease [42, 48]. Recent results suggest that the 5' to 3' exonucleolytic activity of RNase J homologues may also be involved in RNA processing in Euryarchaea and Crenarchaea [49, 50].
Our detection of crRNAs is the first demonstration of the existence of an active CRISPR defense system in P. abyssi. The CRISPR loci 1 and 4 appear to be expressed independently of growth phase in our experimental conditions. In agreement with the findings in P. furiosus, no antisense transcription of these CRISPR loci was detected (data not shown) excluding the possibility that double strand RNA intermediates are formed. This contrasts with the report that RNA is transcribed from the complementary spacer strand in S. acidocaldaricus. This difference could reflect a specific feature of the Sulfolobales or a technical issue regarding the sensitivity of the Northern blots analysis used in the studies of P. furiosus and P. abyssi. Based on our data, we propose a new annotation for the four P. abyssi CRISPR arrays that differs from earlier proposals (UCSC archaeal genome browser and ). The presence of a typical consensual promoter in the leader sequences of CRISPR 1 and CRISPR 4 correlates with the detection of CRISPR-derived RNAs from these loci. It has been reported that CRISPR cassettes expressing crRNAs provide acquired immunity [28, 53, 54]. It should be noted that short regions of spacers 7 and 19 of CRISPR 1 (14 bp and 17 bp, respectively) are complementary to two coding regions of PAV1, a virus that infects P. abyssi (data not shown) . No other similarities were observed between CRISPR 1 and 4 spacer sequences and known P. abyssi mobile genetic elements [56, 57]. The silent CRISPR 2 and 3 loci differ from the expressed cassettes in several ways. First, the CRISPR 2 locus has divergent direct repeat and leader sequences. The CRIPSR 3 locus is atypical by its leader sequence, the reduced number of repeats, and the sequence degeneracy of its repeats, suggesting that it is a relic of an active CRISPR cassette as mentioned in [27, 28, 58]. In general, CRISPR cassettes are physically linked to a cohort of conserved cas or cmr genes, in varying orientation and order, that encode CRISPR-associated proteins (reviewed in [24, 27, 31]). It is interesting to note that the only cas gene (Cas6, PAB1613) in the vicinity of P. abyssi CRISPR cassettes is located upstream of CRISPR 3 (Figure 1A). In the P. abyssi genome, all the other genes encoding Cas/Cmr core proteins are grouped into an operon located far from any CRISPR locus (Figure 1A). Further studies will be required to elucidate the mode of action and the dynamics of the CRISPR/Cas system in P. abyssi.
5' UTR-derived ncRNAs
Distinct RNA species are not always indicative of independently synthesized RNAs. As revealed by RACE experiments and analysis of genomic context, several ncRNAs (Table 1) seem to derive from mRNA leaders. This class of ncRNAs appears to originate from maturation of longer transcripts, as we demonstrated for sRkB and propose for sRk28, sRk33 and sRk61. Stable RNAs derived from RNA leaders were initially identified in E. coli and some appear to correspond to riboswitch elements such as RFN and THI . It is now well established that 5'UTRs can encompass transcription-termination signals involving riboswitch structures [60, 61]. In Archaea, it has recently been suggested that riboswitches may also exist. Based on comparative genomics, crcB has been found in Archaea and Bacteria, and TPP and flp in Euryarchaea [60, 62]. The crcB element was found within our set of 82 candidates (Additional file 2, Table S1), but ncRNA corresponding to this element was not detected by Northern blotting under the conditions employed in this study. To date, no biological evidence showing gene regulation by riboswitches has been reported in Archaea. Nevertheless, the conserved 70 nt sRk28 shares several feature with the bacterial preQ1 riboswitch [63, 64]. PAB 1992, adjacent to sRk28, is predicted to encode an ATPase with a queosine synthesis-like protein domain (QueC). The growth-regulated sRkB RNA is transcribed from a locus adjacent to PAB0571, which encodes a putative ATPase. In bacteria, the expression of several ATPase-like proteins is controlled by cis-mechanisms involving the SAM-I riboswitch. It is interesting to note in Listeria monocytogenes that this type of riboswitch also acts in trans as a small ncRNA to regulate expression of the virulence factor PrfA . Moreover, our data strongly suggest that sRkB has the propensity to form K-turn structures, which are found to be functional RNA motifs involved in the assembly of ribonucleoprotein complexes and in the orientation of pseudoknot interactions in the aptamer structures of the SAM-I and lysine riboswitches [36, 37]. Therefore, sRkB could have a specific function in the P. abyssi cell that involves the regulation of the expression of the putative ATPase encoded by PAB0571. We speculate that sRk28 and sRkB may be representative of a novel family of archaeal ncRNAs produced by transcription attenuation or RNA processing. Finally, the promoter-associated sRk33 and sRk61 could be related to the recently discovered ncRNAs associated with transcription starts sites reported in higher eukaryotes [65, 66], in yeast [67, 68] and in Salmonella enterica serovar typhi . This class of ncRNAs is thought to interfere with transcription of the downstream promoter by cis-mediated occlusion or by a trans-mechanism.
In this study, sRk48, sRk52, sRkB and sRkC ncRNAs, which are transcribed from repetitive loci, were shown to be related to a CDS annotated as an orphan OrfB element of the IS607 and IS200/605 family, . For example, in the case of sRk48 and sRk52, similar sequences correspond to 3' end coding regions annotated as partial OrfB-like IS607 transposase in P. furiosus and T. sibiricus. In bacteria, it has been shown that IS10 transposition events are controlled by a complementary cis-encoded antisense ncRNA , but this seems unrelated to our observations since the ncRNA sequences matched the sense sequence of the 3' end of the annotated ORFs. However, it is pertinent to note that a link between ncRNAs and transposase ORFs was already observed in Sulfolobus solfataricus[16, 17] and in Salmonella . Moreover, in Sulfolobus solfataricus it was shown that several small RNAs linked to annotated transposase ORFs could bind the multifunctional ribosomal protein L7Ae through recognition of an RNA kink turn as is the case for sRkB, sRkC and sRk52.
Two novel P. abyssi ncRNA families, the CRIPSR and the 5' UTR-related ncRNAs are described in this study. We certainly missed other P. abyssi ncRNA families because the biased composition screen only identifies GC rich ncRNAs. This is the case for the box C/D guide RNAs with low GC content. Nevertheless, the high number of known ncRNAs including the box H/ACA RNAs found by this approach confirms its relevance in searching for ncRNAs in AT-rich genomes. A limitation to this search, which is a bias due to misannotated ORFs and therefore intergenic regions, could not be excluded since sequence signals for archaeal transcription and translation are not as well defined as their counterparts in Bacteria. The discovery of novel ncRNAs by our combined computational approaches emphasizes the potential diversity of ncRNAs in Archaea, which could be enlarged by global RNomic approaches such as RNAseq technology. This would provide a deeper insight into the P. abyssi ncRNA world and help improve our knowledge of the specific roles of P. abyssi ncRNAs and their relationship to the 5'UTRs described in this study.
Genome sequences and related annotation files (gbk, .fna and .ptt extensions) were downloaded from the NCBI ftp database for P. abyssi, P. furiosus, P. horikoshii and T. kodakaraensis genomes. For these four genomes, a comparative analysis of all intergenic sequences was realized using RNAsim as described in . RNAsim searches for conserved sequences and structured regions between different genomes by using Wu-blast 2.0 to select pairwise alignments including conserved sequences (here with W = 7, E < 0.0001) and QRNA  to identify in pairwise alignments base substitution patterns that could correspond to a conserved RNA secondary structure. In a final step, RNAsim combines this information to predict loci that are conserved in multiple genomes. In this analysis, all the regions encoding annotated ncRNAs except for the box H/ACA RNA genes were excluded from the set of intergenic regions (46 tRNAs, 1 rRNA operon, 2 5S rRNAs, 59 C/D guide RNAs, 1 RNase P RNA and 1 7S/SRP RNA). GC-rich regions were predicted in P. abyssi by using the same segmentation approach as previously used by Klein and colleagues to predict new ncRNAs in P. furiosus. However, our approach differed in that the transition probabilities were adjusted to enrich the number of candidate regions. Additional sequence comparisons were performed on putative ncRNA candidates with Blastn (default values and W = 7) to add to this initial set other P. abyssi and archaeal homolog sequences. BRE/TATA, consensus boxes 5'_RNNANNTTTAWA_3' and 5'_RAAANNTWWWWA_3', and TATA consensus box 5'_TTWWWWA_3', K-turn and K-loop consensus motifs, inferred from the literature [34, 38, 39, 73], were identified using Patscan . Regions of favorable free energy were computed by setting the free energy threshold such that all tRNAs were found in sliding windows of 150 bp (here threshold=-32.3). Only regions of free energy below the threshold and longer than 50 bp were displayed in ApolloRNA, an extension of the annotation environment Apollo , developed to support ncRNA analysis. Highly structured hairpins with minimal hairpin stems of 6 bp (including G-U and U-G pairs) were searched with Patscan. RNA secondary structures were proposed on the basis of multiple alignments of similar regions within the Pyrococcales using Multalin  that were improved manually by combining RNAfold predictions and covariations. All data including predictions, motifs and annotations were integrated and visualized in ApolloRNA.
Oligonucleotides used in this work
Additional file 7 Table S2 list the primers used for Northern blot detection, CR-RT-PCR, primer extension, and the preparation of transcription templates by PCR.
Preparation of total cellular RNA
P. abyssi strain GE5 cells were grown as described in  at 95°C in Vent Sulfothermophiles Medium. P. abyssi cells were stopped at three different growth phases: exponential, end of exponential and stationary phases. The exponential and stationary phase cell paste was provided by A. Hecker (Institut de Génétique et de Microbiologie, Paris Sud-Orsay). Entry into stationary phase cell paste was purchased (Reims University, France). Total RNA was prepared from P. abyssi cell paste by Trizol extraction followed by treatments with DNase RQ1 (RNA-free, Promega), proteinase K and phenol/chloroform extraction followed by ethanol precipitation.
Northern blotting analysis
Total RNA (10 μg)extracted from cells in exponential growth phase (E), entry into stationary phase (ES) or stationary phase (S), and a 5' [32P]-end-labeled denatured PhiX174/HinfI DNA ladder were separated on a denaturing 6% polyacrylamide gel (8M urea, 1 × TBE buffer). Gels were transferred onto Hybond N+ nylon membrane using a Transphor Power Lid (Hoefer) apparatus in 0.5 × TBE buffer. The RNAs were UV cross-linked to the membrane (1200 J/cm²). Prehybridization was carried out for 30 minutes at 50°C in hybridization buffer (5 × SSC, 1 × Denhartd's solution, 1% SDS, 0.05 mg/mL sperm DNA herring). DNA oligonucleotides were designed using Primer designer or Vector NTI and 5' end labeled with [γ-32P] ATP and polynucleotide kinase. Hybridizations were carried out at 50°C for 16 h followed by two washes in 0.1 × SSC, 0.1% SDS buffer at room temperature for 10 min. The blots were analyzed by phosphorimaging (Molecular Dynamics) or autoradiography using MR or MS film (Kodak).
Primer extension and Circular RACE (CR-RT-PCR)
Primer extension analysis was performed using 10 μg of total RNA prepared from P. abyssi cells in entry into stationary phase. Total RNA was reverse transcribed at 42°C by AMV reverse transcriptase (Promega) using a 5' end labeled specific primer. CR-RT-PCR was performed with 20 μg of total RNA prepared from P. abyssi cells in entry into stationary phase, treated with (+) or without (-) 25U of Tobacco Acid Pyrophosphatase (TAP) according to manufacturer's protocol (Epicentre Biotechnologies). RNAs were extracted with phenol/chloroform then precipitated with ethanol. RNA (1 μg) +/- TAP was circularized with 40U of T4 RNA ligase (New England Biolabs), extracted with phenol/chloroform, ethanol precipitated and reverse transcribed by AMV reverse transcriptase using specific primers. After ethanol precipitation, the reverse transcripts were PCR amplified using Crimson Taq (New England Biolabs) and appropriate primers. The products were separated on a 6% native polyacrylamide gel (1% glycerol, 0.5 × TBE), treated with TAP, cloned in pCR®II-TOPO® with TOPO-TA Cloning® Kit according to the manufacturer's instructions (Invitrogen) and sequenced by Eurofins MWG Operon. About 100 sequences were analyzed for each CR-RT-PCR experiment.
In vitro transcription
A portion of the intergenic region corresponding to sRk52, sRkB and sRkC, respectively, was amplified by PCR from P. abyssi genomic DNA using specific primers (Additional file 2, Table S1). PCR products served as templates for in vitro transcription using T7 RNA polymerase as previously described .
Electrophoretic mobility shift assay (EMSA) and enzymatic structural probing
EMSA was performed using E.coli tRNA as a non-specific competitor as previously described . RNA and RNP complexes were separated on a native 6% (19:1) polyacrylamide gel containing 0.5 × TBE and 5% glycerol. Electrophoresis was performed at room temperature at 250 V in 0.5 × TBE running buffer containing 5% glycerol. The gels were dried and visualized using a Fuji-Bas 1000 phosphorImager.
We thank A. Hecker and P. Forterre (UMR CNRS 8621, Orsay, France) for kindly providing P. abyssi exponential and stationary cell paste, members of the Carpousis group for helpful discussions and critical comments on the manuscript, and M.J. Cros for ApolloRNA support. Our research is supported by the Centre National de la Recherche Scientifique (CNRS) with additional funding from the Agence Nationale de la Recherche (ANR) (grant BLAN08-1_329396). KP was supported by a predoctoral fellowship from the Ministère de l'Enseignement Supérieur et de la Recherche (MESR).
- Waters LS, Storz G: Regulatory RNAs in bacteria. Cell. 2009, 136: 615-628. 10.1016/j.cell.2009.01.043.PubMed CentralPubMedView Article
- Geissmann T, Marzi S, Romby P: The role of mRNA structure in translational control in bacteria. RNA Biol. 2009, 6: 153-160. 10.4161/rna.6.2.8047.PubMedView Article
- Ghildiyal M, Zamore PD: Small silencing RNAs: an expanding universe. Nat Rev Genet. 2009, 10: 94-108. 10.1038/nrg2504.PubMed CentralPubMedView Article
- Vogel J, Wagner EG: Target identification of small noncoding RNAs in bacteria. Curr Opin Microbiol. 2007, 10: 262-270. 10.1016/j.mib.2007.06.001.PubMedView Article
- Backofen R, Hess WR: Computational prediction of sRNAs and their targets in bacteria. RNA Biol. 2010, 7: 33-42. 10.4161/rna.7.1.10655.PubMedView Article
- Majdalani N, Vanderpool CK, Gottesman S: Bacterial small RNA regulators. Crit Rev Biochem Mol Biol. 2005, 40: 93-113. 10.1080/10409230590918702.PubMedView Article
- Loh E, Dussurget O, Gripenland J, Vaitkevicius K, Tiensuu T, Mandin P, Repoila F, Buchrieser C, Cossart P, Johansson J: A trans-acting riboswitch controls expression of the virulence regulator PrfA in Listeria monocytogenes. Cell. 2009, 139: 770-779. 10.1016/j.cell.2009.08.046.PubMedView Article
- Papenfort K, Vogel J: Regulatory RNA in bacterial pathogens. Cell Host Microbe. 2010, 8: 116-127. 10.1016/j.chom.2010.06.008.PubMedView Article
- Roth A, Breaker RR: The structural and functional diversity of metabolite-binding riboswitches. Annu Rev Biochem. 2009, 78: 305-334. 10.1146/annurev.biochem.78.070507.135656.PubMedView Article
- Klinkert B, Narberhaus F: Microbial thermosensors. Cell Mol Life Sci. 2009, 66: 2661-2676. 10.1007/s00018-009-0041-3.PubMedView Article
- Omer AD, Lowe TM, Russell AG, Ebhardt H, Eddy SR, Dennis PP: Homologs of small nucleolar RNAs in Archaea. Science. 2000, 288: 517-522. 10.1126/science.288.5465.517.PubMedView Article
- Gaspin C, Cavaille J, Erauso G, Bachellerie JP: Archaeal homologs of eukaryotic methylation guide small nucleolar RNAs: lessons from the Pyrococcus genomes. J Mol Biol. 2000, 297: 895-906. 10.1006/jmbi.2000.3593.PubMedView Article
- Muller S, Leclerc F, Behm-Ansmant I, Fourmann JB, Charpentier B, Branlant C: Combined in silico and experimental identification of the Pyrococcus abyssi H/ACA sRNAs and their target sites in ribosomal RNAs. Nucleic Acids Res. 2008, 36: 2459-2475. 10.1093/nar/gkn077.PubMed CentralPubMedView Article
- Grosjean H, Gaspin C, Marck C, Decatur WA, de Crecy-Lagard V: RNomics and Modomics in the halophilic archaea Haloferax volcanii: identification of RNA modification genes. BMC Genomics. 2008, 9: 470-10.1186/1471-2164-9-470.PubMed CentralPubMedView Article
- Tang TH, Bachellerie JP, Rozhdestvensky T, Bortolin ML, Huber H, Drungowski M, Elge T, Brosius J, Huttenhofer A: Identification of 86 candidates for small non-messenger RNAs from the archaeon Archaeoglobus fulgidus. Proc Natl Acad Sci USA. 2002, 99: 7536-7541. 10.1073/pnas.112047299.PubMed CentralPubMedView Article
- Tang TH, Polacek N, Zywicki M, Huber H, Brugger K, Garrett R, Bachellerie JP, Huttenhofer A: Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus. Mol Microbiol. 2005, 55: 469-481.PubMedView Article
- Zago MA, Dennis PP, Omer AD: The expanding world of small RNAs in the hyperthermophilic archaeon Sulfolobus solfataricus. Mol Microbiol. 2005, 55: 1812-1828. 10.1111/j.1365-2958.2005.04505.x.PubMedView Article
- Schattner P: Searching for RNA genes using base-composition statistics. Nucleic Acids Res. 2002, 30: 2076-2082. 10.1093/nar/30.9.2076.PubMed CentralPubMedView Article
- Klein RJ, Misulovin Z, Eddy SR: Noncoding RNA genes identified in AT-rich hyperthermophiles. Proc Natl Acad Sci USA. 2002, 99: 7542-7547. 10.1073/pnas.112063799.PubMed CentralPubMedView Article
- Straub J, Brenneis M, Jellen-Ritter A, Heyer R, Soppa J, Marchfelder A: Small RNAs in haloarchaea: identification, differential expression and biological function. RNA Biol. 2009, 6: 281-292. 10.4161/rna.6.3.8357.PubMedView Article
- Soppa J, Straub J, Brenneis M, Jellen-Ritter A, Heyer R, Fischer S, Granzow M, Voss B, Hess WR, Tjaden B, Marchfelder A: Small RNAs of the halophilic archaeon Haloferax volcanii. Biochem Soc Trans. 2009, 37: 133-136. 10.1042/BST0370133.PubMedView Article
- Jager D, Sharma CM, Thomsen J, Ehlers C, Vogel J, Schmitz RA: Deep sequencing analysis of the Methanosarcina mazei Go1 transcriptome in response to nitrogen availability. Proc Natl Acad Sci USA. 2009, 106: 21878-21882. 10.1073/pnas.0909051106.PubMed CentralPubMedView Article
- Lillestol RK, Redder P, Garrett RA, Brugger K: A putative viral defence mechanism in archaeal cells. Archaea. 2006, 2: 59-72. 10.1155/2006/542818.PubMed CentralPubMedView Article
- Hale C, Kleppe K, Terns RM, Terns MP: Prokaryotic silencing (psi)RNAs in Pyrococcus furiosus. RNA. 2008, 14: 2572-2579. 10.1261/rna.1246808.PubMed CentralPubMedView Article
- Makarova KS, Grishin NV, Shabalina SA, Wolf YI, Koonin EV: A putative RNA-interference-based immune system in prokaryotes: computational analysis of the predicted enzymatic machinery, functional analogies with eukaryotic RNAi, and hypothetical mechanisms of action. Biol Direct. 2006, 1: 7-10.1186/1745-6150-1-7.PubMed CentralPubMedView Article
- Kunin V, Sorek R, Hugenholtz P: Evolutionary conservation of sequence and secondary structures in CRISPR repeats. Genome Biol. 2007, 8: R61-10.1186/gb-2007-8-4-r61.PubMed CentralPubMedView Article
- Deveau H, Garneau JE, Moineau S: CRISPR/Cas System and Its Role in Phage-Bacteria Interactions. Annu Rev Microbiol. 2010
- Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, Romero DA, Horvath P: CRISPR provides acquired resistance against viruses in prokaryotes. Science. 2007, 315: 1709-1712. 10.1126/science.1138140.PubMedView Article
- Horvath P, Barrangou R: CRISPR/Cas, the immune system of bacteria and archaea. Science. 2010, 327: 167-170. 10.1126/science.1179555.PubMedView Article
- Mojica FJ, Diez-Villasenor C, Garcia-Martinez J, Soria E: Intervening sequences of regularly spaced prokaryotic repeats derive from foreign genetic elements. J Mol Evol. 2005, 60: 174-182. 10.1007/s00239-004-0046-3.PubMedView Article
- Karginov FV, Hannon GJ: The CRISPR system: small RNA-guided defense in bacteria and archaea. Mol Cell. 2010, 37: 7-19. 10.1016/j.molcel.2009.12.033.PubMed CentralPubMedView Article
- Rozhdestvensky TS, Tang TH, Tchirkova IV, Brosius J, Bachellerie JP, Huttenhofer A: Binding of L7Ae protein to the K-turn of archaeal snoRNAs: a shared RNA binding motif for C/D and H/ACA box snoRNAs in Archaea. Nucleic Acids Res. 2003, 31: 869-877. 10.1093/nar/gkg175.PubMed CentralPubMedView Article
- Cho IM, Lai LB, Susanti D, Mukhopadhyay B, Gopalan V: Ribosomal protein L7Ae is a subunit of archaeal RNase P. Proc Natl Acad Sci USA. 2010, 107: 14573-14578. 10.1073/pnas.1005556107.PubMed CentralPubMedView Article
- Nolivos S, Carpousis AJ, Clouet-d'Orval B: The K-loop, a general feature of the Pyrococcus C/D guide RNAs, is an RNA structural motif related to the K-turn. Nucleic Acids Res. 2005, 33: 6507-6514. 10.1093/nar/gki962.PubMed CentralPubMedView Article
- Schroeder KT, McPhee SA, Ouellet J, Lilley DM: A structural database for k-turn motifs in RNA. RNA. 2010, 16: 1463-1468. 10.1261/rna.2207910.PubMed CentralPubMedView Article
- Montange RK, Batey RT: Structure of the S-adenosylmethionine riboswitch regulatory mRNA element. Nature. 2006, 441: 1172-1175. 10.1038/nature04819.PubMedView Article
- Blouin S, Lafontaine DA: A loop loop interaction and a K-turn motif located in the lysine aptamer domain are important for the riboswitch gene regulation control. RNA. 2007, 13: 1256-1267. 10.1261/rna.560307.PubMed CentralPubMedView Article
- Cohen GN, Barbe V, Flament D, Galperin M, Heilig R, Lecompte O, Poch O, Prieur D, Querellou J, Ripp R, et al: An integrated analysis of the genome of the hyperthermophilic archaeon Pyrococcus abyssi. Mol Microbiol. 2003, 47: 1495-1512. 10.1046/j.1365-2958.2003.03381.x.PubMedView Article
- Leontis NB, Lescoute A, Westhof E: The building blocks and motifs of RNA architecture. Curr Opin Struct Biol. 2006, 16: 279-287. 10.1016/j.sbi.2006.05.009.PubMedView Article
- Portillo MC, Gonzalez JM: CRISPR elements in the Thermococcales: evidence for associated horizontal gene transfer in Pyrococcus furiosus. J Appl Genet. 2009, 50: 421-430. 10.1007/BF03195703.PubMedView Article
- Schneider KL, Pollard KS, Baertsch R, Pohl A, Lowe TM: The UCSC Archaeal Genome Browser. Nucleic Acids Res. 2006, 34: D407-410. 10.1093/nar/gkj134.PubMed CentralPubMedView Article
- Carte J, Wang R, Li H, Terns RM, Terns MP: Cas6 is an endoribonuclease that generates guide RNAs for invader defense in prokaryotes. Genes Dev. 2008, 22: 3489-3496. 10.1101/gad.1742908.PubMed CentralPubMedView Article
- Filee J, Siguier P, Chandler M: Insertion sequence diversity in archaea. Microbiol Mol Biol Rev. 2007, 71: 121-157. 10.1128/MMBR.00031-06.PubMed CentralPubMedView Article
- Hobbs EC, Fontaine F, Yin X, Storz G: An expanding universe of small proteins. Curr Opin Microbiol. 2011, 14: 167-173. 10.1016/j.mib.2011.01.007.PubMed CentralPubMedView Article
- Wurtzel O, Sapra R, Chen F, Zhu Y, Simmons BA, Sorek R: A single-base resolution map of an archaeal transcriptome. Genome Res. 2010, 20: 133-141. 10.1101/gr.100396.109.PubMed CentralPubMedView Article
- Ramos CR, Oliveira CL, Torriani IL, Oliveira CC: The Pyrococcus exosome complex: structural and functional characterization. J Biol Chem. 2006, 281: 6751-6759. 10.1074/jbc.M512495200.PubMedView Article
- Tsai HY, Pulukkunat DK, Woznick WK, Gopalan V: Functional reconstitution and characterization of Pyrococcus furiosus RNase P. Proc Natl Acad Sci USA. 2006, 103: 16147-16152. 10.1073/pnas.0608000103.PubMed CentralPubMedView Article
- Evguenieva-Hackenberg E, Klug G: RNA degradation in Archaea and Gram-negative bacteria different from Escherichia coli. Prog Mol Biol Transl Sci. 2009, 85: 275-317.PubMedView Article
- Clouet-d'Orval B, Rinaldi D, Quentin Y, Carpousis AJ: Euryarchaeal beta-CASP proteins with homology to bacterial RNase J Have 5'- to 3'-exoribonuclease activity. J Biol Chem. 2010, 285: 17574-17583. 10.1074/jbc.M109.095117.PubMed CentralPubMedView Article
- Hasenohrl D, Konrat R, Blasi U: Identification of an RNase J ortholog in Sulfolobus solfataricus: implications for 5'-to-3' directional decay and 5'-end protection of mRNA in Crenarchaeota. RNA. 2011, 17: 99-107. 10.1261/rna.2418211.PubMed CentralPubMedView Article
- Hale CR, Zhao P, Olson S, Duff MO, Graveley BR, Wells L, Terns RM, Terns MP: RNA-guided RNA cleavage by a CRISPR RNA-Cas protein complex. Cell. 2009, 139: 945-956. 10.1016/j.cell.2009.07.040.PubMed CentralPubMedView Article
- Lillestol RK, Shah SA, Brugger K, Redder P, Phan H, Christiansen J, Garrett RA: CRISPR families of the crenarchaeal genus Sulfolobus: bidirectional transcription and dynamic properties. Mol Microbiol. 2009, 72: 259-272. 10.1111/j.1365-2958.2009.06641.x.PubMedView Article
- Marraffini LA, Sontheimer EJ: CRISPR interference limits horizontal gene transfer in staphylococci by targeting DNA. Science. 2008, 322: 1843-1845. 10.1126/science.1165771.PubMed CentralPubMedView Article
- Brouns SJ, Jore MM, Lundgren M, Westra ER, Slijkhuis RJ, Snijders AP, Dickman MJ, Makarova KS, Koonin EV, van der Oost J: Small CRISPR RNAs guide antiviral defense in prokaryotes. Science. 2008, 321: 960-964. 10.1126/science.1159689.PubMedView Article
- Geslin C, Gaillard M, Flament D, Rouault K, Le Romancer M, Prieur D, Erauso G: Analysis of the first genome of a hyperthermophilic marine virus-like particle, PAV1, isolated from Pyrococcus abyssi. J Bacteriol. 2007, 189: 4510-4519. 10.1128/JB.01896-06.PubMed CentralPubMedView Article
- Erauso G, Marsin S, Benbouzid-Rollet N, Baucher MF, Barbeyron T, Zivanovic Y, Prieur D, Forterre P: Sequence of plasmid pGT5 from the archaeon Pyrococcus abyssi: evidence for rolling-circle replication in a hyperthermophile. J Bacteriol. 1996, 178: 3232-3237.PubMed CentralPubMed
- Soler N, Marguet E, Cortez D, Desnoues N, Keller J, van Tilbeurgh H, Sezonov G, Forterre P: Two novel families of plasmids from hyperthermophilic archaea encoding new families of replication proteins. Nucleic Acids Res. 2010, 38: 5088-5104. 10.1093/nar/gkq236.PubMed CentralPubMedView Article
- Stern A, Keren L, Wurtzel O, Amitai G, Sorek R: Self-targeting by CRISPR: gene regulation or autoimmunity?. Trends Genet. 2010, 26: 335-340. 10.1016/j.tig.2010.05.008.PubMed CentralPubMedView Article
- Vogel J, Bartels V, Tang TH, Churakov G, Slagter-Jager JG, Huttenhofer A, Wagner EG: RNomics in Escherichia coli detects new sRNA species and indicates parallel transcriptional output in bacteria. Nucleic Acids Res. 2003, 31: 6435-6443. 10.1093/nar/gkg867.PubMed CentralPubMedView Article
- Barrick JE, Breaker RR: The distributions, mechanisms, and structures of metabolite-binding riboswitches. Genome Biol. 2007, 8: R239-10.1186/gb-2007-8-11-r239.PubMed CentralPubMedView Article
- Blouin S, Mulhbacher J, Penedo JC, Lafontaine DA: Riboswitches: ancient and promising genetic regulators. Chembiochem. 2009, 10: 400-416. 10.1002/cbic.200800593.PubMedView Article
- Weinberg Z, Wang JX, Bogue J, Yang J, Corbino K, Moy RH, Breaker RR: Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea, and their metagenomes. Genome Biol. 2010, 11: R31-10.1186/gb-2010-11-3-r31.PubMed CentralPubMedView Article
- Roth A, Winkler WC, Regulski EE, Lee BW, Lim J, Jona I, Barrick JE, Ritwik A, Kim JN, Welz R, et al: A riboswitch selective for the queuosine precursor preQ1 contains an unusually small aptamer domain. Nat Struct Mol Biol. 2007, 14: 308-317. 10.1038/nsmb1224.PubMedView Article
- Kang M, Peterson R, Feigon J: Structural Insights into riboswitch control of the biosynthesis of queuosine, a modified nucleotide found in the anticodon of tRNA. Mol Cell. 2009, 33: 784-790. 10.1016/j.molcel.2009.02.019.PubMedView Article
- Kapranov P, Cheng J, Dike S, Nix DA, Duttagupta R, Willingham AT, Stadler PF, Hertel J, Hackermuller J, Hofacker IL, et al: RNA maps reveal new RNA classes and a possible function for pervasive transcription. Science. 2007, 316: 1484-1488. 10.1126/science.1138341.PubMedView Article
- Taft RJ, Glazov EA, Cloonan N, Simons C, Stephen S, Faulkner GJ, Lassmann T, Forrest AR, Grimmond SM, Schroder K, et al: Tiny RNAs associated with transcription start sites in animals. Nat Genet. 2009, 41: 572-578. 10.1038/ng.312.PubMedView Article
- Hirota K, Miyoshi T, Kugou K, Hoffman CS, Shibata T, Ohta K: Stepwise chromatin remodelling by a cascade of transcription initiation of non-coding RNAs. Nature. 2008, 456: 130-134. 10.1038/nature07348.PubMedView Article
- Martens JA, Laprade L, Winston F: Intergenic transcription is required to repress the Saccharomyces cerevisiae SER3 gene. Nature. 2004, 429: 571-574. 10.1038/nature02538.PubMedView Article
- Chinni SV, Raabe CA, Zakaria R, Randau G, Hoe CH, Zemann A, Brosius J, Tang TH, Rozhdestvensky TS: Experimental identification and characterization of 97 novel npcRNA candidates in Salmonella enterica serovar Typhi. Nucleic Acids Res. 2010, 38: 5893-5908. 10.1093/nar/gkq281.PubMed CentralPubMedView Article
- Ma C, Simons RW: The IS10 antisense RNA blocks ribosome binding at the transposase translation initiation site. EMBO J. 1990, 9: 1267-1274.PubMed CentralPubMed
- Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, Binnewies TT, Hinton JC, Vogel J: Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet. 2008, 4: e1000163-10.1371/journal.pgen.1000163.PubMed CentralPubMedView Article
- Rivas E, Eddy SR: Noncoding RNA gene detection using comparative sequence analysis. BMC Bioinformatics. 2001, 2: 8-10.1186/1471-2105-2-8.PubMed CentralPubMedView Article
- Littlefield O, Korkhin Y, Sigler PB: The structural basis for the oriented assembly of a TBP/TFB/promoter complex. Proc Natl Acad Sci USA. 1999, 96: 13668-13673. 10.1073/pnas.96.24.13668.PubMed CentralPubMedView Article
- Dsouza M, Larsen N, Overbeek R: Searching for patterns in genomic data. Trends Genet. 1997, 13: 497-498.PubMedView Article
- Lewis SE, Searle SM, Harris N, Gibson M, Lyer V, Richter J, Wiel C, Bayraktaroglir L, Birney E, Crosby MA, et al: Apollo: a sequence annotation editor. Genome Biol. 2002, 3: RESEARCH0082-PubMed CentralPubMedView Article
- Corpet F: Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 1988, 16: 10881-10890. 10.1093/nar/16.22.10881.PubMed CentralPubMedView Article
- Charbonnier F, Erauso G, Barbeyron T, Prieur D, Forterre P: Evidence that a plasmid from a hyperthermophilic archaebacterium is relaxed at physiological temperatures. J Bacteriol. 1992, 174: 6103-6108.PubMed CentralPubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.