Characterization of a Haloferax volcanii member of the enolase superfamily: deletion mutant construction, expression analysis, and transcriptome comparison

C(P)XCG Proteins of Haloferax volcanii with Predicted Zinc Finger Domains: The Majority Bind Zinc, but Several Do Not

In recent years, interest in very small proteins (µ-proteins) has increased significantly, and they were found to fulfill important functions in all prokaryotic and eukaryotic species. The halophilic archaeon Haloferax volcanii encodes about 400 µ-proteins of less than 70 amino acids, 49 of which contain at least two C(P)XCG motifs and are, thus, predicted zinc finger proteins. The determination of the NMR solution structure of HVO_2753 revealed that only one of two predicted zinc fingers actually bound zinc, while a second one was metal-free. Therefore, the aim of the current study was the homologous production of additional C(P)XCG proteins and the quantification of their zinc content. Attempts to produce 31 proteins failed, underscoring the particular difficulties of working with µ-proteins. In total, 14 proteins could be produced and purified, and the zinc content was determined. Only nine proteins complexed zinc, while five proteins were zinc-free. Three of the latter could be analyzed using ESI-MS and were found to contain another metal, most likely cobalt or nickel. Therefore, at least in haloarchaea, the variability of predicted C(P)XCG zinc finger motifs is higher than anticipated, and they can be metal-free, bind zinc, or bind another metal. Notably, AlphaFold2 cannot correctly predict whether or not the four cysteines have the tetrahedral configuration that is a prerequisite for metal binding.

Identification and characterization of a novel type of ketohexokinase from the haloarchaeon Haloferax volcanii

Ketohexokinase (KHK) catalyzes the ATP-dependent phosphorylation of fructose, forming fructose-1-phosphate and ADP. The enzyme is well studied in Eukarya, in particular in humans and other vertebrates, but homologs have not been identified in Bacteria and Archaea. Here we report the identification of a novel type of KHK from the haloarchaeon Haloferax volcanii (HvKHK). The encoding gene khk was identified as HVO_1812. The gene was expressed as a 90-kDa homodimeric protein, catalyzing the phosphorylation of fructose with a Vmax value of 59 U/mg and apparent KM values for ATP and fructose of 0.47 and 1.29 mM, respectively. Homologs of HvKHK were only identified in a few haloarchaea and halophilic Bacteria. The protein showed low sequence identity to characterized KHKs from Eukarya and phylogenetic analyses indicate that haloarchaeal KHKs are largely separated from eukaryal KHKs. This is the first report of the identification of KHKs in prokaryotes that form a novel cluster of sugar kinases within the ribokinase/pfkB superfamily.

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.

Characterization of Non-selected Intermolecular Gene Conversion in the Polyploid Haloarchaeon Haloferax volcanii

Gene conversion is defined as the non-reciprocal transfer of genetic information from one site to a homologous, but not identical site of the genome. In prokaryotes, gene conversion can increase the variance of sequences, like in antigenic variation, but can also lead to a homogenization of sequences, like in the concerted evolution of multigene families. In contrast to these intramolecular mechanisms, the intermolecular gene conversion in polyploid prokaryotes, which leads to the equalization of the multiple genome copies, has hardly been studied. We have previously shown the intermolecular gene conversion in halophilic and methanogenic archaea is so efficient that it can be studied without selecting for conversion events. Here, we have established an approach to characterize unselected intermolecular gene conversion in Haloferax volcanii making use of two genes that encode enzymes involved in carotenoid biosynthesis. Heterozygous strains were generated by protoplast fusion, and gene conversion was quantified by phenotype analysis or/and PCR. It was verified that unselected gene conversion is extremely efficient and it was shown that gene conversion tracts are much longer than in antigenic variation or concerted evolution in bacteria. Two sites were nearly always co-converted when they were 600 bp apart, and more than 30% co-conversion even occurred when two sites were 5 kbp apart. The gene conversion frequency was independent from the extent of genome differences, and even a one nucleotide difference triggered conversion.

Glucose Metabolism and Acetate Switch in Archaea: the Enzymes in Haloferax volcanii

The halophilic archaeon Haloferax volcanii has been proposed to degrade glucose via the semiphosphorylative Entner-Doudoroff (spED) pathway. Following our previous studies on key enzymes of this pathway, we now focus on the characterization of enzymes involved in 3-phosphoglycerate conversion to pyruvate, in anaplerosis, and in acetyl coenzyme A (acetyl-CoA) formation from pyruvate. These enzymes include phosphoglycerate mutase, enolase, pyruvate kinase, phosphoenolpyruvate carboxylase, and pyruvate-ferredoxin oxidoreductase. The essential function of these enzymes were shown by transcript analyses and growth experiments with respective deletion mutants. Furthermore, we show that H. volcanii-during aerobic growth on glucose-excreted significant amounts of acetate, which was consumed in the stationary phase (acetate switch). The enzyme catalyzing the conversion of acetyl-CoA to acetate as part of the acetate overflow mechanism, an ADP-forming acetyl-CoA synthetase (ACD), was characterized. The functional involvement of ACD in acetate formation and of AMP-forming acetyl-CoA synthetases (ACSs) in activation of excreted acetate was proven by using respective deletion mutants. Together, the data provide a comprehensive analysis of enzymes of the spED pathway and of anaplerosis and report the first genetic evidence of the functional involvement of enzymes of the acetate switch in archaea.IMPORTANCE In this work, we provide a comprehensive analysis of glucose degradation via the semiphosphorylative Entner-Doudoroff pathway in the haloarchaeal model organism Haloferax volcanii The study includes transcriptional analyses, growth experiments with deletion mutants. and characterization of all enzymes involved in the conversion of 3-phosphoglycerate to acetyl coenzyme A (acetyl-CoA) and in anaplerosis. Phylogenetic analyses of several enzymes indicate various lateral gene transfer events from bacteria to haloarchaea. Furthermore, we analyzed the key players involved in the acetate switch, i.e., in the formation (overflow) and subsequent consumption of acetate during aerobic growth on glucose. Together, the data provide novel aspects of glucose degradation, anaplerosis, and acetate switch in H. volcanii and thus expand our understanding of the unusual sugar metabolism in archaea.

Several One-Domain Zinc Finger µ-Proteins of Haloferax Volcanii Are Important for Stress Adaptation, Biofilm Formation, and Swarming

Zinc finger domains are highly structured and can mediate interactions to DNA, RNA, proteins, lipids, and small molecules. Accordingly, zinc finger proteins are very versatile and involved in many biological functions. Eukaryotes contain a wealth of zinc finger proteins, but zinc finger proteins have also been found in archaea and bacteria. Large zinc finger proteins have been well studied, however, in stark contrast, single domain zinc finger µ-proteins of less than 70 amino acids have not been studied at all, with one single exception. Therefore, 16 zinc finger µ-proteins of the haloarchaeon Haloferax volcanii were chosen and in frame deletion mutants of the cognate genes were generated. The phenotypes of mutants and wild-type were compared under eight different conditions, which were chosen to represent various pathways and involve many genes. None of the mutants differed from the wild-type under optimal or near-optimal conditions. However, 12 of the 16 mutants exhibited a phenotypic difference under at least one of the four following conditions: Growth in synthetic medium with glycerol, growth in the presence of bile acids, biofilm formation, and swarming. In total, 16 loss of function and 11 gain of function phenotypes were observed. Five mutants indicated counter-regulation of a sessile versus a motile life style in H. volcanii. In conclusion, the generation and analysis of a set of deletion mutants demonstrated the high importance of zinc finger µ-proteins for various biological functions, and it will be the basis for future mechanistic insight.

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.

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.

Key Enzymes of the Semiphosphorylative Entner-Doudoroff Pathway in the Haloarchaeon Haloferax volcanii: Characterization of Glucose Dehydrogenase, Gluconate Dehydratase, and 2-Keto-3-Deoxy-6-Phosphogluconate Aldolase

The halophilic archaeon Haloferax volcanii has been proposed to degrade glucose via the semiphosphorylative Entner-Doudoroff (spED) pathway. So far, the key enzymes of this pathway, glucose dehydrogenase (GDH), gluconate dehydratase (GAD), and 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase (KDPGA), have not been characterized, and their functional involvement in glucose degradation has not been demonstrated. Here we report that the genes HVO_1083 and HVO_0950 encode GDH and KDPGA, respectively. The recombinant enzymes show high specificity for glucose and KDPG and did not convert the corresponding C4 epimers galactose and 2-keto-3-deoxy-6-phosphogalactonate at significant rates. Growth studies of knockout mutants indicate the functional involvement of both GDH and KDPGA in glucose degradation. GAD was purified from H. volcanii, and the encoding gene, gad, was identified as HVO_1488. GAD catalyzed the specific dehydration of gluconate and did not utilize galactonate at significant rates. A knockout mutant of GAD lost the ability to grow on glucose, indicating the essential involvement of GAD in glucose degradation. However, following a prolonged incubation period, growth of the Δgad mutant on glucose was recovered. Evidence is presented that under these conditions, GAD was functionally replaced by xylonate dehydratase (XAD), which uses both xylonate and gluconate as substrates. Together, the characterization of key enzymes and analyses of the respective knockout mutants present conclusive evidence for the in vivo operation of the spED pathway for glucose degradation in H. volcanii

IMPORTANCE

The work presented here describes the identification and characterization of the key enzymes glucose dehydrogenase, gluconate dehydratase, and 2-keto-3-deoxy-6-phosphogluconate aldolase and their encoding genes of the proposed semiphosphorylative Entner-Doudoroff pathway in the haloarchaeon Haloferax volcanii The functional involvement of the three enzymes was proven by analyses of the corresponding knockout mutants. These results provide evidence for the in vivo operation of the semiphosphorylative Entner-Doudoroff pathway in haloarchaea and thus expand our understanding of the unusual sugar degradation pathways in the domain Archaea.

The oxidative pentose phosphate pathway in the haloarchaeon Haloferax volcanii involves a novel type of glucose-6-phosphate dehydrogenase--The archaeal Zwischenferment

The oxidative pentose phosphate pathway (OPPP), catalyzing the oxidation of glucose-6-phosphate to ribulose-5-phosphate is ubiquitous in eukarya and bacteria but has not yet been reported in archaea. In haloarchaea a putative 6-phosphogluconate dehydrogenase (6PGDH) is annotated, whereas a gene coding for glucose-6-phosphate dehydrogenase (Glc6PDH) could not be identified. Here we report the purification and characterization of a novel type of Glc6PDH in Haloferax volcanii that is not related to bacterial and eukaryal Glc6PDHs and the encoding gene is designated as azf (archaeal zwischenferment). Further, recombinant H. volcanii 6PGDH was characterized. Deletion mutant analyses indicate that both, Glc6PDH and 6PGDH, are functionally involved in pentose phosphate formation in vivo. This is the first report on the operation of the OPPP in the domain of archaea.

A comprehensive analysis of the importance of translation initiation factors for Haloferax volcanii applying deletion and conditional depletion mutants

Translation is an important step in gene expression. The initiation of translation is phylogenetically diverse, since currently five different initiation mechanisms are known. For bacteria the three initiation factors IF1 – IF3 are described in contrast to archaea and eukaryotes, which contain a considerably higher number of initiation factor genes. As eukaryotes and archaea use a non-overlapping set of initiation mechanisms, orthologous proteins of both domains do not necessarily fulfill the same function. The genome of Haloferax volcanii contains 14 annotated genes that encode (subunits of) initiation factors. To gain a comprehensive overview of the importance of these genes, it was attempted to construct single gene deletion mutants of all genes. In 9 cases single deletion mutants were successfully constructed, showing that the respective genes are not essential. In contrast, the genes encoding initiation factors aIF1, aIF2γ, aIF5A, aIF5B, and aIF6 were found to be essential. Factors aIF1A and aIF2β are encoded by two orthologous genes in H. volcanii. Attempts to generate double mutants failed in both cases, indicating that also these factors are essential. A translatome analysis of one of the single aIF2β deletion mutants revealed that the translational efficiency of the second ortholog was enhanced tenfold and thus the two proteins can replace one another. The phenotypes of the single deletion mutants also revealed that the two aIF1As and aIF2βs have redundant but not identical functions. Remarkably, the gene encoding aIF2α, a subunit of aIF2 involved in initiator tRNA binding, could be deleted. However, the mutant had a severe growth defect under all tested conditions. Conditional depletion mutants were generated for the five essential genes. The phenotypes of deletion mutants and conditional depletion mutants were compared to that of the wild-type under various conditions, and growth characteristics are discussed.

Fructose degradation in the haloarchaeon Haloferax volcanii involves a bacterial type phosphoenolpyruvate-dependent phosphotransferase system, fructose-1-phosphate kinase, and class II fructose-1,6-bisphosphate aldolase

The halophilic archaeon Haloferax volcanii utilizes fructose as a sole carbon and energy source. Genes and enzymes involved in fructose uptake and degradation were identified by transcriptional analyses, deletion mutant experiments, and enzyme characterization. During growth on fructose, the gene cluster HVO_1495 to HVO_1499, encoding homologs of the five bacterial phosphotransferase system (PTS) components enzyme IIB (EIIB), enzyme I (EI), histidine protein (HPr), EIIA, and EIIC, was highly upregulated as a cotranscript. The in-frame deletion of HVO_1499, designated ptfC (ptf stands for phosphotransferase system for fructose) and encoding the putative fructose-specific membrane component EIIC, resulted in a loss of growth on fructose, which could be recovered by complementation in trans. Transcripts of HVO_1500 (pfkB) and HVO_1494 (fba), encoding putative fructose-1-phosphate kinase (1-PFK) and fructose-1,6-bisphosphate aldolase (FBA), respectively, as well as 1-PFK and FBA activities were specifically upregulated in fructose-grown cells. pfkB and fba knockout mutants did not grow on fructose, whereas growth on glucose was not inhibited, indicating the functional involvement of both enzymes in fructose catabolism. Recombinant 1-PFK and FBA obtained after homologous overexpression were characterized as having kinetic properties indicative of functional 1-PFK and a class II type FBA. From these data, we conclude that fructose uptake in H. volcanii involves a fructose-specific PTS generating fructose-1-phosphate, which is further converted via fructose-1,6-bisphosphate to triose phosphates by 1-PFK and FBA. This is the first report of the functional involvement of a bacterial-like PTS and of class II FBA in the sugar metabolism of archaea.

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.

Hot Transcriptomics

DNA microarray technology allows for a quick and easy comparison of complete transcriptomes, resulting in improved molecular insight in fluctuations of gene expression. After emergence of the microarray technology about a decade ago, the technique has now matured and has become routine in many molecular biology laboratories. Numerous studies have been performed that have provided global transcription patterns of many organisms under a wide range of conditions. Initially, implementation of this high-throughput technology has lead to high expectations for ground breaking discoveries. Here an evaluation is performed of the insight that transcriptome analysis has brought about in the field of hyperthermophilic archaea. The examples that will be discussed have been selected on the basis of their impact, in terms of either biological insight or technological progress.

Genome information management and integrated data analysis with HaloLex

HaloLex is a software system for the central management, integration, curation, and web-based visualization of genomic and other -omics data for any given microorganism. The system has been employed for the manual curation of three haloarchaeal genomes, namely Halobacterium salinarum (strain R1), Natronomonas pharaonis, and Haloquadratum walsbyi. HaloLex, in particular, enables the integrated analysis of genome-wide proteomic results with the underlying genomic data. This has proven indispensable to generate reliable gene predictions for GC-rich genomes, which, due to their characteristically low abundance of stop codons, are known to be hard targets for standard gene finders, especially concerning start codon assignment. The proteomic identification of more than 600 N-terminal peptides has greatly increased the reliability of the start codon assignment for Halobacterium salinarum. Application of homology-based methods to the published genome of Haloarcula marismortui allowed to detect 47 previously unidentified genes (a problem that is particularly serious for short protein sequences) and to correct more than 300 start codon misassignments.

Genomics and functional genomics with haloarchaea

The first haloarchaeal genome was published in 2000 and today five genome sequences are available. Transcriptome and proteome analyses have been established for two and three haloarchaeal species, respectively, and more than 20 studies using these functional genomic approaches have been published in the last two years. These studies gave global overviews of metabolic regulation (aerobic and anaerobic respiration, phototrophy, carbon source usage), stress response (UV, X-rays, transition metals, osmotic and temperature stress), cell cycle-dependent transcript level regulation, and transcript half-lives. The only translatome analysis available for any prokaryotic species revealed that 10 and 20% of all transcripts are translationally regulated in Haloferax volcanii and Halobacterium salinarum, respectively. Very effective methods for the construction of in frame deletion mutants have been established recently for haloarchaea and are intensively used to unravel the biological roles of genes in this group. Bioinformatic analyses include both cross-genome comparisons as well as integration of genomic data with experimental results. The first systems biology approaches have been performed that used experimental data to construct predictive models of gene expression and metabolism, respectively. In this contribution the current status of genomics, functional genomics, and molecular genetics of haloarchaea is summarized and selected examples are discussed.