Phylogenetic depth of Thermotoga maritima inferred from analysis of the fus gene: amino acid sequence of elongation factor G and organization of the Thermotoga str operon

Effect of Oxygen and Redox Potential on Glucose Fermentation in Thermotoga maritima under Controlled Physicochemical Conditions

Batch cultures of Thermotoga maritima were performed in a bioreactor equipped with instruments adapted for experiments performed at 80°C to mimic the fluctuating oxidative conditions in the hot ecosystems it inhabits. When grown anaerobically on glucose, T. maritima was shown to significantly decrease the redox potential (Eh) of the culture medium down to about −480 mV, as long as glucose was available. Addition of oxygen into T. maritima cultures during the stationary growth phase led to a drastic reduction in glucose consumption rate. However, although oxygen was toxic, our experiment unambiguously proved that T. maritima was able to consume it during a 12-hour exposure period. Furthermore, a shift in glucose metabolism towards lactate production was observed under oxidative conditions.

Characterization of cDNA encoding a L37a ribosomal protein from Taenia crassiceps and its potential use in phylogenetic reconstructions

In this study, we characterized for the first time the complete sequence of a L37a cDNA from a cestode specie: Taenia crassiceps. A phylogenetic analysis of L37a ribosomal proteins from distant animal species is presented and the potential use of such proteins in molecule-based phylogeny is discussed.

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.

Protein phylogenies and signature sequences: A reappraisal of evolutionary relationships among archaebacteria, eubacteria, and eukaryotes

The presence of shared conserved insertion or deletions (indels) in protein sequences is a special type of signature sequence that shows considerable promise for phylogenetic inference. An alternative model of microbial evolution based on the use of indels of conserved proteins and the morphological features of prokaryotic organisms is proposed. In this model, extant archaebacteria and gram-positive bacteria, which have a simple, single-layered cell wall structure, are termed monoderm prokaryotes. They are believed to be descended from the most primitive organisms. Evidence from indels supports the view that the archaebacteria probably evolved from gram-positive bacteria, and I suggest that this evolution occurred in response to antibiotic selection pressures. Evidence is presented that diderm prokaryotes (i.e., gram-negative bacteria), which have a bilayered cell wall, are derived from monoderm prokaryotes. Signature sequences in different proteins provide a means to define a number of different taxa within prokaryotes (namely, low G+C and high G+C gram-positive, Deinococcus-Thermus, cyanobacteria, chlamydia-cytophaga related, and two different groups of Proteobacteria) and to indicate how they evolved from a common ancestor. Based on phylogenetic information from indels in different protein sequences, it is hypothesized that all eukaryotes, including amitochondriate and aplastidic organisms, received major gene contributions from both an archaebacterium and a gram-negative eubacterium. In this model, the ancestral eukaryotic cell is a chimera that resulted from a unique fusion event between the two separate groups of prokaryotes followed by integration of their genomes.

Reverse gyrase from the hyperthermophilic bacterium Thermotoga maritima: properties and gene structure

The hyperthermophilic bacterium Thermotoga maritima MSB8 possesses a reverse gyrase whose enzymatic properties are very similar to those of archaeal reverse gyrases. It catalyzes the positive supercoiling of the DNA in an Mg2+- and ATP-dependent process. Its optimal temperature of activity is around 90 degrees C, and it is highly thermostable. We have cloned and DNA sequenced the corresponding gene (T. maritima topR). This is the first report describing the analysis of a gene encoding a reverse gyrase in bacteria. The T. maritima topR gene codes for a protein of 1,104 amino acids with a deduced molecular weight of 128,259, a value in agreement with that estimated from the denaturing gel electrophoresis of the purified enzyme. Like its archaeal homologs, the T. maritima reverse gyrase exhibits helicase and topoisomerase domains, and its sequence matches very well the consensus sequence for six reverse gyrases now available. Phylogenetic analysis shows that all reverse gyrases, including the T. maritima enzyme, form a very homogeneous group, distinct from the type I 5' topoisomerases of the TopA subfamily, for which we have previously isolated a representative gene in T. maritima (topA). The coexistence of these two distinct genes, coding for a reverse gyrase and an omega-like topoisomerase, respectively, together with the recent description of a gyrase in T. maritima (O. Guipaud, E. Marguet, K. M. Noll, C. Bouthier de la Tour, and P. Forterre, Proc. Natl. Acad. Sci. USA 94:10606-10611, 1977) addresses the question of the control of the supercoiling in this organism.

Archaea and the prokaryote-to-eukaryote transition

Since the late 1970s, determining the phylogenetic relationships among the contemporary domains of life, the Archaea (archaebacteria), Bacteria (eubacteria), and Eucarya (eukaryotes), has been central to the study of early cellular evolution. The two salient issues surrounding the universal tree of life are whether all three domains are monophyletic (i.e., all equivalent in taxanomic rank) and where the root of the universal tree lies. Evaluation of the status of the Archaea has become key to answering these questions. This review considers our cumulative knowledge about the Archaea in relationship to the Bacteria and Eucarya. Particular attention is paid to the recent use of molecular phylogenetic approaches to reconstructing the tree of life. In this regard, the phylogenetic analyses of more than 60 proteins are reviewed and presented in the context of their participation in major biochemical pathways. Although many gene trees are incongruent, the majority do suggest a sisterhood between Archaea and Eucarya. Altering this general pattern of gene evolution are two kinds of potential interdomain gene transferrals. One horizontal gene exchange might have involved the gram-positive Bacteria and the Archaea, while the other might have occurred between proteobacteria and eukaryotes and might have been mediated by endosymbiosis.

A chimeric disposition of the elongation factor genes in Rickettsia prowazekii

An exceptional disposition of the elongation factor genes is observed in Rickettsia prowazekii, in which there is only one tuf gene, which is distant from the lone fus gene. In contrast, the closely related bacterium Agrobacterium tumefaciens has the normal bacterial arrangement of two tuf genes, of which one is tightly linked to the fus gene. Analysis of the flanking sequences of the single tuf gene in R. prowazekii shows that it is preceded by two of the four tRNA genes located in the 5' region of the Escherichia coli tufB gene and that it is followed by rpsJ as well as associated ribosomal protein genes, which in E. coli are located downstream of the tufA gene. The fus gene is located within the str operon and is followed by one tRNA gene as well as by the genes secE and nusG, which are located in the 3' region of tufB in E. coli. This atypical disposition of genes suggests that intrachromosomal recombination between duplicated tuf genes has contributed to the evolution of the unique genomic architecture of R. prowazekii.

Sensitivity of ribosomes of the hyperthermophilic bacterium Aquifex pyrophilus to aminoglycoside antibiotics

A poly(U)-programmed cell-free system from the hyperthermophilic bacterium Aquifex pyrophilus has been developed, and the susceptibility of Aquifex ribosomes to the miscoding-inducing and inhibitory actions of all known classes of aminoglycoside antibiotics has been assayed at temperatures (75 to 80 degrees C) close to the physiological optimum for cell growth. Unlike Thermotoga maritima ribosomes, which are systematically refractory to all known classes of aminoglycoside compounds (P. Londei, S. Altamura, R. Huber, K. O. Stetter, and P. Cammarano, J. offteriol. 170-4353-4360, 1988), Aquifex ribosomes are susceptible to all of the aminoglycosides tested (disubstituted 2-deoxystreptamines, monosubstituted 2-deoxystreptamines, sand streptidine compounds). The significance of this result in light of the Aquifex and Thermotoga placements in phylogenetic trees of molecular sequences is discussed.

Arrangement and nucleotide sequence of the gene (fus) encoding elongation factor G (EF-G) from the hyperthermophilic bacterium Aquifex pyrophilus: phylogenetic depth of hyperthermophilic bacteria inferred from analysis of the EF-G/fus sequences

The gene fus (for EF-G) of the hyperthermophilic bacterium Aquifex pyrophilus was cloned and sequenced. Unlike the other bacteria, which display the streptomycin-operon arrangement of EF genes (5'-rps12-rps7-fus-tuf-3'), the Aquifex fus gene (700 codons) is not preceded by the two small ribosomal subunit genes although it is still followed by a tuf gene (for EF-Tu). The opposite strand upstream from the EF-G coding locus revealed an open reading frame (ORF) encoding a polypeptide having 52.5% identity with an E. coli protein (the pdxJ gene product) involved in pyridoxine condensation. The Aquifex EF-G was aligned with available homologs representative of Deinococci, high G+C Gram positives, Proteobacteria, cyanobacteria, and several Archaea. Outgroup-rooted phylogenies were constructed from both the amino acid and the DNA sequences using first and second codon positions in the alignments except sites containing synonymous changes. Both datasets and alternative tree-making methods gave a consistent topology, with Aquifex and Thermotoga maritima (a hyperthermophile) as the first and the second deepest offshoots, respectively. However, the robustness of the inferred phylogenies is not impressive. The branching of Aquifex more deeply than Thermotoga and the branching of Thermotoga more deeply than the other taxa examined are given at bootstrap values between 65 and 70% in the fus-based phylogenies, while the EF-G(2)-based phylogenies do not provide a statistically significant level of support (< or = 50% bootstrap confirmation) for the emergence of Thermotoga between Aquifex and the successive offshoot (Thermus genus). At present, therefore, the placement of Aquifex at the root of the bacterial tree, albeit reproducible, can be asserted only with reservation, while the emergence of Thermotoga between the Aquificales and the Deinococci remains (statistically) indeterminate.

Chromosomal organization and nucleotide sequence of the genes for elongation factors EF-1 alpha and EF-2 and ribosomal proteins S7 and S10 of the hyperthermophilic archaeum Desulfurococcus mobilis

The Desulfurococcus mobilis genes fus (encoding EF-2) and tuf (for EF-1 alpha) were cloned and sequenced together with genes for ribosomal proteins S10 (rps10) and S7 (rps7). Unlike Methanococcus, which displays the bacterial-like fus and tuf gene context 5'-rps12-rps7-fus-tuf-3', and similar to Sulfolobus and Pyrococcus, the Desulfurococcus fus gene (734 codons) has a distinct chromosomal location. Moreover, tuf (441 codons) is the promoter-proximal unit of a three-gene cluster comprising the genes rps10 (98 codons) and tRNA(Ser); the arrangement of the cluster is 5'-tuf-91 bp spacer -rps10-138 bp spacer -tRNA(Ser)-3' and the tuf gene is preceded by a canonical archaeal promoter. The D. mobilis gene rps7 (198 codons) is located further upstream from tuf (535 bp 'silent' intergenic spacing) and no rps12 homolog occurs in its immediate vicinity. Also, judging from putative promoter and transcription termination sequences, rps7 appears to be separately transcribed. Analysis of the predicted fus and tuf gene products revealed the three consensus motifs characteristic of GTP-binding proteins, and the fus-encoded EF-2 protein also displayed the consensus sequence required for ADP-ribosylation by Diphtheria toxin. Both EF sequences were definitely crenarchaeal by comparison with available homologs from other Archaea. Outgroup-rooted phylogenies derived from the sequences of ribosomal proteins S10 and S7 yielded the Sulfolobus-Desulfurococcus association at a high bootstrap confidence level.

Chromosomal organization and nucleotide sequence of the fus-gene encoding elongation factor 2 (EF-2) of the hyperthermophilic archaeum Pyrococcus woesei

A Pyrococcus woesei EcoRI DNA fragment (3400 bp) harbouring the gene fus for elongation factor 2 (EF-2) was cloned and almost completely sequenced. Unlike Methanococcus vannielii (which displays the 'str operon'-like fus and tuf gene context, 5'-rps12-rps7-fus-tuf-3'), and similar to Sulfolobus acidocaldarius and Desulfurococcus mobilis, the Pyrococcus fus gene (732 codons) is unlinked to the rps and tuf genes, and is immediately followed (57 bp intergenic spacing) by an ORF of 106 codons. Both ORFs are preceded by potential archaeal promoters located 52 bp (for fus) and 37 bp (for ORF106) upstream of the putative start codons. The Pyrococcus EF-2(G) equivalent factor is somewhat closer to the eukaryal than to the bacterial homolog, and also shares with the former the C-terminal sequence required for ADP ribosylation of EF-2 by Diphtheria toxin.

Phylogenetic depth of S10 and spc operons: cloning and sequencing of a ribosomal protein gene cluster from the extremely thermophilic bacterium Thermotoga maritima

A segment of Thermotoga maritima DNA spanning 6,613 bp downstream from the gene tuf for elongation factor Tu was sequenced by use of a chromosome walking strategy. The sequenced region comprised a string of 14 tightly linked open reading frames (ORFs) starting 50 bp downstream from tuf. The first 11 ORFs were identified as homologs of ribosomal protein genes rps10, rpl3, rpl4, rpl23, rpl2, rps19, rpl22, rps3, rpl16, rpl29, and rps17 (which in Escherichia coli constitute the S10 operon, in that order); the last three ORFs were homologous to genes rpl14, rpl24, and rpl5 (which in E. coli constitute the three promoter-proximal genes of the spectinomycin operon). The 14-gene string was preceded by putative -35 and -10 promoter sequences situated 5' to gene rps10, within the 50-bp spacing between genes tuf and rps10; the same region exhibited a potential transcription termination signal for the upstream gene cluster (having tuf as the last gene) but displayed also the potential for formation of a hairpin loop hindering the terminator; this suggests that transcription of rps10 and downstream genes may start farther upstream. The similar organization of the sequenced rp genes in the deepest-branching bacterial phyla (T. maritima) and among Archaea has been interpreted as indicating that the S10-spc gene arrangement existed in the (last) common ancestor. The phylogenetic depth of the Thermotoga lineage was probed by use of r proteins as marker molecules: in all except one case (S3), Proteobacteria or the gram-positive bacteria, and not the genus Thermotoga, were the deepest-branching lineage; in only two cases, however, was the inferred branching order substantiated by bootstrap analysis.

Synechocystis sp. PCC6803 fusB gene, located outside of the str operon, encodes a polypeptide related to protein synthesis factor EF-G

Synechocystis sp. PCC6803, a cyanobacterium, possesses an unusual gene (fusB) which encodes a protein with strong homology to protein synthesis elongation factor G (EF-G), although it is not linked to the classical str operon. The fusB gene is redundant, since a Synechocystis gene similar to str operon-encoded fusA genes of other bacteria is also present (based on PCR and hybridization results). There is no evidence for the presence of a fusB homologue in other bacteria. The Synechocystis fusB gene encodes unusual amino acids at some positions that are highly conserved in fusA genes of other prokaryotes.

Evolutionary relationships of bacterial and archaeal glutamine synthetase genes

Glutamine synthetase (GS), an essential enzyme in ammonia assimilation and glutamine biosynthesis, has three distinctive types: GSI, GSII and GSIII. Genes for GSI have been found only in bacteria (eubacteria) and archaea (archaebacteria), while GSII genes only occur in eukaryotes and a few soil-dwelling bacteria. GSIII genes have been found in only a few bacterial species. Recently, it has been suggested that several lateral gene transfers of archaeal GSI genes to bacteria may have occurred. In order to study the evolution of GS, we cloned and sequenced GSI genes from two divergent archaeal species: the extreme thermophile Pyrococcus furiosus and the extreme halophile Haloferax volcanii. Our phylogenetic analysis, which included most available GS sequences, revealed two significant prokaryotic GSI subdivisions: GSI-alpha and GSI-beta. GSI-alpha-genes are found in the thermophilic bacterium, Thermotoga maritima, the low G+C Gram-positive bacteria, and the Euryarchaeota (includes methanogens, halophiles, and some thermophiles). GSI-beta-type genes occur in all other bacteria. GSI-alpha- and GSI-beta-type genes also differ with respect to a specific 25-amino-acid insertion and adenylylation control of GS enzyme activity, both absent in the former but present in the latter. Cyanobacterial genes lack adenylylation regulation of GS and may have secondarily lost it. The GSI gene of Sulfolobus solfataricus, a member of the Crenarchaeota (extreme thermophiles), is exceptional and could not be definitely placed in either subdivision.

Evolution of translational elongation factor (EF) sequences: reliability of global phylogenies inferred from EF-1 alpha(Tu) and EF-2(G) proteins

The EF-2 coding genes of the Archaea Pyrococcus woesei and Desulfurococcus mobilis were cloned and sequenced. Global phylogenies were inferred by alternative tree-making methods from available EF-2(G) sequence data and contrasted with phylogenies constructed from the more conserved but shorter EF-1 alpha(Tu) sequences. Both the monophyly (sensu Hennig) of Archaea and their subdivision into the kingdoms Crenarchaeota and Euryarchaeota are consistently inferred by analysis of EF-2(G) sequences, usually at a high bootstrap confidence level. In contrast, EF-1 alpha(Tu) phylogenies tend to be inconsistent with one another and show low bootstrap confidence levels. While evolutionary distance and DNA maximum parsimony analyses of EF-1 alpha(Tu) sequences do show archaeal monophyly, protein parsimony and DNA maximum-likelihood analyses of these data do not. In no case, however, do any of the tree topologies inferred from EF-1 alpha(Tu) sequence analyses receive significant bootstrap support.

Molecular phylogenies based on ribosomal protein L11, L1, L10, and L12 sequences

Available sequences that correspond to the E. coli ribosomal proteins L11, L1, L10, and L12 from eubacteria, archaebacteria, and eukaryotes have been aligned. The alignments were analyzed qualitatively for shared structural features and for conservation of deletions or insertions. The alignments were further subjected to quantitative phylogenetic analysis, and the amino acid identity between selected pairs of sequences was calculated. In general, eubacteria, archaebacteria, and eukaryotes each form coherent and well-resolved nonoverlapping phylogenetic domains. The degree of diversity of the four proteins between the three groups is not uniform. For L11, the eubacterial and archaebacterial proteins are very similar whereas the eukaryotic L11 is clearly less similar. In contrast, in the case of the L12 proteins and to a lesser extent the L10 proteins, the archaebacterial and eukaryotic proteins are similar whereas the eubacterial proteins are different. The eukaryotic L1 equivalent protein has yet to be identified. If the root of the universal tree is near or within the eubacterial domain, our ribosomal protein-based phylogenies indicate that archaebacteria are monophyletic. The eukaryotic lineage appears to originate either near or within the archaebacterial domain.

Cloning and sequencing of the gene encoding thermostable elongation factor 2 in Sulfolobus solfataricus

The gene (aEF-2) coding for the translation elongation factor 2 (aEF-2) in the thermoacidophilic archaebacterium, Sulfolobus solfataricus, has been cloned and sequenced. The deduced primary structure of aEF-2 is composed of 735 amino acids (aa), excluding the Met start residue. There are no Cys residues and the calculated M(r) is 81,699. In the coding region of aEF-2, the high A + T content greatly influences the codon usage. From the alignment of the primary structure of aEF-2 with that of the analogous factors from the three kingdoms, aa identities were derived. The greatest identity (82%) was found with EF-2 from Sulfolobus acidocaldarius; lower values were observed with other archaebacterial EF-2 (45-47%), eukaryotic EF-2 (38-40%) and with the functional eubacterial analogue EF-G (28-31%). aEF-2 possesses the consensus sequences required for a GTP-binding protein and the four regions which are supposed to be involved in the functional regulation of EF-2/EF-G. These data should have phylogenetic implications.

The DNA-dependent RNA-polymerase of Thermotoga maritima; characterisation of the enzyme and the DNA-sequence of the genes for the large subunits

An improved purification procedure for Thermotoga maritima RNA-polymerase holoenzyme was developed. The enzyme is highly active with poly dAT or T7 phage DNA as template. DNA gyrase was found to be a side product of this RNA-polymerase purification. The genes for the large subunits beta and beta' of RNA-polymerase were cloned and sequenced. The phylogenetic position of T.maritima within the bacterial domain was determined by various methods. It is the lowest bacterial offspring but slightly higher than the chloroplasts.

Early evolutionary relationships among known life forms inferred from elongation factor EF-2/EF-G sequences: phylogenetic coherence and structure of the archaeal domain

Phylogenies were inferred from both the gene and the protein sequences of the translational elongation factor termed EF-2 (for Archaea and Eukarya) and EF-G (for Bacteria). All treeing methods used (distance-matrix, maximum likelihood, and parsimony), including evolutionary parsimony, support the archaeal tree and disprove the "eocyte tree" (i.e., the polyphyly and paraphyly of the Archaea). Distance-matrix trees derived from both the amino acid and the DNA sequence alignments (first and second codon positions) showed the Archaea to be a monophyletic-holophyletic grouping whose deepest bifurcation divides a Sulfolobus branch from a branch comprising Methanococcus, Halobacterium, and Thermoplasma. Bootstrapped distance-matrix treeing confirmed the monophyly-holophyly of Archaea in 100% of the samples and supported the bifurcation of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 97% of the samples. Similar phylogenies were inferred by maximum likelihood and by maximum (protein and DNA) parsimony. DNA parsimony trees essentially identical to those inferred from first and second codon positions were derived from alternative DNA data sets comprising either the first or the second position of each codon. Bootstrapped DNA parsimony supported the monophyly-holophyly of Archaea in 100% of the bootstrap samples and confirmed the division of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 93% of the bootstrap samples. Distance-matrix and maximum likelihood treeing under the constraint that branch lengths must be consistent with a molecular clock placed the root of the universal tree between the Bacteria and the bifurcation of Archaea and Eukarya. The results support the division of Archaea into the kingdoms Crenarchaeota (corresponding to the Sulfolobus branch and Euryarchaeota). This division was not confirmed by evolutionary parsimony, which identified Halobacterium rather than Sulfolobus as the deepest offspring within the Archaea.

Nucleotide sequence of a DNA region comprising the gene for elongation factor 1 alpha (EF-1 alpha) from the ultrathermophilic archaeote Pyrococcus woesei: phylogenetic implications

The gene encoding elongation factor 1 alpha (EF-1 alpha, 1290 bp) of the ultrathermophilic, sulfur-reducing archaeote Pyrococcus woesei was localized within a Bg/II fragment of chromosomal DNA. Sequence analysis showed that the EF-1 alpha gene is the upstream unit of a three-gene cluster comprising the genes for ribosomal protein S10 (306 bp) and transfer RNAser (GGA). The three genes follow each other immediately in the order EF-1 alpha.S10.tRNA(ser) after a putative promoter located 55 bp upstream of the EF-1 alpha gene. Alignment of the derived EF-1 alpha sequence with the corresponding sequences from Eukarya, Bacteria/organelles, and with available archaeal sequences (Sulfolobus, Thermococcus, Methanococcus, Halobacterium) showed that Pyrococcus EF-1 alpha is highly homologous (89% identity) to Thermococcus celer EF-1 alpha, both being strikingly more similar to eukaryotic EF-1 alpha than to bacterial EF-Tu. Unrooted dendrograms computed from aligned sequences by distance matrix and DNA parsimony methods, including evolutionary parsimony, showed the Archaea to be a monophyletic-holophyletic cluster closer to Eukarya than to Bacteria. Both distance matrix and DNA parsimony--although not evolutionary parsimony--support the partition of the known archaeal lineages between the kingdoms Crenarchaeota and Euryarchaeota, and the affiliation of the Pyrococcus-Thermococcus lineage to the Euryarchaeota, of which it is the most primitive offspring. A closer relation of Pyrococcus to Euryarchaeota than to Crenarchaeota was also inferred from sequence analysis of S10 ribosomal proteins.