RNA-Seq analyses reveal the order of tRNA processing events and the maturation of C/D box and CRISPR RNAs in the hyperthermophile Methanopyrus kandleri

Genetic and Functional Analyses of Archaeal ATP-Dependent RNA Ligase in C/D Box sRNA Circularization and Ribosomal RNA Processing

RNA ligases play important roles in repairing and circularizing RNAs post-transcriptionally. In this study, we generated an allelic knockout of ATP-dependent RNA ligase (Rnl) in the hyperthermophilic archaeon Thermococcus kodakarensis to identify its biological targets. A comparative analysis of circular RNA reveals that the Rnl-knockout strain represses circularization of C/D box sRNAs without affecting the circularization of tRNA and rRNA processing intermediates. Recombinant archaeal Rnl could circularize C/D box sRNAs with a mutation in the conserved C/D box sequence element but not when the terminal stem structures were disrupted, suggesting that proximity of the two ends could be critical for intramolecular ligation. Furthermore, T. kodakarensis accumulates aberrant RNA fragments derived from ribosomal RNA in the absence of Rnl. These results suggest that Rnl is responsible for C/D box sRNA circularization and may also play a role in ribosomal RNA processing.

Conservation of Archaeal C/D Box sRNA-Guided RNA Modifications

Post-transcriptional modifications fulfill many important roles during ribosomal RNA maturation in all three domains of life. Ribose 2'-O-methylations constitute the most abundant chemical rRNA modification and are, for example, involved in RNA folding and stabilization. In archaea, these modification sites are determined by variable sets of C/D box sRNAs that guide the activity of the rRNA 2'-O-methyltransferase fibrillarin. Each C/D box sRNA contains two guide sequences that can act in coordination to bridge rRNA sequences. Here, we will review the landscape of archaeal C/D box sRNA genes and their target sites. One focus is placed on the apparent accelerated evolution of guide sequences and the varied pairing of the two individual guides, which results in different rRNA modification patterns and RNA chaperone activities.

Heavily Armed Ancestors: CRISPR Immunity and Applications in Archaea with a Comparative Analysis of CRISPR Types in Sulfolobales

Prokaryotes are constantly coping with attacks by viruses in their natural environments and therefore have evolved an impressive array of defense systems. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) is an adaptive immune system found in the majority of archaea and about half of bacteria which stores pieces of infecting viral DNA as spacers in genomic CRISPR arrays to reuse them for specific virus destruction upon a second wave of infection. In detail, small CRISPR RNAs (crRNAs) are transcribed from CRISPR arrays and incorporated into type-specific CRISPR effector complexes which further degrade foreign nucleic acids complementary to the crRNA. This review gives an overview of CRISPR immunity to newcomers in the field and an update on CRISPR literature in archaea by comparing the functional mechanisms and abundances of the diverse CRISPR types. A bigger fraction is dedicated to the versatile and prevalent CRISPR type III systems, as tremendous progress has been made recently using archaeal models in discerning the controlled molecular mechanisms of their unique tripartite mode of action including RNA interference, DNA interference and the unique cyclic-oligoadenylate signaling that induces promiscuous RNA shredding by CARF-domain ribonucleases. The second half of the review spotlights CRISPR in archaea outlining seminal in vivo and in vitro studies in model organisms of the euryarchaeal and crenarchaeal phyla, including the application of CRISPR-Cas for genome editing and gene silencing. In the last section, a special focus is laid on members of the crenarchaeal hyperthermophilic order Sulfolobales by presenting a thorough comparative analysis about the distribution and abundance of CRISPR-Cas systems, including arrays and spacers as well as CRISPR-accessory proteins in all 53 genomes available to date. Interestingly, we find that CRISPR type III and the DNA-degrading CRISPR type I complexes co-exist in more than two thirds of these genomes. Furthermore, we identified ring nuclease candidates in all but two genomes and found that they generally co-exist with the above-mentioned CARF domain ribonucleases Csx1/Csm6. These observations, together with published literature allowed us to draft a working model of how CRISPR-Cas systems and accessory proteins cross talk to establish native CRISPR anti-virus immunity in a Sulfolobales cell.

Recent Development of Extremophilic Bacteria and Their Application in Biorefinery

The biorefining technology for biofuels and chemicals from lignocellulosic biomass has made great progress in the world. However, mobilization of laboratory research toward industrial setup needs to meet a series of criteria, including the selection of appropriate pretreatment technology, breakthrough in enzyme screening, pathway optimization, and production technology, etc. Extremophiles play an important role in biorefinery by providing novel metabolic pathways and catalytically stable/robust enzymes that are able to act as biocatalysts under harsh industrial conditions on their own. This review summarizes the potential application of thermophilic, psychrophilic alkaliphilic, acidophilic, and halophilic bacteria and extremozymes in the pretreatment, saccharification, fermentation, and lignin valorization process. Besides, the latest studies on the engineering bacteria of extremophiles using metabolic engineering and synthetic biology technologies for high-efficiency biofuel production are also introduced. Furthermore, this review explores the comprehensive application potential of extremophiles and extremozymes in biorefinery, which is partly due to their specificity and efficiency, and points out the necessity of accelerating the commercialization of extremozymes.

Approaches to study CRISPR RNA biogenesis and the key players involved

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins provide an inheritable and adaptive immune system against phages and foreign genetic elements in many bacteria and archaea. The three stages of CRISPR-Cas immunity comprise adaptation, CRISPR RNA (crRNA) biogenesis and interference. The maturation of the pre-crRNA into mature crRNAs, short guide RNAs that target invading nucleic acids, is crucial for the functionality of CRISPR-Cas defense systems. Mature crRNAs assemble with Cas proteins into the ribonucleoprotein (RNP) effector complex and guide the Cas nucleases to the cognate foreign DNA or RNA target. Experimental approaches to characterize these crRNAs, the specific steps toward their maturation and the involved factors, include RNA-seq analyses, enzyme assays, methods such as cryo-electron microscopy, the crystallization of proteins, or UV-induced protein-RNA crosslinking coupled to mass spectrometry analysis. Complex and multiple interactions exist between CRISPR-cas-encoded specific riboendonucleases such as Cas6, Cas5d and Csf5, endonucleases with dual functions in maturation and interference such as the enzymes of the Cas12 and Cas13 families, and nucleases belonging to the cell's degradosome such as RNase E, PNPase and RNase J, both in the maturation as well as in interference. The results of these studies have yielded a picture of unprecedented diversity of sequences, enzymes and biochemical mechanisms.

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

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

Next Generation DNA-Seq and Differential RNA-Seq Allow Re-annotation of the Pyrococcus furiosus DSM 3638 Genome and Provide Insights Into Archaeal Antisense Transcription

Pyrococcus furiosus DSM 3638 is a model organism for hyperthermophilic archaea with an optimal growth temperature near 100°C. The genome was sequenced about 18 years ago. However, some publications suggest that in contrast to other Pyrococcus species, the genome of P. furiosus DSM 3638 is prone to genomic rearrangements. Therefore, we re-sequenced the genome using third generation sequencing techniques. The new de novo assembled genome is 1,889,914 bp in size and exhibits high sequence identity to the published sequence. However, two major deviations were detected: (1) The genome is 18,342 bp smaller than the NCBI reference genome due to a recently described deletion. (2) The region between PF0349 and PF0388 is inverted most likely due an assembly problem for the original sequence. In addition, numerous minor variations, ranging from single nucleotide exchanges, deletions or insertions were identified. The total number of insertion sequence (IS) elements is also reduced from 30 to 24 in the new sequence. Re-sequencing of a 2-year-old “lab culture” using Nanopore sequencing confirmed the overall stability of the P. furiosus DSM 3638 genome even under normal lab conditions without taking any special care. To improve genome annotation, the updated DNA sequence was combined with an RNA sequencing approach. Here, RNAs from eight different growth conditions were pooled to increase the number of detected transcripts. Furthermore, a differential RNA-Seq approach was employed for the identification of transcription start sites (TSSs). In total, 2515 TSSs were detected and classified into 834 primary (pTSS), 797 antisense (aTSS), 739 internal and 145 secondary TSSs. Our analysis of the upstream regions revealed a well conserved archaeal promoter structure. Interrogation of the distances between pTSSs and aTSSs revealed a significant number of antisense transcripts, which are a result of bidirectional transcription from the same TATA box. This mechanism of antisense transcript production could be further confirmed by in vitro transcription experiments. We assume that bidirectional transcription gives rise to non-functional antisense RNAs and that this is a widespread phenomenon in archaea due to the architecture of the TATA element and the symmetric structure of the TATA-binding protein.

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.

RNA ribose methylation (2'-O-methylation): Occurrence, biosynthesis and biological functions

Methylation of riboses at 2'-OH group is one of the most common RNA modifications found in number of cellular RNAs from almost any species which belong to all three life domains. This modification was extensively studied for decades in rRNAs and tRNAs, but recent data revealed the presence of 2'-O-methyl groups also in low abundant RNAs, like mRNAs. Ribose methylation is formed in RNA by two alternative enzymatic mechanisms: either by stand-alone protein enzymes or by complex assembly of proteins associated with snoRNA guides (sno(s)RNPs). In that case one catalytic subunit acts at various RNA sites, the specificity is provided by base pairing of the sno(s)RNA guide with the target RNA. In this review we compile available information on 2'-OH ribose methylation in different RNAs, enzymatic machineries involved in their biosynthesis and dynamics, as well as on the physiological functions of these modified residues.

Transfer RNA Modification Enzymes from Thermophiles and Their Modified Nucleosides in tRNA

To date, numerous modified nucleosides in tRNA as well as tRNA modification enzymes have been identified not only in thermophiles but also in mesophiles. Because most modified nucleosides in tRNA from thermophiles are common to those in tRNA from mesophiles, they are considered to work essentially in steps of protein synthesis at high temperatures. At high temperatures, the structure of unmodified tRNA will be disrupted. Therefore, thermophiles must possess strategies to stabilize tRNA structures. To this end, several thermophile-specific modified nucleosides in tRNA have been identified. Other factors such as RNA-binding proteins and polyamines contribute to the stability of tRNA at high temperatures. Thermus thermophilus, which is an extreme-thermophilic eubacterium, can adapt its protein synthesis system in response to temperature changes via the network of modified nucleosides in tRNA and tRNA modification enzymes. Notably, tRNA modification enzymes from thermophiles are very stable. Therefore, they have been utilized for biochemical and structural studies. In the future, thermostable tRNA modification enzymes may be useful as biotechnology tools and may be utilized for medical science.

The RNA-splicing endonuclease from the euryarchaeaon Methanopyrus kandleri is a heterotetramer with constrained substrate specificity

Four different types (α4, α'2, (αβ)2 and ϵ2) of RNA-splicing endonucleases (EndAs) for RNA processing are known to exist in the Archaea. Only the (αβ)2 and ϵ2 types can cleave non-canonical introns in precursor (pre)-tRNA. Both enzyme types possess an insert associated with a specific loop, allowing broad substrate specificity in the catalytic α units. Here, the hyperthermophilic euryarchaeon Methanopyrus kandleri (MKA) was predicted to harbor an (αβ)2-type EndA lacking the specific loop. To characterize MKA EndA enzymatic activity, we constructed a fusion protein derived from MKA α and β subunits (fMKA EndA). In vitro assessment demonstrated complete removal of the canonical bulge-helix-bulge (BHB) intron structure from MKA pre-tRNAAsn. However, removal of the relaxed BHB structure in MKA pre-tRNAGlu was inefficient compared to crenarchaeal (αβ)2 EndA, and the ability to process the relaxed intron within mini-helix RNA was not detected. fMKA EndA X-ray structure revealed a shape similar to that of other EndA types, with no specific loop. Mapping of EndA types and their specific loops and the tRNA gene diversity among various Archaea suggest that MKA EndA is evolutionarily related to other (αβ)2-type EndAs found in the Thaumarchaeota, Crenarchaeota and Aigarchaeota but uniquely represents constrained substrate specificity.

General response of Salmonella enterica serovar Typhimurium to desiccation: A new role for the virulence factors sopD and sseD in survival

Salmonella can survive for long periods under extreme desiccation conditions. This stress tolerance poses a risk for food safety, but relatively little is known about the molecular and cellular regulation of this adaptation mechanism. To determine the genetic components involved in Salmonella’s cellular response to desiccation, we performed a global transcriptomic analysis comparing S. enterica serovar Typhimurium cells equilibrated to low water activity (a_w 0.11) and cells equilibrated to high water activity (a_w 1.0). The analysis revealed that 719 genes were differentially regulated between the two conditions, of which 290 genes were up-regulated at a_w 0.11. Most of these genes were involved in metabolic pathways, transporter regulation, DNA replication/repair, transcription and translation, and, more importantly, virulence genes. Among these, we decided to focus on the role of sopD and sseD. Deletion mutants were created and their ability to survive desiccation and exposure to a_w 0.11 was compared to the wild-type strain and to an E. coli O157:H7 strain. The sopD and sseD mutants exhibited significant cell viability reductions of 2.5 and 1.3 Log (CFU/g), respectively, compared to the wild-type after desiccation for 4 days on glass beads. Additional viability differences of the mutants were observed after exposure to a_w 0.11 for 7 days. E. coli O157:H7 lost viability similarly to the mutants. Scanning electron microscopy showed that both mutants displayed a different morphology compared to the wild-type and differences in production of the extracellular matrix under the same conditions. These findings suggested that sopD and sseD are required for Salmonella’s survival during desiccation.

Biology and clinical relevance of noncoding sno/scaRNAs

Small nucleolar RNAs (snoRNAs) are a group of noncoding RNAs that perform various biological functions, including biochemical modifications of other RNAs, precursors of miRNA, splicing, and telomerase activity. The small Cajal body-associated RNAs (scaRNAs) are a subset of the snoRNA family and collect in the Cajal body where they perform their canonical function to biochemically modify spliceosomal RNAs prior to maturation. Failure of sno/scaRNAs have been implicated in pathology such as congenital heart anomalies, neuromuscular disorders, and various malignancies. Thus, understanding of sno/scaRNAs demonstrates the clinical value.

Using RNA sequencing to identify putative competing endogenous RNAs (ceRNAs) potentially regulating fat metabolism in bovine liver

RNA sequencing has been extensively used to study specific gene expression patterns to discover potential key genes related to complex traits of interest in animals. Of note, a new regulatory mechanism builds a large-scale regulatory network among transcriptome, where lncRNAs act as competing endogenous RNAs (ceRNAs) to sponge miRNAs to regulate the expression of miRNA target genes post-transcriptionally. In this study, we sequenced the cDNA and sRNA libraries of nine liver samples from three Holstein cows during dry period, early lactation, and peak of lactation with HiSeq platform. As a result, we identified 665 genes, 57 miRNAs and 33 lncRNAs that displayed differential expression patterns across periods. Subsequently, a total of 41ceRNA pairs (lncRNA-mRNA) sharing 11 miRNAs were constructed including 30 differentially expressed genes. Importantly, 12 among them were presented in our large metabolic networks, and predicted to influence the lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, thus, these genes were considered as the most promising candidates for milk fat formation. To our knowledge, this is first investigation to profile the ceRNA regulatory networks of liver transcriptome that could affect milk fat synthesis in bovine, providing a new view of the regulatory mechanism of RNAs.

Unsupervised detection of regulatory gene expression information in different genomic regions enables gene expression ranking

Background

The regulation of all gene expression steps (e.g., Transcription, RNA processing, Translation, and mRNA Degradation) is known to be primarily encoded in different parts of genes and in genomic regions in proximity to genes (e.g., promoters, untranslated regions, coding regions, introns, etc.). However, the entire gene expression codes and the genomic regions where they are encoded are still unknown.

Results

Here, we employ an unsupervised approach to estimate the concentration of gene expression codes in different non-coding parts of genes and transcripts, such as introns and untranslated regions, focusing on three model organisms (Escherichia coli, Saccharomyces cerevisiae, and Schizosaccharomyces pombe). Our analyses support the conjecture that regions adjacent to the beginning and end of ORFs and the beginning and end of introns tend to include higher concentration of gene expression information relatively to regions further away. In addition, we report the exact regions with elevated concentration of gene expression codes. Furthermore, we demonstrate that the concentration of these codes in different genetic regions is correlated with the expression levels of the corresponding genes, and with splicing efficiency measurements and meiotic stage gene expression measurements in S. cerevisiae.

Conclusion

We suggest that these discoveries improve our understanding of gene expression regulation and evolution; they can also be used for developing improved models of genome/gene evolution and for engineering gene expression in various biotechnological and synthetic biology applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1497-z) contains supplementary material, which is available to authorized users.

Unsupervised statistical discovery of spaced motifs in prokaryotic genomes

BACKGROUND

DNA sequences contain repetitive motifs which have various functions in the physiology of the organism. A number of methods have been developed for discovery of such sequence motifs with a primary focus on detection of regulatory motifs and particularly transcription factor binding sites. Most motif-finding methods apply probabilistic models to detect motifs characterized by unusually high number of copies of the motif in the analyzed sequences.

RESULTS

We present a novel method for detection of pairs of motifs separated by spacers of variable nucleotide sequence but conserved length. Unlike existing methods for motif discovery, the motifs themselves are not required to occur at unusually high frequency but only to exhibit a significant preference to occur at a specific distance from each other. In the present implementation of the method, motifs are represented by pentamers and all pairs of pentamers are evaluated for statistically significant preference for a specific distance. An important step of the algorithm eliminates motif pairs where the spacers separating the two motifs exhibit a high degree of sequence similarity; such motif pairs likely arise from duplications of the whole segment including the motifs and the spacer rather than due to selective constraints indicative of a functional importance of the motif pair. The method was used to scan 569 complete prokaryotic genomes for novel sequence motifs. Some motifs detected were previously known but other motifs found in the search appear to be novel. Selected motif pairs were subjected to further investigation and in some cases their possible biological functions were proposed.

CONCLUSIONS

We present a new motif-finding technique that is applicable to scanning complete genomes for sequence motifs. The results from analysis of 569 genomes suggest that the method detects previously known motifs that are expected to be found as well as new motifs that are unlikely to be discovered by traditional motif-finding methods. We conclude that our approach to detection of significant motif pairs can complement existing motif-finding techniques in discovery of novel functional sequence motifs in complete genomes.

Plasticity of archaeal C/D box sRNA biogenesis

Archaeal and eukaryotic organisms contain sets of C/D box s(no)RNAs with guide sequences that determine ribose 2'-O-methylation sites of target RNAs. The composition of these C/D box sRNA sets is highly variable between organisms and results in varying RNA modification patterns which are important for ribosomal RNA folding and stability. Little is known about the genomic organization of C/D box sRNA genes in archaea. Here, we aimed to obtain first insights into the biogenesis of these archaeal C/D box sRNAs and analyzed the genetic context of more than 300 archaeal sRNA genes. We found that the majority of these genes do not possess independent promoters but are rather located at positions that allow for co-transcription with neighboring genes and their start or stop codons were frequently incorporated into the conserved boxC and D motifs. The biogenesis of plasmid-encoded C/D box sRNA variants was analyzed in vivo in Sulfolobus acidocaldarius. It was found that C/D box sRNA maturation occurs independent of their genetic context and relies solely on the presence of intact RNA kink-turn structures. The observed plasticity of C/D box sRNA biogenesis is suggested to enable their accelerated evolution and, consequently, allow for adjustments of the RNA modification landscape.

Genome-wide identification of transcriptional start sites in the haloarchaeon Haloferax volcanii based on differential RNA-Seq (dRNA-Seq)

BACKGROUND

Differential RNA-Seq (dRNA-Seq) is a recently developed method of performing primary transcriptome analyses that allows for the genome-wide mapping of transcriptional start sites (TSSs) and the identification of novel transcripts. Although the transcriptomes of diverse bacterial species have been characterized by dRNA-Seq, the transcriptome analysis of archaeal species is still rather limited. Therefore, we used dRNA-Seq to characterize the primary transcriptome of the model archaeon Haloferax volcanii.

RESULTS

Three independent cultures of Hfx. volcanii grown under optimal conditions to the mid-exponential growth phase were used to determine the primary transcriptome and map the 5'-ends of the transcripts. In total, 4749 potential TSSs were detected. A position weight matrix (PWM) was derived for the promoter predictions, and the results showed that 64 % of the TSSs were preceded by stringent or relaxed basal promoters. Of the identified TSSs, 1851 belonged to protein-coding genes. Thus, fewer than half (46 %) of the 4040 protein-coding genes were expressed under optimal growth conditions. Seventy-two percent of all protein-coding transcripts were leaderless, which emphasized that this pathway is the major pathway for translation initiation in haloarchaea. A total of 2898 of the TSSs belonged to potential non-coding RNAs, which accounted for an unexpectedly high fraction (61 %) of all transcripts. Most of the non-coding TSSs had not been previously described (2792) and represented novel sequences (59 % of all TSSs). A large fraction of the potential novel non-coding transcripts were cis-antisense RNAs (1244 aTSSs). A strong negative correlation between the levels of antisense transcripts and cognate sense mRNAs was found, which suggested that the negative regulation of gene expression via antisense RNAs may play an important role in haloarchaea. The other types of novel non-coding transcripts corresponded to internal transcripts overlapping with mRNAs (1153 iTSSs) and intergenic small RNA (sRNA) candidates (395 TSSs).

CONCLUSION

This study provides a comprehensive map of the primary transcriptome of Hfx. volcanii grown under optimal conditions. Fewer than half of all protein-coding genes have been transcribed under these conditions. Unexpectedly, more than half of the detected TSSs belonged to several classes of non-coding RNAs. Thus, RNA-based regulation appears to play a more important role in haloarchaea than previously anticipated.

Coevolution Theory of the Genetic Code at Age Forty: Pathway to Translation and Synthetic Life

The origins of the components of genetic coding are examined in the present study. Genetic information arose from replicator induction by metabolite in accordance with the metabolic expansion law. Messenger RNA and transfer RNA stemmed from a template for binding the aminoacyl-RNA synthetase ribozymes employed to synthesize peptide prosthetic groups on RNAs in the Peptidated RNA World. Coevolution of the genetic code with amino acid biosynthesis generated tRNA paralogs that identify a last universal common ancestor (LUCA) of extant life close to Methanopyrus, which in turn points to archaeal tRNA introns as the most primitive introns and the anticodon usage of Methanopyrus as an ancient mode of wobble. The prediction of the coevolution theory of the genetic code that the code should be a mutable code has led to the isolation of optional and mandatory synthetic life forms with altered protein alphabets.

Interference activity of a minimal Type I CRISPR-Cas system from Shewanella putrefaciens

Type I CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas (CRISPR-associated) systems exist in bacterial and archaeal organisms and provide immunity against foreign DNA. The Cas protein content of the DNA interference complexes (termed Cascade) varies between different CRISPR-Cas subtypes. A minimal variant of the Type I-F system was identified in proteobacterial species including Shewanella putrefaciens CN-32. This variant lacks a large subunit (Csy1), Csy2 and Csy3 and contains two unclassified cas genes. The genome of S. putrefaciens CN-32 contains only five Cas proteins (Cas1, Cas3, Cas6f, Cas1821 and Cas1822) and a single CRISPR array with 81 spacers. RNA-Seq analyses revealed the transcription of this array and the maturation of crRNAs (CRISPR RNAs). Interference assays based on plasmid conjugation demonstrated that this CRISPR-Cas system is active in vivo and that activity is dependent on the recognition of the dinucleotide GG PAM (Protospacer Adjacent Motif) sequence and crRNA abundance. The deletion of cas1821 and cas1822 reduced the cellular crRNA pool. Recombinant Cas1821 was shown to form helical filaments bound to RNA molecules, which suggests its role as the Cascade backbone protein. A Cascade complex was isolated which contained multiple Cas1821 copies, Cas1822, Cas6f and mature crRNAs.

A Member of the Arabidopsis Mitochondrial Transcription Termination Factor Family Is Required for Maturation of Chloroplast Transfer RNAIle(GAU)

Plastid gene expression is crucial for organelle function, but the factors that control it are still largely unclear. Members of the so-called mitochondrial transcription termination factor (mTERF) family are found in metazoans and plants and regulate organellar gene expression at different levels. Arabidopsis (Arabidopsis thaliana) mTERF6 is localized in chloroplasts and mitochondria, and its knockout perturbs plastid development and results in seedling lethality. In the leaky mterf6-1 mutant, a defect in photosynthesis is associated with reduced levels of photosystem subunits, although corresponding messenger RNA levels are unaffected, whereas translational capacity and maturation of chloroplast ribosomal RNAs (rRNAs) are perturbed in mterf6-1 mutants. Bacterial one-hybrid screening, electrophoretic mobility shift assays, and coimmunoprecipitation experiments reveal a specific interaction between mTERF6 and an RNA sequence in the chloroplast isoleucine transfer RNA gene (trnI.2) located in the rRNA operon. In vitro, recombinant mTERF6 bound to its plastid DNA target site can terminate transcription. At present, it is unclear whether disturbed rRNA maturation is a primary or secondary defect. However, it is clear that mTERF6 is required for the maturation of trnI.2. This points to an additional function of mTERFs.

C/D box sRNA-guided 2'-O-methylation patterns of archaeal rRNA molecules

Background

In archaea and eukaryotes, ribonucleoprotein complexes containing small C/D box s(no)RNAs use base pair complementarity to target specific sites within ribosomal RNA for 2'-O-ribose methylation. These modifications aid in the folding and stabilization of nascent rRNA molecules and their assembly into ribosomal particles. The genomes of hyperthermophilic archaea encode large numbers of C/D box sRNA genes, suggesting an increased necessity for rRNA stabilization at extreme growth temperatures.

Results

We have identified the complete sets of C/D box sRNAs from seven archaea using RNA-Seq methodology. In total, 489 C/D box sRNAs were identified, each containing two guide regions. A combination of computational and manual analyses predicts 719 guide interactions with 16S and 23S rRNA molecules. This first pan-archaeal description of guide sequences identifies (i) modified rRNA nucleotides that are frequently conserved between species and (ii) regions within rRNA that are hotspots for 2'-O-methylation. Gene duplication, rearrangement, mutational drift and convergent evolution of sRNA genes and guide sequences were observed. In addition, several C/D box sRNAs were identified that use their two guides to target locations distant in the rRNA sequence but close in the secondary and tertiary structure. We propose that they act as RNA chaperones and facilitate complex folding events between distant sequences.

Conclusions

This pan-archaeal analysis of C/D box sRNA guide regions identified conserved patterns of rRNA 2'-O-methylation in archaea. The interaction between the sRNP complexes and the nascent rRNA facilitates proper folding and the methyl modifications stabilize higher order rRNA structure within the assembled ribosome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1839-z) contains supplementary material, which is available to authorized users.

Biogenesis pathways of RNA guides in archaeal and bacterial CRISPR-Cas adaptive immunity

CRISPR-Cas is an RNA-mediated adaptive immune system that defends bacteria and archaea against mobile genetic elements. Short mature CRISPR RNAs (crRNAs) are key elements in the interference step of the immune pathway. A CRISPR array composed of a series of repeats interspaced by spacer sequences acquired from invading mobile genomes is transcribed as a precursor crRNA (pre-crRNA) molecule. This pre-crRNA undergoes one or two maturation steps to generate the mature crRNAs that guide CRISPR-associated (Cas) protein(s) to cognate invading genomes for their destruction. Different types of CRISPR-Cas systems have evolved distinct crRNA biogenesis pathways that implicate highly sophisticated processing mechanisms. In Types I and III CRISPR-Cas systems, a specific endoribonuclease of the Cas6 family, either standalone or in a complex with other Cas proteins, cleaves the pre-crRNA within the repeat regions. In Type II systems, the trans-acting small RNA (tracrRNA) base pairs with each repeat of the pre-crRNA to form a dual-RNA that is cleaved by the housekeeping RNase III in the presence of the protein Cas9. In this review, we present a detailed comparative analysis of pre-crRNA recognition and cleavage mechanisms involved in the biogenesis of guide crRNAs in the three CRISPR-Cas types.

This review presents a detailed comparative analysis of pre-crRNA recognition and cleavage mechanisms involved in the biogenesis of guide crRNAs in the different bacterial and archaeal CRISPR-Cas immune systems.

Inosine modifications in human tRNAs are incorporated at the precursor tRNA level

Transfer RNAs (tRNAs) are key adaptor molecules of the genetic code that are heavily modified post-transcriptionally. Inosine at the first residue of the anticodon (position 34; I34) is an essential widespread tRNA modification that has been poorly studied thus far. The modification in eukaryotes results from a deamination reaction of adenine that is catalyzed by the heterodimeric enzyme adenosine deaminase acting on tRNA (hetADAT), composed of two subunits: ADAT2 and ADAT3. Using high-throughput small RNA sequencing (RNAseq), we show that this modification is incorporated to human tRNAs at the precursor tRNA level and during maturation. We also functionally validated the human genes encoding for hetADAT and show that the subunits of this enzyme co-localize in nucleus in an ADAT2-dependent manner. Finally, by knocking down HsADAT2, we demonstrate that variations in the cellular levels of hetADAT will result in changes in the levels of I34 modification in all its potential substrates. Altogether, we present RNAseq as a powerful tool to study post-transcriptional tRNA modifications at the precursor tRNA level and give the first insights on the biology of I34 tRNA modification in metazoans.

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.

Investigating CRISPR RNA Biogenesis and Function Using RNA-seq

The development of deep sequencing technology has greatly facilitated transcriptome analyses of both prokaryotes and eukaryotes. RNA-sequencing (RNA-seq), which is based on massively parallel sequencing of cDNAs, has been used to annotate transcript boundaries and revealed widespread antisense transcription as well as a wealth of novel noncoding transcripts in many bacteria. Moreover, RNA-seq is nowadays widely used for gene expression profiling and about to replace hybridization-based approaches such as microarrays. RNA-seq has also informed about the biogenesis and function of CRISPR RNAs (crRNAs) of different types of bacterial RNA-based CRISPR-Cas immune systems. Here we describe several studies that employed RNA-seq for crRNA analyses, with a particular focus on a differential RNA-seq (dRNA-seq) approach, which can distinguish between primary and processed transcripts and allows for a genome-wide annotation of transcriptional start sites. This approach helped to identify a new crRNA biogenesis pathway of Type II CRISPR-Cas systems that involves a trans-encoded small RNA, tracrRNA, and the host factor RNase III.

Transcription start site associated RNAs (TSSaRNAs) are ubiquitous in all domains of life

A plethora of non-coding RNAs has been discovered using high-resolution transcriptomics tools, indicating that transcriptional and post-transcriptional regulation is much more complex than previously appreciated. Small RNAs associated with transcription start sites of annotated coding regions (TSSaRNAs) are pervasive in both eukaryotes and bacteria. Here, we provide evidence for existence of TSSaRNAs in several archaeal transcriptomes including: Halobacterium salinarum, Pyrococcus furiosus, Methanococcus maripaludis, and Sulfolobus solfataricus. We validated TSSaRNAs from the model archaeon Halobacterium salinarum NRC-1 by deep sequencing two independent small-RNA enriched (RNA-seq) and a primary-transcript enriched (dRNA-seq) strand-specific libraries. We identified 652 transcripts, of which 179 were shown to be primary transcripts (∼7% of the annotated genome). Distinct growth-associated expression patterns between TSSaRNAs and their cognate genes were observed, indicating a possible role in environmental responses that may result from RNA polymerase with varying pausing rhythms. This work shows that TSSaRNAs are ubiquitous across all domains of life.

RNA-Seq analyses reveal CRISPR RNA processing and regulation patterns

In bacteria and archaea, RNA-Seq deep sequencing methodology allows for the detection of abundance and processing sites of the small RNAs that comprise a CRISPR (clustered regularly interspaced short palindromic repeats) RNome. Comparative analyses of these CRISPR RNome sets highlight conserved patterns that include the gradual decline of CRISPR RNA abundance from the leader-proximal to the leader-distal end. In the present review, we discuss exceptions to these patterns that indicate the extensive impact of individual spacer sequences on CRISPR array transcription and RNA maturation. Spacer sequences can contain promoter and terminator elements and can promote the formation of CRISPR RNA-anti-CRISPR RNA duplexes. In addition, potential RNA duplex formation with host tRNA was observed. These factors can influence the functionality of CRISPR-Cas (CRISPR-associated) systems and need to be considered in the design of synthetic CRISPR arrays.

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

Small regulatory RNAs in Archaea

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

In vitro assembly and activity of an archaeal CRISPR-Cas type I-A Cascade interference complex

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR-associated (Cas) systems of type I use a Cas ribonucleoprotein complex for antiviral defense (Cascade) to mediate the targeting and degradation of foreign DNA. To address molecular features of the archaeal type I-A Cascade interference mechanism, we established the in vitro assembly of the Thermoproteus tenax Cascade from six recombinant Cas proteins, synthetic CRISPR RNAs (crRNAs) and target DNA fragments. RNA-Seq analyses revealed the processing pattern of crRNAs from seven T. tenax CRISPR arrays. Synthetic crRNA transcripts were matured by hammerhead ribozyme cleavage. The assembly of type I-A Cascade indicates that Cas3′ and Cas3′′ are an integral part of the complex, and the interference activity was shown to be dependent on the crRNA and the matching target DNA. The reconstituted Cascade was used to identify sequence motifs that are required for efficient DNA degradation and to investigate the role of the subunits Cas7 and Cas3′′ in the interplay with other Cascade subunits.

Structure and RNA-binding properties of the type III-A CRISPR-associated protein Csm3

The prokaryotic adaptive immune system is based on the incorporation of genome fragments of invading viral genetic elements into clusters of regulatory interspaced short palindromic repeats (CRISPRs). The CRISPR loci are transcribed and processed into crRNAs, which are then used to target the invading nucleic acid for degradation. The large family of CRISPR-associated (Cas) proteins mediates this interference response. We have characterized Methanopyrus kandleri Csm3, a protein of the type III-A CRISPR-Cas complex. The 2.4 Å resolution crystal structure shows an elaborate four-domain fold organized around a core RRM-like domain. The overall architecture highlights the structural homology to Cas7, the Cas protein that forms the backbone of type I interference complexes. Csm3 binds unstructured RNAs in a sequence non-specific manner, suggesting that it interacts with the variable spacer sequence of the crRNA. The structural and biochemical data provide insights into the similarities and differences in this group of Cas proteins.