- Research article
- Open Access
The core and unique proteins of haloarchaea
BMC Genomics volume 13, Article number: 39 (2012)
Since the first genome of a halophilic archaeon was sequenced in 2000, biologists have been advancing the understanding of genomic characteristics that allow for survival in the harsh natural environments of these organisms. An increase in protein acidity and GC-bias in the genome have been implicated as factors in tolerance to extreme salinity, desiccation, and high solar radiation. However, few previous attempts have been made to identify novel genes that would permit survival in such extreme conditions.
With the recent release of several new complete haloarchaeal genome sequences, we have conducted a comprehensive comparative genomic analysis focusing on the identification of unique haloarchaeal conserved proteins that likely play key roles in environmental adaptation. Using bioinformatic methods, we have clustered 31,312 predicted proteins from nine haloarchaeal genomes into 4,455 haloarchaeal orthologous groups (HOGs). We assigned likely functions by association with established COG and KOG databases in NCBI. After identifying homologs in four additional haloarchaeal genomes, we determined that there were 784 core haloarchaeal protein clusters (cHOGs), of which 83 clusters were found primarily in haloarchaea. Further analysis found that 55 clusters were truly unique (tucHOGs) to haloarchaea and qualify as signature proteins while 28 were nearly unique (nucHOGs), the vast majority of which were coded for on the haloarchaeal chromosomes. Of the signature proteins, only one example with any predicted function, Ral, involved in desiccation/radiation tolerance in Halobacterium sp. NRC-1, was identified. Among the core clusters, 33% was predicted to function in metabolism, 25% in information transfer and storage, 10% in cell processes and signaling, and 22% belong to poorly characterized or general function groups.
Our studies have established conserved groups of nearly 800 protein clusters present in all haloarchaea, with a subset of 55 which are predicted to be accessory proteins that may be critical or essential for success in an extreme environment. These studies support core and signature genes and proteins as valuable concepts for understanding phylogenetic and phenotypic characteristics of coherent groups of organisms.
Extremely halophilic Archaea (haloarchaea) have adapted to thrive in environments of high salinity, desiccation, and intense solar radiation. These microorganisms require at least 1.5 - 2.5 M NaCl for viability and typically display optimal growth in NaCl concentrations at or above 3.5 M. Haloarchaea commonly inhabit hypersaline environments, e.g. salt lakes, salterns, and heavily salted hides, meats, fish, and sauces [1–3]. Additionally, haloarchaea have been shown to survive space conditions  and viable cells have been reported from ancient deep underground salt deposits [5, 6]. Unlike most other extremophilic and archaeal microorganisms, haloarchaea form a monophyletic and coherent taxonomic group, the family Haloarchaeaceae .
The Halobacterium sp. NRC-1 genome sequence gave researchers the first opportunity, at the genome level, to probe the mechanisms of adaptation to hypersaline brine [8, 9]. Characterization of the 2 Mbp chromosome and two large megaplasmids showed that the overwhelming majority of predicted proteins were highly acidic, with a pI mode of 4.2, and very few neutral or basic proteins [10, 11]. In contrast, predicted proteins from most other non-haloarchaeal and bacterial organisms had equal fractions of acidic and basic components. The negatively charged residues in haloarchaeal proteins were predominantly found at the protein surface and predicted to function in enhancing their solubility and stability in high salt concentrations. A few individual haloarchaeal proteins have been crystallized, e.g. malate dehydrogenase, dihydrofolate reductase, and DNA sliding clamp (PCNA), and they all display markedly more acidic residues than non-haloarchaeal homologs. They also possess clusters of negative charges on the surface [12–14]. The high prevalence of negatively charged surface residues produces tightly bound hydration shells with salt ions bound at the protein surface [16, 17].
Several previous studies have examined the gene content in haloarchaea, including one aimed at identifying information transfer genes and another concerning metabolic genes [18, 19]. While a significant degree of conservation was found among the essential components of DNA replication, repair, and recombination, transcription, and translation, the study of metabolic genes showed substantially more diversity. Indeed, this diversity was illustrated by the recent identification of genes for a new pathway in central carbon metabolism, the methylaspartate cycle, in several haloarchaea . An additional characteristic observed in most haloarchaeal genomes is the presence of large megaplasmids or minichromosomes which often harbor important or essential genes . Gene content in these large extrachromosomal elements was compared and resulted in the finding of expanded gene families for replication and transcription initiation, e.g. orc and tfb , as well as the presence of a variety of genes needed for cell survival, e.g. an amino-acyl tRNA synthetase , resistance to arsenic , and production of buoyant gas vesicles .
In the current study, we present a comprehensive analysis of haloarchaeal genomes aimed at identifying the core haloarchaeal proteins and uniquely haloarchaeal groups. Halophilic Archaea representing thirteen different genera were included, all within the Haloarchaeaceae family. These microorganisms represent both geographic and phylogenetic diversity, including isolates from all 7 continents (Figure 1) and almost half of the genera in this tight clade of the Euryarchaea . The genome-wide analysis produced nearly 800 protein clusters that are completely conserved among sequenced haloarchaea and a subset of 55 protein families that are unique to this family of extremophilic microbes.
Haloarchaeal orthologous groups (HOGs)
Using the best reciprocal hit method, 31,312 predicted proteins from nine complete haloarchaeal genomes (Halobacterium sp. NRC-1, Haloarcula marismortui, Natronomonas pharaonis, Haloquadratum walsbyi, Halorubrum lacusprofundi, Halogeometricum borinquense, Halomicrobium mukohataei, Halorhabdus utahensis, and Haloferax volcanii) were initially compared to form 4,455 haloarchaeal orthologous groups (HOGs) (see Table 1 and Table 2; Figure 1 and 2) [23, 24]. Our results showed that the overwhelming majority of predicted haloarchaeal proteins were members of HOGs, ranging from a high of 82.8% for Halobacterium sp. NRC-1 to a low of 73.9% for H. utahensis. These results underscored the close relationship of these haloarchaeal species.
Core HOGs (cHOGs)
We examined the abundance of the haloarchaeal proteins present in these 4,455 HOGs and found a bimodal distribution (Figure 3). The largest number of protein clusters were found in either 2 or 3 haloarchaea (1358 or 716, respectively) or all 9 members (799 protein clusters), and the protein clusters with an intermediate (4 - 8) number of haloarchaea were less abundant (250 - 442). The 799 clusters conserved in all nine organisms were designated as core haloarchaeal orthologous groups (cHOGs) (see Additional file 1) and represented proteins that are known or expected to be important or essential in all of the haloarchaea (see below). Taking into account that several HOGs correspond to more than a single COG and KOG, comparison of the cHOGs to the COG and KOG databases in NCBI showed that of the 799 cHOGs, 422 corresponded to both COGs and KOGs and 288 corresponded to COGs only, with 89 novel clusters unique to haloarchaea.
Uniquely haloarchaeal orthologous groups (ucHOGs, tucHOGs, and nucHOGs)
Of the 799 cHOGs present in all nine haloarchaea, 89 (11%) appeared to be unique to haloarchaea based on their absence in both the COG and KOG databases. These unique core HOGs (ucHOGs) were candidates for being 'signature' proteins for this clade, based on their ubiquity among haloarchaea and absence in non-haloarchaeal clades (Figure 2). However, since the members of these protein clusters were quite diverse, with the percent identity varying widely (between 22% and 85%), we re-appraised the statistical significance of group members by carrying out pairwise alignments of the proteins within each cluster, including randomized global alignments for statistical analysis using the Needleman and Wunsch algorithm [25, 26]. Using this approach, we were able to establish a 99.9999% confidence level for pairs of sequences among proteins within each cluster.
With the rapid sequencing of new haloarchaeal genomes, we further scrutinized the 89 ucHOGs using a sequential multi-step approach: (1) protein sequences were BLASTed against four recently available complete haloarchaeal genome sequences (Haloterrigena turkmenica, Natrialba magadii, Halalkalicoccus jeotgali, and Halopiger xanaduensis) to find conserved haloarchaeal homologs, (2) protein sequences were BLASTed against the NCBI non-redundant database to find non-haloarchaeal hits, and (3) any non-haloarchaeal hits identified were aligned with each member of the cHOG cluster. Of the 89 clusters with no associated COGs or KOGs, all members of 55 ucHOG clusters were found to be truly unique core haloarchaeal orthologous groups and named tucHOGs (Figure 2). Of the remaining 34 clusters, 6 were absent in one or more of the four newly sequenced genomes, and 29 had one or more members with at least one hit to a non-haloarchaeal peptide. Proteins from six clusters had hits to over a dozen different non-haloarchaeal proteins and proteins from the remaining 23 clusters had fewer hits, ranging from 1 - 10 per cluster. The significance of hits was evaluated by base composition-preserved randomized alignments. This analysis showed that the 28 cHOG clusters with hits to non-haloarchaeal proteins were not entirely unique to the haloarchaea with a 99.0% or higher level of confidence, and were named nearly unique core haloarchaeal orthologous groups, or nucHOGs (Figure 2).
Genomic locations and functions of ucHOGs
Consistent with a critical role in the biology of haloarchaea, ucHOGs were found to be encoded overwhelmingly on the main chromosomes of haloarchaeal organisms. Indeed, in five, Halobacterium sp. NRC-1, H. marismortui, N. pharaonis, H. utahensis, and H. walsby i, all of the tucHOG polypeptides were chromosomally encoded and dispersed relatively evenly over the entire chromosome (Figure 4). Only seven tucHOG protein genes did not map to large chromosomal replicons, with two on the small chromosome in H. lacusprofundi, one each on the pHB200 and pHB500 megaplasmids in H. borinquense, two on the pHV4 megaplasmid in H. volcanii, and one on the pHM61 megaplasmid in H. mukohataei (see Additional file 1). Similarly, all of the nucHOG proteins mapped to the large chromosomes of N. pharaonis, H. walsbyi, H. borinquense, H. mukohataei, H. utahensis, and H. lacusprofundi. A single nucHOG protein is coded on both the smaller chromosome II and pNG600 in H. marismortui, two nucHOGs are coded on pHV4 in H. volcanii, and one nucHOGs is found on the common inverted repeats of pNRC100 and 200 of Halobacterium sp. NRC-1 (see Additional file 1).
The function of only a single uniquely conserved haloarchaeal orthologous protein gene, vng2163 (cluster tucHOG0456), has so far been investigated in any detail . In Halobacterium sp. NRC-1, the gene coding for this protein was annotated as ral (rfa-linked) due to its transcriptional linkage to two genes, rfa3 and rfa8, which encode eukaryotic replication protein A (RPA)-like single-stranded DNA binding protein subunits . The genes around ral showed a significant degree of synteny among the haloarchaeal genomes (Figure 5), consistent with a conserved function in haloarchaea.
Functional classification of HOGs
Biological functional categories were assigned to HOGs by membership of Halobacterium sp. NRC-1 HOG proteins in COGs, where possible (Table 3; Figure 3). However, the majority (86%) of accessory HOGs (aHOGs), protein clusters with peptide sequences from eight or fewer haloarchaea, were not members of any COGs or KOGs, or were members of poorly characterized COGs and could not be assigned to a functional class. Of the aHOGs that could be assigned functions based on COG-association, 3% were classified as being involved in information transfer and storage, or in cellular processing and signaling, and 8% were classified as being involved in metabolism.
In contrast, the great majority (89%) of cHOGs was associated with one or more COGs and KOGs, and a large fraction, 69%, was assigned to a functional class based on this criterion.
Among cHOGs, we classified 25% of the protein clusters as being involved in information transfer and storage . Half of the proteins in these cHOGs were involved in translation, ribosomal structure, and biogenesis, including 25 50S ribosomal subunit clusters (cHOG0202, cHOG0218, cHOG0230, cHOG0241, cHOG0248, cHOG0414, cHOG0415, cHOG0438, cHOG0478, cHOG0485, cHOG0512, cHOG0543, cHOG0560, cHOG0572, cHOG0579, cHOG0690, cHOG0700, cHOG0703, cHOG0737, cHOG0743, cHOG0745, cHOG0752, cHOG0753, cHOG0757, and cHOG0772), 21 30S ribosomal subunit clusters (cHOG0154, cHOG0271, cHOG0274, cHOG0379, cHOG0396, cHOG0539, cHOG0564, cHOG0655, cHOG0660, cHOG0675, cHOG0680, cHOG0692, cHOG0709, cHOG0726, cHOG0740, cHOG0750, cHOG0758, cHOG0760, cHOG0770, cHOG0771, and cHOG0774), and 13 amino-acyl tRNA synthetase clusters (cHOG0160, cHOG0184, cHOG0199, cHOG0250, cHOG0289, cHOG0306, cHOG0435, cHOG0468, cHOG0484, cHOG0487, cHOG0514, cHOG0536, and cHOG0672). In addition, we identified 11 cHOGs as containing RNA polymerase II-like enzymes (cHOG0165, cHOG0338, cHOG0407, cHOG0412, cHOG0492, cHOG0507, cHOG0679, cHOG0722, cHOG0741, cHOG0773, and cHOG0779), two containing origin recognition complex homologs (cHOG0234 and cHOG0244), three containing histone acetyltransferases (cHOG0049, cHOG0352, and cHOG0398), two containing transcription initiation factor IIB homologs (cHOG0004 and cHOG0018), and one containing transcription initiation factor IID homologs (cHOG0044).
An additional 10% of cHOG protein clusters was classified as being involved in cellular processing and signaling. Half of the proteins in these cHOGs were involved in posttranslational modification, protein turnover, or assembly, including two proteasome subunit clusters (cHOG0058 and cHOG0127), four heat shock protein clusters (cHOG0150, cHOG0156, cHOG0458, and cHOG0678), and two thermosome subunit clusters (cHOG0320 and cHOG0344). Three categories of COGs, nuclear structure, cytoskeleton, and extracellular structure, were not represented in any of the HOGs.
The largest number of cHOGs (33%) was classified as being involved in metabolism. Unlike cHOGs involved in information transfer and storage and cellular processes and signaling, there was no single category of metabolism that was overwhelmingly abundant. Four categories, energy production and conversion, amino acid transport and metabolism, nucleotide transport and metabolism, and coenzyme transport and metabolism, each contained over 40 cHOGs and accounted for 5% or more of the core clusters. Included in these cHOGs were nine ATP synthase subunit clusters (cHOG0124, cHOG0195, cHOG0233, cHOG0293, cHOG0302, cHOG0527, cHOG0600, cHOG0616, and cHOG0909), and ten NADH:ubiquinone oxidoreductase subunit clusters (cHOG0036, cHOG0126, cHOG0132, cHOG0187, cHOG0381, cHOG0453, cHOG0496, cHOG0599, cHOG0775, and cHOG0777).
The number of cHOGs associated with a cellular process did not necessarily correlate with the degree of conservation of that process. In particular, while there was a smaller number of cHOGs associated with information transfer and storage than metabolism, the proteins involved in information transfer and storage were more conserved in haloarchaea than those of metabolism or cellular processing and signaling. A large majority (65%) of HOGs associated with information transfer and storage was conserved in all nine genomes, whereas only 46% and 38% of the metabolism and cellular processing and signaling HOGs, respectively, were conserved in all of the genomes.
Newly sequenced haloarchaeal genomes
We also used BLAST analysis to determine if the cHOG proteins were conserved in four recently completed genomes (Table 4; see also Additional file 1). Homologs of the overwhelming majority of the cHOGs (784 out of 799) were identified in the recently completed genomes of Haloterrigena turkmenica, Halopiger xanaduensis, Natrialba magadii, and Halalkalicoccus jeotgali, with only six, two, five, and six clusters absent in these species, respectively (Table 4). Among the unique genes, five out of 60 tucHOGs and one of the 29 ucHOGs in the nine original genomes analyzed were absent in one or more of the four newer haloarchaeal genomes (Table 4).
Our current study has established core and unique haloarchaeal proteins and assigned likely functions to these conserved haloarchaeal proteins among sequenced haloarchaea. The core haloarchaeal orthologous groups (cHOGs) contained nearly 800 protein clusters that accounted for 21 - 33% of each predicted haloarchaeal proteome. The majority (89%) of the core proteins could be assigned specific or general functions based on association with NCBI KOGs and/or COGs, while the remainder (11%) were novel and could not be correlated to any previously known protein clusters. Based on further analysis of four recently sequenced haloarchaeal genomes and statistical analysis of alignments with non-haloarchaeal homologs, 55 protein clusters (named tucHOGs) were identified as haloarchaeal signature proteins.
The precise functions of the signature proteins are not clear because of their unique nature and the dearth of experimental studies focused on these genes. Only a single example among the truly unique haloarchaeal orthologous groups, Ral (tucHOG0456), was examined in any previous experimental work and was suggested to function in double-stranded DNA break repair and desiccation/radiation tolerance in the model haloarchaeon, Halobacterium sp. NRC-1 . Transcriptome analysis of both UV irradiated Halobacterium sp. NRC-1 and its highly ionizing radiation resistant mutants showed an up-regulation of the rfa3-rfa8-ral operon, consistent with their involvement in DNA repair and protection . Due to the transcriptional linkage of the three genes, and the presence of oligonucleotide binding (OB) folds in rfa3 and rfa8, the ral gene was also hypothesized to function as part of the eukaryotic-type single-stranded DNA binding RPA complex. However, analysis of the amino acid sequence of Ral did not reveal an OB fold domain, and it is not clear whether it serves as the third subunit of the RPA complex, replacing the RPA14 subunit found in eukaryotic organisms. While additional experimental studies are still required to determine the precise function of Ral, the possibility that it, as well as those of the other uniquely conserved haloarchaeal proteins, functions in adaptation of these organisms to their naturally extreme environments is an attractive hypothesis.
A somewhat larger (83) group of protein clusters, unique core haloarchaeal proteins (ucHOGs), includes 28 members which are nearly unique to haloarchaea (nucHOGs) and 55 which are truly unique to haloarchaea (tucHOGs). Our bioinformatic analysis of the ucHOGs suggested that they are quite typical of haloarchaeal proteins in pI, molecular weight, and GC-composition of their genes. The average pI of the ucHOGs is 4.7, similar to other haloarchaeal proteins (see Additional file 2). Similarly, the average G + C content of the ucHOGs are typical for each haloarchaeal chromosome (ranging from 68.5% for Halobacterium sp NRC-1 to 48.0% for H. walsbyi) (see Additional file 3). Their average molecular weight, 19.8 kDa, is somewhat smaller than predicted haloarchaeal proteins in general, 31 kDa (see Additional file 4). Their smaller size is consistent with their role as accessories to protein complexes, as suggested for the Ral protein in single-stranded DNA binding and DNA repair and protection. For example, as a group, ucHOGs may improve the activity or function of complexes in the cytoplasm with essentially saturating concentrations of KCl . The great majority of ucHOGs appear to be soluble proteins (unpublished data).
The genomic distribution of ucHOG protein genes was examined and they were found to map overwhelmingly on the chromosomes in all of the haloarchaeal microorganisms (see Additional file 1). In the case of Halobacterium sp. NRC-1, all of the tucHOGs and all but one of the nucHOGs were located on the chromosome (Figure 4). The haloarchaea do not contain more than one or at most two of these proteins on megaplasmids. These findings suggest that the ucHOG proteins serve integral functions in these microorganisms and are likely important and possibly critical for survival. In addition to their likely important function, the ucHOGs, and especially the signature proteins (tucHOGs) and their genes, will also be useful as markers for the presence of members of the Haloarchaeaceae family in the environment.
Of the 83 ucHOGs, 28 were not completely unique to haloarchaea, with one or a few homologs present in non-haloarchaea (see Additional file 1). A large fraction (46%) of the hits were to methanogenic Archaea belonging to the Methanosarcinaceae, Methanosaetaceae, and Methanocellaceae families, which are relatively close to haloarchaea based on phylogenetic analysis of 16S sequences and include some moderate halophiles . There were also a number of hits to halophilic bacteria, e.g., Salinibacter ruber, which may be the result of lateral gene transfer between species in a common environment . Of the clusters determined to not be uniquely haloarchaeal, 14 were associated with archaeal COGs (arCOGs) containing non-haloarchaeal homologs, consistent with their presence in more than a single family of Archaea  (see Additional file 1). This may reflect the distinct and common ancestry of the Archaea.
Prior to our study, an analysis of conserved proteins in the Archaea was first completed on eight archaeal genomes which did not include any haloarchaeal genomes . In this early study, 351 signature proteins present in at least two of the archaeal genomes were identified. In a subsequent study, 11 archaeal genomes were compared, including two haloarchaeal genomes . The number of signature proteins shared by all 11 genomes decreased to only six and an additional 30 were identified in the majority of archaeal genomes. In an analysis of four haloarchaeal genomes, 127 haloarchaeal-specific proteins were reported . Of these, we classified 51 as signature proteins or tucHOGs, 13 as nucHOGs, while the remaining 63 were either missing in one or more of the 13 haloarchaeal genomes or were associated with a COG (see Additional file 5). In another report, ten haloarchaeal genomes were recently compared and 112 'signature' clusters were reported , of which we found that 50 were similar to tucHOGs and 11 are like nucHOGs (see Additional file 6).
Several studies aimed at identifying signature proteins in other taxonomic groups have been conducted for organisms from other domains of life. Among bacteria, an analysis of actinobacterial genomes found 29 signature proteins present in the majority of genomes and an additional 204 that are found in some, but not all of the genomes . In another study , five Chlamydial genomes and one Parachlamydial genome were compared, and 59 proteins were conserved in all six genomes, coded by hypothetical genes with no known functions. Two subsequent studies of α-proteobacterial genomes reported signature proteins [33, 34]. Initially three genomes were compared and six signature proteins were identified in the majority of α-proteobacterial genomes and an additional 47 proteins were identified in some but not all subgroups . With the increase to 12 α-proteobacterial genomes, further work showed that only four of the original six signature proteins were present in all of the genomes . Among eukaryotes, 300 conserved signature proteins were identified in sequenced genomes, including the deeply branching Giardia lamblia species [35–37].
The entire set of genes within a given species or group of organisms, in essence, the combination of the core and all dispensable genes, is sometimes referred to as the "pan-genome" . With this approach, as more whole genomes become available, the size of the pan-genome increases due to an increase in the number of accessory genes, while the size of the core-genome asymptotically reaches a minimum. While there are numerous studies of species level pan-genomes, there are only a few published studies at the genus or family level. A study of 26 genomes from the Streptococcus genus found that the core-genome contains 611 orthologous groups, which constituted 26 - 30% of any one genome . Analysis of 11 genomes from the Vibrionaceae family found the core-genome of 1,882 orthologous groups constituted 32 - 50% of these genomes . Analysis of six genomes from the Enterobacteriaceae family identified 2,125 core orthologous groups that accounted for 43 - 88% of these genomes .
Our result from this study of the Haloarchaeaceae family showed that 21 - 33% of each genome constituted the core-genome and was similar to the results reported in earlier studies on other groups. Moreover, the great majority of core orthologous groups identified in the first nine haloarchaea were conserved in the subsequent four sequenced species. Our preliminary results with analysis of the pan-genome of haloarchaea show an expanding number of dispensable genes among members of this group (data not shown). The sequencing of additional haloarchaeal genomes and metagenomes and further bioinformatic analysis are likely to yield additional insights into the genetic composition of this interesting group of extremophilic microorganisms .
The signature and core genes and proteins are valuable concepts for understanding phylogenetic and phenotypic characteristics of coherent groups of organisms. Our analysis of 13 haloarchaea from different genera has established that the haloarchaeal proteome consists of 4,455 orthologous groups (HOGs), 784 of which form the core proteome (cHOGs), and 55 of which constitute haloarchaeal signature proteins (tucHOGs). The conservation of the cHOG and tucHOG clusters suggests that they may be essential or vital for survival. An attractive hypothesis, similar to what has been suggested for Ral, the only tucHOG with a predicted function, is that these small, chromosomally encoded proteins may act as accessory proteins enhancing macromolecular function in extreme conditions.
Sources of nucleotide and protein sequences
Nucleotide and protein sequences were obtained for completed haloarchaeal genomes from NCBI: Halobacterium sp. NRC-1 ATCC 700922 (NRC-1) , Haloarcula marismortui ATCC 43049 (Hma) , Natronomonas pharaonis DSM 2160 (Nph) , Haloquadratum walsbyi DSM 16790 (Hwa) , Halorubrum lacusprofundi ATCC 49239 (Hla) , Halogeometricum borinquense DSM 11551 (Hbo) , Halomicrobium mukohataei DSM 12286 (Hmu) , Halorhabdus utahensis DSM 12940 (Hut) , Haloferax volcanii DS2 (Hvo) , Haloterrigena turkmenica DSM 5511 (Htu) , Natrialba magadii ATCC 43099 (Nma) , Halalkalicoccus jeotgali B3 (Hje) , and Halopiger xanaduensis SH-6 (Hxa) .
Construction of protein clusters
For the initial nine genomes, we used the method of Tatusov [23, 24] to determine best reciprocal hits and the program MUSCLE for multiple sequence alignments . Conserved protein clusters were used to construct orthologous groups using in-house Perl scripts and manual navigation of data stored in a MySQL database and served on our Linux-Apache servers (HaloWeb - http://halo4.umbi.umd.edu) . Subsequently, we analyzed four additional sequences using our HOGnitor, via BLAST analysis. Similar non-haloarchaeal proteins were identified with BLAST analysis using HOG proteins as query sequences against the NCBI non-redundant database (June 5, 2011 version).
Statistical analysis of proteins clusters
Significance of protein assignment to clusters was established by base composition-preserved randomized pairwise global alignments using the method of Needleman and Wunsch [26, 57]. Scores of paired alignments were compared to scores and standard deviation for 50 randomized sequences with base composition-preserved. Protein families displaying greater than 99.9999% confidence were grouped into haloarchaeal orthologous groups (HOGs), and families with similar non-haloarchaeal proteins displaying greater than 99.0% confidence were grouped into nearly unique haloarchaeal orthologous groups (nucHOGs) [25, 58].
Correlation with COGs, KOGs, and arCOGs of haloarchaeal orthologous groups and functional classification
Haloarchaeal orthologous groups or HOGs were correlated with prokaryotic (COGs) and eukaryotic (KOGs) orthologous groups at NCBI using one of three methods: (1) HOGs were correlated to COGs using the Halobacterium sp. NRC-1 COGs as reference [23, 24]. (2) COGs and KOGs were correlated based on the Saccharomyces cerevisiae predicted proteins. (3) HOGs associated KOGs were also identified using the KOGnitor tool . HOGs were correlated with the clusters of archaeal orthologous groups (arCOGs) based on Halobacterium sp. NRC-1 proteins .
Genomic and protein analysis
Genomic analysis was conducted using tools available on our HaloWeb servers . Protein analysis was carried out using either stand-alone Perl scripts or Perl scripts running the Wisconsin Package protein analysis programs . Chromosome maps were generated using either our HaloWeb servers or GenomeVx software [56, 60].
Grant WD, Kamekura M, McGenity TJ, Ventosa A: Class III. Halobacteria class. nov. Bergey's Manual of Systematic Bacteriology. Edited by: Boone D, Castenholz R. 2001, New York: Springer, 2:
DasSarma S, DasSarma P: Halophiles. Encyclopedia of Life Sciences. 2012, John Wiley & Sons, Ltd
DasSarma P, Coker J, Huse V, DasSarma S: Halophiles, Biotechnology. Encyclopedia of Industrial Biotechnology, Bioprocess, Bioseparation, and Cell Technology. Edited by: Flickinger M. 2010, Wiley, John & Sons, Incorporated, 2769-2777.
Horneck G, Klaus DM, Mancinelli RL: Space microbiology. Microbiol Mol Biol Rev. 2010, 74: 121-156. 10.1128/MMBR.00016-09.
Vreeland RH, Jones J, Monson A, Rosenzweig WD, Lowenstein TK, Timofeeff M, Satterfield C, Cho BC, Park JS, Wallace A, Grant WD: Isolation of Live Cretaceous (121-112 Million Years Old) Halophilic Archaea from Primary Salt Crystals. Geomicrobiology Journal. 2007, 24: 275-282. 10.1080/01490450701456917.
Pääbo S, Poinar H, Serre D, Jaenicke-Després V, Hebler J, Rohland N, Kuch M, Krause J, Vigilant L, Hofreiter M: Genetic analysis from ancient DNA. Annual Review of Genetics. 2004, 38:
DasSarma P, DasSarma S: On the origin of prokaryotic "species": the taxonomy of halophilic Archaea. Saline Systems. 2008, 4: 5-10.1186/1746-1448-4-5.
Ng WV, Kennedy SP, Mahairas GG, Berquist B, Pan M, Shukla HD, Lasky SR, Baliga NS, Thorsson V, Sbrogna J, Swartzell S, Weir D, Hall J, Dahl TA, Welti R, Goo YA, Leithauser B, Keller K, Cruz R, Danson MJ, Hough DW, Maddocks DG, Jablonski PE, Krebs MP, Angevine CM, Dale H, Isenbarger TA, Peck RF, Pohlschroder M, Spudich JL, Jung KW, Alam M, Freitas T, Hou S, Daniels CJ, Dennis PP, Omer AD, Ebhardt H, Lowe TM, Liang P, Riley M, Hood L, DasSarma S: Genome sequence of Halobacterium species NRC-1. Proc Natl Acad Sci USA. 2000, 97: 12176-12181. 10.1073/pnas.190337797.
DasSarma S: Genome sequence of an extremely halophilic archaeon. Microbial Genomes. Edited by: Fraser CM, Read TD, Nelson KE. 2004, Totowa, NJ: Humana Press, Inc, 383-399.
Kennedy SP, Ng WV, Salzberg SL, Hood L, DasSarma S: Understanding the adaptation of Halobacterium species NRC-1 to its extreme environment through computational analysis of its genome sequence. Genome Res. 2001, 11: 1641-1650. 10.1101/gr.190201.
Paul S, Bag SK, Das S, Harvill ET, Dutta C: Molecular signature of hypersaline adaptation: insights from genome and proteome composition of halophilic prokaryotes. Genome Biol. 2008, 9: R70-10.1186/gb-2008-9-4-r70.
Dym O, Mevarech M, Sussman JL: Structural features that stabilize halophilic malate dehydrogenase from an Archaebacterium. Science. 1995, 267: 1344-1346. 10.1126/science.267.5202.1344.
Richard SB, Madern D, Garcin E, Zaccai G: Halophilic adaptation: novel solvent protein interactions observed in the 2.9 and 2.6 Å resolution structures of the wild type and a mutant of malate dehydrogenase from Haloarcula marismortui. Biochemistry. 2000, 39: 992-1000. 10.1021/bi991001a.
Pieper U, Kapadia G, Mevarech M, Herzberg O: Structural features of halophilicity derived from the crystal structure of dihydrofolate reductase from the Dead Sea halophilic archaeon, Haloferax volcanii. Structure. 1998, 6: 75-88. 10.1016/S0969-2126(98)00009-4.
Winter JA, Christofi P, Morroll S, Bunting KA: The crystal structure of Haloferax volcanii proliferating cell nuclear antigen reveals unique surface charge characteristics due to halophilic adaptation. BMC Struct Biol. 2009, 9: 55-10.1186/1472-6807-9-55.
Arakawa T, Tokunaga M: Electrostatic and hydrophobic interactions play a major role in the stability and refolding of halophilic proteins. Protein Pept Lett. 2004, 11: 125-132. 10.2174/0929866043478220.
Britton KL, Baker PJ, Fisher M, Ruzheinikov S, Gilmour DJ, Bonete MJ, Ferrer J, Pire C, Esclapez J, Rice DW: Analysis of protein solvent interactions in glucose dehydrogenase from the extreme halophile Haloferax mediterranei. Proc Natl Acad Sci USA. 2006, 103: 4846-4851. 10.1073/pnas.0508854103.
Capes MD, Coker JA, Gessler R, Grinblat-Huse V, DasSarma SL, Jacob CG, Kim JM, DasSarma P, DasSarma S: The information transfer system of halophilic archaea. Plasmid. 2011, 65: 77-101. 10.1016/j.plasmid.2010.11.005.
Anderson I, Scheuner C, Goker M, Mavromatis K, Hooper SD, Porat I, Klenk HP, Ivanova N, Kyrpides N: Novel insights into the diversity of catabolic metabolism from ten haloarchaeal genomes. PLoS One. 2011, 6: e20237-10.1371/journal.pone.0020237.
Khomyakova M, Bukmez O, Thomas LK, Erb TJ, Berg IA: A methylaspartate cycle in haloarchaea. Science. 2011, 331: 334-337. 10.1126/science.1196544.
DasSarma S, Capes M, DasSarma P: Haloarchaeal megaplasmids. Megaplasmids. Edited by: Schwartz E. 2008, Berlin: Springer-Verlag Berlin and Heidelberg GmbH & Co. KG, 3-30.
Wang G, Kennedy SP, Fasiludeen S, Rensing C, DasSarma S: Arsenic resistance in Halobacterium sp. strain NRC-1 examined by using an improved gene knockout system. J Bacteriol. 2004, 186: 3187-3194. 10.1128/JB.186.10.3187-3194.2004.
Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science. 1997, 278: 631-637. 10.1126/science.278.5338.631.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41-10.1186/1471-2105-4-41.
Doolittle RF: OF URFS AND ORFS A Primer on How to Analyze Derived Amino Acid Sequences. 1987, Mill Valley: University Science Books
Lipman DJ, Wilbur WJ, Smith TF, Waterman MS: On the statistical significance of nucleic acid similarities. Nucleic Acids Res. 1984, 12: 215-226. 10.1093/nar/12.1Part1.215.
DeVeaux LC, Müller JA, Smith J, Petrisko J, Wells DP, DasSarma S: Extremely radiation-resistant mutants of a halophilic archaeon with increased single-stranded DNA-binding protein (RPA) gene expression. Radiat Res. 2007, 168: 507-514. 10.1667/RR0935.1.
Makarova KS, Sorokin AV, Novichkov PS, Wolf YI, Koonin EV: Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea. Biol Direct. 2007, 2: 33-10.1186/1745-6150-2-33.
Graham DE, Overbeek R, Olsen GJ, Woese CR: An archaeal genomic signature. Proc Natl Acad Sci USA. 2000, 97: 3304-3308. 10.1073/pnas.050564797.
Gao B, Gupta RS: Phylogenomic analysis of proteins that are distinctive of Archaea and its main subgroups and the origin of methanogenesis. BMC Genomics. 2007, 8: 86-10.1186/1471-2164-8-86.
Gao B, Paramanathan R, Gupta RS: Signature proteins that are distinctive characteristics of Actinobacteria and their subgroups. Antonie Van Leeuwenhoek. 2006, 90: 69-91. 10.1007/s10482-006-9061-2.
Griffiths E, Ventresca MS, Gupta RS: BLAST screening of chlamydial genomes to identify signature proteins that are unique for the Chlamydiales, Chlamydiaceae, Chlamydophila and Chlamydia groups of species. BMC Genomics. 2006, 7: 14-10.1186/1471-2164-7-14.
Kainth P, Gupta RS: Signature proteins that are distinctive of alpha proteobacteria. BMC Genomics. 2005, 6: 94-10.1186/1471-2164-6-94.
Gupta RS, Mok A: Phylogenomics and signature proteins for the alpha proteobacteria and its main groups. BMC Microbiol. 2007, 7: 106-10.1186/1471-2180-7-106.
Hartman H, Fedorov A: The origin of the eukaryotic cell: a genomic investigation. Proc Natl Acad Sci USA. 2002, 99: 1420-1425. 10.1073/pnas.032658599.
Kurland CG, Collins LJ, Penny D: Genomics and the irreducible nature of eukaryote cells. Science. 2006, 312: 1011-1014. 10.1126/science.1121674.
Kersey PJ, Duarte J, Williams A, Karavidopoulou Y, Birney E, Apweiler R: The International Protein Index: an integrated database for proteomics experiments. Proteomics. 2004, 4: 1985-1988. 10.1002/pmic.200300721.
Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, Angiuoli SV, Crabtree J, Jones AL, Durkin AS, Deboy RT, Davidsen TM, Mora M, Scarselli M, Margarit y, Ros I, Peterson JD, Hauser CR, Sundaram JP, Nelson WC, Madupu R, Brinkac LM, Dodson RJ, Rosovitz MJ, Sullivan SA, Daugherty SC, Haft DH, Selengut J, Gwinn ML, Zhou L, Zafar N, Khouri H, Radune D, Dimitrov G, Watkins K, O'Connor KJ, Smith S, Utterback TR, White O, Rubens CE, Grandi G, Madoff LC, Kasper DL, Telford JL, Wessels MR, Rappuoli R, Fraser CM: Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome". Proc Natl Acad Sci USA. 2005, 102: 13950-13955. 10.1073/pnas.0506758102.
Lefébure T, Stanhope M: Evolution of the core and pan-genome of Streptococcus: positive selection, recombination, and genome composition. Genome Biology. 2007, 8: R71-10.1186/gb-2007-8-5-r71.
Lilburn T, Gu J, A Cai H, A Wang Y: Comparative genomics of the family Vibrionaceae reveals the wide distribution of genes encoding virulence-associated proteins. BMC Genomics. 2010, 11: 369-10.1186/1471-2164-11-369.
Uchiyama I: Multiple genome alignment for identifying the core structure among moderately related microbial genomes. BMC Genomics. 2008, 9: 515-10.1186/1471-2164-9-515.
Narasingarao P, Podell S, Ugalde JA, Brochier-Armanet C, Emerson JB, Brocks JJ, Heidelberg KB, Banfield JF, Allen EE: De novo metagenomic assembly reveals abundant novel major lineage of Archaea in hypersaline microbial communities. The ISME Journal. 2011
Baliga NS, Bonneau R, Facciotti MT, Pan M, Glusman G, Deutsch EW, Shannon P, Chiu Y, Weng RS, Gan RR, Hung P, Date SV, Marcotte E, Hood L, Ng V: Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea. Genome Res. 2004, 14: 2221-2234. 10.1101/gr.2700304.
Falb M, Pfeiffer F, Palm P, Rodewald K, Hickmann V, Tittor J, Oesterhelt D: Living with two extremes: conclusions from the genome sequence of Natronomonas pharaonis. Genome Res. 2005, 15: 1336-1343. 10.1101/gr.3952905.
Bolhuis H, Palm P, Wende A, Falb M, Rampp M, Rodriguez-Valera F, Pfeiffer F, Oesterhelt D: The genome of the square archaeon Haloquadratum walsbyi : life at the limits of water activity. BMC Genomics. 2006, 7: 169-10.1186/1471-2164-7-169.
Franzmann PD, Stackebrandt E, Sanderson K, Volkman JK, Cameron DE, Stevenson PL, McMeekin TA, Burton HR: Halobacterium lacusprofundi sp. nov., a halophilic bacterium isolated from Deep Lake, Antarctica. Syst Appl Microbiol. 1988, 11: 20-27.
Malfatti S, Tindall B, Schneider S, Fähnrich R, Lapidus A, Labutti K, Copeland A, Glavina del Rio T, Nolan M, Chen F, Lucas S, Tice H, Cheng J-F, Bruce D, Goodwin L, Pitluck S, Anderson I, Pati A, Ivanova N, Mavrommatis K, Chen A, Palaniappan K, D'Haeseleer P, Göker M, Bristow J, Eisen J, Markowitz V, Hugenholtz P, Kyrpides N, Klenk H, Chain P: Complete genome sequence of Halogeometricum borinquense type strain (PR3T). Standards in Genomic Sciences. 2009, 1:
Tindall BJ, Schneider S, Lapidus A, Copeland A, Rio TGD, Nolan M, Lucas S, Chen F, Tice H, Cheng J-F, Saunders E, Bruce D, Goodwin L, Pitluck S, Mikhailova N, Pati A, Ivanova N, Mavrommatis K, Chen A, Palaniappan K, Chain P, Land M, Hauser L, Chang Y-J, Jeffries CD, Brettin T, Han C, Rohde M, Göker M, Bristow J, Eisen JA, Markowitz V, Hugenholtz P, Klenk H-P, Kyrpides NC, Detter JC: Complete genome sequence of Halomicrobium mukohataei type strain (arg-2T). Standards in Genomic Sciences. 2009, 1:
Bakke P, Carney N, Deloache W, Gearing M, Ingvorsen K, Lotz M, McNair J, Penumetcha P, Simpson S, Voss L, Win M, Heyer LJ, Campbell AM: Evaluation of three automated genome annotations for Halorhabdus utahensis. PLoS One. 2009, 4: e6291-10.1371/journal.pone.0006291.
Hartman AL, Norais C, Badger JH, Delmas S, Haldenby S, Madupu R, Robinson J, Khouri H, Ren Q, Lowe TM, Maupin-Furlow J, Pohlschroder M, Daniels C, Pfeiffer F, Allers T, Eisen JA: The complete genome sequence of Haloferax volcanii DS2, a model archaeon. PLoS One. 2010, 5: e9605-10.1371/journal.pone.0009605.
Saunders E, Tindall BJ, Fähnrich R, Lapidus A, Copeland A, Rio TGD, Lucas S, Chen F, Tice H, Cheng J-F, Han C, Detter JC, Bruce D, Goodwin L, Chain P, Pitluck S, Pati A, Ivanova N, Mavromatis K, Chen A, Palaniappan K, Land M, Hauser L, Chang Y-J, Jeffries CD, Brettin T, Rohde M, Göker M, Bristow J, Eisen JA, Markowitz V, Hugenholtz P, Klenk H-P, Kyrpides NC: Complete genome sequence of Haloterrigena turkmenica type strain (4kT). Standards in Genomic Sciences. 2010, 2: 107-116. 10.4056/sigs.681272.
Kamekura M, Dyall-Smith ML, Upasani V, Ventosa A, Kates M: Diversity of alkaliphilic halobacteria: proposals for transfer of Natronobacterium vacuolatum, Natronobacterium magadii, and Natronobacterium pharaonis to Halorubrum, Natrialba, and Natronomonas gen. nov., respectively, as Halorubrum vacuolatum comb. nov., Natrialba magadii comb. nov., and Natronomonas pharaonis comb. nov., respectively. Int J Syst Bacteriol. 1997, 47: 853-857. 10.1099/00207713-47-3-853.
Roh SW, Nam YD, Nam SH, Choi SH, Park HS, Bae JW: Complete genome sequence of Halalkalicoccus jeotgali B3T, an extremely halophilic archaeon. J Bacteriol. 2010, 192: 4528-4529. 10.1128/JB.00663-10.
Gutiérrez MC, Castillo AM, Kamekura M, Xue Y, Ma Y, Cowan DA, Jones BE, Grant WD, Ventosa A: Halopiger xanaduensis gen. nov., sp. nov., an extremely halophilic archaeon isolated from saline Lake Shangmatala in Inner Mongolia, China. Int J Syst Evol Microbiol. 2007, 57: 1402-1407. 10.1099/ijs.0.65001-0.
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.
DasSarma SL, Capes MD, DasSarma P, DasSarma S: HaloWeb: the haloarchaeal genomes database. Saline Systems. 2010, 6: 12-10.1186/1746-1448-6-12.
Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970, 48: 443-453. 10.1016/0022-2836(70)90057-4.
Lipman DJ, Wilbur WJ: Interaction of silent and replacement changes in eukaryotic coding sequences. J Mol Evol. 1984, 21: 161-167.
Devereux J, Haeberli P, Smithies O: A comprehensive set of sequence analysis programs for the VAX. Nucleic Acids Res. 1984, 12: 387-395. 10.1093/nar/12.1Part1.387.
Conant GC, Wolfe KH: GenomeVx: simple web-based creation of editable circular chromosome maps. Bioinformatics. 2008, 24: 861-862. 10.1093/bioinformatics/btm598.
Acknowledgements and funding
We thank Satyajit L. DasSarma for assistance with databases and servers and R. Tatusov for the initial round of reciprocal best hit analysis. This work was supported by grants from the Henry M. Jackson Foundation grant HU0001-09-1-0002-660883 and the National Aeronautics and Space Administration grant NNX10AP47G.
The authors declare that they have no competing interests.
MDC, PD, and SD contributed to the bioinformatic analysis and writing the manuscript. All authors read and approved the final manuscript.