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

Proteomic Sample Preparation and Data Analysis in Line with the Archaeal Proteome Project

Despite the ecological, evolutionary and economical significance of archaea, key aspects of their cell biology, metabolic pathways, and adaptations to a wide spectrum of environmental conditions, remain to be elucidated. Proteomics allows for the system-wide analysis of proteins, their changes in abundance between different conditions, as well as their post-translational modifications, providing detailed insights into the function of proteins and archaeal cell biology. In this chapter, we describe a sample preparation and mass spectrometric analysis workflow that has been designed for Haloferax volcanii but can be applied to a broad range of archaeal species. Furthermore, proteomics experiments provide a wealth of data that is invaluable to various disciplines. Therefore, we previously initiated the Archaeal Proteome Project (ArcPP), a community project that combines the analysis of multiple datasets with expert knowledge in various fields of archaeal research. The corresponding bioinformatic analysis, allowing for the integration of new proteomics data into the ArcPP, as well as the interactive exploration of ArcPP results is also presented here. In combination, these protocols facilitate an optimized, detailed and collaborative approach to archaeal proteomics.

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.

Open Issues for Protein Function Assignment in Haloferax volcanii and Other Halophilic Archaea

Background: Annotation ambiguities and annotation errors are a general challenge in genomics. While a reliable protein function assignment can be obtained by experimental characterization, this is expensive and time-consuming, and the number of such Gold Standard Proteins (GSP) with experimental support remains very low compared to proteins annotated by sequence homology, usually through automated pipelines. Even a GSP may give a misleading assignment when used as a reference: the homolog may be close enough to support isofunctionality, but the substrate of the GSP is absent from the species being annotated. In such cases, the enzymes cannot be isofunctional. Here, we examined a variety of such issues in halophilic archaea (class Halobacteria), with a strong focus on the model haloarchaeon Haloferax volcanii. Results: Annotated proteins of Hfx. volcanii were identified for which public databases tend to assign a function that is probably incorrect. In some cases, an alternative, probably correct, function can be predicted or inferred from the available evidence, but this has not been adopted by public databases because experimental validation is lacking. In other cases, a probably invalid specific function is predicted by homology, and while there is evidence that this assigned function is unlikely, the true function remains elusive. We listed 50 of those cases, each with detailed background information, so that a conclusion about the most likely biological function can be drawn. For reasons of brevity and comprehension, only the key aspects are listed in the main text, with detailed information being provided in a corresponding section of the Supplementary Materials. Conclusions: Compiling, describing and summarizing these open annotation issues and functional predictions will benefit the scientific community in the general effort to improve the evaluation of protein function assignments and more thoroughly detail them. By highlighting the gaps and likely annotation errors currently in the databases, we hope this study will provide a framework for experimentalists to systematically confirm (or disprove) our function predictions or to uncover yet more unexpected functions.

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.

Transcriptional regulation of dimethyl sulfoxide respiration in a haloarchaeon, Haloferax volcanii

The halophilic euryarchaeon Haloferax volcanii can grow anaerobically by DMSO respiration. DMSO reductase was induced by DMSO respiration not only under anaerobic growth conditions but also in denitrifying cells of H. volcanii. Deletion of the dmsR gene, encoding a putative regulator for the DMSO reductase, resulted in the loss of anaerobic growth by DMSO respiration. Reporter experiments revealed that only the anaerobic condition was essential for transcription of the dmsEABCD genes encoding DMSO reductase and that transcription was enhanced threefold by supplementation of DMSO. In the ∆dmsR mutant, transcription of the dmsEABCD genes induced by the anaerobic condition was not enhanced by DMSO, suggesting that DmsR is a DMSO-responsive regulator. Transcriptions of the dmsR and mgd genes for Mo-bisMGD biosynthesis were regulated in the same manner as the dmsEABCD genes. These results suggest that the genetic regulation of DMSO respiration in H. volcanii is controlled by at least two systems: one is the DMSO-responsive DmsR, and the other is an unknown anaerobic regulator.

Transcriptional profiles of Haloferax mediterranei based on nitrogen availability

The haloarchaeon Haloferax mediterranei is able to grow in the presence of different inorganic and organic nitrogen sources by means of the assimilatory pathway under aerobic conditions. In order to identify genes of potential importance in nitrogen metabolism and its regulation in the halophilic microorganism, we have analysed its global gene expression in three culture media with different nitrogen sources: (a) cells were grown stationary and exponentially in ammonium, (b) cells were grown exponentially in nitrate, and (c) cells were shifted to nitrogen starvation conditions. The main differences in the transcriptional profiles have been identified between the cultures with ammonium as nitrogen source and the cultures with nitrate or nitrogen starvation, supporting previous results which indicate the absence of ammonium as the factor responsible for the expression of genes involved in nitrate assimilation pathway. The results have also permitted the identification of transcriptional regulators and changes in metabolic pathways related to the catabolism and anabolism of amino acids or nucleotides. The microarray data was validated by real-time quantitative PCR on 4 selected genes involved in nitrogen metabolism. This work represents the first transcriptional profiles study related to nitrogen assimilation metabolism in extreme halophilic microorganisms using microarray technology.

Phylogenetically driven sequencing of extremely halophilic archaea reveals strategies for static and dynamic osmo-response

Organisms across the tree of life use a variety of mechanisms to respond to stress-inducing fluctuations in osmotic conditions. Cellular response mechanisms and phenotypes associated with osmoadaptation also play important roles in bacterial virulence, human health, agricultural production and many other biological systems. To improve understanding of osmoadaptive strategies, we have generated 59 high-quality draft genomes for the haloarchaea (a euryarchaeal clade whose members thrive in hypersaline environments and routinely experience drastic changes in environmental salinity) and analyzed these new genomes in combination with those from 21 previously sequenced haloarchaeal isolates. We propose a generalized model for haloarchaeal management of cytoplasmic osmolarity in response to osmotic shifts, where potassium accumulation and sodium expulsion during osmotic upshock are accomplished via secondary transport using the proton gradient as an energy source, and potassium loss during downshock is via a combination of secondary transport and non-specific ion loss through mechanosensitive channels. We also propose new mechanisms for magnesium and chloride accumulation. We describe the expansion and differentiation of haloarchaeal general transcription factor families, including two novel expansions of the TATA-binding protein family, and discuss their potential for enabling rapid adaptation to environmental fluxes. We challenge a recent high-profile proposal regarding the evolutionary origins of the haloarchaea by showing that inclusion of additional genomes significantly reduces support for a proposed large-scale horizontal gene transfer into the ancestral haloarchaeon from the bacterial domain. The combination of broad (17 genera) and deep (≥5 species in four genera) sampling of a phenotypically unified clade has enabled us to uncover both highly conserved and specialized features of osmoadaptation. Finally, we demonstrate the broad utility of such datasets, for metagenomics, improvements to automated gene annotation and investigations of evolutionary processes.

The ability to adjust to changing osmotic conditions (osmoadaptation) is crucial to the survival of organisms across the tree of life. However, significant gaps still exist in our understanding of this important phenomenon. To help fill some of these gaps, we have produced high-quality draft genomes for 59 osmoadaptation “experts” (extreme halophiles of the euryarchaeal family Halobacteriaceae). We describe the dispersal of osmoadaptive protein families across the haloarchaeal evolutionary tree. We use this data to suggest a generalized model for haloarchaeal ion transport in response to changing osmotic conditions, including proposed new mechanisms for magnesium and chloride accumulation. We describe the evolutionary expansion and differentiation of haloarchaeal general transcription factor families and discuss their potential for enabling rapid adaptation to environmental fluxes. Lastly, we challenge a recent high-profile proposal regarding the evolutionary origins of the haloarchaea by showing that inclusion of additional genomes significantly reduces support for a proposed large-scale horizontal gene transfer into the ancestral haloarchaeon from the bacterial domain. This result highlights the power of our dataset for making evolutionary inferences, a feature which will make it useful to the broader evolutionary community. We distribute our genomic dataset through a user-friendly graphical interface.

A simple technique to improve the resolution of membrane acidic proteins of the haloarchaeon Haloferax volcanii by 2D electrophoresis

Proteins present in the archaeal cell envelope play key roles in a variety of processes necessary for survival in extreme environments. The haloarchaeon Haloferax volcanii is a good model for membrane proteomic studies because its genome sequence is known, it can be genetically manipulated, and a number of studies at the "omics" level have been performed in this organism. This work reports an easy strategy to improve the resolution of acidic membrane proteins from H. volcanii by 2DE. The method is based on the solubilization, delipidation, and salt removal from membrane proteins. Due to the abundance of the S-layer glycoprotein (SLG) in membrane protein extracts, other proteins from the envelope are consequently underrepresented. Thus, a protocol to reduce the amount of the SLG by EDTA treatment was applied and 11 cm narrow range pH (3.9-5.1) IPG strips were used to fractionate the remaining proteins. Using this method, horizontal streaking was substantially decreased and at least 75 defined spots (20% of the predicted membrane proteome within this pI/Mw range) were reproducibly detected. Two of these spots were identified as thermosome subunit 1 and NADH dehydrogenase from H. volcanii, confirming that proteins from the membrane fraction were enriched. Removal of the SLG from membrane protein extracts can be applied to increase protein load for 2DE as well as for other proteomic methods.

The effects of extremes of pH on the growth and transcriptomic profiles of three haloarchaea

The halophilic archaea (haloarchaea) live in saline environments, which are found across the globe.  In addition to salinity, these niches can be quite dynamic and experience extreme conditions such as low oxygen content, radiation (gamma and UV), pH and temperature.  However, of all the naturally occurring stresses faced by the haloarchaea, only one, pH, has not been previously investigated in regard to the changes induced in the transcriptome. Therefore, we endeavored to determine the responses in three haloarchaea: Halorubrum lacusprofundi (Hla), Haloferax volcanii (Hvo), and Halobacterium sp. NRC-1 (NRC-1) to growth under acidic and alkaline pH. Our observations showed that the transcriptomes of Hvo and NRC-1 regulated stress, motility, and ABC transporters in a similar manner, which is in line with previous reports from other prokaryotes when grown in an acidic environment.  However, the pattern for Hla was more species specific. For alkaline stress, all three haloarchaea responded in a manner similar to well-studied archaea and bacteria showing the haloarchaeal response was general to prokaryotes. Additionally, we performed an analysis on the changes in the transcriptomes of the three haloarchaea when shifting from one pH extreme to the other. The results showed that the transcriptomes of all three haloarchaea respond more similarly when moving from alkaline to acidic conditions compared to a shift in the opposite direction. Interestingly, our studies also showed that individual genes of multiple paralogous gene families ( tbp, tfb, orc/ cdc6, etc.) found in the haloarchaea were regulated under specific stresses thereby providing evidence that they modulate the response to various environmental stresses. The studies described here are the first to catalog the changes in the haloarchaeal transcriptomes under growth in extreme pH and help us understand how life is able to thrive under all conditions present on Earth and, if present, on extraterrestrial bodies as well.

Functional implications of ribosomal RNA methylation in response to environmental stress

The study of post-transcriptional RNA modifications has long been focused on the roles these chemical modifications play in maintaining ribosomal function. The field of ribosomal RNA modification has reached a milestone in recent years with the confirmation of the final unknown ribosomal RNA methyltransferase in Escherichia coli in 2012. Furthermore, the last 10 years have brought numerous discoveries in non-coding RNAs and the roles that post-transcriptional modification play in their functions. These observations indicate the need for a revitalization of this field of research to understand the role modifications play in maintaining cellular health in a dynamic environment. With the advent of high-throughput sequencing technologies, the time is ripe for leaps and bounds forward. This review discusses ribosomal RNA methyltransferases and their role in responding to external stress in Escherichia coli, with a specific focus on knockout studies and on analysis of transcriptome data with respect to rRNA methyltransferases.

Overview of the genetic tools in the Archaea

This section provides an overview of the genetic systems developed in the Archaea. Genetic manipulation is possible in many members of the halophiles, methanogens, Sulfolobus, and Thermococcales. We describe the selection/counterselection principles utilized in each of these groups, which consist of antibiotics and their resistance markers, and auxotrophic host strains and complementary markers. The latter strategy utilizes techniques similar to those developed in yeast. However, Archaea are resistant to many of the antibiotics routinely used for selection in the Bacteria, and a number of strategies specific to the Archaea have been developed. In addition, examples utilizing the genetic systems developed for each group will be briefly described.

High throughput sequencing reveals a plethora of small RNAs including tRNA derived fragments in Haloferax volcanii

To define the complete sRNA population of the halophilic archaeon Haloferax volcanii, we employed high throughput sequencing. cDNAs were generated from RNA ranging in size from 17 to 500 nucleotides isolated from cells grown at three different conditions to exponential and stationary phase, respectively. Altogether, 145 intergenic and 45 antisense sRNAs were identified. Comparison of the expression profile showed different numbers of reads at the six different conditions for the majority of sRNAs. A striking difference in the number of sRNA reads was observed between cells grown under standard vs. low salt conditions. Furthermore, the six highest numbers of reads were found for low salt conditions. In contrast, only slight differences between sRNA reads at different growth temperatures were detected. Attempts to delete four sRNA genes revealed that one sRNA gene is essential. The three viable sRNA gene deletion mutants possessed distinct phenotypes. According to microarray analyses, the removal of the sRNA gene resulted in a profound change of the transcriptome when compared with the wild type. High throughput sequencing also showed the presence of high concentrations of tRNA derived fragments in H. volcanii. These tRF molecules were shown to have different amounts of reads at the six conditions analyzed. Northern analysis was used to confirm the presence of the tRNA-derived fragments.