Regulation of translation in haloarchaea: 5'- and 3'-UTRs are essential and have to functionally interact in vivo

5'-untranslated region sequences enhance plasmid-based protein production in Sulfolobus acidocaldarius

Sulfolobus acidocaldarius, a thermoacidophilic archaeon of the phylum Thermoproteota (former Crenarchaeota), is a widely used model organism for gene deletion studies and recombinant protein production. Previous research has demonstrated the efficacy of the saci_2122 promoter (Para), providing low basal activity and high pentose-dependent induction. However, the available expression vector does not include a 5'-terminal untranslated region (5'-UTR), a typical element found in bacterial expression vectors that usually enhances protein production in bacteria. To establish S. acidocaldarius as a production strain in biotechnology in the long term, it is intrinsically relevant to optimize its tools and capacities to increase production efficiencies. Here we show that protein production is increased by the integration of S. acidocaldarius 5'-UTRs into Para expression plasmids. Using the esterase Saci_1116 as a reporter protein, we observed a four-fold increase in soluble and active protein yield upon insertion of the saci_1322 (alba) 5'-UTR. Screening of four additional 5'-UTRs from other highly abundant proteins (thα, slaA, slaB, saci_0330) revealed a consistent enhancement in target protein production. Additionally, site-directed mutagenesis of the Shine-Dalgarno (SD) motif within the alba 5'-UTR revealed its significance for protein synthesis. Ultimately, the alba 5'-UTR optimized expression vector improved the expression of various proteins, including six glycosyltransferases and one hydroxyacyl-CoA dehydratase from S. acidocaldarius, and a malto-oligosyltrehalose trehalohydrolase from Saccharolobus solfataricus, demonstrating its applicability. Our results show that the integration of SD-motif containing 5'-UTRs significantly enhanced plasmid-based protein production in S. acidocaldarius. This advancement in recombinant expression not only broadens the utility of S. acidocaldarius as an archaeal expression platform but also marks an important step toward potential biotechnological applications.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

Tuning Gene Expression by Phosphate in the Methanogenic Archaeon Methanococcus maripaludis

Methanococcus maripaludis is a rapidly growing, hydrogenotrophic, and genetically tractable methanogen with unique capabilities to convert formate and CO₂ to CH₄. The existence of genome-scale metabolic models and an established, robust system for both large-scale and continuous cultivation make it amenable for industrial applications. However, the lack of molecular tools for differential gene expression has hindered its application as a microbial cell factory to produce biocatalysts and biochemicals. In this study, a library of differentially regulated promoters was designed and characterized based on the pst promoter, which responds to the inorganic phosphate concentration in the growth medium. Gene expression increases by 4- to 6-fold when the medium phosphate drops to growth-limiting concentrations. Hence, this regulated system decouples growth from heterologous gene expression without the need for adding an inducer. The minimal pst promoter is identified and contains a conserved AT-rich region, a factor B recognition element, and a TATA box for phosphate-dependent regulation. Rational changes to the factor B recognition element and start codon had no significant impact on expression; however, changes to the transcription start site and the 5' untranslated region resulted in the differential protein production with regulation remaining intact. Compared to a previous expression system based upon the histone promoter, this regulated expression system resulted in significant improvements in the expression of a key methanogenic enzyme complex, methyl-coenzyme M reductase, and the potentially toxic arginine methyltransferase MmpX.

Elucidation of the Translation Initiation Factor Interaction Network of Haloferax volcanii Reveals Coupling of Transcription and Translation in Haloarchaea

Translation is an important step in gene expression. Initiation of translation is rate-limiting, and it is phylogenetically more diverse than elongation or termination. Bacteria contain only three initiation factors. In stark contrast, eukaryotes contain more than 10 (subunits of) initiation factors (eIFs). The genomes of archaea contain many genes that are annotated to encode archaeal homologs of eukaryotic initiation factors (aIFs). However, experimental characterization of aIFs is scarce and mostly restricted to very few species. To broaden the view, the protein-protein interaction network of aIFs in the halophilic archaeon Haloferax volcanii has been characterized. To this end, tagged versions of 14 aIFs were overproduced, affinity isolated, and the co-isolated binding partners were identified by peptide mass fingerprinting and MS/MS analyses. The aIF-aIF interaction network was resolved, and it was found to contain two interaction hubs, (1) the universally conserved factor aIF5B, and (2) a protein that has been annotated as the enzyme ribose-1,5-bisphosphate isomerase, which we propose to rename to aIF2Bα. Affinity isolation of aIFs also led to the co-isolation of many ribosomal proteins, but also transcription factors and subunits of the RNA polymerase (Rpo). To analyze a possible coupling of transcription and translation, seven tagged Rpo subunits were overproduced, affinity isolated, and co-isolated proteins were identified. The Rpo interaction network contained many transcription factors, but also many ribosomal proteins as well as the initiation factors aIF5B and aIF2Bα. These results showed that transcription and translation are coupled in haloarchaea, like in Escherichia coli. It seems that aIF5B and aIF2Bα are not only interaction hubs in the translation initiation network, but also key players in the transcription-translation coupling.

Halobacterium salinarum and Haloferax volcanii Comparative Transcriptomics Reveals Conserved Transcriptional Processing Sites

Post-transcriptional processing of messenger RNA is an important regulatory strategy that allows relatively fast responses to changes in environmental conditions. In halophile systems biology, the protein perspective of this problem (i.e., ribonucleases which implement the cleavages) is generally more studied than the RNA perspective (i.e., processing sites). In the present in silico work, we mapped genome-wide transcriptional processing sites (TPS) in two halophilic model organisms, Halobacterium salinarum NRC-1 and Haloferax volcanii DS2. TPS were established by reanalysis of publicly available differential RNA-seq (dRNA-seq) data, searching for non-primary (monophosphorylated RNAs) enrichment. We found 2093 TPS in 43% of H. salinarum genes and 3515 TPS in 49% of H. volcanii chromosomal genes. Of the 244 conserved TPS sites found, the majority were located around start and stop codons of orthologous genes. Specific genes are highlighted when discussing antisense, ribosome and insertion sequence associated TPS. Examples include the cell division gene ftsZ2, whose differential processing signal along growth was detected and correlated with post-transcriptional regulation, and biogenesis of sense overlapping transcripts associated with IS200/IS605. We hereby present the comparative, transcriptomics-based processing site maps with a companion browsing interface.

A Synthetic Riboswitch to Regulate Haloarchaeal Gene Expression

In recent years, synthetic riboswitches have become increasingly important to construct genetic circuits in all three domains of life. In bacteria, synthetic translational riboswitches are often employed that modulate gene expression by masking the Shine-Dalgarno (SD) sequence in the absence or presence of a cognate ligand. For (halo-)archaeal translation, a SD sequence is not strictly required. The application of synthetic riboswitches in haloarchaea is therefore limited so far, also because of the molar intracellular salt concentrations found in these microbes. In this study, we applied synthetic theophylline-dependent translational riboswitches in the archaeon Haloferax volcanii. The riboswitch variants A through E and E^∗ were chosen since they not only mask the SD sequence but also the AUG start codon by forming a secondary structure in the absence of the ligand theophylline. Upon addition of the ligand, the ribosomal binding site and start codon become accessible for translation initiation. Riboswitch E mediated a dose-dependent, up to threefold activation of the bgaH reporter gene expression. Raising the salt concentration of the culture media from 3 to 4 M NaCl resulted in a 12-fold increase in the switching capacity of riboswitch E, and switching activity increased up to 26-fold when the cultivating temperature was reduced from 45 to 30°C. To construct a genetic circuit, riboswitch E was applied to regulate the synthesis of the transcriptional activator GvpE allowing a dose-dependent activation of the mgfp6 reporter gene under P _pA promoter control.

Influence of the spacer region between the Shine-Dalgarno box and the start codon for fine-tuning of the translation efficiency in Escherichia coli

Translation efficiency contributes several orders of magnitude difference in the overall yield of exogenous gene expression in bacteria. In diverse bacteria, the translation initiation site, whose sequence is the primary determinant of the translation performance, is comprised of the start codon and the Shine-Dalgarno box located upstream. Here, we have examined how the sequence of a spacer between these main components of the translation initiation site contributes to the yield of synthesized protein. We have created a library of reporter constructs with the randomized spacer region, performed fluorescently activated cell sorting and applied next-generation sequencing analysis (the FlowSeq protocol). As a result, we have identified sequence motifs for the spacer region between the Shine-Dalgarno box and AUG start codon that may modulate the translation efficiency in a 100-fold range.

Ribosome profiling in archaea reveals leaderless translation, novel translational initiation sites, and ribosome pausing at single codon resolution

High-throughput methods, such as ribosome profiling, have revealed the complexity of translation regulation in Bacteria and Eukarya with large-scale effects on cellular functions. In contrast, the translational landscape in Archaea remains mostly unexplored. Here, we developed ribosome profiling in a model archaeon, Haloferax volcanii, elucidating, for the first time, the translational landscape of a representative of the third domain of life. We determined the ribosome footprint of H. volcanii to be comparable in size to that of the Eukarya. We linked footprint lengths to initiating and elongating states of the ribosome on leadered transcripts, operons, and on leaderless transcripts, the latter representing 70% of H. volcanii transcriptome. We manipulated ribosome activity with translation inhibitors to reveal ribosome pausing at specific codons. Lastly, we found that the drug harringtonine arrested ribosomes at initiation sites in this archaeon. This drug treatment allowed us to confirm known translation initiation sites and also reveal putative novel initiation sites in intergenic regions and within genes. Ribosome profiling revealed an uncharacterized complexity of translation in this archaeon with bacteria-like, eukarya-like, and potentially novel translation mechanisms. These mechanisms are likely to be functionally essential and to contribute to an expanded proteome with regulatory roles in gene expression.

Translational coupling via termination-reinitiation in archaea and bacteria

The genomes of many prokaryotes contain substantial fractions of gene pairs with overlapping stop and start codons (ATGA or TGATG). A potential benefit of overlapping gene pairs is translational coupling. In 720 genomes of archaea and bacteria representing all major phyla, we identify substantial, albeit highly variable, fractions of co-directed overlapping gene pairs. Various patterns are observed for the utilization of the SD motif for de novo initiation at upstream genes versus reinitiation at overlapping gene pairs. We experimentally test the predicted coupling in 9 gene pairs from the archaeon Haloferax volcanii and 5 gene pairs from the bacterium Escherichia coli. In 13 of 14 cases, translation of both genes is strictly coupled. Mutational analysis of SD motifs located upstream of the downstream genes indicate that the contribution of the SD to translational coupling widely varies from gene to gene. The nearly universal, abundant occurrence of overlapping gene pairs suggests that tight translational coupling is widespread in archaea and bacteria.

Archaea and bacteria often have gene pairs with overlapping stop and start codons, suggesting translational coupling. Here, Huber et al. analyse overlapping gene pairs from 720 genomes, and validate translational coupling via termination-reinitiation for 14 gene pairs in Haloferax volcanii and Escherichia coli.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

Unique Archaeal Small RNAs

Advances in genome-wide sequence technologies allow for detailed insights into the complexity of RNA landscapes of organisms from all three domains of life. Recent analyses of archaeal transcriptomes identified interaction and regulation networks of noncoding RNAs in this understudied domain. Here, we review current knowledge of small, noncoding RNAs with important functions for the archaeal lifestyle, which often requires adaptation to extreme environments. One focus is RNA metabolism at elevated temperatures in hyperthermophilic archaea, which reveals elevated amounts of RNA-guided RNA modification and virus defense strategies. Genome rearrangement events result in unique fragmentation patterns of noncoding RNA genes that require elaborate maturation pathways to yield functional transcripts. RNA-binding proteins, e.g., L7Ae and LSm, are important for many posttranscriptional control functions of RNA molecules in archaeal cells. We also discuss recent insights into the regulatory potential of their noncoding RNA partners.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Genome-wide gene expression and RNA half-life measurements allow predictions of regulation and metabolic behavior in Methanosarcina acetivorans

BACKGROUND

While a few studies on the variations in mRNA expression and half-lives measured under different growth conditions have been used to predict patterns of regulation in bacterial organisms, the extent to which this information can also play a role in defining metabolic phenotypes has yet to be examined systematically. Here we present the first comprehensive study for a model methanogen.

RESULTS

We use expression and half-life data for the methanogen Methanosarcina acetivorans growing on fast- and slow-growth substrates to examine the regulation of its genes. Unlike Escherichia coli where only small shifts in half-lives were observed, we found that most mRNA have significantly longer half-lives for slow growth on acetate compared to fast growth on methanol or trimethylamine. Interestingly, half-life shifts are not uniform across functional classes of enzymes, suggesting the existence of a selective stabilization mechanism for mRNAs. Using the transcriptomics data we determined whether transcription or degradation rate controls the change in transcript abundance. Degradation was found to control abundance for about half of the metabolic genes underscoring its role in regulating metabolism. Genes involved in half of the metabolic reactions were found to be differentially expressed among the substrates suggesting the existence of drastically different metabolic phenotypes that extend beyond just the methanogenesis pathways. By integrating expression data with an updated metabolic model of the organism (iST807) significant differences in pathway flux and production of metabolites were predicted for the three growth substrates.

CONCLUSIONS

This study provides the first global picture of differential expression and half-lives for a class II methanogen, as well as provides the first evidence in a single organism that drastic genome-wide shifts in RNA half-lives can be modulated by growth substrate. We determined which genes in each metabolic pathway control the flux and classified them as regulated by transcription (e.g. transcription factor) or degradation (e.g. post-transcriptional modification). We found that more than half of genes in metabolism were controlled by degradation. Our results suggest that M. acetivorans employs extensive post-transcriptional regulation to optimize key metabolic steps, and more generally that degradation could play a much greater role in optimizing an organism's metabolism than previously thought.

In silico analysis of 5'-UTRs highlights the prevalence of Shine-Dalgarno and leaderless-dependent mechanisms of translation initiation in bacteria and archaea, respectively

In prokaryotes, a heterogeneous set of protein translation initiation mechanisms such as Shine-Dalgarno (SD) sequence-dependent, SD sequence-independent or ribosomal protein S1 mediated and leaderless transcript-dependent exists. To estimate the distribution of coding sequences employing a particular translation initiation mechanism, a total of 107 prokaryotic genomes were analysed using in silico approaches. Analysis of 5'-untranslated regions (UTRs) of genes reveals the existence of three types of mRNAs described as transcripts with and without SD motif and leaderless transcripts. Our results indicate that although all the three types of translation initiation mechanisms are widespread among prokaryotes, the number of SD-dependent genes in bacteria is higher than that of archaea. In contrast, archaea contain a significantly higher number of leaderless genes than SD-led genes. The correlation analysis between genome size and SD-led & leaderless genes suggests that the SD-led genes are decreasing (increasing) with genome size in bacteria (archaea). However, the leaderless genes are increasing (decreasing) in bacteria (archaea) with genome size. Moreover, an analysis of the start-codon biasness confirms that among ATG, GTG and TTG codons, ATG is indeed the most preferred codon at the translation initiation site in most of the coding sequences. In leaderless genes, however, the codons GTG and TTG are also observed at the translation initiation site in some species contradicting earlier studies which suggested the usage of only ATG codon. Henceforth, the conventional mechanism of translation initiation cannot be generalized as an exclusive way of initiating the process of protein biosynthesis in prokaryotes.

Enhanced Salt Tolerance Conferred by the Complete 2.3 kb cDNA of the Rice Vacuolar Na(+)/H(+) Antiporter Gene Compared to 1.9 kb Coding Region with 5' UTR in Transgenic Lines of Rice

Soil salinity is one of the most challenging problems that restricts the normal growth and production of rice worldwide. It has therefore become very important to produce more saline tolerant rice varieties. This study shows constitutive over-expression of the vacuolar Na(+)/H(+) antiporter gene (OsNHX1) from the rice landrace (Pokkali) and attainment of enhanced level of salinity tolerance in transgenic rice plants. It also shows that inclusion of the complete un-translated regions (UTRs) of the alternatively spliced OsNHX1 gene provides a higher level of tolerance to the transgenic rice. Two separate transformation events of the OsNHX1 gene, one with 1.9 kb region containing the 5' UTR with CDS and the other of 2.3 kb, including 5' UTR, CDS, and the 3' UTR regions were performed. The transgenic plants with these two different constructs were advanced to the T3 generation and physiological and molecular screening of homozygous plants was conducted at seedling and reproductive stages under salinity (NaCl) stress. Both transgenic lines were observed to be tolerant compared to WT plants at both physiological stages. However, the transgenic lines containing the CDS with both the 5' and 3' UTR were significantly more tolerant compared to the transgenic lines containing OsNHX1 gene without the 3' UTR. At the seedling stage at 12 dS/m stress, the chlorophyll content was significantly higher (P < 0.05) and the electrolyte leakage significantly lower (P < 0.05) in the order 2.3 kb > 1.9 kb > and WT lines. Yield in g/plant in the best line from the 2.3 kb plants was significantly more (P < 0.01) compared, respectively, to the best 1.9 kb line and WT plants at stress of 6 dS/m. Transformation with the complete transcripts rather than the CDS may therefore provide more durable level of tolerance.

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.

Haloferax volcanii, a prokaryotic species that does not use the Shine Dalgarno mechanism for translation initiation at 5'-UTRs

It was long assumed that translation initiation in prokaryotes generally occurs via the so-called Shine Dalgarno (SD) mechanism. Recently, it became clear that translation initiation in prokaryotes is more heterogeneous. In the haloarchaeon Haloferax volcanii, the majority of transcripts is leaderless and most transcripts with a 5′-UTR lack a SD motif. Nevertheless, a bioinformatic analysis predicted that 20–30% of all genes are preceded by a SD motif in haloarchaea. To analyze the importance of the SD mechanism for translation initiation in haloarchaea experimentally the monocistronic sod gene was chosen, which contains a 5′-UTR with an extensive SD motif of seven nucleotides and a length of 19 nt, the average length of 5′UTRs in this organism. A translational fusion of part of the sod gene with the dhfr reporter gene was constructed. A mutant series was generated that matched the SD motif from zero to eight positions, respectively. Surprisingly, there was no correlation between the base pairing ability between transcripts and 16S rRNA and translational efficiency in vivo under several different growth conditions. Furthermore, complete replacement of the SD motif by three unrelated sequences did not reduce translational efficiency. The results indicate that H. volcanii does not make use of the SD mechanism for translation initiation in 5′-UTRs. A genome analysis revealed that while the number of SD motifs in 5′-UTRs is rare, their fraction within open reading frames is high. Possible biological functions for intragenic SD motifs are discussed, including re-initiation of translation at distal genes in operons.

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.

Contribution of transcriptomics to systems-level understanding of methanogenic Archaea

Methane-producing Archaea are of interest due to their contribution to atmospheric change and for their roles in technological applications including waste treatment and biofuel production. Although restricted to anaerobic environments, methanogens are found in a wide variety of habitats, where they commonly live in syntrophic relationships with bacterial partners. Owing to tight thermodynamic constraints of methanogenesis alone or in syntrophic metabolism, methanogens must carefully regulate their catabolic pathways including the regulation of RNA transcripts. The transcriptome is a dynamic and important control point in microbial systems. This paper assesses the impact of mRNA (transcriptome) studies on the understanding of methanogenesis with special consideration given to how methanogenesis is regulated to cope with nutrient limitation, environmental variability, and interactions with syntrophic partners. In comparison with traditional microarray-based transcriptome analyses, next-generation high-throughput RNA sequencing is greatly advantageous in assessing transcription start sites, the extent of 5′ untranslated regions, operonic structure, and the presence of small RNAs. We are still in the early stages of understanding RNA regulation but it is already clear that determinants beyond transcript abundance are highly relevant to the lifestyles of methanogens, requiring further study.

Regulatory RNAs in archaea: first target identification in Methanoarchaea

sRNAs (small non-coding RNAs) representing important players in many cellular and regulatory processes have been identified in all three domains of life. In Eukarya and Bacteria, functions have been assigned for many sRNAs, whereas the sRNA populations in Archaea are considerably less well characterized. Recent analyses on a genome-wide scale particularly using high-throughput sequencing techniques demonstrated the presence of high numbers of sRNA candidates in several archaea. However, elucidation of the molecular mechanism of sRNA action, as well as understanding their physiological roles, is in general still challenging, particularly in Archaea, since efficient genetic tools are missing. The identification of cellular targets of identified archaeal sRNAs by experimental approaches or computational prediction programs has begun only recently. At present, targets have been identified for one archaeal sRNA, sRNA162 in Methanosarcina mazei, which interacts with the 5′ region of its targets, a cis-encoded and a trans-encoded target, blurring the paradigm of a border between cis- and trans-encoded sRNAs. Besides, the first experimental implications have been obtained in Haloarchaea and Pyrobaculum that archaeal sRNAs also target 3′ regions of mRNAs. The present review summarizes our current knowledge on archaeal sRNAs and their biological functions and targets.

An archaeal sRNA targeting cis- and trans-encoded mRNAs via two distinct domains

We report on the characterization and target analysis of the small (s)RNA₁₆₂ in the methanoarchaeon Methanosarcina mazei. Using a combination of genetic approaches, transcriptome analysis and computational predictions, the bicistronic MM2441-MM2440 mRNA encoding the transcription factor MM2441 and a protein of unknown function was identified as a potential target of this sRNA, which due to processing accumulates as three stabile 5′ fragments in late exponential growth. Mobility shift assays using various mutants verified that the non-structured single-stranded linker region of sRNA₁₆₂ (SLR) base-pairs with the MM2440-MM2441 mRNA internally, thereby masking the predicted ribosome binding site of MM2441. This most likely leads to translational repression of the second cistron resulting in dis-coordinated operon expression. Analysis of mutant RNAs in vivo confirmed that the SLR of sRNA₁₆₂ is crucial for target interactions. Furthermore, our results indicate that sRNA₁₆₂-controlled MM2441 is involved in regulating the metabolic switch between the carbon sources methanol and methylamine. Moreover, biochemical studies demonstrated that the 5′ end of sRNA₁₆₂ targets the 5′-untranslated region of the cis-encoded MM2442 mRNA. Overall, this first study of archaeal sRNA/mRNA-target interactions unraveled that sRNA₁₆₂ acts as an antisense (as)RNA on cis- and trans-encoded mRNAs via two distinct domains, indicating that cis-encoded asRNAs can have larger target regulons than previously anticipated.

Functional genomic and advanced genetic studies reveal novel insights into the metabolism, regulation, and biology of Haloferax volcanii

The genome sequence of Haloferax volcanii is available and several comparative genomic in silico studies were performed that yielded novel insight for example into protein export, RNA modifications, small non-coding RNAs, and ubiquitin-like Small Archaeal Modifier Proteins. The full range of functional genomic methods has been established and results from transcriptomic, proteomic and metabolomic studies are discussed. Notably, Hfx. volcanii is together with Halobacterium salinarum the only prokaryotic species for which a translatome analysis has been performed. The results revealed that the fraction of translationally-regulated genes in haloarchaea is as high as in eukaryotes. A highly efficient genetic system has been established that enables the application of libraries as well as the parallel generation of genomic deletion mutants. Facile mutant generation is complemented by the possibility to culture Hfx. volcanii in microtiter plates, allowing the phenotyping of mutant collections. Genetic approaches are currently used to study diverse biological questions–from replication to posttranslational modification—and selected results are discussed. Taken together, the wealth of functional genomic and genetic tools make Hfx. volcanii a bona fide archaeal model species, which has enabled the generation of important results in recent years and will most likely generate further breakthroughs in the future.

Translation initiation: variations in the mechanism can be anticipated

Translation initiation is a critical step in protein synthesis. Previously, two major mechanisms of initiation were considered as essential: prokaryotic, based on SD interaction; and eukaryotic, requiring cap structure and ribosomal scanning. Although discovered decades ago, cap-independent translation has recently been acknowledged as a widely spread mechanism in viruses, which may take place in some cellular mRNA translations. Moreover, it has become evident that translation can be initiated on the leaderless mRNA in all three domains of life. New findings demonstrate that other distinguishable types of initiation exist, including SD-independent in Bacteria and Archaea, and various modifications of 5' end-dependent and internal initiation mechanisms in Eukarya. Since translation initiation has developed through the loss, acquisition, and modification of functional elements, all of which have been elevated by competition with viral translation in a large number of organisms of different complexity, more variation in initiation mechanisms can be anticipated.

The complete genome sequence of Haloferax volcanii DS2, a model archaeon

BACKGROUND

Haloferax volcanii is an easily culturable moderate halophile that grows on simple defined media, is readily transformable, and has a relatively stable genome. This, in combination with its biochemical and genetic tractability, has made Hfx. volcanii a key model organism, not only for the study of halophilicity, but also for archaeal biology in general.

METHODOLOGY/PRINCIPAL FINDINGS

We report here the sequencing and analysis of the genome of Hfx. volcanii DS2, the type strain of this species. The genome contains a main 2.848 Mb chromosome, three smaller chromosomes pHV1, 3, 4 (85, 438, 636 kb, respectively) and the pHV2 plasmid (6.4 kb).

CONCLUSIONS/SIGNIFICANCE

The completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

Deep sequencing analysis of the Methanosarcina mazei Gö1 transcriptome in response to nitrogen availability

Methanosarcina mazei and related mesophilic archaea are the only organisms fermenting acetate, methylamines, and methanol to methane and carbon dioxide, contributing significantly to greenhouse gas production. The biochemistry of these metabolic processes is well studied, and genome sequences are available, yet little is known about the overall transcriptional organization and the noncoding regions representing 25% of the 4.01-Mb genome of M. mazei. We present a genome-wide analysis of transcription start sites (TSS) in M. mazei grown under different nitrogen availabilities. Pyrosequencing-based differential analysis of primary vs. processed 5' ends of transcripts discovered 876 TSS across the M. mazei genome. Unlike in other archaea, in which leaderless mRNAs are prevalent, the majority of the detected mRNAs in M. mazei carry long untranslated 5' regions. Our experimental data predict a total of 208 small RNA (sRNA) candidates, mostly from intergenic regions but also antisense to 5' and 3' regions of mRNAs. In addition, 40 new small mRNAs with ORFs of < or = 30 aa were identified, some of which might have dual functions as mRNA and regulatory sRNA. We confirmed differential expression of several sRNA genes in response to nitrogen availability. Inspection of their promoter regions revealed a unique conserved sequence motif associated with nitrogen-responsive regulation, which might serve as a regulator binding site upstream of the common IIB recognition element. Strikingly, several sRNAs antisense to mRNAs encoding transposases indicate nitrogen-dependent transposition events. This global TSS map in archaea will facilitate a better understanding of transcriptional and posttranscriptional control in the third domain of life.

Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity

Transcriptome-wide studies in eukaryotes have been instrumental in the characterization of fundamental regulatory mechanisms for more than a decade. By contrast, in prokaryotes (bacteria and archaea) whole-transcriptome studies have not been performed until recently owing to the general view that microbial gene structures are simple, as well as technical difficulties in enriching for mRNAs that lack poly(A) tails. Deep RNA sequencing and tiling array studies are now revolutionizing our understanding of the complexity, plasticity and regulation of microbial transcriptomes.

Proteins participating to the post-transcriptional regulation of the mitochondrial cytochrome c oxidase subunit IV via elements located in the 3'UTR

In developing rat brain cytochrome c oxidase subunit IV (COXIV) expression is also regulated at post-transcriptional level and two 3'UTR-COXIV RNA-binding factors have been identified. Here, we report the enrichment and identification of the factors from just born rat brains by affinity chromatography of biotinylated 3'UTR-COXIV RNA-protein complexes on streptavidin-conjugated paramagnetic particles. We successfully isolated two main proteins of about 86 and 42kDa, whose sequences were highly attributable to Hsp90 and Actin. The purified proteins maintain RNA-binding ability and specificity for COXIV messenger and, interacting with the 3'UTR, then could negatively modulate mRNA translation. We also studied the content of Hsp90 and Actin during postnatal brain development and demonstrated that in just born rat brain, when the COXIV protein appears at low level, Hsp90 was not phosphorylated. Vice versa in the adult tissue, when COXIV accumulates, Hsp90 appeared phosphorylated in serine therefore it could be unable to bind COXIV messenger, suggesting that the phosphorylation event could provoke the loss of Hsp90 binding to mRNA. We hypothesize a new post-transcriptional mechanism regulating a messenger encoded by nuclear genome for a mitochondrial protein and that Hsp90 and Actin, could represent key players in COXIV translation.

Begin at the beginning: evolution of translational initiation

Initiation of protein synthesis, entailing ribosomal recognition of the mRNA start codon and setting of the correct reading frame, is the rate-limiting step in translation and the main target of translation regulation in all modern cells. As efficient selection of the translation start site is vital for survival of extant cells, a mechanism for ensuring this may already have been in existence in the last universal common ancestor of present-day cells. This article reviews known features of the molecular machinery for initiation in the primary domains of life, Bacteria, Archaea and Eukarya, and attempts to identify conserved features that may be useful for reconstructing a model of the ancestral initiation apparatus.