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

Metabolic and evolutionary relationships among Pyrococcus Species: genetic exchange within a hydrothermal vent environment

Pyrococcus furiosus and Pyrococcus woesei grow optimally at temperatures near 100 degrees C and were isolated from the same shallow marine volcanic vent system. Hybridization of genomic DNA from P. woesei to a DNA microarray containing all 2,065 open reading frames (ORFs) annotated in the P. furiosus genome, in combination with PCR analysis, indicated that homologs of 105 ORFs present in P. furiosus are absent from the uncharacterized genome of P. woesei. Pulsed-field electrophoresis indicated that the sizes of the two genomes are comparable, and the results were consistent with the hypothesis that P. woesei lacks the 105 ORFs found in P. furiosus. The missing ORFs are present in P. furiosus mainly in clusters. These clusters include one cluster (Mal I, PF1737 to PF1751) involved in maltose metabolism and another cluster (PF0691 to PF0695) whose products are thought to remove toxic reactive nitrogen species. Accordingly, it was found that P. woesei, in contrast to P. furiosus, is unable to utilize maltose as a carbon source for growth, and the growth of P. woesei on starch was inhibited by addition of a nitric oxide generator. In P. furiosus the ORF clusters not present in P. woesei are bracketed by or are in the vicinity of insertion sequences or long clusters of tandem repeats (LCTRs). While the role of LCTRs in lateral gene transfer is not known, the Mal I cluster in P. furiosus is a composite transposon that undergoes replicative transposition. The same locus in P. woesei lacks any evidence of insertion activity, indicating that P. woesei is a sister or even the parent of P. furiosus. P. woesei may have acquired by lateral gene transfer more than 100 ORFs from other organisms living in the same thermophilic environment to produce the type strain of P. furiosus.

GRS: a graphic tool for genome retrieval and segment analysis

GRS is a graphic tool for retrieval and visualization of genome segments from partially or completely sequenced genomes. To facilitate visual identification of conserved genomic motifs, genes are color-coded according to their presumed functional roles. Aligned genes can be rapidly screened for potential homology by automatic retrieval and alignment of the corresponding protein sequences. Furthermore, the map location of any genome segment can be visually compared to the position of the same segment in other genomes or to the position of other segments within the same genome. The gene string analysis option of GRS allows the identification of genes that are identically arranged in any pairwise set of genomes. Finally, the program allows the user to create new gene table format files to enable comparisons of gene order structures in recently determined sequence data to the patterns of genes in already existing microbial and organellar databases. With the help of GRS, the genomic contexts of genes for which no identifiable homologues exist can be analyzed to provide an additional source of information for sequence annotations. We illustrate the use of GRS by analyzing the structure and distribution of phylogenetically conserved motifs in closely as well as more distantly related microbial genomes.

Molecular biology of hyperthermophilic Archaea

The sequences of a number of archaeal genomes have recently been completed, and many more are expected shortly. Consequently, the research of Archaea in general and hyperthermophiles in particular has entered a new phase, with many exciting discoveries to be expected. The wealth of sequence information has already led, and will continue to lead to the identification of many enzymes with unique properties, some of which have potential for industrial applications. Subsequent functional genomics will help reveal fundamental matters such as details concerning the genetic, biochemical and physiological adaptation of extremophiles, and hence give insight into their genomic evolution, polypeptide structure-function relations, and metabolic regulation. In order to optimally exploit many unique features that are now emerging, the development of genetic systems for hyperthermophilic Archaea is an absolute requirement. Such systems would allow the application of this class of Archaea as so-called "cell factories": (i) expression of certain archaeal enzymes for which no suitable conventional (mesophilic bacterial or eukaryal) systems are available, (ii) selection for thermostable variants of potentially interesting enzymes from mesophilic origin, and (iii) the development of in vivo production systems by metabolic engineering. An overview is given of recent insight in the molecular biology of hyperthermophilic Archaea, as well as of a number of promising developments that should result in the generation of suitable genetic systems in the near future.

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.

Direct linkage of str-, S10- and spc-related gene clusters in Thermus thermophilus HB8, and sequences of ribosomal proteins L4 and S10

We have analyzed the genomic region harboring str-S10-spc-related gene clusters in the thermophilic bacterium, Thermus thermophilus (Tt) HB8. This study was initiated for the purpose of isolating the gene encoding ribosomal (r-) protein S10 which is assumed to be involved in the antitermination of transcription at the rRNA-encoding genes in Bacteria. The S10-related gene cluster encodes the same set of r-proteins as in Escherichia coli. However, the gene coding for elongation factor Tu (the last gene of the str operon in E. coli) is separated by only eight nucleotides (nt) from the gene encoding r-protein S10 (the first gene of the S10 operon in E. coli), and the genes encoding r-protein S17 (the last gene of the S10 operon in E. coli) and L14 (the first gene of the spc operon in E. coli) overlap. This suggests that all three gene clusters are cotranscribed from a single promoter preceding the str-related operon. In addition, we determined the complete nt sequences of the Tt genes encoding r-proteins L4 and S10. Tt L4 shows the lowest degree of conservation among the known L4 r-proteins from Bacteria. Tt S10 has the highest proportion of basic amino acids (aa) and the lowest number of acidic aa, as compared with its homologues from Bacteria and Archaea, which might be related to its possible role in binding to the boxA RNA of nascent rRNA transcripts at high temperatures.

Genomic evolution drives the evolution of the translation system

Our thesis is that the characteristics of the translational machinery and its organization are selected in part by evolutionary pressure on genomic traits have nothing to do with translation per se. These genomic traits include size, composition, and architecture. To illustrate this point, we draw parallels between the structure of different genomes that have adapted to intracellular niches independently of each other. Our starting point is the general observation that the evolutionary history of organellar and parasitic bacteria have favored bantam genomes. Furthermore, we suggest that the constraints of the reductive mode of genomic evolution account for the divergence of the genetic code in mitochondria and the genetic organization of the translational system observed in parasitic bacteria. In particular, we associate codon reassignments in animal mitochondria with greatly simplified tRNA populations. Likewise, we relate the organization of translational genes in the obligate intracellular parasite Rickettsia prowazekii to the processes supporting the reductive mode of genomic evolution. Such findings provide strong support for the hypothesis that genomes of organelles and of parasitic bacteria have arisen from the much larger genomes of ancestral bacteria that have been reduced by intrachromosomal recombination and deletion events. A consequence of the reductive mode of genomic evolution is that the resulting translation systems may deviate markedly from conventional systems.

Calculating the probability of multitaxon evolutionary trees: bootstrappers Gambit

The reconstruction of multitaxon trees from molecular sequences is confounded by the variety of algorithms and criteria used to evaluate trees, making it difficult to compare the results of different analyses. A global method of multitaxon phylogenetic reconstruction described here, Bootstrappers Gambit, can be used with any four-taxon algorithm, including distance, maximum likelihood, and parsimony methods. It incorporates a Bayesian-Jeffreys'-bootstrap analysis to provide a uniform probability-based criterion for comparing the results from diverse algorithms. To examine the usefulness of the method, the origin of the eukaryotes has been investigated by the analysis of ribosomal small subunit RNA sequences. Three common algorithms (paralinear distances, Jukes-Cantor distances, and Kimura distances) support the eocyte topology, whereas one (maximum parsimony) supports the archaebacterial topology, suggesting that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.

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.

Archaebacterial genomes: eubacterial form and eukaryotic content

Since the recognition of the uniqueness and coherence of the archaebacteria (sometimes called Archaea), our perception of their role in early evolution has been modified repeatedly. The deluge of sequence data and rapidly improving molecular systematic methods have combined with a better understanding of archaebacterial molecular biology to describe a group that in some ways appears to be very similar to the eubacteria, though in others is more like the eukaryotes. The structure and contents of archaebacterial genomes are examined here, with an eye to their meaning in terms of the evolution of cell structure and function.

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.

Bacterial phylogeny based on 16S and 23S rRNA sequence analysis

Molecular phylogeny increasingly supports the understanding of organismal relationships and provides the basis for the classification of microorganisms according to their natural affiliations. Comparative sequence analysis of ribosomal RNAs or the corresponding genes currently is the most widely used approach for the reconstruction of microbial phylogeny. The highly and less conserved primary and higher order structure elements of rRNAs document the history of microbial evolution and are informative for definite phylogenetic levels. An optimal alignment of the primary structures and a careful data selection are prerequisites for reliable phylogenetic conclusions. rRNA based phylogenetic trees can be reconstructed and the significance of their topologies evaluated by applying distance, maximum parsimony and maximum likelihood methods of phylogeny inference in comparison, and by fortuitous or directed resampling of the data set. Phylogenetic trees based on almost equivalent data sets of bacterial 23S and 16S rRNAs are in good agreement and their overall topologies are supported by alternative phylogenetic markers such as elongation factors and ATPase subunits. Besides their phylogenetic information content, the differently conserved primary structure regions of rRNAs provide target sites for specific hybridization probes which have been proven to be powerful tools for the identification of microbes on the basis of their phylogenetic relationships.

The TATA-binding protein: a general transcription factor in eukaryotes and archaebacteria

The TATA-binding protein TBP appears to be essential for all transcription in eukaryotic cell nuclei, which suggests that its function was established early in evolution. Archaebacteria constitute a kingdom of organisms distinct from eukaryotes and eubacteria. Archaebacterial gene regulatory sequences often map to TATA box-like motifs. Here it is shown that the archaebacterium Pyrococcus woesei expresses a protein with structural and functional similarity to eukaryotic TBP molecules. This suggests that TBP's role in transcription was established before the archaebacterial and eukaryotic lineages diverged and that the transcription systems of archaebacteria and eukaryotes are fundamentally homologous.

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.

The nucleotide sequence of the gene coding for the elongation factor 1 alpha in Sulfolobus solfataricus. Homology of the product with related proteins

The cloning and sequencing of the gene coding for the archaebacterial elongation factor 1 alpha (aEF-1 alpha) was performed by screening a Sulfolobus solfataricus genomic library using a probe constructed from the eptapeptide KNMITGA that is conserved in all the EF-1 alpha/EF-Tu known so far. The isolated recombinant phage contained the part of the aEF-1 alpha gene from amino acids 1 to 171. The other part (amino acids 162-435) was obtained through the amplification of the S. solfataricus DNA by PCR. The codon usage by the aEF-1 alpha gene showed a preference for triplets ending in A and/or T. This behavior was almost identical to that of the S. acidocaldarius EF-1 alpha gene but differed greatly from that of EF-1 alpha/EF-Tu genes in other archaebacteria eukaryotes and eubacteria. The translated protein is made of 435 amino acid residues and contains sequence motifs for the binding of GTP, tRNA and ribosome. Alignments of aEF-1 alpha with several EF-1 alpha/EF-Tu revealed that aEF-1 alpha is more similar to its eukaryotic than to its eubacterial counterparts.

DNA-dependent RNA polymerase subunit B as a tool for phylogenetic reconstructions: branching topology of the archaeal domain

The branching topology of the archaeal (archaebacterial) domain was inferred from sequence comparisons of the largest subunit (B) of DNA-dependent RNA polymerases (RNAP). Both the nucleic acid sequences of the genes coding for RNAP subunit B and the amino acid sequences of the derived gene products were used for phylogenetic reconstructions. Individual analysis of the three nucleotide positions of codons revealed significant inequalities with respect to guanosine and cytosine (GC) content and evolutionary rates. Only the nucleotides at the second codon positions were found to be unbiased by varied GC contents and sufficiently conserved for reliable phylogenetic reconstructions. A decision matrix was used for the combination of the results of distance matrix, maximum parsimony, and maximum likelihood methods. For this purpose the original results (sums of squares, steps, and logarithms of likelihoods) were transformed into comparable effective values and analyzed with methods known from the theory of statistical decisions. Phylogenetic invariants and statistical analysis with resampling techniques (bootstrap and jackknife) confirmed the preferred branching topology, which is significantly different from the topology known from phylogenetic trees based on 16S rRNA sequences. The preferred topology reconstructed by this analysis shows a common stem for the Methanococcales and Methanobacteriales and a separation of the thermophilic sulfur archaea from the methanogens and halophiles. The latter coincides with a unique phylogenetic location of a characteristic splitting event replacing the largest RNAP subunit of thermophilic sulfur archaea by two fragments in methanogens and halophiles. This topology is in good agreement with physiological and structural differences between the various archaea and demonstrates RNAP to be a suitable phylogenetic marker molecule.

Reconstructing evolutionary trees from DNA and protein sequences: paralinear distances

The reconstruction of phylogenetic trees from DNA and protein sequences is confounded by unequal rate effects. These effects can group rapidly evolving taxa with other rapidly evolving taxa, whether or not they are genealogically related. All algorithms are sensitive to these effects whenever the assumptions on which they are based are not met. The algorithm presented here, called paralinear distances, is valid for a much broader class of substitution processes than previous algorithms and is accordingly less affected by unequal rate effects. It may be used with all nucleic acid, protein, or other sequences, provided that their evolution may be modeled as a succession of Markov processes. The properties of the method have been proven both analytically and by computer simulations. Like all other methods, paralinear distances can fail when sequences are misaligned or when site-to-site sequence variation of rates is extensive. To examine the usefulness of paralinear distances, the "origin of the eukaryotes" has been investigated by the analysis of elongation factor Tu sequences with a variety of sequence alignments. It has been found that the order in which sequences are pairwise aligned strongly determines the topology which is reconstructed by paralinear distances (as it does for all other reconstruction methods tested). When the parts of the alignment that are unaffected by alignment order are analyzed, paralinear distances strongly select the eocyte topology. This provides evidence that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.

Phylogenetic relationships of Bacteria based on comparative sequence analysis of elongation factor Tu and ATP-synthase beta-subunit genes

Comparative sequence analyses were performed on 14 genes encoding bacterial elongation factors EF-Tu and 7 genes encoding the beta-subunit of bacterial F1F0 type ATP-synthases. The corresponding predicted amino acid sequences were compared with published primary structures of homologous molecules. Phylogenetic trees were reconstructed from both data sets of aligned protein sequences and from an equivalent selection of 16S rRNA sequences by applying distance matrix and maximum parsimony methods. The EF-Tu data were in very good agreement with the rRNA data, although the resolution within the EF-Tu tree was reduced at certain phylogenetic levels. The resolution power of the ATPase beta-subunit sequence data were more reduced than those of the EF-Tu data. In comparison with the 16S rRNA tree there are minor differences in the order of adjacent branchings within the ATPase beta-subunit tree.

Evidence for a chimeric nature of nuclear genomes: eubacterial origin of eukaryotic glyceraldehyde-3-phosphate dehydrogenase genes

Higher plants process two distinct, nuclear gene-encoded glyceraldehyde-3-phosphate dehydrogenase (GAPDH) proteins, a Calvin-cycle enzyme active within chloroplasts and a glycolytic enzyme active within the cytosol. The gene for the chloroplast enzyme was previously suggested to be of endosymbiotic origin. Since the ancestors of plastids were related to cyanobacteria, we have studied GAPDH genes in the cyanobacterium Anabaena variabilis. Our results confirm that the nuclear gene for higher plant chloroplast GAPDH indeed derives from the genome of a cyanobacterium-like endosymbiont. But two additional GAPDH genes were found in the Anabaena genome and, surprisingly, one of these sequences is very similar to nuclear genes encoding the GAPDH enzyme of glycolysis in plants, animals, and fungi. Evidence that the eukaryotic nuclear genes for glycolytic GAPDH, as well as the Calvin-cycle genes, are of eubacterial origin suggests that eukaryotic genomes are more highly chimeric than previously assumed.

Nucleotide sequence and molecular evolution of the gene coding for glyceraldehyde-3-phosphate dehydrogenase in the thermoacidophilic archaebacterium Sulfolobus solfataricus

A Sulfolobus solfataricus genomic library cloned in the EMBL3 phage was screened using as probes synthetic oligonucleotides designed from the known amino acid sequence of a peptide obtained from the purified glyceraldehyde-3-phosphate dehydrogenase (aGAPD) protein. The screening led to the isolation of six recombinant phages (lambda G1-lambda G6) and one of them (lambda G4) contained the entire GAPD gene. The deduced amino acid sequence accounts for a protein made of 341 amino acids and the initial methionine is encoded by a GTG triplet. Alignment of the S. solfataricus aGAPD sequence versus GAPD from archaea, eukarya, and bacteria showed that aGAPD is very similar to other archaebacterial but not to eukaryotic or eubacterial GAPD. For known archaebacterial GAPD sequences, the rate of nucleotide substitutions per site per year showed that these sequences are homologous not only at the amino acid but also at the nucleotide level. The evolutionary rates are nearly similar to those reported for other eukaryotic genes.

Cloning and sequencing of the gene encoding glutamine synthetase I from the archaeum Pyrococcus woesei: anomalous phylogenies inferred from analysis of archaeal and bacterial glutamine synthetase I sequences

The gene glnA encoding glutamine synthetase I (GSI) from the archaeum Pyrococcus woesei was cloned and sequenced with the Sulfolobus solfataricus glnA gene as the probe. An operon reading frame of 448 amino acids was identified within a DNA segment of 1,528 bp. The encoded protein was 49% identical with the GSI of Methanococcus voltae and exhibited conserved regions characteristic of the GSI family. The P. woesei GSI was aligned with available homologs from other archaea (S. solfataricus, M. voltae) and with representative sequences from cyanobacteria, proteobacteria, and gram-positive bacteria. Phylogenetic trees were constructed from both the amino acid and the nucleotide sequence alignments. In accordance with the sequence similarities, archaeal and bacterial sequences did not segregate on a phylogeny. On the basis of sequence signatures, the GSI trees could be subdivided into two ensembles. One encompassed the GSI of cyanobacteria and proteobacteria, but also that of the high-G + C gram-positive bacterium Streptomyces coelicolor (all of which are regulated by the reversible adenylylation of the enzyme subunits); the other embraced the GSI of the three archaea as well as that of the low-G + C gram-positive bacteria (Clostridium acetobutilycum, Bacillus subtilis) and Thermotoga maritima (none of which are regulated by subunit adenylylation). The GSIs of the Thermotoga and the Bacillus-Clostridium lineages shared a direct common ancestor with that of P. woesei and the methanogens and were unrelated to their homologs from cyanobacteria, proteobacteria, and S. coelicolor. The possibility is presented that the GSI gene arose among the archaea and was then laterally transferred from some early methanogen to a Thermotoga-like organism. However, the relationship of the cyanobacterial-proteobacterial GSIs to the Thermotoga GSI and the GSI of low-G+C gram-positive bacteria remains unexplained.

Evolution of glutamate dehydrogenase genes: evidence for two paralogous protein families and unusual branching patterns of the archaebacteria in the universal tree of life

The existence of two families of genes coding for hexameric glutamate dehydrogenases has been deduced from the alignment of 21 primary sequences and the determination of the percentages of similarity between each pair of proteins. Each family could also be characterized by specific motifs. One family (Family I) was composed of gdh genes from six eubacteria and six lower eukaryotes (the primitive protozoan Giardia lamblia, the green alga Chlorella sorokiniana, and several fungi and yeasts). The other one (Family II) was composed of gdh genes from two eubacteria, two archaebacteria, and five higher eukaryotes (vertebrates). Reconstruction of phylogenetic trees using several parsimony and distance methods confirmed the existence of these two families. Therefore, these results reinforced our previously proposed hypothesis that two close but already different gdh genes were present in the last common ancestor to the three Ur-kingdoms (eubacteria, archaebacteria, and eukaryotes). The branching order of the different species of Family I was found to be the same whatever the method of tree reconstruction although it varied slightly according the region analyzed. Similarly, the topological positions of eubacteria and eukaryotes of Family II were independent of the method used. However, the branching of the two archaebacteria in Family II appeared to be unexpected: (1) the thermoacidophilic Sulfolobus solfataricus was found clustered with the two eubacteria of this family both in parsimony and distance trees, a situation not predicted by either one of the contradictory trees recently proposed; and (2) the branching of the halophilic Halobacterium salinarium varied according to the method of tree construction: it was closer to the eubacteria in the maximum parsimony tree and to eukaryotes in distance trees. Therefore, whatever the actual position of the halophilic species, archaebacteria did not appear to be monophyletic in these gdh gene trees. This result questions the firmness of the presently accepted interpretation of previous protein trees which were supposed to root unambiguously the universal tree of life and place the archaebacteria in this tree.

A single amino acid substitution in elongation factor Tu disrupts interaction between the ternary complex and the ribosome

Elongation factor Tu (EF-Tu).GTP has the primary function of promoting the efficient and correct interaction of aminoacyl-tRNA with the ribosome. Very little is known about the elements in EF-Tu involved in this interaction. We describe a mutant form of EF-Tu, isolated in Salmonella typhimurium, that causes a severe defect in the interaction of the ternary complex with the ribosome. The mutation causes the substitution of Val for Gly-280 in domain II of EF-Tu. The in vivo growth and translation phenotypes of strains harboring this mutation are indistinguishable from those of strains in which the same tuf gene is insertionally inactivated. Viable cells are not obtained when the other tuf gene is inactivated, showing that the mutant EF-Tu alone cannot support cell growth. We have confirmed, by partial protein sequencing, that the mutant EF-Tu is present in the cells. In vitro analysis of the natural mixture of wild-type and mutant EF-Tu allows us to identify the major defect of this mutant. Our data shows that the EF-Tu is homogeneous and competent with respect to guanine nucleotide binding and exchange, stimulation of nucleotide exchange by EF-Ts, and ternary complex formation with aminoacyl-tRNA. However various measures of translational efficiency show a significant reduction, which is associated with a defective interaction between the ribosome and the mutant EF-Tu.GTP.aminoacyl-tRNA complex. In addition, the antibiotic kirromycin, which blocks translation by binding EF-Tu on the ribosome, fails to do so with this mutant EF-Tu, although it does form a complex with EF-Tu. Our results suggest that this region of domain II in EF-Tu has an important function and influences the binding of the ternary complex to the codon-programmed ribosome during protein synthesis. Models involving either a direct or an indirect effect of the mutation are discussed.

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.

The nature of the last universal ancestor and the root of the tree of life, still open questions

The nature of the last universal ancestor to all extent cellular organisms and the rooting of the universal tree of life are fundamental questions which can now be addressed by molecular evolutionists. Several scenarios have been proposed during the last years, based on the phylogenies of ribosomal RNA and of duplicated proteins, which suggest that the last universal ancestor was either an RNA progenote or an hyperthermophilic prokaryote. We discuss these hypotheses in the light of new data on the evolution of DNA metabolizing enzymes and of contradictions between different protein phylogenies. We conclude that the last universal ancestor was a member of the DNA world already containing several DNA polymerases and DNA topoisomerases. Furthermore, we criticize current data which suggest that the rooting of the universal tree of life is located in the eubacterial branch and we conclude that both rooting the universal tree and the nature of the last universal ancestor are still open questions.

Comparative biochemistry of Archaea and Bacteria

This review compares exemplary molecular and metabolic features of Archaea and Bacteria in terms of phylogenetic aspects. The results of the comparison confirm the coherence of the Archaea as postulated by Woese. Archaea and Bacteria share many basic features of their genetic machinery and their central metabolism. Similarities and distinctions allow projections regarding the nature of the common ancestor and the process of lineage diversification.