The complete genome sequence of Nab. magadii consisted of four replicons (total size 4,443,643 bp). Three of these elements had a GC content of ~61% whereas pNMAG02 had a GC content of 56.82%. A comparison of some of the relevant features of these four elements is shown in Table 1. A BLASTN analysis of pNMAG03 on the NCBI database revealed 99% identity to halovirus φCh1 (58,498 bp; GenBank accession number NC_004084), a bacteriophage-like element isolated from Nab. magadii. Since halovirus φCh1 has already been described elsewhere 303132, the analysis of pNMAG03 was excluded from the scope of the current work. The large chromosome of Nab. magadii contained two genes encoding putative replication factor C-like proteins (Nmag_1868 and 1910). The large chromosome, pNMAG01, and pNMAG02 were predicted to replicate using a conserved archaeal mechanism 33, since each of these replicons contained at least one gene encoding an Orc1/Cdc6 family replication initiation protein. For the large chromosome, the major replication origin was predicted to be at ca 1.9 Mb, located between Orc1 (Nmag_1930) on the forward strand and a three-gene operon on the reverse strand (Nmag_1927-1929). This set of four highly conserved genes was found adjacent to the replication origin in almost all halophilic archaea.

Archaeal genomes can contain a large number of transposable elements and the variety of archaeal insertion sequences is thought to approximate that of bacteria 34. However, most archaeal genomes lack prophage elements 35. Manual curation indicated that the genome of Nab. magadii contained ~36 full-length or truncated genes encoding putative transposases. These insertion sequence elements were scattered throughout the chromosomes and about 20 of these belong to the IS605 OrfB family. The IS605 OrfB transposase (also called IS1341-type transposase) genes were highly diverse, as is typical of halophilic archaea 36. A single IS605 OrfA (Nmag_4105, also called IS200-type transposase) was identified in the genome. Other transposase genes in Nab. magadii include 7 of the broad category IS4 (3 IS7-type and 4 IS9-type), a single IS240-type, and 4 related to ISSod10. The small number of transposase genes and their heterogeneity may indicate that Nab. magadii is only minimally affected by these elements. The genome also contained several genes related to bacteriophage elements (e.g., PhiH1 repressor protein, phage tail proteins, and phage protein D) and a vgr-like gene related to recombination hot spot elements. In addition, there were 13 genes encoding integrase/recombinase-like proteins (Additional file 1: Table S1).

Archaeal genomes generally have 1–4 rRNA operons consisting of the 16 S, 23 S, and 5 S rRNA genes with a tRNA^Ala gene located in the internal transcribed spacer 37. The large chromosome of Nab. magadii contained two copies of 16 S rRNA-tRNA^Ala-23 S rRNA-5 S rRNA sequences, one each on the plus and minus strands, as well as two genes encoding components of the RNA guide machinery (Nmag_0693-0694) with fibrillarin-like RNA methyltransferase as the catalytic component. The small chromosome pNMAG01 contained a copy of 16 S rRNA-tRNA^Ala-23 S rRNA-5 S rRNA sequence on the minus strand and a copy of 23 S rRNA-5 S rRNA sequence on the plus strand. The three 16 S rRNA-tRNA^Ala-23 S rRNA-5 S rRNA sequences of Nab. magadii had 99% nucleotide identity to each other. The small chromosome pNMAG01 also contained an orphan 5 S rRNA sequence that had 89% nucleotide identity to the other four 5 S rRNA genes of Nab. magadii. Since pNMAG02 lacked rRNA operons and had a lesser GC content than the large and small chromosomes, this self replicating element could be considered a large plasmid. The heterogeneity of the rRNA operons within Nab. magadii is not a unique feature and the occurrence of such rRNA operons among halobacterial genomes is thought to be due to recombination between rRNA genes of different strains or species 38. The 16 S rRNA genes of Nab. magadii were closely related to those of Nab. asiatica (97% identity), Htg. turkmenica (96% identity), and Nmn. pharaonis (90% identity). Furthermore, the genome of Nab. magadii was compared to 17 complete haloarchaeal genomes available in the public databases (Additional file 2: Table S2). Based on this analysis, Htg. turkmenica contained the highest number of orthologs (2601 symmetrical hits), followed by Halopiger xanaduensis strain SH-6 (2533 symmetrical hits). There were lesser number of orthologs (1805 symmetrical hits) in Nmn. pharaonis, which has a relatively smaller genome. However, when the data for the percentage of proteins having a bidirectional best blast hit in Nab. magadii was computed (Additional file 2: Table S2), Nmn. pharaonis was the top (63% of the proteins having a bidirectional best blast pair), followed by Hpg. xanaduensis (60%) and Htg. turkmenica (51%). Results from Nmn. pharaonis and Htg. turkmenica are emphasized in this paper since the former was the only other haloalkaliphilic archaeon with a complete genome sequence and the latter contained the highest number of orthologs.

The combined size of the complete genome of Nab. magadii was 1.7 Mb larger than the complete genome of Nmn. pharaonis, which consists of three replicons (total size 2,749,696 bp). However, Nab. magadii genome was 1 Mb smaller than the complete genome of Htg. turkmenica, which consists of seven replicons (total size 5,440,782 bp). The GC content (61.42%) of the large chromosome of Nab. magadii was slightly lesser than that of the large chromosomes of Htg. turkmenica (~66% GC) and Nmn. pharaonis (~63% GC). Alignment of the large chromosome of Nmn. pharaonis (2,595,221 bp) with that of Nab. magadii using MAUVE showed the presence of very few short syntenic regions (data not shown), whereas a similar alignment using the large chromosome of Htg. turkmenica (3,889,038 bp) showed the presence of numerous short syntenic regions (Figure 1A). To further dissect this co-linearity, a BLASTN comparison of the large chromosomes of Nab. magadii and Htg. turkmenica was performed. This analysis revealed the presence of 400 homologous regions (226 plus/plus and 174 plus/minus, sequence range >300 bp, but <2,000 bp) with an average nucleic acid identity of 89% (the identity range was 84%-98%, E-value = 0). The plus and minus strand matches among the chromosomes of Htg. turkmenica and Nab. magadii generated by NUCmer are shown in Figure 1B.

A three-way comparison of all predicted protein-coding genes of Nab. magadii using the TaxPlot tool of NCBI revealed that Htg. turkmenica contained 2387 orthologs, whereas Nmn. pharaonis contained only 426 orthologs (symmetrical hits). These analyses further confirmed that Nab. magadii was more closely related to Htg. turkmenica than to Nmn. pharaonis. In Nab. magadii, 1518 genes could not be assigned a function based on BLAST analysis and were annotated as encoding hypothetical or conserved hypothetical proteins. The isoelectric point (pI) of most of the predicted proteins of Nab. magadii was in the 3–5 range, indicating that the general proteome is acidic, which is typical of most halophilic archaea. A two-way comparison of the large chromosomes revealed that Nab. magadii contained ~945 putative protein-coding genes that had no homologs in Htg. turkmenica. A vast majority (~75%) of these Nab. magadii-specific genes encoded hypothetical proteins. Other genome-specific genes in Nab. magadii encoded ABC-type transporters, ATPases, kinases, phosphatases, proteases, and oxidoreductases. The genome of Nab. magadii also contained a variety of simple sequence repeats encoding characteristic peptide repeat patterns.

In addition to maintaining an acidic proteome and a cell wall composed of acidic glycoproteins, haloalkaliphilic species appear to have evolved several other mechanisms of adaptation to their niche 394041. These include, but are not limited to, intracellular accumulation of inorganic cations and/or neutral organic compounds 424344. Halophilic archaea maintain the necessary water balance and osmotic pressure even when the extracellular Na⁺ concentration exceeds 5 M by pumping Na⁺ out and K⁺ into the cell using a variety of cation/proton antiporters 124546. The genome of Nab. magadii contained an operon of nine genes encoding a putative pH adaptation K⁺ efflux system (Nmag_3445-3453). Genes related to this operon were present in several halophilic archaea, indicating that they may not encode a specialized system involved in stress response to alkaline growth conditions. Apart from this operon, the genome contained three other genes encoding putative cation/proton antiporters and a gene encoding a putative OsmC family protein (Additional file 1: Table S1).

Low molecular weight organic compounds such as amino acids, polyols, and sugars facilitate cellular adaptation to high-osmolarity and are referred to as osmoprotectants or compatible solutes 47. Halophilic species also accumulate neutral organic compounds as a means of adaptation to their niche 144344. The large chromosome of Nab. magadii contained a locus encoding a putative trehalose-phosphate synthase and a trehalose-phosphatase, which may be involved in the biosynthesis of the osmoprotectant 2-sulfotrehalose. The orthologs of these genes were found in few other halophilic archaea and the osmolyte has been detected by nuclear magnetic resonance spectroscopic analysis in Nab. magadii14. Nab. magadii also contained genes encoding the biosynthesis of spermine as well as transporters for the uptake of choline/carnitine/betaine and spermidine/putrescine, which may also provide protection at high-osmolarity (Additional file 1: Table S1). Therefore, it appeared that Nab. magadii had multiple mechanisms (e.g., intracellular accumulation of inorganic cations as well as neutral and/or charged organic compounds) for osmotic adaptation.

The scarcity of molecular oxygen in a hypersaline milieu could be a growth-limiting factor for aerobic chemoorganotrophic prokaryotes 484950. It has been proposed that some archaeal species accumulate intracellular gas vesicles that help them float on the surface and perform oxidative respiration in their native saturated saltwater habitats 12. The large chromosome of Nab. magadii contained a cluster of 11 genes encoding putative gas vesicle synthesis family proteins, which were related to the gas vesicle family proteins of Hbt. salinarum. However, Nab. magadii lacked genes related to those encoding the minor gas vesicle protein (GvpC) and the regulators (GvpD and GvpE). The gas vesicle gene clusters of Nab. magadii and Htg. turkmenica were highly similar to each other and appeared to contain a distant homolog of the Hbt. salinarum gvpI gene (Nmag_0338 and Htur_2370, respectively). Nevertheless, these genes encode putative proteins with an N-terminal extension of more than 200 residues not found in GvpI of Hbt. salinarum. Furthermore, the gas vesicle clusters of Nab. magadii and Htg. turkmenica contained an additional gene (Nmag_0337 and Htur_2371, respectively) that was absent in Hbt. salinarum. Nab. magadii also contained a gene encoding a hemAT-type aerotactic transducer with a putative globin-coupled sensor protein comprising of a globin fold domain and a methyl-accepting chemotaxis transducer domain (Additional file 1: Table S1). These traits, in addition to the osmotic adaptation mechanisms discussed above, may play a role in the survival of Nab. magadii in its natural environment.

Other features likely to facilitate the adaptation of Nab. magadii to its niche included genes encoding putative mechanosensitive ion channels (MscS) that afford protection against hypoosmotic shock, chaperone proteins DnaJ and DnaK, a thermosome, and heat shock proteins that may participate in protein quality control and cellular response to stress. Nab. magadii also contained 47 genes (26 on the large chromosome, 18 on pNMAG01, and 3 on pNMAG02) encoding putative proteins of various sizes with a universal stress protein (UspA) domain. One of these genes (Nmag_1302) appeared to form an operon with a gene encoding a putative GCN5-related N-acetyltransferase (GNAT, Nmag_1303), and a similar gene pair was found in Htg. turkmenica (Htur_3429-3430, 80% identity, E-value = 2e-82 to 1e-52) and Nmn. pharaonis (NP1710A-NP1712A, 49% identity, E-value = 2e-41 to 4e-28). It is possible that the GNAT is involved in the acetylation of the linked universal stress protein in these species. In addition, the large chromosome of Nab. magadii contained genes encoding a superoxide dismutase (sodA), two catalases (katG, which is common to most halophiles, and katE, which is closely related to katE of methanogens, also present in Htg. turkmenica and Hpg. xanaduensis, but not in other halophiles), two alkyl-hydroperoxidase-like proteins, a carbonic anhydrase, and methionine sulfoxide reductases (msrA and msrB; Additional file 1: Table S1). It is possible that these enzymes have a role in the adaptation of this haloarchaeon to various oxidative stresses associated with energy metabolism.

Furthermore, Nab. magadii contained genes encoding metal transport proteins and a putative copper resistance protein. Nab. magadii copper resistance protein appears to contain fused CopC-CopD domains and a distant homolog of this protein occurs in Nmn. pharaonis (NP4610A, 39% identity, E-value = 2e-105), but not in other archaea. These genes may be involved in metal homeostasis in the hypersaline environment that Nab. magadii inhabits. Nab. magadii also encoded DNA methylases, DNA damage repair excinuclease ABC subunits, DNA mismatch repair proteins, and DNA repair/recombination proteins RadA and RadB (Additional file 1: Table S1). Homologs of these genes are found in several other archaea and they are predicted to be involved in stress response and maintaining genetic integrity.

At least 83 genes encoding various types of peptidases/proteases were identified in the genome of Nab. magadii by manual curation (~1.8% of all protein-coding genes; Additional file 3: Table S3). Interestingly, Nab. magadii appears to encode a larger set of proteolytic enzymes compared to most halophilic archaea, including Nmn. pharaonis, Hfx. volcanii and Hbt. salinarum. This suggests that the natural environment inhabited by Nab. magadii contains an ample supply of protein debris, which could be used as a major carbon and nitrogen source. The closest homologs of the vast majority of Nab. magadii genes (~50%) encoding putative peptidases/proteases were found in Htg. turkmenica (Additional file 3: Table S3). Most of the Nab. magadii predicted proteases belong to the catalytic type of metallo- and serine proteases. Other proteases include various amino- and carboxypeptidases, oligopeptidases, signal peptidases, ATP-dependent proteases, and intramembrane cleaving proteases (I-CLiPs).

Subtilases (COG1404; subtilisin-like serine proteases) are a large superfamily of functionally diverse endo- and exo-peptidases that occur in prokaryotes and eukaryotes 51. Nab. magadii contained nine genes encoding putative S8 and S53 subtilisin kexin sedolisins (Additional file 3: Table S3). Although the predicted subtilisins of Nab. magadii had diverse sizes (ranging from 402 to 1710 aa), the amino acid motifs containing the catalytic triad (Asp-His-Ser) were conserved in all of them. Six of the predicted subtilisins (Nmag_0073, 0714, 0715, 1249, 1874, and 3633) of Nab. magadii contained putative targeting signals for translocation through the twin-arginine transport (Tat) pathway, suggesting that these proteases are most likely exported out of the cell. Within this group, Nmag_0715 has been biochemically characterized and designated as the Natrialba extracellular protease (Nep) 16. Nep was demonstrated to be alkali-resistant, a feature that correlates with the conditions that predominate in the natural environment of Nab. magadii16. Interestingly, the C-terminal domain of Nep contains an acidic patch composed of 12 amino acid residues that is absent in the subtilases of neutrophilic organisms 52. This distinctive feature of Nep may be involved in its stability at high salt and/or high pH. In addition, pNMAG01 contained a gene encoding a putative microcystin LR degradation protein (Nmag_3774, MlrC-like-protein). MlrC peptidases, initially isolated from the bacterium Sphingomonas, are a specialized group of metalloproteases assigned to M81 family and they participate in the last step of the degradation pathway of microcystin LR 53. These enzymes rarely occur in the archaeal domain and the homologs of Nmag_3774 were not found in Nmn. pharaonis and Htg. turkmenica (Additional file 3: Table S3).

All archaeal genomes studied to date are predicted to encode self-compartmentalized proteases (20 S proteasomes and Lon-type proteases) likely to function in energy-dependent proteolysis and an ubiquitin-type mechanism for targeting proteins to proteasomes termed sampylation 5455. In archaea, 20 S proteasomes of α- and β-type subunits are thought to function with AAA-ATPases such as the proteasome-activating nucleotidase (PAN) in degrading folded proteins 54. In addition, ubiquitin-like small archaeal modifier proteins (SAMPs) appear to be conjugated to protein targets by an E1-like enzyme termed ubiquitin-like conjugating enzyme of archaea or UbaA (based on study of Hfx. volcanii56). The genome of Nab. magadii contained an operon encoding putative 20 S proteasome α and β subunits (Nmag_0515-0514, respectively). Apart from this operon, the genome contained separate genes encoding 20 S proteasome α and β subunit homologs (Nmag_3313 and Nmag_3351, respectively). Nab. magadii was also predicted to encode homologs of PAN (Nmag_1362 and 2440) and ubiquitin-like small archaeal modifier proteins (SAMPs; Nmag_0567, 1914, 2668, and 2971). The genome of Nab. magadi contained two genes encoding putative ubiquitin-like activating enzymes of archaea (UbaA; Nmag_1394 and 3812). Furthermore, it also encoded a distant homolog of UbaA (Nmag_0356) containing a C-terminal JAB1/MPN/Mov34 metalloenzyme (JAMM) domain that was predicted to remove SAMPs from target proteins. In contrast, Hfx. volcanii encodes only a single UbaA-type protein that functions in both protein conjugation (sampylation) and sulfur mobilization 56. Nab. magadii also encoded an archaeal-type LonB protease (Nmag_2822), which was demonstrated in its cell membranes 57. While LonB homologs are conserved and likely act as key energy-dependent proteases in archaea, the physiological significance of these enzymes has not been addressed.

The tetrahedral aminopeptidase (TET protease) is an energy-independent protein complex (with a peptidase domain of the clan MH, family M42, according to the MEROPS database) that was isolated from the neutrophilic haloarchaeon Har. marismortui58. It has been suggested that TET degrades oligopeptides released by ATP-dependent proteases such as the proteasome and LonB. Nab. magadii encodes a homolog of TET (Nmag_1335, peptidase M42 family protein), which, in combination with the energy-dependent proteases, may participate in the intracellular protein turnover in this extremophile. Furthermore, similar to the majority of haloarchaea, Nab. magadii appears to encode homologs of the three families of membrane-embedded regulatory proteases denoted as I-CLiPs. These include sppA-type signal peptide peptidases (SPPs, Nmag_2612 and 2635), site-2 protease class of zinc metalloproteases that cleave transmembrane domains (S2P type peptidases, Nmag_1508, 1514, 2136, and 3494), and rhomboids (Rho, Nmag_1128, 1636, 2518, and 3579). Furthermore, Nab. magadii contained genes encoding type I signal peptidases (sec11-type, Nmag_1326, 1932, 1944, 3375, 3743, and 4175) and a type IV prepilin peptidase (Nmag_1752). The type I signal peptidases and the type IV prepilin peptidase are predicted to be involved in the processing of N-terminal signal peptides of exported proteins and flagellin precursors, respectively.

Cellular protease activity is frequently controlled by endogenous protease inhibitors 59. Genes encoding putative homologs of protease inhibitors of the serpin (Nmag_2110) and phosphatidylethanolamine-binding protein (Nmag_0329) types were present in Nab. magadii. A subtilisin protease inhibitor from this archaeon, denoted NSI, was previously purified and biochemically characterized 60. This protease inhibitor remains to be investigated at the molecular level and the availability of its gene sequence could facilitate cloning and expression of the recombinant protein for further analysis. A representation of the major proteolytic systems (predicted and/or validated by detecting mRNA and/or assaying protein activity) of Nab. magadii is presented in Figure 2. Although this depiction assumes that the proteolytic systems of Nab. magadii are independent of each other, their synergistic action in vivo cannot be ruled out.

Protein translocation across the cell membrane in prokaryotes is facilitated by at least three mechanisms including the general secretion (Sec) system, the specialized Tat system, and the highly ornate, substrate-specific secretion systems for delivering effector proteins to target sites 616263. Nab. magadii contained genes that encoded putative components of the Sec system (secYEGDF, Nmag_0233, 1140, 1564, 2707, and 2708; srp19, Nmag_3604; srp54, Nmag_1802; ftsY, Nmag_0182) and of the Sec-independent Tat protein translocase complex (tatC, Nmag_2050-2051; tatA, Nmag_3135). While the Tat pathway is commonly used for a small subset of exported proteins in bacteria, it is a dominant export route in halophilic archaea. Many of the exported proteins are subsequently attached to the cell membrane by a lipid anchor and Nab. magadii has 119 genes encoding lipid-modified Tat target proteins, as detected by TatLipo analysis 64. Furthermore, Nab. magadii contained genes encoding putative components of a type II secretion system (Nmag_3137-3138) and an archaeosortase (Nmag_2750) for which 17 targets with PGF_CTERM motif were identified 65.

N-glycosylation in archaea and eukaryotes uses dolichol phosphate as the lipid base for the assembly of oligosaccharides 6667. Glycosyltransferases (GTs) are key components of N-glycosylation in all three domains of life, and the genome of Nab. magadii contained 23 genes encoding putative GTs. Based on BLASTP analysis on the NCBI database and the presence of conserved domains, these genes were assigned into the GTA (Nmag_0916, 0926, 1132, 1200, 2046, 2620, 2830, 3015, 3273, 3275, 3807, and 4148) and GTB (Nmag_0132, 0135, 0432, 0925, 2541, 3017, 3285, 3431, 3512, 3832, and 3843) superfamilies. One of these genes (Nmag_3015) is in an operon with Nmag_3011 (hexapeptide repeat-containing transferase, 192 aa), Nmag_3012 (aminotransferase, 481 aa), Nmag_3013 (oxidoreductase, 340 aa), and Nmag_3014 (nucleotide sugar dehydrogenase, 544 aa). Nab. magadii also contained genes encoding a putative oligosaccharyltransferase subunit (Nmag_0927, aglB homolog) and a dolichol kinase-like protein (Nmag_1986). Therefore, Nab. magadii appears to have the genetic potential for N-glycosylation.

Several species of halophilic archaea are known to produce copious amounts of extracellular polysaccharides 68. Although transmission electron microscopic (TEM) images show the presence of an exopolysaccharide-like material around Nab. magadii cells (Figure 3A), purification and biochemical analyses of this material are yet to be accomplished. Nab. magadii contained six genes encoding putative polysaccharide biosynthesis proteins (Nmag_0147, 0922, 2457, 3122, 3272, and 3437). Other genes in the genome that encoded putative enzymes involved in polysaccharide biosynthesis included six polysaccharide deacetylases (Nmag_1899, 2045, 2647, 3024, 3271, and 3278), two polyprenyl glycosyl-phosphotransferases (Nmag_0111 and 1184, 65% identity at the predicted protein level), an O-antigen polymerase (Nmag_0143), two UDP-N-acetylglucosamine 2-epimerases (Nmag_0149 and Nmag_0676), an acylneuraminate cytidylyltransferase (Nmag_0148), an O-acetyltransferase (Nmag_0150), a N-acylneuraminate-9-phosphate synthase (Nmag_0151), and two capsule synthesis proteins (Nmag_1511 and 3999). It is possible that some of these genes are involved in the biosynthesis of Nab. magadii exopolysaccharide- or capsule-like material identified in the TEM images.

Archaeal metabolic pathways are unique and diverse, in comparison to those of eubacteria 24. However, the biosynthesis of several coenzymes involved in archaeal metabolism has only been partially understood. Recent advances in this area include the discovery of a new heme biosynthetic pathway 69, further understanding of the pathway of cobalamin biosynthesis, and the reconstruction of a complete pathway for coenzyme F₄₂₀ biosynthesis in haloarchaea. Analyses of genes putatively involved in coenzyme biosynthesis in Nab. magadii was performed in light of these new discoveries. This section describes the genes related to the biosynthesis of thiamine, riboflavin, NAD, coenzyme F₄₂₀, folate, heme, and cobalamin.

Vitamin B₁ (thiamine pyrophosphate) is involved in several microbial metabolic functions 70. Prokaryotes have evolved elaborate mechanisms to either synthesize this important co-factor de novo or acquire it from their niche 71. Thiamine biosynthetic pathways among prokaryotes are very diverse 70. Thiamine biosynthesis is accomplished by joining two intermediate molecules that are synthesized separately. One of these molecules is hydroxymethylpyrimidine pyrophospate (HMP-PP), which is made from aminoimidazole ribotide (AIR, an intermediate of purine biosynthesis) using ThiC and ThiD. The other molecule is hydroxyethylthiazole phosphate (HET-P), which in bacteria is generated by ThiGH and TenI and involves the sulfur carrier ThiS. The sulfur carrier is activated for thiolylation via C-terminal adenylation catalyzed by the N-terminal E1-like domain of ThiI. Nab. magadii contained a gene (Nmag_3460) encoding a putative ThiI and several ubiquitin-like β-grasp fold proteins (Nmag_0567, 1914, 2668, and 2971). However, β-grasp fold proteins have multiple functions in halophilic archaea, being involved in sulfur chemistry as well as ubiquitin-like protein modification by SAMPylation 56. The specific β-grasp fold protein likely to participate in thiamine biosynthesis in Nab. magadii remains unidentified.

Furthermore, Nab. magadii lacked homologs of ThiG, ThiH, and TenI involved in HET-P biosynthesis in bacteria. Interestingly, Nab. magadii encoded a homolog of the yeast HET-P synthase THI4 (Nmag_2419). However, Nmag_2419 is currently annotated as ribose-1,5-bisphosphate isomerase based on the functional characterization of the ortholog MJ0601 from Methanocaldococcus jannaschii. In contrast, the ortholog of Nmag_2419 in Pyrococcus kodokaraensis (TK0434) does not have ribose-1,5-bisphosphate isomerase activity 72. Biochemical characterization is required to ascertain the potential role of Nmag_2419 in thiamine biosynthesis. The genome of Nab. magadii contained purM (Nmag_1281) and thiC (Nmag_2593) homologs, which were predicted to be involved in AIR and HMP biosynthesis, respectively (Table 2).

Two distinct proteins, ThiE and ThiN, are known to join HMP-PP and HET-P to generate thiamine phosphate. Nab. magadii and other halophilic archaea contain both enzymes. Whereas ThiE (Nmag_1811) is a monofunctional protein, ThiN (Nmag_1282) exists as a C-terminal domain in a ThiDN fusion protein. At the last step, thiamine phosphate is predicted to be further phosphorylated to thiamine pyrophosphate by ThiL (Nmag_1515). Therefore, the conversion of AIR to HMP-PP in Nab. magadii appears to be similar to the bacterial pathway and may involve ThiC and ThiD, whereas HET-P biosynthesis in this haloarchaeon appears to be similar to the eukaryotic pathway and may involve Nmag_2419. In addition, Nab. magadii contained genes encoding a HET kinase (ThiM, Nmag_1810, predicted to be involved in thiamine salvage) and a thiamine transporter (ThiBPQ, Nmag_460-462; ThiB2, Nmag_1940).

Vitamin B₂ (riboflavin) is the precursor of coenzymes flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD), which are cofactors for several biochemical reactions 73. Most bacteria, fungi, and plants can synthesize riboflavin de novo using one molecule of GTP and two molecules of ribulose 5-phosphate as substrates 74. Riboflavin biosynthesis has been described in M jannaschii75. In general, reduction precedes deamination in the archaeal riboflavin biosynthesis pathway, which appears to be similar to the fungal pathway. In M. jannaschii, GTP cyclohydrolase III, the first enzyme of the riboflavin biosynthesis pathway, produces an archaeal-specific formylated intermediate that requires a subsequent deformylation step 75. However, the haloarchaeal homolog of GTP cyclohydrolase III has not been identified thus far. Furthermore, riboflavin kinases of halophilic archaea are homologous to those of bacteria, but are unrelated to M. jannaschii riboflavin kinase. Conversely, riboflavin synthases of halophilic archaea are related to those of M. jannaschii, but are unrelated to bacterial riboflavin synthases. Overall, six genes encoding putative enzymes of the riboflavin biosynthesis pathway were identified in Nab. magadii and only two of these (Nmag_0941 and 0942) were clustered together (Table 2).

Vitamin B₃ (nicotinic acid) is the central component of coenzymes NAD and NADP, which are essential redox cofactors in metabolism. Like most bacteria, halophilic archaea synthesize NAD from aspartate via quinolinate. Nab. magadii contained 7 genes (Nmag_2920-2922 form an operon and encode nadABC, Nmag_2823 encodes nadM, and Nmag_1544, 2475, and 2846 encode nadE) that were predicted to be involved in NAD biosynthesis.

Coenzyme F₄₂₀ is involved in methanogenesis 7677 and other metabolic pathways (e.g., aflatoxin reduction in mycobacteria 7879) that require hydride transfer from the low-potential reduced deazaflavin F₄₂₀ to substrates with electron-deficient ring systems. Furthermore, 5-amino-6-(D-ribitylamino)uracil is an intermediate of coenzyme F₄₂₀ and riboflavin biosynthesis pathways 7677. Although coenyzme F₄₂₀ has been detected in some halophilic archaea 80, and coenyzme F₄₂₀ biosynthesis genes have been identified among the methanogens 81, its precise function in halophilic archaea is unknown. Nab. magadii and other halophilic archaea contained several genes encoding putative enzymes of the coenyzme F₄₂₀ biosynthesis pathway. These genes were identified based on the presence of their homologs among the methanogens. Furthermore, using SIMBAL analysis 82, coenyzme F₄₂₀ dependent enzymes have been predicted in halophilic archaea (e.g., NP1902A, TIGRFam 04024, D. Haft, personal communication). Nab. magadii predicted proteins with an assigned F₄₂₀-related TIGRFam are shown in Table 2. A possible function of coenzyme F₄₂₀ in the respiratory chain of Nab. magadii is discussed later in this paper.

Tetrahydrofolate participates in a number of biochemical reactions and reduced folate cofactors are required for the biosynthesis of a variety of molecules in both prokaryotes and eukaryotes 8384. The production of folate involves several enzymes catalyzing the pterin and para-aminobenzoic acid branches of the pathway 85. Nab. magadii genes putatively involved in folate biosynthesis were generally similar to those described in Nmn. pharaonis, including the archaeal-type GTP cyclohydrolase (Nmag_2853). Nab. magadii contained a pabABC operon (Nmag_2792-2794) and the homologs of these genes were predicted to be involved in para-aminobenzoate biosynthesis in Nmn, pharaonis24. Nab. magadii also contained genes encoding a dihydropteroate synthase, a fused dihydropteroate-dihydrofolate synthase, a dihydrofolate reductase, and a methenyltetrahydrofolate cyclohydrolase (Table 2). The latter enzyme is predicted to participate in the conversion of the C1 metabolite attached to tetrahydrofolate. However, none of the other genes encoding C1-converting enzymes identified in Haloquadratum or Haloarcula were found in Nab. magadii and Nmn. pharaonis. Furthermore, Nab. magadii contained three genes encoding putative enzymes of the later stages of the folate biosynthesis pathway (FolP, Nmag_0002; FolC-Prd-FolP, Nmag_2554; FolA, Nmag_2988). Nevertheless, genes encoding the bacterial homologs of FolQBK, the atypical FolQ described in Desulfovibrio, or the alternative pathway bypassing FolQB (described in Plasmodium and several bacteria) were absent in Nab. magadii and Nmn. pharaonis. Therefore, precise mechanisms of folate biosynthesis in these species remain to be discovered.

Environmental bacteria utilize a variety of redox molecules such as porphyrins and other modified tetrapyrroles like heme, siroheme, and adenosylcobalamin for catalysis, energy transfer, and signal transduction 86. These tetrapyrroles are synthesized de novo using a branched pathway and aminolevulinic acid as the precursor 8788. In most prokaryotes, the conversion of glutamate to aminolevulinic acid is catalyzed by glutamyl-tRNA synthetase, glutamyl-tRNA reductase, and glutamate-1-semialdehyde aminotransferase. Two molecules of aminolevulinic acid are condensed by the action porphobilinogen synthase to form porphobilinogen. Four molecules of porphobilinogen are polymerized by the action of the porphobilinogen deaminase to form the tetrapyrrole hydroxymethylbilane. Uroporphyrinogen III methyltransferase cyclizes hydroxymethylbilane to produce uroporphyrinogen III. Uroporphyrinogen III is converted to precorrin-2 in the biosynthetic pathway of adenosylcobalamin and siroheme, which was recently found to be an intermediate of heme biosynthesis (see below).

The complete pathway for the biosynthesis of adenosylcobalamin from precorrin-2 involves two major branches and several enzymes 8990, some of which are archaea-specific (e.g.cobYcbiZ) 9192. Halophilic archaea (e.g.Nmn. pharaonis) use the “anaerobic” branch, which is characterized by an oxygen-independent ring contraction process (cibG) 93. However, it has been shown that Halobacterium synthesizes cobalamin de novo under aerobic conditions 91. The “anaerobic” branch is also characterized by early cobalt insertion and Nmn. pharaonis has homologs of the ATP-independent early cobalt chelatase (CbiX) from Bacillus halodurans and Archaeoglobus fulgidus9495. In the “anaerobic” branch, seven archaeal enzymes are known to be involved in the conversion of precorrin-2 into cobyrinic acid (cbiLcbiHcbiFcbiGcbiEcbiT, and cbiC), but two pathway gaps (corresponding to cbiD/cobF and cbiJ/cobK) still remain. A set of 11 genes is known to be involved in conversion of cobyrinic acid into adenosylcobalamin (cbiAcobA, pduO, cbiPcbiZcbiBcobYcobScobCcobD, and cobT).

Based on genome analyses, it appeared that Nab. magadii was incapable of de novo cobalamin biosynthesis since it lacked the genes encoding enzymes for conversion of precorrin-2 into cobyrinic acid. This is in contrast to Htg. turkmenica, which was predicted to be capable of de novo cobalamin biosynthesis since it contained the corresponding genes. However, Nab. magadii was predicted to be capable of corrinoid salvage since it contained a gene encoding a putative corrinoid transporter. Nab. magadii also contained a set of genes that were predicted to be involved in the conversion of cobyrinic acid into adenosylcobalamin, including a gene (cbiZ) that is specific to the archaeal corrinoid salvage pathway (Table 2).

The heme biosynthesis pathway in archaea involving uroporphyrinogen III, precorrin-2, and siroheme appears to be similar to that of Desulfovibrio69. Conversion of uroporphyrinogen III into siroheme requires three functions (methylation by SirA, iron chelation, and oxidation by SirC). The enzyme catalyzing iron chelation is unknown since the haloarchaeal precorrin-2 dehydrogenase might be monofunctional (as in Bacillus megaterium SirC) or might also be a ferrochelatase (as in trifunctional Escherichia coli CysG or the bifunctional yeast MET8). From comparison of Nab. magadii with other halophilic archaea, another possibility emerges: iron chelation may be performed by one of the proteins annotated as CbiX-type cobalt chelatase. Nmn. pharaonis has three cbiX paralogs (cbiX1, NP1108A; cbiX2, NP1588A; cbiX3, NP0734A), two of which have closely related orthologs in nearly all other halophilic archaea. Htg. turkmenika, probably capable of de novo cobalamin biosynthesis, has orthologs of both cbiX1 and cbiX2. However, Nab. magadii, which is predicted to be incapable of de novo cobalamin biosynthesis, and therefore expected to lack these proposed early cobalt chelatases, surprisingly contained a cbiX2 ortholog (Nmag_3212). It is possible that cbiX2 functions as a ferrochelatase during siroheme/heme biosynthesis rather than as a cobaltochelatase during de novo cobalamin biosynthesis.

Didecarboxysiroheme, a common intermediate of heme and heme d1 biosynthesis, is generated by the decarboxylation of siroheme on the C12 and C18 acetyl groups 69. Siroheme decarboxylase activity is attributed to the nirDLGH gene set, which is represented by a pair of two-domain proteins in halophilic archaea (AhbA/NirDL, Nmag_2894; AhbB/NirGH, Nmag_1221). Heme d1 is a coenzyme of dissimilatory nitrite reductase (also called cytochrome cd₁) and is not required by organisms lacking this enzyme. The last steps of heme biosynthesis include the removal of acetyl side chains of Fe-coproporphyrin by AhbC (encoded by Mbar_A1793 in Methanosarcina barkeri) and the oxidative decarboxylation of heme by AhbD (encoded by DVU_0855 in Desulfovibrio) 69 [gene assignments: M. Warren, personal communication, including the information that Supplemental Figure 3 shows the sequence of ahbC and not ahbD. Orthologs encoding AhbC and AhbD were present in Nmn. pharaonis (NP1542A and NP1546A) and Htg. turkmenica (Htur_1726 and 1728), but not in Nab. magadii. The presence of ahbC and ahbD in some halophilic archaea but not in others is believed to be due to metabolic heterogeneity rather than incomplete heme biosynthesis. Conversion of heme (also called heme B) into heme A in Nab. magadii was predicted to be catalyzed by CtaA and CtaB homologs (ctaA, Nmag_0636; ctaB, Nmag_2302).

Vitamin H, commonly known as biotin, acts as a coenzyme in several enzyme-catalyzed carboxylation and decarboxylation reactions 96. Most bacteria can synthesize biotin de novo using pimelic acid as a precursor, and some others have evolved mechanisms for importing this essential cofactor from their natural environments 9798. Whereas Nab. magadii is a biotin auxotroph, Nmn. pharaonis is a biotin prototroph and the genome of this haloalkaliphilic archaeon has been shown to contain at least three genes putatively involved in the biosynthesis of biotin 24. The absence of genes for the biosynthesis of biotin in Nab. magadii was apparent from the analyses of its genome sequence. However, the large chromosome of Nab. magadii contained a locus encoding a putative biotin transporter (BioYMN, Nmag_0886-0888), which may facilitate the uptake of biotin from the environment.

Nutritional requirements of halophilic archaea in the laboratory are as diverse as their observed phenotypes, suggesting that the metabolic pathways in these organisms are quite intricate 2499. The analysis of the genome sequence provided an unprecedented opportunity to comprehend the metabolic versatility of Nab. magadii. Additional file 4: Table S4 contains a comprehensive list of genes predicted to be involved in a diverse array of functions. Furthermore, genes encoding putative enzymes for archaeal modified pathways of gluconeogenesis and glycolysis as well as those of ribose metabolism and the tricarboxylic acid cycle were present in Nab. magadii. Genes that encoded putative enzymes for glycerol utilization, aromatic amino acid catabolism, ureagenesis, and urea degradation were also identified in Nab. magadii. Other putative metabolic features of Nab. magadii included xylose isomerases, an alpha amylase, a methylglyoxal synthase, sulfatases, a chlorite dismutase, sarcosine oxidases, and aldehyde dehydrogenases (Additional file 4: Table S4).

Molybdenum cofactor (MoCo) is essential for the functioning of molybdoenzymes such as dimethylsulfoxide and trimethylamine-N-oxide reductases, formate dehydrogenases, and nitrate reductases 24100. Molybdopterin is the dithiolene-containing tricyclic moiety found within MoCo of all molybdoenzymes except nitrogenases 101. In bacteria, genes of the moa, mob, mod, moe, and mog loci have been implicated in the biosynthesis of MoCo using GTP as the substrate 102. The large chromosome of Nab. magadii contained 9 genes encoding MoCo biosynthesis functions (Additional file 4: Table S4). It is uncertain if this subset of genes is sufficient for MoCo biosynthesis in Nab. magadii and biochemical studies are required to test whether this haloarchaeon is molybdenum-dependent.

The haloarchaeon Haloarcula marismortui converts acetyl-CoA to glyoxylate via the key intermediate methylaspartate 103. Glyoxylate is condensed with a second molecule of acetyl-CoA to form malate, which is an intermediate of the tricarboxylic acid cycle. Malate can subsequently be converted to oxaloacetate, which is used by phosphoenolpyruvate carboxykinase for gluconeogenesis. In Nab. magadii, activities of the enzymes of the methylaspartate cycle, but not those of the key enzymes of the glyoxylate cycle, were detected 103. An operon (Nmag_3333-3338) encoding putative homologs of the methylaspartate cycle and a gene encoding a putative phosphoenolpyruvate carboxykinase (Nmag_3507) were present in Nab. magadii. The square archaeon Haloquadratum walsbyi contains a gene (HQ2709A) encoding a phosphoenolpyruvate-dependent phosphotransferase system (PTS) that is predicted to be involved in the phosphorylation of dihydroxyacetone 104. Homologs of HQ2709A and genes encoding additional PTS components were present in Htg. turkmenicaHfx. volcanii, and several other haloarchaeal genomes. However, Nab. magadii and Nmn. pharaonis lacked homologs of these genes encoding PTS components.

Running a proton-driven, energy-conserving ATP synthase at high extracellular pH is an obvious challenge. Energy coupling of sodium ions instead of protons was proposed to be an adaptation to alkaliphilic growth conditions and an ATP synthase driven by Na⁺ is the hallmark of such an adaptation. Nab. magadii had a cluster of eight genes that form the atpHIKECFAB operon encoding putative ATP synthase subunits (Nmag_1370-1377) and an unlinked atpD homolog (Nmag_1366). Similar gene clusters were found in several halophilic archaea. Ion specificity of the ATP synthase is determined by the c-ring, which is encoded by the atpK gene (Nmag_1375) for A-type ATP synthases. Nab. magadii may have a proton-driven ATP synthase since its predicted AtpK lacks the sequence signature of Na⁺-dependent ATP synthases 105. Instead, within the ion-determining region of AtpK, the sequence (PETLVIL) is identical to that of the proton-driven ATP synthases from Hfx. volcanii106, Hbt. salinarum107, and Nmn. pharaonis23.

Reduction of oxygen and the associated proton-coupled electron transfer (respiration) is the primary source of energy among aerobic organisms. Respiratory complexes, which include a variety of cytochromes and terminal oxidases, are essential components of this process. Biochemical and comparative genomic analyses of the electron transport chain of Nmn. pharaonis have revealed several novel features, including a gene encoding a type II NADH dehydrogenase (NP3508A) 23108. A homolog of NP3508A in Acidianus ambivalens was proposed to be involved in NADH reoxidation, feeding into the lipid-soluble quinone pool 109. A homolog of NP3508A was also present in Nab. magadii (Nmag_0301) and several other halophilic archaea. Nab. magadii also contained genes encoding a putative nuo complex (Nmag_3245-3255), which was similar to the mitochrondrial NADH dehydrogenase (complex I). Although 13 nuo cluster subunits were conserved among halophilic archaea and E. coli, the nuoEFG subcomplex, which is involved in accepting NADH, was missing in halophilic archaea 23. Furthermore, involvement of a type I complex in NADH reoxidation has been ruled out in Hbt. salinarum110. It is speculated that reduced coenzyme F₄₂₀, which is similar to NADH in its redox potential, may interact with the nuo complex in halophilic archaea. In addition to the NADH dehydrogenases, Nab. magadii and other halophilic archaea are predicted to encode a succinate dehydrogenase that may oxidize succinate and reduce the quinone pool of the electron transport chain.

A number of cytochromes involved in respiratory electron transport have been characterized among the archaea 12111. Terminal oxidases, also known as oxygen reductases, can accept electrons from a variety of donors and reduce dioxygen to water. The large chromosome of Nab. magadii contained loci encoding putative cytochrome c-type terminal oxidase subunits I and II (Nmag_0263-0264) and cytochrome ubiquinol oxidase subunits I and II (Nmag_1036-1037). Furthermore, pNMAG02 contained an operon encoding putative cytochrome ubiquinol oxidase subunits I and II (Nmag_4038-4039) that were related to the proteins encoded by Nmag_1036-1037 (46% identity at the protein level). The homologs of Nmag_0263-0264 and Nmag_1036-1037 were present in Htg. turkmenica (Htur_2248- 2249 and Htur_4570- 4571, respectively), but not in Nmn. pharaonis. Two cbaDBAC operons (plus a monocistronic cbaE) encoding putative cytochrome ba₃ terminal oxidase complexes were identified in pNMAG01 (Nmag_3754-3758 and 3802–3805), and these operons appeared to be related to each other (38-65% identity at the predicted protein level). The homologs of these ORFs were present in Htg. turkmenica (Htur_0462- 0466) and Nmn. pharaonis (NP2960A-NP2968A), the latter of which have been functionally characterized 112. Halocyanins are predicted to act as one-electron carriers to the terminal oxidases in halophilic archaea 23108. This prediction is supported by the observation that cbaD is fused to halocyanin in Hbt. salinarum and Har. marismortui113. Nab. magadii contained several genes encoding putative halocyanin-like proteins (Nmag_2576, 2424, 0446, 0741, 3725, 1878, 3525, 1774, and 3800). Halocyanins are coupled to the reduced quinone pool by the cytochrome bc₁ complex. Although genes encoding a cytochrome bc₁ complex are present in Hbt. salinarumHfx. volcanii, and several other halophilic archaea, they were absent in Htg. turkmenicaNmn. pharaonis, and Nab. magadii. The electron transfer from the reduced lipid-soluble quinone pool to halocyanin remains unresolved in species that lack a cytochrome bc₁ complex. Based on genome comparisons, it appeared that the respiratory chain of Nab. magadii was more similar to that of Htg. turkmenica than to Nmn. pharaonis.

Other genes encoding putative cytochrome-related proteins in the large chromosome of Nab. magadii included Nmag_0636 (cytochrome oxidase assembly protein), Nmag_1972 (cytochrome P450), and Nmag_3057 (cytochrome c biogenesis protein). The small chromosome pNMAG01 contained a gene encoding a putative cytochrome c biogenesis protein (Nmag_3708) that had 64% identity to the protein encoded by Nmag_3057. Nab. magadii also contained genes (Nmag_2430-2432) encoding a putative sulfur utilization factor (SUF) system, which was shown to be important for Fe-S cluster biogenesis during stress in E. coli114. Other genes predicted to participate in bioenergetic conversion in Nab. magadii include those encoding electron transfer flavoprotein subunits (Nmag_0482-0483 and 1388–1389), SCO1/SenC electron transport proteins (Nmag_0793, 3059, and 3710), and a redoxin domain protein (Nmag_3709).

Two-component signal transduction systems consisting of a histidine kinase (HK) and a response regulator (RR) constitute one of the most frequently encountered bacterial and archaeal communication circuits 115. Nab. magadii contained a pair of genes (Nmag_3535-3536) encoding putative HK-RR proteins. Nab. magadii also contained two pairs of genes (Nmag_2042-2043 and Nmag_3740-3741) encoding putative HK-RR proteins with an additional RR domain in the N-terminus of the predicted HK protein. Furthermore, another set of genes (Nmag_1129-1130) encoding a HK-like protein, which was distantly related to CheA of Nmn. pharaonis, and a putative protein containing a RR domain was also predicted in Nab. magadii. Interestingly, Nab. magadii also contained 15 genes encoding putative HK proteins without a cognate RR gene (Nmag_0168, 0323, 0435, 1296, 1909, 2062, 3106, 3188, 3216, 3297, 3954, 4064, 4101, 4114, and 4119) and 11 genes encoding putative RR proteins without a cognate HK gene (Nmag_0951, 1095, 1797, 2147, 2692, 2693, 2877, 3151, 3379, 3389, and 3679).

Halobacterial perception of and response to physical stimulus such as light (phototaxis) is mediated by photoreceptors 116. The genome of Nab. magadii contained a single rhodopsin gene (Nmag_2582, the chloride pump halorhodopsin). Nab. magadii also contained 3 genes (Nmag_1701, 2879, and 2881) encoding distant homologs of rhodopsins that were related to each other (36-49% sequence identity). They were also related to NP1758A from Nmn. pharaonis (33%-42% sequence identity) and predicted to encode distant rhodopsin homologs that lack the Lys residue involved in covalent retinal attachment. Whether these rhodopsin homologs interact with retinal noncovalently, or if they interact with retinal at all, is unknown. Two of the three genes encoding retinal homologs in Nab. magadii (Nmag_2879 and 2881) were located adjacent to genes (Nmag_2880 and 2882) encoding putative methyl-accepting chemotactic transducer proteins. In Nmn. pharaonis, NP1758A and NP3132A, which are homologs of Nmag_2879 and 2881, were also found adjacent to genes (NP1756A and NP3134A, respectively) encoding putative methyl-accepting chemotactic transducer proteins 23. Therefore, this group of distant rhodopsin homologs may be involved in perception of external stimuli, although it remains to be determined if they are involved in light perception. Although Htg. turkmenica lacked genes encoding bacteriorhodopsin and halorhodopsin, it contained a single gene (Htur_3663) that appeared to be a distant homolog of Nmag_2879 and 2881. Furthermore, similar to Nmn. pharaonisNab. magadii lacked a gene encoding a proton pump bacteriorhodopsin. However, Nab. magadii contained a locus encoding putative phytoene desaturase, UbiA prenyltransferase, carotene biosynthesis protein, and phytoene synthase (Nmag_1001-1004) as well as a gene encoding a putative squalene/phytoene synthase (Nmag_2309, unrelated to Nmag_1004).

Two unrelated enzyme families are used for the cleavage of β-carotene into retinal in halophilic archaea. Distant paralogs, which belong to one of these β-carotenase enzyme families, have been identified in Hbt. salinarum and designated Brp and Blh 117. Although homologs of brp and blh were present in Nmn. pharaonis (NP0650A and NP0206A, respectively), they were absent in Nab. magadii. However, Nab. magadii contained a homolog (Nmag_4083) of Hqr. walsbyi HQ2020A, which was predicted to encode a distinct β-carotenase unrelated to those mentioned above 104. Interestingly, Htg. turkmenica lacked homologs of brp and blh as well as HQ2020A, which is consistent with the absence of all canonical rhodopsins in this organism.

Microbial response to chemical stimulus (chemotaxis, movement toward nutrients or away from stressors) is mediated by chemoreceptors 118. The large chromosome of Nab. magadii contained two loci encoding putative motility and signal transduction functions. One of them contained exclusively “che genes” in a cheYBACCDR operon (Nmag_3145-3151), which is preceded by two divergently transcribed and distantly related cheW genes (Nmag_3152-3153). The cheYBACCDR operon encodes a very long signal transduction histidine kinase (CheA, 1576 amino acids), a response regulator receiver protein (CheY), a CheR-type MCP methyltransferase, and a response regulator receiver-modulated methylesterase (CheB). While CheD has been reported to function as glutamine deamidase in some organisms or as methylesterase in others, CheB functions as both 119. A similar locus was also present in Htg. turkmenica (Htur_0954-0962) and based on predicted protein homology, it appeared that the two loci were evolutionarily very closely linked. Highly similar gene clusters were also found in the genomes of Halopiger xanaduensisNatrinema pellirubrum, and Natronobacterium gregoryi. The second locus (Nmag_2859-2889) contained “che genes” along with “fla genes” encoding flagellin biosynthesis and assembly functions (see below). The “che genes” in this locus encode putative CheA, CheY, CheR, CheB, CheD, a two-domain CheC, and two CheW proteins. Nab. magadii contained two cheF genes within this locus and homologs of these genes were shown to be involved in chemotaxis in Halobacterium120. This gene cluster also encoded three methyl-accepting chemotaxis sensory transducers (Nmag_2880, 2882, 2884), two of which were adjacent to genes encoding distant rhodopsin homologs (Nmag_2879 and 2881). Other genes encoding putative methyl-accepting chemotaxis sensory transducers in Nab. magadii include Nmag_0478, 0937, 1253, 1386, 1542, 2639, 3325, 3638, and 3856). Among these, two (Nmag_1253, 1386) were adjacent to genes encoding periplasmic ligand-binding proteins.

Archaeal flagella are very different in composition and assembly in comparison to bacterial flagella 121. In contrast to the bacterial flagellar motor, which is driven by an ion gradient, the archaeal flagellar motor is driven by ATP, as shown in Hbt. salinarum122. Within the second motility and signal transduction gene cluster (Nmag_2859-2889) of Nab. magadii is a region (14,766 bp, 58.66% GC) with 13 predicted ORFs encoding putative flagellin biosynthesis and assembly proteins (Nmag_2862-2874). Except Nmag_2871, which encoded a protein of unknown function, all other ORFs were located on the plus strand. This region contains 4 flagellin genes (Nmag_2862-2865), which encode the flagella proteins previously identified 123 and characterized 17. Furthermore, Nab. magadii contained homologs of flaF (Nmag_2866), flaG (Nmag_2867), flaH (Nmag_2873), flaI (Nmag_2874), and flaJ (Nmag_2869). The latter two genes (flaI and flaJ) encode putative proteins homologous to the type II secretion system proteins E and F, respectively. In several archaea, FlaI has been shown to be involved in flagellin assembly 124125, and was recently proposed as a motor component 126. The motility gene clusters of halophilic archaea are generally polymorphic, probably due to divergence of genome organization and deletion/duplication of the accessory genes 23. Nevertheless, flaHflaI, and flaJ represent a core set of highly conserved genes presumably crucial for archaeal motility.

Since previous electron microscopic analyses have demonstrated that Nab. magadii contains distinctive flagella 127128129, and structures resembling flagella are also visible in the TEM images in Figure 3, it is likely that the fla locus (Nmag_2862-2874) of Nab. magadii was involved in flagellin biosynthesis and motility. In addition, the large chromosome of Nab. magadii contained genes encoding a putative full-length PilT protein (Nmag_1543) and a prepilin peptidase (Nmag_1752), whose homologs were found in Nmn. pharaonis (NP0198A and NP1276A) and Htg. turkmenica (Htur_3514 and Htur_0098).

Archaeal basal transcription machinery has many similarities to the eukaryotic RNA polymerase II apparatus. However, the mechanisms of transcription regulation and the transcriptional regulators among archaea are distinct from those of eukaryotes 130131. Nab. magadii contained 90 genes (~2.2% of all protein-coding genes) encoding putative transcriptional regulators (TRs). BLASTP analyses indicated that most of these predicted proteins were related to bacterial TRs. These TRs were categorized into the following families based on their helix-turn-helix (HTH) motifs (numbers in parenthesis indicate the number of proteins in each family): AsnC (16), IclR (13), PadR (10), ArsR/GntR/Fur (11), XRE (6), TrmB (8), AbrB (3), TetR (3), MarR (3), CopG/Arc/MetJ/NikR (10), HxlR (2), TenA (1), ModE (1), and unassigned (3). Apart from these TRs, Nab. magadii also contained 27 genes encoding TRs with an HTH-10 domain, which was also found in bacterio-opsin activators . In addition, Nab. magadii encoded a two-domain archaeal histone (Nmag_2205), a single TATA-binding transcription initiation factor (Nmag_1863, 93% protein sequence identity to TbpE from Hbt. salinarum), a single transcription initiation factor TFE (Nmag_3157), and a set of seven transcription initiation factors TFB (Nmag_0308, 0527, 2197, 2805, 3037, 3548, and 4179).