BACKGROUND
archaea are widely regarded as one of the three main domains of life  <cit> , although their origin is a subject of debate  <cit> .
archaeal species were earlier believed to inhabit only extreme environments such as extremely hot, or hot and acidic, extremely saline, or very acidic or alkaline conditions  <cit> .
however, recent studies provide evidence that they are widespread in different environments  <cit> .
the archaea also include methanogens, which grow under strictly anaerobic and often thermophilic conditions, and are the only organisms that derive all of their metabolic energy by reduction of co <dig> by hydrogen to produce methane  <cit> .
the archaeal species branch distinctly from all other organisms in phylogenetic trees based on 16s rrna and many other gene/protein sequences  <cit> .
in addition, many morphological or physiological characteristics such as the presence of branched-chain ether-linked lipids in their cell membrane, lack of peptidoglycan in their cell wall, characteristic subunit pattern of rna polymerase, presence of modified bases in trna, presence of a unique form of dna polymerase, have been previously indicated as defining characteristics of archaea  <cit> .
however, as noted by walsh and doolittle  <cit> , many of these features are either not shared by all archaea or they are also present in various eukaryotes or some thermophilic bacteria, indicating that they do not constitute distinctive characteristics of all archaea.
the phylogenetic analyses of archaea have led to their division into two major groups or phyla designated as crenarchaeota and euryarchaeota  <cit> .
the crenarchaeota species have also been referred to as 'eocytes' by lake and coworkers  <cit> .
the species from both these groups, particularly euryarchaeota, are highly diverse in terms of their metabolism and physiology.
based on their metabolic and physiological characteristics and other unique features, five functionally distinct groups within euryarchaeota are currently recognized: methanogens, sulfate reducers, extreme halophiles, cell wall-less archaea, and extremely thermophilic sulfur metabolizing archaea  <cit> .
some of these groups, such as methanogens, are polyphyletic in different phylogenetic trees  <cit> .
however, the sets of genes or proteins that are unique to these different functional groups and distinguish them from all others remain to be identified.
in the past  <dig> years, complete genomes of many archaeal species  covering all major divisions within the archaea have been sequenced .
comparative analyses of these sequences provide a valuable resource for identifying different genes/proteins that are distinctive characteristics of various taxonomic and functional groups within archaea  <cit> .
whole proteins that are uniquely present in particular groups or subgroups of organisms but not found anywhere else provide valuable molecular markers for taxonomic, phylogenetic and biochemical studies.
these proteins, which we refer to as signature proteins in our work, and others have called them as orfans or conserved hypothetical proteins, are present at different phylogenetic depths, such as genus, family, order or even phylum  <cit> .
in our recent work, a large number of such proteins that are distinctive characteristics of several groups within bacteria , and also their subgroups, were identified  <cit> .
these proteins provide not only valuable molecular markers for identifying and circumscribing species belonging to these major groups  in molecular terms, but their species distribution pattern also provides useful information about the branching order within these groups.
as archaea constitute a very diverse group, identification of sets of proteins that are specific for its main groups and subgroups should prove useful in terms of identifying molecular characteristics that are unique to them.
additionally, this information should also be helpful in understanding the evolutionary relationships among different groups.
comparative studies on limited numbers of archaeal genomes have been carried out by a number of investigators using different criteria.
graham et al.
<cit>  analyzed  <dig> archaeal genomes to identify signature proteins that function uniquely within the archaea.
their definition of an archaeal signature protein required it to be present in only two different euryarchaeal species and they identified  <dig> archaeal signature proteins.
makarova and koonin  <cit>  have analyzed archaeal genomes to identify core sets of genes, which are present in all archaeal species, but which are not restricted to the archaeal species.
recently, walsh and doolittle have analyzed prokaryotic genomes to measure dissimilarity between archaea and bacteria  <cit> .
although it was reported that 28% of the proteins from archaeal genomes are restricted to the archaea, specific proteins that were present in different groups of archaea were not identified.
other comparative studies using different criteria have been conducted on smaller groups within archaea such as pyrococcus, sulfolobus and thermoacidophilic organisms .
however, thus far no comprehensive phylogenomics study on different archaeal genomes has been carried out using the same standard criteria to identify proteins or orfs that are shared by all archaea or its different major lineages.
in this study we have carried out comparative analyses of archaeal genomes using uniform criteria to identify proteins that are uniquely present in archaeal species at different phylogenetic depths  representing all major groups within the archaea.
RESULTS
a. phylogenetic analyses of archaeal species
prior to undertaking comparative studies on archaeal genomes, phylogenetic analysis of sequenced archaeal species was carried out so that the results of phylogenomics analyses could be compared with those obtained by traditional phylogenetic approaches.
phylogenetic trees for the archaeal species based on 16s rrna as well as concatenated sequences of translation and transcription-related proteins have been published by other investigators  <cit> .
in the present work, we have constructed phylogenetic trees for  <dig> archaeal species  using a set of  <dig> universally distributed proteins that are involved in a broad range of functions  <cit> .
the sequence of haloquadratum walsbyi dsm  <dig>  which became available afterward, was not included in these studies.
phylogenetic trees based on a concatenated sequence alignment of these proteins were constructed using the neighbour-joining , maximum-likelihood  and maximum-parsimony  methods.
the results of these analyses are presented in fig.
<dig>  all three methods gave very similar tree topologies except for the branching positions of m. kandleri and methanospirillum hungatei, which were found to be variable.
except for this, the branching pattern of the archaeal species based on our dataset is very similar to that reported by gribaldo et al.
<cit>  based on concatenated sequences of translation and transcription-related proteins.
in the tree shown, the crenarchaeota and euryarchaeota, the two major phyla within archaea were clearly distinguished from each other.
the phylogenetic affinity of nanoarchaeum, which has a long-branch length, was not resolved in this or various other trees  <cit> .
within crenarchaeota, pyrobaculum was indicated to be a deeper branch, and aeropyrum branched in between the pyrobaculum and sulfolobus.
within euryarchaeota, the clades corresponding to halobacteria, thermococci and thermoplasmata were resolved with high bootstrap scores, but the methanogens were split into 2– <dig> clusters.
one of these clusters that has low bootstrap score consisted of methanobacteriales and methanococcales with m. kandleri  branching in its vicinity  <cit> .
the second cluster, with higher bootstrap score, showed a grouping of methanomicrobiales and methanosarcinales.
these two clusters, which are separated by thermoplasmata, archaeoglobi and halobacteria, have been referred to as class i and class ii methanogens by bapteste et al.
<cit> .
b. phylogenomic analyses of archaeal genomes
to search for proteins , which are uniquely present in either all archaea or various subgroups of them, blast searches were performed on each open reading frame  from a total of  <dig> archaeal genomes .
these genomes included  <dig> crenarchaeota   <cit>  and  <dig> divergent euryarchaeota species covering all main functional and phylogenetic groups .
the euryarchaeota genomes analyzed included: pyrococcus abyssi from extremely thermophilic sulfur metabolizing archaea  <cit> , methanococcus maripaludis  <cit>  from methanococcales, halobacterium sp.
nrc- <dig> and h. walsbyi from extreme halophiles  <cit> , thermoplasma acidophilum and picrophilus torridus belonging to the cell wall-less archaea  <cit> , methanococcoides burtonii from methanosarcinales and methanopyrus kandleri from the methanopyrales order  <cit> .
the chosen genomes should provide information regarding all archaeal proteins that are shared at a taxonomic level higher than a genus.
the analysis of the remainder of the genomes, which was expected to provide information regarding proteins that are only unique to a given species, was not carried out.
each orf from these genomes was examined by means of blastp and psi-blast searches against all available sequences from different organisms to identify proteins that are specific for only archaeal lineages.
the methods and the criteria that we have used to identify proteins that are specific for either all or various subgroups of archaea are described in the methods section.
generally, a protein was considered to be specific for a given archaeal lineage if all significant hits or alignments in the blastp and psi-blast searches with the query protein were from the indicated group of archaeal species.
in a few cases, where 1– <dig> isolated species from other groups also exhibited significant similarity, such proteins were retained as they provide interesting examples of lateral gene transfer  from archaea to other groups.
our analyses have identified  <dig> proteins that are unique to different groups of archaea and for which no homologues are generally found in any bacterial or eukaryotic species.
based on their specificity for different taxonomic groups, these proteins have been divided into a number of different groups .
a brief description of the different subsets of archaeal-specific proteins and functional information regarding them, where known, is given below.
in the description of these proteins that follows, the 'ape', 'hq', 'mbu', 'mk', 'mmp', 'pab', 'pae', 'pto', 'saci', 'ta', 'vng', and 'neq' part of the descriptors in proteins indicate that the original query protein sequence was from the genome of a. pernix k <dig>  h. walsbyi dsm  <dig>  m. burtonii dsm  <dig>  m. kandleri av <dig>  m. maripaludis s <dig>  p. abyssi ge <dig>  p. aerophilum str.
im <dig>  p. torridus dsm  <dig>  s. acidocaldarius dsm  <dig>  t. acidophilum dsm  <dig>  halobacterium sp.
nrc- <dig> and n. equitans, respectively.
proteins that are specific for all archaea
of the proteins that are uniquely present in all archaea, pab <dig> corresponds to trna nucleotidyltransferase , which builds and repairs the 3' end of trna  <cit> .
functionally similar enzymes are also present in bacteria and eukaryotes , but their sequences share very little homology with the archaeal cca-adding enzyme , which explains why no homologs were detected in any bacteria or eukaryotes in blast searches.
the main mechanistic difference between class i and class ii enzymes is that the trna substrate is required to fully define the nucleotide binding site in class i enzyme, whereas class ii has a preformed nucleotide binding site that recognizes ctp and atp in the absence of trna  <cit> .
another protein pab <dig> is assigned as archaeal type dna primase, which also has its synonymous counterparts in bacterial and eukaryotic species, but shows very little homology to them  <cit> .
in the same way, protein pab <dig> is annotated as a pilt family atpase, which showed very little similarity to bacterial atpases involved in type iv pili biogenesis  <cit> .
further studies of this protein could provide insights into novel aspects of the archaeal flagellar system.
a number of other proteins viz.
pab <dig>  pab0018a, pab <dig>  pab <dig> and pab <dig>  have also been assigned putative functions based on sequence analysis, but their exact roles in archaeal cells remains to be determined.
interestingly, for protein pab <dig>  two gene copies with acceptable e-values are also present in the genomes of dehalococcoides ethenogenes  <dig>  dehalococcoides sp.
cbdb <dig> and dehalococcoides sp.
bav <dig>  which belong to chloroflexi  <cit> .
because no homologue of pab <dig> is present in other bacteria, it is likely that this protein was transferred from archaea to the common ancestor of dehalococcoides followed by a gene duplication event.
proteins that are specific for crenarchaeota
as mentioned in the introduction, the archaea are divided into  <dig> main groups, crenarchaeota and euryarchaeota, based on 16s rrna trees as well many other gene trees and characteristics.
the crenarchaeota are also indicated to differ from euryarchaeota in terms of their ribosome structure  <cit> .
in comparison to euryarchaeota, which contain physiologically and metabolically diverse groups of organisms, the crenarchaeota were thought to be a pure collection of extreme thermophiles and most members metabolize sulfur.
however, recent studies indicate that crenarchaeota are much more diverse in their physiology and ecology than was previously believed  <cit> .
many species living in the cold ocean also belong to this group based on their branching pattern in 16s rrna trees, although most of them have not been cultivated  <cit> .
currently, this phylum is comprised of one single class thermoprotei containing three orders: thermoproteales, desulfurococcales and sulfolobales.
fortunately, every order has a completely sequenced representative  <cit> , which provide a platform to explore the characteristics that are unique to crenarchaeal species.
comparative genomic surveys have revealed some molecular features that are shared by crenarchaea but not euryarchaea, such as the lack of histones, absence of the ftsz-mincde system and distinctive rrna operon organization  <cit> .
lake et al.
have also identified distinctive differences in ribosome structure and an insert in elongation factor ef-g and ef-tu, which can be used to distinguish crenarchaeota from euryarchaeota  <cit> .
however, these features are not unique characteristics of the crenarchaeota.
blast searches on each orf from the genomes of a. pernix and s. acidocaldarius dsm  <dig>  <cit>  have identified  <dig> proteins which are shared by all five crenarchaeal species, but whose homologs are not found in other archaea, or any bacteria or eukaryotes with only  <dig> exceptions ).
the genes for these proteins likely evolved in a common ancestor of the crenarchaeota and they provide potential molecular markers for species from this phylum.
additionally,  <dig> proteins that are listed in table  <dig> are only found in a. pernix and three sulfolobus genomes.
these proteins suggest that aeropyrum and sulfolobus may have shared a common ancestor exclusive of pyrobaculum.
however, we have also come across  <dig> proteins that are shared by aeropyrum and pyrobaculum ) and  <dig> proteins that are exclusively present in the  <dig> sulfolobus species and pyrobaculum ).
hence, based upon the species distributions of these proteins, the relationships among the aeropyrum, sulfolobales and pyrobaculum are not entirely clear .
in phylogenetic trees thermoproteales  branches consistently earlier than desulfurococcales  and sulfolobales   <cit> .
this observation in conjunction with the fact that aeropyrum and sulfolobus share larger numbers of proteins in common with each other suggests that these two groups likely shared a common ancestor exclusive of pyrobaculum .
the proteins that are only found in aeropyrum and pyrobaculum, or in sulfolobus and pyrobaculum, most likely evolved in a common ancestor of the crenarchaea, but were subsequently lost in either the sulfolobales or a. pernix lineages.
in addition to these proteins that are uniquely present in either all sequenced crenarchaeota genomes or different groups of crenarchaeota species, these analyses have also identified  <dig> proteins that are unique for the sulfolobales species .
of these,  <dig> proteins are present in all  <dig> sequenced sulfolobus genomes, whereas the remaining  <dig> are present in at least two of the three sulfolobus genomes.
in this work, since blast analyses were not carried out on all three sulfolobus genomes, it is likely that the numbers of genes or proteins that are uniquely shared by only two sulfolobus genomes is much higher than indicated here.
chen et al.
<cit>  have previously analyzed the genome of s. acidocaldarius dsm  <dig> and indicated the presence of  <dig> genes that were specific for crenarchaeota and  <dig> genes that were specific to sulfolobus genus.
however, in the present work, relatively few genes that are uniquely shared by various crenarchaeota species were identified.
this difference could be due to more stringent criteria that we have employed for identification of proteins that are specific to different groups.
the genome of thermofilum pendens hrk  <dig>  which belongs to thermoproteales, has also been partially sequenced and information for large numbers of genes/proteins from this species is available in the ncbi database.
by carrying out blast searches on each orf from p. aerophilum genome  <cit> , we have identified  <dig> proteins that are only found in the above  <dig> thermoproteales species .
the numbers of proteins shared by these two species will likely increase once complete genome of t. pendens becomes available.
many of these proteins are expected to provide markers for the thermoproteales order.
proteins that are specific for euryarchaeota
the euryarchaeota, which comprise a majority of the cultured and sequenced archaea, is a morphologically, metabolically and physiologically diverse collection of species as evidenced by the presence in this group of various methanogens, extreme halophiles, cell wall-less archaea and sulfate reducing microbes  <cit> .
no unique biochemical or molecular characteristic that is commonly shared by all of the different lineages is known.
the present study has identified  <dig> proteins that are only found in euryarchaeota species with  <dig> exceptions .
in this table, the first  <dig> proteins ) are present in most euryarchaeota species.
of these proteins, pab <dig> and pab <dig> were found in all sequenced euryarchaeota species.
pab <dig> was also present in n. equitans, supporting its placement within the euryarchaeota  <cit> .
the protein pab <dig> is annotated as archaeosine trna-ribosyltransferase , which catalyzes the exchange of guanine with a free 7-cyano-7-deazaguanine  base, as the first step in the biosynthesis of an archaea-specific modified base, archaeosine   <cit> .
it should be mentioned that there is another protein pab <dig> in the same genome, which is also annotated and experimentally confirmed as arctgt  <cit> .
the latter belongs to a family of proteins that are highly conserved in all archaea species  and some bacteria.
it seems that pab <dig> might be involved in rna modification since it possesses a pua domain , but its function is likely different from pab <dig>  the protein pab <dig>  which is annotated as dna polymerase ii large subunit, is highly conserved within euryarchaeota, but is not found anywhere else except in nanoarchaeum.
this enzyme is the major dna replicase in euryarchaeota and also a distinctive molecular marker for this group  <cit> .
the genes for the above proteins likely evolved in a common ancestor of euryarchaeota  and they provide molecular markers for this diverse group of organisms.
another  <dig> proteins listed in table  <dig> are found in almost all euryarchaeota, but they are missing in thermoplasmata.
their distribution suggests that either thermoplasmata is a deep branching lineage within euryarchaeota or that the genes for these proteins have been selectively lost from thermoplasmata  <cit> .
of these proteins, pab <dig> is also present in n. equitans supporting its placement with euryarchaeota.
five other proteins from the first two columns in table  <dig>  are also not found in the  <dig> thermococci species.
these results can again be explained by either selective loss of these genes from these particular groups or deeper branching of these lineages within the euryarchaeota species.
on the basis of proteins listed in table  <dig>  although one can infer that thermoplasmata and thermococci are deeper branching lineages within euryarchaeota in comparison to methanogens, their relative branching order cannot be resolved.
proteins that are specific for different main groups within euryarchaeota
proteins specific for methanogenic archaea and their various subgroups
currently, the methanogens form the largest group within the euryarchaeota.
they are distinguished from all other prokaryotes by their ability to obtain all or most of their energy via the reduction of co <dig> to methane or by the process of methanogenesis.
in the bergey's manual  <cit> , the methanogenes are divided into  <dig> distinct orders .
some studies have suggested that these organisms possess a set of unique enzymes which are responsible for methanogenesis, such as coenzyme m, factor  <dig> and methanopterin  <cit> .
however, no systematic study has been carried out thus far to identify proteins that are uniquely present in different methanogens.
our blast searches of proteins from different methanogens have led to identification of  <dig> proteins, which are uniquely found in various methanogenic archaea.
twenty of these  <dig> proteins are present in all sequenced methanogens, while  <dig> proteins are missing only in m. stadtmanae, which is a human intestinal inhabitant .
this archaeon generate methane by reduction of methanol with h <dig> and lacks many proteins present in the genomes of other methanogens  <cit> .
thus, it is highly likely that the  <dig> proteins missing in m. stadtmanae were selectively lost from this species.
therefore, it is very likely that the genes for these  <dig> proteins that are commonly shared by virtually all methanogens ) evolved in a common ancestor of all methanogens.
these analyses have also identified  <dig> proteins that are uniquely shared by various methanogens as well as a. fulgidus ).
the genes for these proteins likely evolved in a common ancestor of a. fulgidus and various methanogenic archaea and they point to a close relationship between these two groups of organisms .
ten additional proteins are present in a. fulgidus as well as various methanosarcinales and m. hungatei  ).
it is likely that the genes for these proteins also evolved in a common ancestor of a. fulgidus and various methanogenic archaea, but they were selectively lost in other methanogens.
of the proteins that are commonly shared by a. fulgidus and various methanogenic archaea, mmp <dig> is reported to be a novel repressor of nif and glna genes, which are involved in nitrogen assimilation  <cit> .
interestingly,  <dig> homologs of this protein are also found in  <dig> dehalococcoides species, but nowhere else, which are very likely due to lgt.
protein mmp <dig> is the ε-subunit of carbon-monoxide dehydrogenase complex, which is made up of five subunits in different methanogens  <cit> .
the epsilon subunits are required for the reversible oxidation of co to co <dig>  <cit> .
all of the other components could be found in a few bacterial species, while the ε-subunit is restricted to methanogenic archaea and a. fulgidus  <cit> .
protein mmp <dig> is identified as a transcriptional regulator with a helix-turn-helix  motif, but its exact role has not been reported.
among the genes that are uniquely shared by various methanogenic archaea , two large gene clusters responsible for methanogenesis are found.
the proteins mmp <dig>  mmp1560–mmp <dig> and mmp1566–mmp <dig>  are parts of an eight-component complex, coenzyme m methyltransferase , which catalyzes an energy-conserving, sodium-ion-translocating step in methanogenesis from h <dig> and co <dig>  <cit> .
m. maripaludis contains all of the known mtr subunits, but the gene coding for mtrf is fused into the n-terminal region of mtra  <cit> .
all other methanogenic archaeal genomes contain complete set of mtr genes.
it is of interest to note that for the protein mmp <dig> , homologues with low e-values are also found in two desulfitobacterium hafniense strains as well as in three rhizobiales species  .
these three rhizobiae species can use methyl halides as a sole source of carbon and energy, and all of them possess a set of cmu genes which are essential for methyl chloride degradation  <cit> .
in particular, the cmub protein which is homologous to mmp <dig> transfers a methyl group to methylcobalamin:h <dig> folate , which is analogous to the reverse of the reaction catalyzed by mtrh in archaea  <cit> .
in view of the sequence and functional similarity between mtrh and cmub proteins, it is likely that the mtrh gene was laterally transferred from a methanogenic archaeon to the common ancestor of the above three rhizobiae species to serve the new functional role.
the function of the laterally transferred mtrh related gene in d. hafniense is not known at present.
the proteins mmp1555–mmp <dig> in table  <dig> form another gene cluster, encoding the subunits of methyl-coenzyme m reductase .
this complex catalyzes the final reaction of the energy conserving pathway in which methylcoenzyme m and coenzyme b are converted to methane and the heterodisulfide com-s-s-cob  <cit> .
except for these proteins, the other proteins listed in table  <dig> are of putative or unknown functions.
it is likely that these proteins are involved in some aspects of methanogenesis or other unknown pathways unique to methanogenic archaea.
these proteins provide molecular markers for methanogens, which can be used for identification of new archaeal species capable of methane production.
the blast searches of the m. maripaludis  <cit>  and m. kandleri  <cit>  genomes have identified  <dig> proteins that are uniquely shared by all of the following species belonging to the orders methanobacteriales , methanococcales  and methanopyrales  ).
of these, only  <dig> proteins are present in m. stadtmanae, which is also a methanobacteriales that has lost most of its genes due to its adaptation to the human intestine  <cit> .
the genes for these  <dig> proteins likely evolved in a common ancestor of the above groups of methanogens , which corresponds to the cluster of methanogenic archaea referred to as "class i methanogens"  <cit> .
interestingly, these studies have also identified  <dig> proteins that are uniquely shared by these methanogenic orders and m. hungatei ), which branches distantly in phylogenetic trees  <cit> .
the unique presence of these proteins in these methanogens suggests that species from these groups shared a common ancestor exclusive of other methanogenic archaea .
fifteen additional proteins discovered in this work ) are uniquely present in m. kandleri and various methanobacteriales indicating that these two groups are more closely related to each other than the methanococcales .
we have also come across  <dig> proteins that are uniquely shared by methanococcales and methanobacteriales ), and  <dig> proteins that are only present in methanococcales and methanopyrales ).
the most likely explanation to account for the species distributions of these latter proteins is that their genes also originated in a common ancestor of the above three groups of methanogens, but were selectively lost in either the methanobacteriales or methanopyrales lineages.
these analyses have also identified  <dig> additional proteins that are uniquely present in all  <dig> methanosarcinales species ), as well as  <dig> proteins that are only found in various methanosarcinales and m. hungatei ).
lastly, these studies have also identified  <dig> proteins that are uniquely present in m. maripaludis and m. jannaschii ) and  <dig> proteins that are only present in m. burtonii and  <dig> methanosarcina species, all belonging to the methanosarcinaceae family )  indicating that they are likely distinctive characteristics of species from these groups.
of the proteins that are uniquely found in methanococcales, methanobacteriales, methanopyrales and methanomicrobiales,  <dig> proteins viz.
mmp1448–mmp <dig>  mmp <dig>  mmp1458–mmp <dig> and mmp <dig> are from a big gene cluster eha, which encodes the multisubunit membrane-bound  hydrogenase  <cit> .
two of these proteins, mmp <dig> and mmp <dig>  are only found in methanococcales ).
the whole eha operon is composed of  <dig> orfs in the genome of m. thermoautotrophicus and of these only these  <dig> proteins are restricted to these methanogens while the other subunits have counterparts in bacteria.
the precise roles of these  <dig> proteins, which are predicted to be integral membrane proteins in the hydrogenase complex, have not been determined  <cit> .
among the other proteins that are specific for these groups of methanogens, mmp <dig> and mmp <dig> are hmd homologs, which catalyze the reversible dehydrogenation of n <dig>  n10-methylenetetrahydromethanopterin  <cit> .
in the proteins that are specific for the methanococcales ), one large gene cluster  is found, but no information is available concerning its possible function.
except for these proteins, all other proteins that are specific for these methanogenic archaea are of unknown or putative function.
proteins that are specific for thermococci
thermococci are obligately thermophilic, strictly anaerobic cocci, which are able to convert elemental sulfur to hydrogen sulfide.
thus, they are so called "extremely thermophilic sulfur metabolizer", which comprise one of the main functional groups within euryarchaeota.
according to the bergey's manual  <cit> , the class thermococci contains only one family, thermococcaceae, consisting of  <dig> genera: thermococcus and pyrococcus.
currently,  <dig> species from this family have been completely sequenced   <cit> .
the blast searches on each protein from p. abyssi have identified  <dig> proteins that are shared by all  <dig> of these species ).
all of these proteins show high degree of conservation within thermococci and they do not have homologs in any other prokaryotes or eukaryotes except one possible lgt event .
the genes for these proteins have likely evolved in a common ancestor of the thermococci .
of these proteins, pab <dig> is annotated as tbp-interacting protein , which forms complex with tbp  to regulate transcription  <cit> .
it is known that the archaeal transcription machinery is strikingly similar to that in eukaryotes  <cit> , but no tbp-binding component was found in archaeal species until the discovery of the tip in t. kodakaraensis  <cit> .
most other themococci-specific proteins are of unknown function, although in a few cases limited similarity to domains in known protein families have been noted.
a number of proteins  are clustered together in the p. abyssi genome, and it is possible that they may form functional units and are involved in related functions.
cohen et al.
<cit>  have reported a large number of proteins which are restricted to the pyrococcus genus.
however, a number of proteins from their list are also found in t. kodakarensis kod <dig>  <cit> , whose genome was not available when their work was published.
some proteins are not specific for either pyrococcus or thermococci according to our criteria and some of them are only found in one species – p. abyssi.
our analysis of the p. abyssi ge <dig> genome has also identified  <dig> proteins that are unique to the pyrococcus genus ).
again, almost all of these proteins are of unknown function except pab <dig>  which is annotated as rnase p, but this annotation seems arbitrary as it does not show significant sequence similarity to known rnases.
the proteins that are uniquely found in the  <dig> pyrococcus genomes likely evolved in a common ancestor of this genus .
proteins that are specific for halobacteria
extreme halophiles constitute another major class within euryarchaeota.
they require 5– <dig> times the salinity of seawater  for optimal growth  <cit> .
in order to grow in such high salinity environments, they have developed a set of physiological adaptation, such as: high internal concentration of potassium chloride, acidic proteome with low pi value, high gc content with gc bias in the wobble position, unique chloride pumps to maintain osmotic balance, etc.
<cit> .
among archaea, halobacteria also have the unique ability to use solar energy to generate a proton gradient to synthesize atp.
so far, the class halobacteria harbors one family with  <dig> genera and  <dig> species have been completely sequenced, including halobacterium sp.
nrc- <dig>  haloarcula marismortui, h. walsbyi and natronomonas pharaonis  <cit> .
by performing blast searches on each protein in the halobacterium sp.
nrc- <dig> genome, we have identified  <dig> proteins, which are only present in all  <dig> halobacteria species with only  <dig> exceptions .
of the proteins listed in this table, vng0016h, vng1096h, vng2414h and vng2563h are annotated as dna-binding proteins or regulators because of the presence of hth domain, but their exact functions have not been reported.
vng0667g is an atp-binding protein of abc transporter family.
several other proteins, such as vng2089h and vng2628h, have also been assigned possible functions based on weak similarity to known conserved domains in the cdd database  <cit> , but their exact functions remain to be determined.
because of their high degree of conservation and uniqueness to halobacteria, the genes for these proteins likely evolved in a common ancestor of halobacteria  and they are presumably involved in unique physiological functions related to their adaptation to the hypersaline environment.
because of their specificity for halobacteria, these proteins provide useful biomarkers for this group of species.
in addition to these proteins that are specific to all sequenced halobacterial species, we have also identified a large number of proteins either uniquely shared by  <dig> halobacterial species or only found in  <dig> halobacterial species .
surprisingly, these proteins are present in different combinations of halobacterial species.
the four-halobacterial species are from  <dig> different genera within the halobacteriales order and their relationships are unclear at present.
the largest numbers of these proteins  are uniquely shared by the haloarcula, haloquadratum and natronomonas species, followed by  <dig> proteins that are restricted to haloquadratum and haloarcula.
these results suggest that of these three species, haloquadratum and haloarcula are more closely related to each other and that halobacterium might be the deepest branching of the four available halobacterial species .
however, the genome size of these halobacterial species varies and some of these protein sequences are present on plasmids found in these species, which makes it difficult to reliably infer their relationships solely based on the number of shared proteins.
among the proteins that are specific for halobacteria, only few have been assigned possible functions.
protein vng2178h is annotated as phih1-like repressor and vng0584h is assigned as a rieske fe-s protein.
two additional proteins vng1720h and vng2562h have been annotated as iron-binding proteins because of their similarity with fhud and troa_a domains, respectively  <cit> .
all of the other proteins are of unknown function.
proteins that are specific for thermoplasmata
the thermoplasmata group is comprised of cell wall-less archaea, which resemble the bacterial mycoplasma species  <cit> .
generally, they are thermoacidophilic, aerobic or facultative anaerobic, and are able to reduce sulfur to h2s under anaerobic conditions  <cit> .
to date, this class include three families-thermoplasmaceae, picrophilaceae, and ferroplasmaceae, each represented by one genus  <cit> .
three complete genomes from this class  are available at present   <cit>  and ferroplasma acidarmanus fer <dig> genome is draft assembled and sequence information for this is also available in the ncbi database.
our analyses have uncovered  <dig> proteins that are uniquely present in all four species belonging to this class ) .
most of these proteins are present in all four available genomes, but a few are missing in one or two species, which is probably due to gene loss.
besides, we have also identified  <dig> proteins, which are shared only by the two thermoplasma species ) and  <dig> proteins unique to p. torridus and f. acidarmanus ).
the latter proteins indicate that species from picrophilaceae and ferroplasmaceae families are more closely related to each other .
all of these proteins are of unknown or predicted functions.
proteins restricted to several archaeal lineages or showing sporadic distribution
in addition to the above proteins that are restricted to specific lineages of archaea, we have also identified  <dig> proteins, which are shared by several archaeal groups .
the distribution pattern of these proteins could provide useful insights concerning the phylogenetic relationship between different groups.
however, their distribution patterns could also be explained by gene losses in specific lineages or lgt between particular groups.
table  <dig> shows many proteins that are uniquely shared by various methanogenic archaea, archaeoglobus and thermococci.
the first  <dig> proteins in table  <dig>  are present in all of the thermococci and most of the methanogens.
four of these proteins are also present in a. fulgidus.
the next  <dig> proteins in this table are also uniquely found in most of the thermococci as well as a number of methanogens and also in many cases in a. fulgidus.
in addition,  <dig> proteins listed in table  <dig> are only found in various thermococci and a. fulgidus.
these results suggest a closer relationship between the methanogenic archaea, a. fulgidus and thermococci within the euryarchaeota lineage.
in conjunction with our earlier inference that a. fulgidus forms an outgroup of the methanogenic archaea, these results suggest that the above three groups are related in the following manner: thermococci → a. fulgidus → methanogens.
although the relationship suggested above is the most likely explanation for the observed results, we have also come across three proteins  that are uniquely present in various halobacteria, a. fulgidus and different methanogens.
to account for their species distribution, one has to postulate that their genes have been selectively lost from the thermococci.
in addition,  <dig> proteins are only found in various halobacteria as well as methanosarcinales and methanomicrobiales ).
their distribution requires again either selective gene losses from other lineages or lgt from halobacteria to these methanogens.
our analyses have also uncovered  <dig> proteins that are uniquely shared by species from thermoplasmata and sulfolobus ).
among these proteins,  <dig> are present in all thermoplasmata and sulfolobus species for which sequence information is available, while the remainder are missing in  <dig> or more species.
it has been reported that there has been much lateral gene transfer between t. acidophilum and s. solfataricus, both of which inhabit the same environment  <cit> .
however, the shared presence of these proteins in these two groups could also result from a unique shared ancestry of these thermo-acidophilic archaea.
another  <dig> archaea-specific proteins are sporadically present in different archaeal species .
a number of proteins in this group are present in a limited number  of archaeal species belonging to different groups.
there are  <dig> possible explanations that can account for their sporadic distribution: first, it is possible that some of these genes are the remnants of sequences that also originated in an ancestral lineage of archaea but they have been selectively lost in many species because they are not required for growth.
second, the sporadic presence of these genes in a number of archaeal species can also be explained if some of these genes originally evolved in a particular group or species of archaea and then transferred to other archaea by lgt  <cit> .
however, in view of the observed specificity of these genes/proteins for archaea, the lgts in these cases need to be selective and limited to within archaea.
CONCLUSIONS
comparative analyses of sequenced archaeal genomes presented here have led to identification of large numbers of proteins that are distinctive characteristics of either all archaea or its different main groups.
based upon these proteins, all of the main groups within archaea  and their subgroups can now be clearly distinguished in molecular terms.
the species distribution of these signature proteins strongly suggests that their genes have evolved or originated at various stages in the evolution of archaea, but once evolved, they are indicated to be generally stably retained in various descendents of these lineages with minimal gene loss or lgts.
based upon the species distributions of these proteins, the evolutionary stages where the genes for these proteins have likely evolved are shown in fig.
<dig>  the evolutionary relationships among archaea have thus far been mainly inferred on the basis of their branching in phylogenetic trees based on 16s rrna and certain protein sequences  <cit> .
the results of our analyses although they support many inferences reached based on phylogenetic trees    <cit> , they also differ from them in important regards.
in particular, our results shed important light on certain phylogenetic relationships that were very puzzling or were not resolved based on earlier studies.
some of these novel inferences are discussed below.
in phylogenetic trees based on 16s rrna and various proteins sequences, the methanogenic archaea form at least two distinct clusters   <cit> .
in addition, in many of these trees, m. kandleri branches distinctly from all other methanogenic archaea  <cit> .
the methanogenic archaea in these trees are interspersed by other groups of non-methanogenic archaea such as halobacteriales, archaeoglobus, thermoplasmatales and thermococcales   <cit> .
this has led to important questions concerning the origin of methanogenesis i.e.
whether it evolved only once and its absence in the intervening lineages  <cit> .
to account for these results, it has been suggested that methanogenesis evolved once in a common ancestor of the above groups, i.e.
different methanogenic archaea, halobacteriales, archaeoglobus, thermoplasmatales and also possibly thermococcales, comprising virtually all euryarchaeota, but that the various genes involved in this process were subsequently lost from different groups except the methanogens  <cit> .
this scenario, in essence, proposes that the common ancestor of different physiologically and metabolically distinct groups within euryarchaeota was a methanogen and this capability was independently lost in all other lineages.
in contrast to this proposal, our phylogenomics analyses have identified  <dig> proteins that are uniquely present in virtually all methanogens, as well as many proteins that are specifically shared by different subgroups of methanogens.
of these proteins only about 1/ <dig> are indicated to be directly involved in methanogenesis and the cellular functions of others are presently not known.
the unique presence of such large numbers of proteins by nearly all methanogens, but none of the above groups of archaea, strongly indicates that the genes for these proteins evolved in a common ancestor of various methanogens.
these results strongly suggest that all methanogenic archaea form a mononphyletic lineage exclusive of all other groups of archaea .
importantly, these studies have also identified  <dig> proteins that are uniquely shared by all methanogens as well as by a. fulgidus.
in contrast, we have not come across any protein that various methanogenic archaea uniquely share with any of the halobacterales or thermoplasmatales.
these observations are highly significant because they strongly suggest that archaeoglobus and all of the methanogens shared a common ancestor exclusive of all other archaea.
in other words, the ancestral lineage that led to the origin of methanogenesis very likely evolved from the archaeoglobus lineage .
it is also significant that of the proteins that are uniquely shared by archaeoglobus and methanogens, several form part of complexes that are important for nitrogen assimilation and methanogenesis.
these results support the view that these characteristics have their origin within the archaeoglobus lineage.
the present work also provides clarification regarding the phylogenetic position of m. kandleri.
in phylogenetic trees based on 16s rrna or different protein sequences, the branching of this species is highly variable  <cit>  and it often forms the deepest branch within the euryarchaeota.
in the present work, we have identified  <dig> proteins that are uniquely shared by all methanogens including m. kandleri, as well as  <dig> proteins that m. kandleri specifically shares with various methanobacteriales and methanococcales, and  <dig> additional proteins that are only found in m. kandleri and the two methanobacteriales species .
these observations reliably place m. kandleri with other methanogenic archaea with the methanobacteriales as its closest relatives .
our results also suggest a closer relationship of the thermococcales to the archaeoglobus and methanogenic archaea, although this relationship is not as strongly supported as between archaeoglobus and methanogens.
the observed differences in the evolutionary relationships among methanogens based upon phylogenomics analyses versus those by traditional phylogenetic methods can in principle be accounted for by three explanations.
first, it is possible that the branching patterns of various clades in phylogenetic trees are misleading and they have been affected by factors such as long branch attraction effect  <cit> .
second, the polyphyletic branching of methanogens can also be explained  if the genes uniquely shared by all methanogens evolved in an early branching lineage such as m. kandleri, but subsequently they were either completely or partially lost from various non-methanogenic  groups that lie in between the two methanogenic clusters .
third, lateral transfer of these genes from one methanogenic archaea to all others can also explain these results.
of these possibilities, we favour the first explanation, as the last two require extensive gene loss or lgt from  multiple independent lineages.
the present work also supports the placement of n. equitans within the euryarchaeota lineage.
n. equitans has a very small genome , which is at least  <dig> times smaller than any other archaeal genome.
due to its very small size, there are only  <dig> genes that n. equitans uniquely shares with all other archaea.
however, our analysis indicates that whereas n. equitans shares a few genes  with most of the euryarchaeota, it does not share any gene uniquely with most of the crenarchaeota species, indicating its closer affinity for the former lineage.
although our analysis of the n. equitans genome has not revealed any strong signals indicating its specific affinity for any of the euryarchaeota groups, the shared presence of some proteins by n. equitans and thermococci  suggest that it may be related to the thermococci.
however, because of the extensive gene losses that have occurred in this genome, we are not able to draw any reliable inference in this regard.
therefore, although we have depicted n. equitans as a deep branching lineage within euryarchaeota , based upon our analysis, its placement within euryarchaeota is not resolved.
the present work also suggests that thermoplasmatales might be a deeper branching lineage within euryarchaeota in comparison to the thermococcales, halobacteriales, archaoglobous and methanogens.
this inference is suggested by the observation that a number of proteins that are uniquely present in almost all other euryarcheota species are missing in the thermoplasmatales.
although the absence of these proteins in the thermoplasmatales can be explained by specific gene loss, the possibility that the genes for at least some of these proteins have evolved after the branching of thermoplasmatales deserves serious consideration.
the deeper branching of the thermoplasmatales within the euryarchaeota will place it closer to the crenarchaeota.
such a placement could prove helpful in understanding why so many genes  are uniquely shared by various thermoplasmatales and the sulfolobales.
for the archaeal-specific proteins identified in the present work, sequence information at present is available from only a limited number of archaeal species.
hence, it is important to obtain information for these genes/proteins from other archaeal species to confirm whether these proteins are distinctive characteristics of the specified groups or a subgroup of such species.
these proteins in addition to their utility for phylogenetic and taxonomic studies also provide valuable means for understanding archaeal biology  <cit> .
the cellular functions of most of these proteins are not known and further studies in this regard should prove very helpful in the discovery of novel biochemical and physiological characteristics that are unique to either all or different groups of archaea  <cit> .
lastly, the primary sequences of many of these genes/proteins are also highly conserved and they provide novel means for identification of different groups of archaea in various environmental settings by means of pcr amplification and other molecular biological and immunological methods.

 archaea are widely regarded as one of the three main domains of life, although their origin is a subject of debate. archaeal species were earlier believed to inhabit only extreme environments such as extremely hot, or hot and acidic, extremely saline, or very acidic or alkaline conditions. however, recent studies provide evidence that they are widespread in different environments. the archaea also include methanogens, which grow under strictly anaerobic and often thermophilic conditions, and are the only organisms that derive all of their metabolic energy by reduction of co by hydrogen to produce methane. the archaeal species branch distinctly from all other organisms in phylogenetic trees based on 16s rrna and many other gene/protein sequences. in addition, many morphological or physiological characteristics such as the presence of branched-chain ether-linked lipids in their cell membrane, lack of peptidoglycan in their cell wall, characteristic subunit pattern of rna polymerase, presence of modified bases in trna, presence of a unique form of dna polymerase, have been previously indicated as defining characteristics of archaea. however, as noted by walsh and doolittle, many of these features are either not shared by all archaea or they are also present in various eukaryotes or some thermophilic bacteria, indicating that they do not constitute distinctive characteristics of all archaea. the phylogenetic analyses of archaea have led to their division into two major groups or phyla designated as crenarchaeota and euryarchaeota. the crenarchaeota species have also been referred to as'eocytes' by lake and coworkers. the species from both these groups, particularly euryarchaeota, are highly diverse in terms of their metabolism and physiology. based on their metabolic and physiological characteristics and other unique features, five functionally distinct groups within euryarchaeota are currently recognized: methanogens, sulfate reducers, extreme halophiles, cell wall-less archaea, and extremely thermophilic sulfur metabolizing archaea. some of these groups, such as methanogens, are polyphyletic in different phylogenetic trees. however, the sets of genes or proteins that are unique to these different functional groups and distinguish them from all others remain to be identified. in the past years, complete genomes of many archaeal species covering all major divisions within the archaea have been sequenced. comparative analyses of these sequences provide a valuable resource for identifying different genes/proteins that are distinctive characteristics of various taxonomic and functional groups within archaea. whole proteins that are uniquely present in particular groups or subgroups of organisms but not found anywhere else provide valuable molecular markers for taxonomic, phylogenetic and biochemical studies. these proteins, which we refer to as signature proteins in our work, and others have called them as orfans or conserved hypothetical proteins, are present at different phylogenetic depths, such as genus, family, order or even phylum. in our recent work, a large number of such proteins that are distinctive characteristics of several groups within bacteria, and also their subgroups, were identified. these proteins provide not only valuable molecular markers for identifying and circumscribing species belonging to these major groups in molecular terms, but their species distribution pattern also provides useful information about the branching order within these groups. as archaea constitute a very diverse group, identification of sets of proteins that are specific for its main groups and subgroups should prove useful in terms of identifying molecular characteristics that are unique to them. additionally, this information should also be helpful in understanding the evolutionary relationships among different groups. comparative studies on limited numbers of archaeal genomes have been carried out by a number of investigators using different criteria. graham et al. analyzed archaeal genomes to identify signature proteins that function uniquely within the archaea. their definition of an archaeal signature protein required it to be present in only two different euryarchaeal species and they identified archaeal signature proteins. makarova and koonin have analyzed archaeal genomes to identify core sets of genes, which are present in all archaeal species, but which are not restricted to the archaeal species. recently, walsh and doolittle have analyzed prokaryotic genomes to measure dissimilarity between archaea and bacteria. although it was reported that 28% of the proteins from archaeal genomes are restricted to the archaea, specific proteins that were present in different groups of archaea were not identified. other comparative studies using different criteria have been conducted on smaller groups within archaea such as pyrococcus, sulfolobus and thermoacidophilic organisms. however, thus far no comprehensive phylogenomics study on different archaeal genomes has been carried out using the same standard criteria to identify proteins or orfs that are shared by all archaea or its different major lineages. in this study we have carried out comparative analyses of archaeal genomes using uniform criteria to identify proteins that are uniquely present in archaeal species at different phylogenetic depths representing all major groups within the archaea.  a. phylogenetic analyses of archaeal species prior to undertaking comparative studies on archaeal genomes, phylogenetic analysis of sequenced archaeal species was carried out so that the results of phylogenomics analyses could be compared with those obtained by traditional phylogenetic approaches. phylogenetic trees for the archaeal species based on 16s rrna as well as concatenated sequences of translation and transcription-related proteins have been published by other investigators. in the present work, we have constructed phylogenetic trees for archaeal species using a set of universally distributed proteins that are involved in a broad range of functions. the sequence of haloquadratum walsbyi dsm which became available afterward, was not included in these studies. phylogenetic trees based on a concatenated sequence alignment of these proteins were constructed using the neighbour-joining, maximum-likelihood and maximum-parsimony methods. the results of these analyses are presented in fig. all three methods gave very similar tree topologies except for the branching positions of m. kandleri and methanospirillum hungatei, which were found to be variable. except for this, the branching pattern of the archaeal species based on our dataset is very similar to that reported by gribaldo et al. based on concatenated sequences of translation and transcription-related proteins. in the tree shown, the crenarchaeota and euryarchaeota, the two major phyla within archaea were clearly distinguished from each other. the phylogenetic affinity of nanoarchaeum, which has a long-branch length, was not resolved in this or various other trees. within crenarchaeota, pyrobaculum was indicated to be a deeper branch, and aeropyrum branched in between the pyrobaculum and sulfolobus. within euryarchaeota, the clades corresponding to halobacteria, thermococci and thermoplasmata were resolved with high bootstrap scores, but the methanogens were split into 2– clusters. one of these clusters that has low bootstrap score consisted of methanobacteriales and methanococcales with m. kandleri branching in its vicinity. the second cluster, with higher bootstrap score, showed a grouping of methanomicrobiales and methanosarcinales. these two clusters, which are separated by thermoplasmata, archaeoglobi and halobacteria, have been referred to as class i and class ii methanogens by bapteste et al.. b. phylogenomic analyses of archaeal genomes to search for proteins, which are uniquely present in either all archaea or various subgroups of them, blast searches were performed on each open reading frame from a total of archaeal genomes. these genomes included crenarchaeota and divergent euryarchaeota species covering all main functional and phylogenetic groups. the euryarchaeota genomes analyzed included: pyrococcus abyssi from extremely thermophilic sulfur metabolizing archaea, methanococcus maripaludis from methanococcales, halobacterium sp. nrc- and h. walsbyi from extreme halophiles, thermoplasma acidophilum and picrophilus torridus belonging to the cell wall-less archaea, methanococcoides burtonii from methanosarcinales and methanopyrus kandleri from the methanopyrales order. the chosen genomes should provide information regarding all archaeal proteins that are shared at a taxonomic level higher than a genus. the analysis of the remainder of the genomes, which was expected to provide information regarding proteins that are only unique to a given species, was not carried out. each orf from these genomes was examined by means of blastp and psi-blast searches against all available sequences from different organisms to identify proteins that are specific for only archaeal lineages. the methods and the criteria that we have used to identify proteins that are specific for either all or various subgroups of archaea are described in the methods section. generally, a protein was considered to be specific for a given archaeal lineage if all significant hits or alignments in the blastp and psi-blast searches with the query protein were from the indicated group of archaeal species. in a few cases, where 1– isolated species from other groups also exhibited significant similarity, such proteins were retained as they provide interesting examples of lateral gene transfer from archaea to other groups. our analyses have identified proteins that are unique to different groups of archaea and for which no homologues are generally found in any bacterial or eukaryotic species. based on their specificity for different taxonomic groups, these proteins have been divided into a number of different groups. a brief description of the different subsets of archaeal-specific proteins and functional information regarding them, where known, is given below. in the description of these proteins that follows, the'ape','hq','mbu','mk','mmp','pab','pae','pto','saci','ta','vng', and'neq' part of the descriptors in proteins indicate that the original query protein sequence was from the genome of a. pernix k h. walsbyi dsm m. burtonii dsm m. kandleri av m. maripaludis s p. abyssi ge p. aerophilum str. im p. torridus dsm s. acidocaldarius dsm t. acidophilum dsm halobacterium sp. nrc- and n. equitans, respectively. proteins that are specific for all archaea of the proteins that are uniquely present in all archaea, pab corresponds to trna nucleotidyltransferase, which builds and repairs the 3' end of trna. functionally similar enzymes are also present in bacteria and eukaryotes, but their sequences share very little homology with the archaeal cca-adding enzyme, which explains why no homologs were detected in any bacteria or eukaryotes in blast searches. the main mechanistic difference between class i and class ii enzymes is that the trna substrate is required to fully define the nucleotide binding site in class i enzyme, whereas class ii has a preformed nucleotide binding site that recognizes ctp and atp in the absence of trna. another protein pab is assigned as archaeal type dna primase, which also has its synonymous counterparts in bacterial and eukaryotic species, but shows very little homology to them. in the same way, protein pab is annotated as a pilt family atpase, which showed very little similarity to bacterial atpases involved in type iv pili biogenesis. further studies of this protein could provide insights into novel aspects of the archaeal flagellar system. a number of other proteins viz. pab pab0018a, pab pab and pab have also been assigned putative functions based on sequence analysis, but their exact roles in archaeal cells remains to be determined. interestingly, for protein pab two gene copies with acceptable e-values are also present in the genomes of dehalococcoides ethenogenes dehalococcoides sp. cbdb and dehalococcoides sp. bav which belong to chloroflexi. because no homologue of pab is present in other bacteria, it is likely that this protein was transferred from archaea to the common ancestor of dehalococcoides followed by a gene duplication event. proteins that are specific for crenarchaeota as mentioned in the introduction, the archaea are divided into main groups, crenarchaeota and euryarchaeota, based on 16s rrna trees as well many other gene trees and characteristics. the crenarchaeota are also indicated to differ from euryarchaeota in terms of their ribosome structure. in comparison to euryarchaeota, which contain physiologically and metabolically diverse groups of organisms, the crenarchaeota were thought to be a pure collection of extreme thermophiles and most members metabolize sulfur. however, recent studies indicate that crenarchaeota are much more diverse in their physiology and ecology than was previously believed. many species living in the cold ocean also belong to this group based on their branching pattern in 16s rrna trees, although most of them have not been cultivated. currently, this phylum is comprised of one single class thermoprotei containing three orders: thermoproteales, desulfurococcales and sulfolobales. fortunately, every order has a completely sequenced representative, which provide a platform to explore the characteristics that are unique to crenarchaeal species. comparative genomic surveys have revealed some molecular features that are shared by crenarchaea but not euryarchaea, such as the lack of histones, absence of the ftsz-mincde system and distinctive rrna operon organization. lake et al. have also identified distinctive differences in ribosome structure and an insert in elongation factor ef-g and ef-tu, which can be used to distinguish crenarchaeota from euryarchaeota. however, these features are not unique characteristics of the crenarchaeota. blast searches on each orf from the genomes of a. pernix and s. acidocaldarius dsm have identified proteins which are shared by all five crenarchaeal species, but whose homologs are not found in other archaea, or any bacteria or eukaryotes with only exceptions ). the genes for these proteins likely evolved in a common ancestor of the crenarchaeota and they provide potential molecular markers for species from this phylum. additionally, proteins that are listed in table are only found in a. pernix and three sulfolobus genomes. these proteins suggest that aeropyrum and sulfolobus may have shared a common ancestor exclusive of pyrobaculum. however, we have also come across proteins that are shared by aeropyrum and pyrobaculum ) and proteins that are exclusively present in the sulfolobus species and pyrobaculum ). hence, based upon the species distributions of these proteins, the relationships among the aeropyrum, sulfolobales and pyrobaculum are not entirely clear. in phylogenetic trees thermoproteales branches consistently earlier than desulfurococcales and sulfolobales. this observation in conjunction with the fact that aeropyrum and sulfolobus share larger numbers of proteins in common with each other suggests that these two groups likely shared a common ancestor exclusive of pyrobaculum. the proteins that are only found in aeropyrum and pyrobaculum, or in sulfolobus and pyrobaculum, most likely evolved in a common ancestor of the crenarchaea, but were subsequently lost in either the sulfolobales or a. pernix lineages. in addition to these proteins that are uniquely present in either all sequenced crenarchaeota genomes or different groups of crenarchaeota species, these analyses have also identified proteins that are unique for the sulfolobales species. of these, proteins are present in all sequenced sulfolobus genomes, whereas the remaining are present in at least two of the three sulfolobus genomes. in this work, since blast analyses were not carried out on all three sulfolobus genomes, it is likely that the numbers of genes or proteins that are uniquely shared by only two sulfolobus genomes is much higher than indicated here. chen et al. have previously analyzed the genome of s. acidocaldarius dsm and indicated the presence of genes that were specific for crenarchaeota and genes that were specific to sulfolobus genus. however, in the present work, relatively few genes that are uniquely shared by various crenarchaeota species were identified. this difference could be due to more stringent criteria that we have employed for identification of proteins that are specific to different groups. the genome of thermofilum pendens hrk which belongs to thermoproteales, has also been partially sequenced and information for large numbers of genes/proteins from this species is available in the ncbi database. by carrying out blast searches on each orf from p. aerophilum genome, we have identified proteins that are only found in the above thermoproteales species. the numbers of proteins shared by these two species will likely increase once complete genome of t. pendens becomes available. many of these proteins are expected to provide markers for the thermoproteales order. proteins that are specific for euryarchaeota the euryarchaeota, which comprise a majority of the cultured and sequenced archaea, is a morphologically, metabolically and physiologically diverse collection of species as evidenced by the presence in this group of various methanogens, extreme halophiles, cell wall-less archaea and sulfate reducing microbes. no unique biochemical or molecular characteristic that is commonly shared by all of the different lineages is known. the present study has identified proteins that are only found in euryarchaeota species with exceptions. in this table, the first proteins ) are present in most euryarchaeota species. of these proteins, pab and pab were found in all sequenced euryarchaeota species. pab was also present in n. equitans, supporting its placement within the euryarchaeota. the protein pab is annotated as archaeosine trna-ribosyltransferase, which catalyzes the exchange of guanine with a free 7-cyano-7-deazaguanine base, as the first step in the biosynthesis of an archaea-specific modified base, archaeosine. it should be mentioned that there is another protein pab in the same genome, which is also annotated and experimentally confirmed as arctgt. the latter belongs to a family of proteins that are highly conserved in all archaea species and some bacteria. it seems that pab might be involved in rna modification since it possesses a pua domain, but its function is likely different from pab the protein pab which is annotated as dna polymerase ii large subunit, is highly conserved within euryarchaeota, but is not found anywhere else except in nanoarchaeum. this enzyme is the major dna replicase in euryarchaeota and also a distinctive molecular marker for this group. the genes for the above proteins likely evolved in a common ancestor of euryarchaeota and they provide molecular markers for this diverse group of organisms. another proteins listed in table are found in almost all euryarchaeota, but they are missing in thermoplasmata. their distribution suggests that either thermoplasmata is a deep branching lineage within euryarchaeota or that the genes for these proteins have been selectively lost from thermoplasmata. of these proteins, pab is also present in n. equitans supporting its placement with euryarchaeota. five other proteins from the first two columns in table are also not found in the thermococci species. these results can again be explained by either selective loss of these genes from these particular groups or deeper branching of these lineages within the euryarchaeota species. on the basis of proteins listed in table although one can infer that thermoplasmata and thermococci are deeper branching lineages within euryarchaeota in comparison to methanogens, their relative branching order cannot be resolved. proteins that are specific for different main groups within euryarchaeota proteins specific for methanogenic archaea and their various subgroups currently, the methanogens form the largest group within the euryarchaeota. they are distinguished from all other prokaryotes by their ability to obtain all or most of their energy via the reduction of co to methane or by the process of methanogenesis. in the bergey's manual, the methanogenes are divided into distinct orders. some studies have suggested that these organisms possess a set of unique enzymes which are responsible for methanogenesis, such as coenzyme m, factor and methanopterin. however, no systematic study has been carried out thus far to identify proteins that are uniquely present in different methanogens. our blast searches of proteins from different methanogens have led to identification of proteins, which are uniquely found in various methanogenic archaea. twenty of these proteins are present in all sequenced methanogens, while proteins are missing only in m. stadtmanae, which is a human intestinal inhabitant. this archaeon generate methane by reduction of methanol with h and lacks many proteins present in the genomes of other methanogens. thus, it is highly likely that the proteins missing in m. stadtmanae were selectively lost from this species. therefore, it is very likely that the genes for these proteins that are commonly shared by virtually all methanogens ) evolved in a common ancestor of all methanogens. these analyses have also identified proteins that are uniquely shared by various methanogens as well as a. fulgidus ). the genes for these proteins likely evolved in a common ancestor of a. fulgidus and various methanogenic archaea and they point to a close relationship between these two groups of organisms. ten additional proteins are present in a. fulgidus as well as various methanosarcinales and m. hungatei ). it is likely that the genes for these proteins also evolved in a common ancestor of a. fulgidus and various methanogenic archaea, but they were selectively lost in other methanogens. of the proteins that are commonly shared by a. fulgidus and various methanogenic archaea, mmp is reported to be a novel repressor of nif and glna genes, which are involved in nitrogen assimilation. interestingly, homologs of this protein are also found in dehalococcoides species, but nowhere else, which are very likely due to lgt. protein mmp is the ε-subunit of carbon-monoxide dehydrogenase complex, which is made up of five subunits in different methanogens. the epsilon subunits are required for the reversible oxidation of co to co. all of the other components could be found in a few bacterial species, while the ε-subunit is restricted to methanogenic archaea and a. fulgidus. protein mmp is identified as a transcriptional regulator with a helix-turn-helix motif, but its exact role has not been reported. among the genes that are uniquely shared by various methanogenic archaea, two large gene clusters responsible for methanogenesis are found. the proteins mmp mmp1560–mmp and mmp1566–mmp are parts of an eight-component complex, coenzyme m methyltransferase, which catalyzes an energy-conserving, sodium-ion-translocating step in methanogenesis from h and co. m. maripaludis contains all of the known mtr subunits, but the gene coding for mtrf is fused into the n-terminal region of mtra. all other methanogenic archaeal genomes contain complete set of mtr genes. it is of interest to note that for the protein mmp, homologues with low e-values are also found in two desulfitobacterium hafniense strains as well as in three rhizobiales species. these three rhizobiae species can use methyl halides as a sole source of carbon and energy, and all of them possess a set of cmu genes which are essential for methyl chloride degradation. in particular, the cmub protein which is homologous to mmp transfers a methyl group to methylcobalamin:h folate, which is analogous to the reverse of the reaction catalyzed by mtrh in archaea. in view of the sequence and functional similarity between mtrh and cmub proteins, it is likely that the mtrh gene was laterally transferred from a methanogenic archaeon to the common ancestor of the above three rhizobiae species to serve the new functional role. the function of the laterally transferred mtrh related gene in d. hafniense is not known at present. the proteins mmp1555–mmp in table form another gene cluster, encoding the subunits of methyl-coenzyme m reductase. this complex catalyzes the final reaction of the energy conserving pathway in which methylcoenzyme m and coenzyme b are converted to methane and the heterodisulfide com-s-s-cob. except for these proteins, the other proteins listed in table are of putative or unknown functions. it is likely that these proteins are involved in some aspects of methanogenesis or other unknown pathways unique to methanogenic archaea. these proteins provide molecular markers for methanogens, which can be used for identification of new archaeal species capable of methane production. the blast searches of the m. maripaludis and m. kandleri genomes have identified proteins that are uniquely shared by all of the following species belonging to the orders methanobacteriales, methanococcales and methanopyrales ). of these, only proteins are present in m. stadtmanae, which is also a methanobacteriales that has lost most of its genes due to its adaptation to the human intestine. the genes for these proteins likely evolved in a common ancestor of the above groups of methanogens, which corresponds to the cluster of methanogenic archaea referred to as "class i methanogens". interestingly, these studies have also identified proteins that are uniquely shared by these methanogenic orders and m. hungatei ), which branches distantly in phylogenetic trees. the unique presence of these proteins in these methanogens suggests that species from these groups shared a common ancestor exclusive of other methanogenic archaea. fifteen additional proteins discovered in this work ) are uniquely present in m. kandleri and various methanobacteriales indicating that these two groups are more closely related to each other than the methanococcales. we have also come across proteins that are uniquely shared by methanococcales and methanobacteriales ), and proteins that are only present in methanococcales and methanopyrales ). the most likely explanation to account for the species distributions of these latter proteins is that their genes also originated in a common ancestor of the above three groups of methanogens, but were selectively lost in either the methanobacteriales or methanopyrales lineages. these analyses have also identified additional proteins that are uniquely present in all methanosarcinales species ), as well as proteins that are only found in various methanosarcinales and m. hungatei ). lastly, these studies have also identified proteins that are uniquely present in m. maripaludis and m. jannaschii ) and proteins that are only present in m. burtonii and methanosarcina species, all belonging to the methanosarcinaceae family ) indicating that they are likely distinctive characteristics of species from these groups. of the proteins that are uniquely found in methanococcales, methanobacteriales, methanopyrales and methanomicrobiales, proteins viz. mmp1448–mmp mmp mmp1458–mmp and mmp are from a big gene cluster eha, which encodes the multisubunit membrane-bound hydrogenase. two of these proteins, mmp and mmp are only found in methanococcales ). the whole eha operon is composed of orfs in the genome of m. thermoautotrophicus and of these only these proteins are restricted to these methanogens while the other subunits have counterparts in bacteria. the precise roles of these proteins, which are predicted to be integral membrane proteins in the hydrogenase complex, have not been determined. among the other proteins that are specific for these groups of methanogens, mmp and mmp are hmd homologs, which catalyze the reversible dehydrogenation of n n10-methylenetetrahydromethanopterin. in the proteins that are specific for the methanococcales ), one large gene cluster is found, but no information is available concerning its possible function. except for these proteins, all other proteins that are specific for these methanogenic archaea are of unknown or putative function. proteins that are specific for thermococci thermococci are obligately thermophilic, strictly anaerobic cocci, which are able to convert elemental sulfur to hydrogen sulfide. thus, they are so called "extremely thermophilic sulfur metabolizer", which comprise one of the main functional groups within euryarchaeota. according to the bergey's manual, the class thermococci contains only one family, thermococcaceae, consisting of genera: thermococcus and pyrococcus. currently, species from this family have been completely sequenced. the blast searches on each protein from p. abyssi have identified proteins that are shared by all of these species ). all of these proteins show high degree of conservation within thermococci and they do not have homologs in any other prokaryotes or eukaryotes except one possible lgt event. the genes for these proteins have likely evolved in a common ancestor of the thermococci. of these proteins, pab is annotated as tbp-interacting protein, which forms complex with tbp to regulate transcription. it is known that the archaeal transcription machinery is strikingly similar to that in eukaryotes, but no tbp-binding component was found in archaeal species until the discovery of the tip in t. kodakaraensis. most other themococci-specific proteins are of unknown function, although in a few cases limited similarity to domains in known protein families have been noted. a number of proteins are clustered together in the p. abyssi genome, and it is possible that they may form functional units and are involved in related functions. cohen et al. have reported a large number of proteins which are restricted to the pyrococcus genus. however, a number of proteins from their list are also found in t. kodakarensis kod, whose genome was not available when their work was published. some proteins are not specific for either pyrococcus or thermococci according to our criteria and some of them are only found in one species – p. abyssi. our analysis of the p. abyssi ge genome has also identified proteins that are unique to the pyrococcus genus ). again, almost all of these proteins are of unknown function except pab which is annotated as rnase p, but this annotation seems arbitrary as it does not show significant sequence similarity to known rnases. the proteins that are uniquely found in the pyrococcus genomes likely evolved in a common ancestor of this genus. proteins that are specific for halobacteria extreme halophiles constitute another major class within euryarchaeota. they require 5– times the salinity of seawater for optimal growth. in order to grow in such high salinity environments, they have developed a set of physiological adaptation, such as: high internal concentration of potassium chloride, acidic proteome with low pi value, high gc content with gc bias in the wobble position, unique chloride pumps to maintain osmotic balance, etc.. among archaea, halobacteria also have the unique ability to use solar energy to generate a proton gradient to synthesize atp. so far, the class halobacteria harbors one family with genera and species have been completely sequenced, including halobacterium sp. nrc- haloarcula marismortui, h. walsbyi and natronomonas pharaonis. by performing blast searches on each protein in the halobacterium sp. nrc- genome, we have identified proteins, which are only present in all halobacteria species with only exceptions. of the proteins listed in this table, vng0016h, vng1096h, vng2414h and vng2563h are annotated as dna-binding proteins or regulators because of the presence of hth domain, but their exact functions have not been reported. vng0667g is an atp-binding protein of abc transporter family. several other proteins, such as vng2089h and vng2628h, have also been assigned possible functions based on weak similarity to known conserved domains in the cdd database, but their exact functions remain to be determined. because of their high degree of conservation and uniqueness to halobacteria, the genes for these proteins likely evolved in a common ancestor of halobacteria and they are presumably involved in unique physiological functions related to their adaptation to the hypersaline environment. because of their specificity for halobacteria, these proteins provide useful biomarkers for this group of species. in addition to these proteins that are specific to all sequenced halobacterial species, we have also identified a large number of proteins either uniquely shared by halobacterial species or only found in halobacterial species. surprisingly, these proteins are present in different combinations of halobacterial species. the four-halobacterial species are from different genera within the halobacteriales order and their relationships are unclear at present. the largest numbers of these proteins are uniquely shared by the haloarcula, haloquadratum and natronomonas species, followed by proteins that are restricted to haloquadratum and haloarcula. these results suggest that of these three species, haloquadratum and haloarcula are more closely related to each other and that halobacterium might be the deepest branching of the four available halobacterial species. however, the genome size of these halobacterial species varies and some of these protein sequences are present on plasmids found in these species, which makes it difficult to reliably infer their relationships solely based on the number of shared proteins. among the proteins that are specific for halobacteria, only few have been assigned possible functions. protein vng2178h is annotated as phih1-like repressor and vng0584h is assigned as a rieske fe-s protein. two additional proteins vng1720h and vng2562h have been annotated as iron-binding proteins because of their similarity with fhud and troa_a domains, respectively. all of the other proteins are of unknown function. proteins that are specific for thermoplasmata the thermoplasmata group is comprised of cell wall-less archaea, which resemble the bacterial mycoplasma species. generally, they are thermoacidophilic, aerobic or facultative anaerobic, and are able to reduce sulfur to h2s under anaerobic conditions. to date, this class include three families-thermoplasmaceae, picrophilaceae, and ferroplasmaceae, each represented by one genus. three complete genomes from this class are available at present and ferroplasma acidarmanus fer genome is draft assembled and sequence information for this is also available in the ncbi database. our analyses have uncovered proteins that are uniquely present in all four species belonging to this class ). most of these proteins are present in all four available genomes, but a few are missing in one or two species, which is probably due to gene loss. besides, we have also identified proteins, which are shared only by the two thermoplasma species ) and proteins unique to p. torridus and f. acidarmanus ). the latter proteins indicate that species from picrophilaceae and ferroplasmaceae families are more closely related to each other. all of these proteins are of unknown or predicted functions. proteins restricted to several archaeal lineages or showing sporadic distribution in addition to the above proteins that are restricted to specific lineages of archaea, we have also identified proteins, which are shared by several archaeal groups. the distribution pattern of these proteins could provide useful insights concerning the phylogenetic relationship between different groups. however, their distribution patterns could also be explained by gene losses in specific lineages or lgt between particular groups. table shows many proteins that are uniquely shared by various methanogenic archaea, archaeoglobus and thermococci. the first proteins in table are present in all of the thermococci and most of the methanogens. four of these proteins are also present in a. fulgidus. the next proteins in this table are also uniquely found in most of the thermococci as well as a number of methanogens and also in many cases in a. fulgidus. in addition, proteins listed in table are only found in various thermococci and a. fulgidus. these results suggest a closer relationship between the methanogenic archaea, a. fulgidus and thermococci within the euryarchaeota lineage. in conjunction with our earlier inference that a. fulgidus forms an outgroup of the methanogenic archaea, these results suggest that the above three groups are related in the following manner: thermococci → a. fulgidus → methanogens. although the relationship suggested above is the most likely explanation for the observed results, we have also come across three proteins that are uniquely present in various halobacteria, a. fulgidus and different methanogens. to account for their species distribution, one has to postulate that their genes have been selectively lost from the thermococci. in addition, proteins are only found in various halobacteria as well as methanosarcinales and methanomicrobiales ). their distribution requires again either selective gene losses from other lineages or lgt from halobacteria to these methanogens. our analyses have also uncovered proteins that are uniquely shared by species from thermoplasmata and sulfolobus ). among these proteins, are present in all thermoplasmata and sulfolobus species for which sequence information is available, while the remainder are missing in or more species. it has been reported that there has been much lateral gene transfer between t. acidophilum and s. solfataricus, both of which inhabit the same environment. however, the shared presence of these proteins in these two groups could also result from a unique shared ancestry of these thermo-acidophilic archaea. another archaea-specific proteins are sporadically present in different archaeal species. a number of proteins in this group are present in a limited number of archaeal species belonging to different groups. there are possible explanations that can account for their sporadic distribution: first, it is possible that some of these genes are the remnants of sequences that also originated in an ancestral lineage of archaea but they have been selectively lost in many species because they are not required for growth. second, the sporadic presence of these genes in a number of archaeal species can also be explained if some of these genes originally evolved in a particular group or species of archaea and then transferred to other archaea by lgt. however, in view of the observed specificity of these genes/proteins for archaea, the lgts in these cases need to be selective and limited to within archaea.  comparative analyses of sequenced archaeal genomes presented here have led to identification of large numbers of proteins that are distinctive characteristics of either all archaea or its different main groups. based upon these proteins, all of the main groups within archaea and their subgroups can now be clearly distinguished in molecular terms. the species distribution of these signature proteins strongly suggests that their genes have evolved or originated at various stages in the evolution of archaea, but once evolved, they are indicated to be generally stably retained in various descendents of these lineages with minimal gene loss or lgts. based upon the species distributions of these proteins, the evolutionary stages where the genes for these proteins have likely evolved are shown in fig. the evolutionary relationships among archaea have thus far been mainly inferred on the basis of their branching in phylogenetic trees based on 16s rrna and certain protein sequences. the results of our analyses although they support many inferences reached based on phylogenetic trees, they also differ from them in important regards. in particular, our results shed important light on certain phylogenetic relationships that were very puzzling or were not resolved based on earlier studies. some of these novel inferences are discussed below. in phylogenetic trees based on 16s rrna and various proteins sequences, the methanogenic archaea form at least two distinct clusters. in addition, in many of these trees, m. kandleri branches distinctly from all other methanogenic archaea. the methanogenic archaea in these trees are interspersed by other groups of non-methanogenic archaea such as halobacteriales, archaeoglobus, thermoplasmatales and thermococcales. this has led to important questions concerning the origin of methanogenesis i.e. whether it evolved only once and its absence in the intervening lineages. to account for these results, it has been suggested that methanogenesis evolved once in a common ancestor of the above groups, i.e. different methanogenic archaea, halobacteriales, archaeoglobus, thermoplasmatales and also possibly thermococcales, comprising virtually all euryarchaeota, but that the various genes involved in this process were subsequently lost from different groups except the methanogens. this scenario, in essence, proposes that the common ancestor of different physiologically and metabolically distinct groups within euryarchaeota was a methanogen and this capability was independently lost in all other lineages. in contrast to this proposal, our phylogenomics analyses have identified proteins that are uniquely present in virtually all methanogens, as well as many proteins that are specifically shared by different subgroups of methanogens. of these proteins only about 1/ are indicated to be directly involved in methanogenesis and the cellular functions of others are presently not known. the unique presence of such large numbers of proteins by nearly all methanogens, but none of the above groups of archaea, strongly indicates that the genes for these proteins evolved in a common ancestor of various methanogens. these results strongly suggest that all methanogenic archaea form a mononphyletic lineage exclusive of all other groups of archaea. importantly, these studies have also identified proteins that are uniquely shared by all methanogens as well as by a. fulgidus. in contrast, we have not come across any protein that various methanogenic archaea uniquely share with any of the halobacterales or thermoplasmatales. these observations are highly significant because they strongly suggest that archaeoglobus and all of the methanogens shared a common ancestor exclusive of all other archaea. in other words, the ancestral lineage that led to the origin of methanogenesis very likely evolved from the archaeoglobus lineage. it is also significant that of the proteins that are uniquely shared by archaeoglobus and methanogens, several form part of complexes that are important for nitrogen assimilation and methanogenesis. these results support the view that these characteristics have their origin within the archaeoglobus lineage. the present work also provides clarification regarding the phylogenetic position of m. kandleri. in phylogenetic trees based on 16s rrna or different protein sequences, the branching of this species is highly variable and it often forms the deepest branch within the euryarchaeota. in the present work, we have identified proteins that are uniquely shared by all methanogens including m. kandleri, as well as proteins that m. kandleri specifically shares with various methanobacteriales and methanococcales, and additional proteins that are only found in m. kandleri and the two methanobacteriales species. these observations reliably place m. kandleri with other methanogenic archaea with the methanobacteriales as its closest relatives. our results also suggest a closer relationship of the thermococcales to the archaeoglobus and methanogenic archaea, although this relationship is not as strongly supported as between archaeoglobus and methanogens. the observed differences in the evolutionary relationships among methanogens based upon phylogenomics analyses versus those by traditional phylogenetic methods can in principle be accounted for by three explanations. first, it is possible that the branching patterns of various clades in phylogenetic trees are misleading and they have been affected by factors such as long branch attraction effect. second, the polyphyletic branching of methanogens can also be explained if the genes uniquely shared by all methanogens evolved in an early branching lineage such as m. kandleri, but subsequently they were either completely or partially lost from various non-methanogenic groups that lie in between the two methanogenic clusters. third, lateral transfer of these genes from one methanogenic archaea to all others can also explain these results. of these possibilities, we favour the first explanation, as the last two require extensive gene loss or lgt from multiple independent lineages. the present work also supports the placement of n. equitans within the euryarchaeota lineage. n. equitans has a very small genome, which is at least times smaller than any other archaeal genome. due to its very small size, there are only genes that n. equitans uniquely shares with all other archaea. however, our analysis indicates that whereas n. equitans shares a few genes with most of the euryarchaeota, it does not share any gene uniquely with most of the crenarchaeota species, indicating its closer affinity for the former lineage. although our analysis of the n. equitans genome has not revealed any strong signals indicating its specific affinity for any of the euryarchaeota groups, the shared presence of some proteins by n. equitans and thermococci suggest that it may be related to the thermococci. however, because of the extensive gene losses that have occurred in this genome, we are not able to draw any reliable inference in this regard. therefore, although we have depicted n. equitans as a deep branching lineage within euryarchaeota, based upon our analysis, its placement within euryarchaeota is not resolved. the present work also suggests that thermoplasmatales might be a deeper branching lineage within euryarchaeota in comparison to the thermococcales, halobacteriales, archaoglobous and methanogens. this inference is suggested by the observation that a number of proteins that are uniquely present in almost all other euryarcheota species are missing in the thermoplasmatales. although the absence of these proteins in the thermoplasmatales can be explained by specific gene loss, the possibility that the genes for at least some of these proteins have evolved after the branching of thermoplasmatales deserves serious consideration. the deeper branching of the thermoplasmatales within the euryarchaeota will place it closer to the crenarchaeota. such a placement could prove helpful in understanding why so many genes are uniquely shared by various thermoplasmatales and the sulfolobales. for the archaeal-specific proteins identified in the present work, sequence information at present is available from only a limited number of archaeal species. hence, it is important to obtain information for these genes/proteins from other archaeal species to confirm whether these proteins are distinctive characteristics of the specified groups or a subgroup of such species. these proteins in addition to their utility for phylogenetic and taxonomic studies also provide valuable means for understanding archaeal biology. the cellular functions of most of these proteins are not known and further studies in this regard should prove very helpful in the discovery of novel biochemical and physiological characteristics that are unique to either all or different groups of archaea. lastly, the primary sequences of many of these genes/proteins are also highly conserved and they provide novel means for identification of different groups of archaea in various environmental settings by means of pcr amplification and other molecular biological and immunological methods. 
these results strongly indicate that all methanogenic archaea form a monophyletic group exclusive of other archaea and that this lineage likely evolved from archaeoglobus. six of these proteins are unique to all archaea, others are only missing in nanoarchaeum equitans and a large number of other proteins are specific for various main groups within the archaea. the species distributions of these proteins provide novel insights into the evolutionary relationships among different groups within archaea, particularly regarding the origin of methanogenesis. most of these proteins are of unknown function and further studies should lead to discovery of novel biochemical and physiological characteristics that are unique to either all archaea or its different subgroups. in contrast to the phylogenomics studies, a monophyletic grouping of archaea is not supported by phylogenetic analyses based on protein sequences. the evolutionary relationships among different groups within the euryarchaeota branch are also not clearly understood. presently, very few molecular characteristics are known that are uniquely shared by either all archaea or the different main groups within archaea. in addition, proteins that are uniquely shared by m. kandleri and methanobacteriales suggest a close evolutionary relationship between them. of particular importance is the observation that proteins are uniquely present in virtually all methanogens and additional proteins are only found in different methanogens as well as a. fulgidus. the identified archaea-specific proteins provide novel molecular markers or signature proteins that are distinctive characteristics of archaea and all of its major subgroups. we have carried out comprehensive analyses on each open reading frame in the genomes of archaea to search for proteins that are unique to either all archaea or for its main subgroups. these studies have identified proteins or orfs that are distinctive characteristics of archaea and its various subgroups and whose homologues are not found in other organisms.
 the archaea are highly diverse in terms of their physiology, metabolism and ecology. presently, very few molecular characteristics are known that are uniquely shared by either all archaea or the different main groups within archaea. the evolutionary relationships among different groups within the euryarchaeota branch are also not clearly understood.  we have carried out comprehensive analyses on each open reading frame in the genomes of archaea to search for proteins that are unique to either all archaea or for its main subgroups. these studies have identified proteins or orfs that are distinctive characteristics of archaea and its various subgroups and whose homologues are not found in other organisms. six of these proteins are unique to all archaea, others are only missing in nanoarchaeum equitans and a large number of other proteins are specific for various main groups within the archaea. of particular importance is the observation that proteins are uniquely present in virtually all methanogens and additional proteins are only found in different methanogens as well as a. fulgidus. in contrast, no protein was exclusively shared by various methanogen and any of the halobacteriales or thermoplasmatales. these results strongly indicate that all methanogenic archaea form a monophyletic group exclusive of other archaea and that this lineage likely evolved from archaeoglobus. in addition, proteins that are uniquely shared by m. kandleri and methanobacteriales suggest a close evolutionary relationship between them. in contrast to the phylogenomics studies, a monophyletic grouping of archaea is not supported by phylogenetic analyses based on protein sequences.  the identified archaea-specific proteins provide novel molecular markers or signature proteins that are distinctive characteristics of archaea and all of its major subgroups. the species distributions of these proteins provide novel insights into the evolutionary relationships among different groups within archaea, particularly regarding the origin of methanogenesis. most of these proteins are of unknown function and further studies should lead to discovery of novel biochemical and physiological characteristics that are unique to either all archaea or its different subgroups. 