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

Surface protein EF3314 contributes to virulence properties of Enterococcus faecalis

PURPOSE

Identification of putative new virulence factors as additional targets for therapeutic approaches alternative to antibiotic treatment of multi-resistant enterococcal infections.

METHODS

The EF3314 gene, coding for a putative surface-exposed antigen, was identified by the analysis of the Enterococcus faecalis V583 genome for LPXTG-motif cell wall anchor surface protein genes. A non-polar EF3314 gene deletion mutant in the E. faecalis 12030 human clinical isolate was obtained. The wild type and the isogenic mutant strain were investigated for biofilm formation, adherence to Hela cells, survival in human macrophages and a Caenorhabditis elegans infection model. The aminoterminal portion of the EF3314 protein was overexpressed in E. coli to obtain mouse polyclonal antibodies for use in Western blotting and immunolocalization experiments.

RESULTS

The EF3314 gene has an unusually high GC content (46.88% vs. an average of 37.5% in the E. faecalis chromosome) and encodes a protein of 1744 amino acids that presents a series of 14 imperfect repeats of 90 amino acids covering almost the entire length of the protein. Its global organization is similar to the alpha-like protein family of group B streptococci, enterococcal surface protein Esp and biofilm associated protein Bap from S. aureus. The EF3314 gene was always present and specific for E. faecalis strains of human, food and animal origin. Differences in size depended on variable numbers of repeats in the repetitive region.

CONCLUSIONS

EF3314 is a newly described, surface exposed protein that contributes to the virulence properties of E. faecalis.

Molecular and functional properties of the psychrophilic elongation factor G from the Antarctic Eubacterium Pseudoalteromonas haloplanktis TAC 125

The molecular and functional properties of the elongation factor (EF) G from the psychrophilic Antarctic eubacterium Pseudoalteromonas haloplanktis (Ph) were studied. PhEF-G catalyzed protein synthesis in vitro that was inhibited by fusidic acid, an antibiotic specifically acting on EF-G. The EF interacted with GDP only in the presence of P. haloplanktis ribosome and fusidic acid with an affinity similar to that displayed by Escherichia coli EF-G. The psychrophilic translocase elicited a ribosome-dependent GTPase that was competitively inhibited by GDP, the slowly hydrolyzable GTP analog GppNHp, and the protein synthesis inhibitor ppGDP. The temperature dependence of the activity of PhEF-G reached its maximum at least 26 degrees C beyond the growth temperature of P. haloplanktis (4-20 degrees C). The heat inactivation profile of the ribosome-dependent GTPase of PhEF-G gave a temperature for half inactivation (46 degrees C), significantly lower than that for half denaturation measured by either UV- (57 degrees C) or fluorescence-melting (62 degrees C). This finding was attributed to a different effect of the temperature on the catalytic domain with respect to that elicited on the other domains constituting the EF, thus confirming the differential molecular flexibility present in psychrophilic enzymes. A molecular model, based on the 3D coordinates of a thermophilic EF-G, showed differences only in connecting loops.

Archaeal elongation factor 1alpha from Sulfolobus solfataricus interacts with the eubacterial antibiotic GE2270A

The thiazolyl-peptide antibiotic GE2270A, an inhibitor of the elongation factor Tu from Escherichia coli (EcEF-Tu), was used to study the effects produced in the biochemical properties of the archaeal functional analogue elongation factor 1alpha from Sulfolobus solfataricus (SsEF-1alpha). GE2270A did not substantially affect the poly(U)-directed-polyPhe incorporation catalyzed by SsEF-1alpha and the formation of the ternary complex SsEF-1alpha.GTP.Phe-tRNAPhe. On the other hand, the antibiotic was able to increase the GDP/GTP exchange rate of SsEF-1alpha; nevertheless, this improvement was not associated with an increase in the catalytic activity of the enzyme. In fact, GE2270A inhibited both the intrinsic GTPase of SsEF-1alpha (GTPaseNa) and that stimulated by ribosomes. Interestingly, GTPaseNa of both intact and C-terminal-deleted SsEF-1alpha resulted in a greater sensitivity to the antibiotic with respect to SsEF-1alpha lacking both the M- and C-terminal domains. This result suggested that, similar to what is found for EcEF-Tu, the M domain of SsEF-1alpha is the region of the enzyme most responsible for the interaction with GE2270A. The different behavior observed in the inhibition of protein synthesis with respect to EcEF-Tu can be ascribed to the different adaptive structural changes that have occurred in SsEF-1alpha during evolution.

Conserved but nonessential interaction of SRP RNA with translation factor EF-G

4.5S RNA is essential for viability of Escherichia coli, and forms a key component of the signal recognition particle (SRP), a ubiquitous ribonucleoprotein complex responsible for cotranslational targeting of secretory proteins. 4.5S RNA also binds independently to elongation factor G (EF-G), a five-domain GTPase that catalyzes the translocation step during protein biosynthesis on the ribosome. Point mutations in EF-G suppress deleterious effects of 4.5S RNA depletion, as do mutations in the EF-G binding site within ribosomal RNA, suggesting that 4.5S RNA might play a critical role in ribosome function in addition to its role in SRP. Here we show that 4.5S RNA and EF-G form a phylogenetically conserved, low-affinity but highly specific complex involving sequence elements required for 4.5S binding to its cognate SRP protein, Ffh. Mutational analysis indicates that the same molecular structure of 4.5S RNA is recognized in each case. Surprisingly, however, the suppressor mutant forms of EF-G bind very weakly or undetectably to 4.5S RNA, implying that cells can survive 4.5S RNA depletion by decreasing the affinity between 4.5S RNA and the translational machinery. These data suggest that SRP function is the essential role of 4.5S RNA in bacteria.

Cloning, expression and evolution of the gene encoding the elongation factor 1alpha from a low thermophilic Sulfolobus solfataricus strain

The gene encoding the elongation factor 1alpha (EF-1alpha) from the archaeon Sulfolobus solfataricus strain MT3 (optimum growth temperature 75 degrees C) was cloned, sequenced and expressed in Escherichia coli. The structural and biochemical properties of the purified enzyme were compared to those of EF-1alpha isolated from S. solfataricus strain MT4 (optimum growth temperature 87 degrees C). Only one amino acid change (Val15-->Ile) was found. Interestingly, the difference was in the first guanine nucleotide binding consensus sequence G(13)HIDHGK and was responsible for a reduced efficiency in protein synthesis, which was accompanied by an increased affinity for both guanosine diphosphate (GDP) and guanosine triphosphate (GTP), and an increased efficiency in the intrinsic GTPase activity. Despite the different thermophilicities of the two microorganisms, only very marginal effects on the thermal properties of the enzyme were observed. Molecular evolution among EF-1alpha genes from Sulfolobus species showed that the average rate of nucleotide substitution per site per year (0.0312x10(-9)) is lower than that reported for other functional genes.

Heterogeneity of genome and proteome content in bacteria, archaea, and eukaryotes

Our analysis compares bacteria, archaea, and eukaryota with respect to a wide assortment of genome and proteome properties. These properties include ribosomal protein gene distributions, chaperone protein contrasts, major variation of transcription/translation factors, gene encoding pathways of energy metabolism, and predicted protein expression levels. Significant differences within and between the three domains of life include protein lengths, information processing procedures, many metabolic and lipid biosynthesis pathways, cellular controls, and regulatory proteins. Differences among genomes are influenced by lifestyle, habitat, physiology, energy sources, and other factors.

Comparative genomic analysis of archaeal genotypic variants in a single population and in two different oceanic provinces

Planktonic crenarchaeotes are present in high abundance in Antarctic winter surface waters, and they also make up a large proportion of total cell numbers throughout deep ocean waters. To better characterize these uncultivated marine crenarchaeotes, we analyzed large genome fragments from individuals recovered from a single Antarctic picoplankton population and compared them to those from a representative obtained from deeper waters of the temperate North Pacific. Sequencing and analysis of the entire DNA insert from one Antarctic marine archaeon (fosmid 74A4) revealed differences in genome structure and content between Antarctic surface water and temperate deepwater archaea. Analysis of the predicted gene products encoded by the 74A4 sequence and those derived from a temperate, deepwater planktonic crenarchaeote (fosmid 4B7) revealed many typical archaeal proteins but also several proteins that so far have not been detected in archaea. The unique fraction of marine archaeal genes included, among others, those for a predicted RNA-binding protein of the bacterial cold shock family and a eukaryote-type Zn finger protein. Comparison of closely related archaea originating from a single population revealed significant genomic divergence that was not evident from 16S rRNA sequence variation. The data suggest that considerable functional diversity may exist within single populations of coexisting microbial strains, even those with identical 16S rRNA sequences. Our results also demonstrate that genomic approaches can provide high-resolution information relevant to microbial population genetics, ecology, and evolution, even for microbes that have not yet been cultivated.

Long branch attraction effects and the status of "basal eukaryotes": phylogeny and structural analysis of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis

The three taxa emerging at the base of the eukaryotic ribosomal RNA phylogenetic tree (Diplomonadida, Microspora, and Parabasalia) include a wide array of parasitic species. and some free-living organisms that appear to be derived from a parasitic ancestry. The basal position of these taxa, which lack mitochondria, has recently been questioned. I sequenced most of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis and a secondary structure model was reconstructed for the SSU rRNA. I conducted a RASA matrix analysis to identify, independently from tree reconstruction, putative long branch attraction effects in the data matrix. The results show that each of the basal clades and the euglenozoan clade act, indeed, as long branches and may have been engaged in a process of accelerated rate of evolution. A nucleotide signature analysis was conducted in the conserved regions for positions defining the three great domains of life (Eubacteria, Archea, and Eukaryota). For the three basal taxa, this analysis showed the presence of a significant number of different non-eukaryotic nucleotides. A precise study of the nature and location of these nucleotides led to conclusions supporting the results of the RASA analysis. Altogether, these findings suggest that the basal placement of these taxa in the SSU ribosomal RNA phylogenetic tree is artifactual, and flawed by long branch attraction effects.

The role of the codon first letter in the relationship between genomic GC content and protein amino acid composition

Analysis of the statistical distribution of amino acid compositions within 22 protein families shows that a GC bias generally affects proteins with a variety of functions from the extreme thermophile Thermus. This results in evident enrichment in amino acids of the group L, V, A, P, R and G and underrepresentation of amino acids of the group I, M, F, S, T, C and W. The strong amino acid composition biases noted in Thermus proteins are not related to thermoadaptation; they were also found in mesophilic homologues encoded by GC-rich genes. The results of a comparative analysis on large samples of translated sequences from 30 organisms, representing the three major kingdoms of life and including extremophiles, indicate a universal correlation between the usage of particular amino acids and the genomic GC content. It is concluded that the codon first letter plays a dominant role in translating the genomic GC signature into protein amino acid composition and sequences.

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.

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.

Duplicate genes and the root of angiosperms, with an example using phytochrome sequences

The root of the angiosperm tree has not yet been established. Major morphological and molecular differences between angiosperms and other seed plants have introduced ambiguities and possibly spurious results. Because it is unlikely that extant species more closely related to angiosperms will be discovered, and because relevant fossils will almost certainly not yield molecular data, the use of duplicate genes for rooting purposes may provide the best hope of a solution. Simultaneous analysis of the genes resulting from a gene duplication event along the branch subtending angiosperms would yield an unrooted network, wherein two congruent gene trees should be connected by a single branch. In these circumstances the best rooted species tree is the one that corresponds to the two gene trees when the network is rooted along the connecting branch. In general, this approach can be viewed as choosing among rooted species trees by minimizing hypothesized events such as gene duplication, gene loss, lineage sorting, and lateral transfer. Of those gene families that are potentially relevant to the angiosperm problem, phytochrome genes warrant special attention. Phylogenetic analysis of a sample of complete phytochrome (PHY) sequences implies that an initial duplication event preceded (or occurred early within) the radiation of seed plants and that each of the two resulting copies duplicated again. In one of these cases, leading to the PHYA and PHYC lineages, duplication appears to have occurred before the diversification of angiosperms. Duplicate gene trees are congruent in these broad analyses, but the sample of sequences is too limited to provide much insight into the rooting question. Preliminary analyses of partial PHYA and PHYC sequences from several presumably basal angiosperm lineages are promising, but more data are needed to critically evaluate the power of these genes to resolve the angiosperm radiation.

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.

Elongation factor Tu1 of the antibiotic GE2270A producer Planobispora rosea has an unexpected resistance profile against EF-Tu targeted antibiotics

Sensitivity of EF-Tu1 of the GE2270A producer Planobispora rosea towards GE2270A, pulvomycin and kirromycin was determined by band-shift assays for EF-Tu1-antibiotic complex formation and by in vitro translation experiments. EF-Tu1 of P. rosea appeared to be not only totally resistant to GE2270A, but also ten times more resistant to kirromycin than EF-Tu1 of Streptomyces coelicolor. In contrast, P. rosea EF-Tu1 was found to be not resistant to pulvomycin, an antibiotic that just like GE2270A blocks EF-Tu x GTP x aminoacyl-tRNA complex formation. Previous in vivo and in vitro experiments with mixed populations of antibiotic resistant and sensitive EF-Tu species had shown that sensitivity to kirromycin and pulvomycin is dominant over resistance. In the case of GE2270A we observed, however, that sensitivity is recessive to resistance, which again points to a different action mechanism than in the case of pulvomycin. Besides the tuf1 gene encoding the regular elongation factor EF-Tu1 a gene similar to S. coelicolor tuf3 for a specialized EF-Tu was located in the P. rosea genome. The tuf1 gene was isolated and sequenced. The amino acid sequence of EF-Tul of P. rosea not only exhibits an unusual Tyr160 substitution (comparable to those described for kirromycin-resistant EF-Tus), but also shows significant changes of conserved amino acids in domain 2 that may be responsible for GE2270A resistance (the latter do not resemble those leading to pulvomycin resistance). P. rosea EF-Tu1 thus is a first example of a bacterial EF-Tu with resistance against two divergently acting antibiotics.

Bacterial diversity based on type II DNA topoisomerase genes

Type II DNA topoisomerases are essential and ubiquitous DNA metabolic enzymes that alter DNA topology. Eubacteria have two indispensable type II DNA topoisomerases, DNA gyrase encoded by gyrB and gyrA and topoisomerase IV encoded by parE and parC. These genes belong to a single family whose members span both eukaryotes and prokaryotes. The highly conserved motifs in these genes provide a rationale for the design of universal primers used in the polymerase chain reaction in order to systematically generate a data set suitable for bacterial diversity studies at the macro-diversity level, as well as at the micro-diversity level displaying individual species and isolates. This family of genes is the subject of intensive biochemical and genetic analyses, which provide an opportunity for comprehensive understanding of sequence conservation and variability and their relationship to function. These genes are ideally suited for microbial identification and biodiversity analyses.

MICROBIAL DIVERSITY: Domains and Kingdoms

▪ Abstract  With the discovery of the eukaryote nucleus, all living organisms were neatly divided into prokaryotes, which lacked a nucleus, and eukaryotes, which possessed it. As data derived directly from the genome became available, it was clear that prokaryotes were comprised of two groups, Eubacteria and Archaebacteria. These were subsequently renamed at the new taxonomic level of Domain as Bacteria and Archaea, with the eukaryotes named as the Eucarya Domain. The interrelationships of the three Domains are still subject to discussion and evaluation, as is their monophyly. Further data, drawn from various protein sequences, suggest conflicting schemes, and resolution may not be straightforward. Additionally, Bacteria and Archaea as well as Eucarya are largely based on organisms already in culture. Investigation of the potentially enormous quantity of uncultured organisms in nature is likely to have as broad-ranging implications as the exploration of new protein sequences.

The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny

The genes for the protein synthesis elongation factors Tu (EF-Tu) and G (EF-G) are the products of an ancient gene duplication, which appears to predate the divergence of all extant organismal lineages. Thus, it should be possible to root a universal phylogeny based on either protein using the second protein as an outgroup. This approach was originally taken independently with two separate gene duplication pairs, (i) the regulatory and catalytic subunits of the proton ATPases and (ii) the protein synthesis elongation factors EF-Tu and EF-G. Questions about the orthology of the ATPase genes have obscured the former results, and the elongation factor data have been criticized for inadequate taxonomic representation and alignment errors. We have expanded the latter analysis using a broad representation of taxa from all three domains of life. All phylogenetic methods used strongly place the root of the universal tree between two highly distinct groups, the archaeons/eukaryotes and the eubacteria. We also find that a combined data set of EF-Tu and EF-G sequences favors placement of the eukaryotes within the Archaea, as the sister group to the Crenarchaeota. This relationship is supported by bootstrap values of 60-89% with various distance and maximum likelihood methods, while unweighted parsimony gives 58% support for archaeal monophyly.

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.

Characterization of uncultivated prokaryotes: isolation and analysis of a 40-kilobase-pair genome fragment from a planktonic marine archaeon

One potential approach for characterizing uncultivated prokaryotes from natural assemblages involves genomic analysis of DNA fragments retrieved directly from naturally occurring microbial biomass. In this study, we sought to isolate large genomic fragments from a widely distributed and relatively abundant but as yet uncultivated group of prokaryotes, the planktonic marine Archaea. A fosmid DNA library was prepared from a marine picoplankton assemblage collected at a depth of 200 m in the eastern North Pacific. We identified a 38.5-kbp recombinant fosmid clone which contained an archaeal small subunit ribosomal DNA gene. Phylogenetic analyses of the small subunit rRNA sequence demonstrated it close relationship to that of previously described planktonic archaea, which form a coherent group rooted deeply within the Crenarchaeota branch of the domain Archaea. Random shotgun sequencing of subcloned fragments of the archaeal fosmid clone revealed several genes which bore highest similarity to archaeal homologs, including large subunit ribosomal DNA and translation elongation factor 2 (EF2). Analyses of the inferred amino acid sequence of archaeoplankton EF2 supported its affiliation with the Crenarchaeote subdivision of Archaea. Two gene fragments encoding proteins not previously found in Archaea were also identified: RNA helicase, responsible for the ATP-dependent alteration of RNA secondary structure, and glutamate semialdehyde aminotransferase, an enzyme involved in initial steps of heme biosynthesis. In total, our results indicate that genomic analysis of large DNA fragments retrieved from mixed microbial assemblages can provide useful perspective on the physiological potential of abundant but as yet uncultivated prokaryotes.

Structure-based sequence alignment of elongation factors Tu and G with related GTPases involved in translation

The G domain and domain II in the crystal structure of Thermus thermophilus elongation factor G (EF-G) were compared with the homologous domains in Thermus aquaticus elongation factor Tu (EF-Tu). Sequence alignment derived from the structural superposition was used to define conserved sequence elements in domain II. These elements and previously known conserved sequence elements in the G domain were used to guide the alignment of the sequences of Sulfolobus acidocaldarius elongation factor 2, human elongation factor 2, and Escherichia coli initiation factor 2 and release factor 3 to the aligned sequences of EF-G and EF-Tu. This alignment, which deviates from previously published alignments, has evolutionary implications and leads to alternative interpretations of biochemical data concerning the interaction of elongation factors with the alpha-sarcin/ricin region of the ribosome. A single conserved sequence motif in domain II was identified and used to further characterize the GTPase subfamily of translation factors and related proteins. It was shown that the motif is found in most if not all the members of the family. Apparently, the common characteristic of these GTPases is an extensive consensus structural unit that possibly accounts for a similar interaction with the ribosome and is composed of two domains homologous to the G domain and domain II in EF-Tu and EF-G.

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.

Root of the universal tree of life based on ancient aminoacyl-tRNA synthetase gene duplications

Universal trees based on sequences of single gene homologs cannot be rooted. Iwabe et al. [Iwabe, N., Kuma, K.-I., Hasegawa, M., Osawa, S. & Miyata, T. (1989) Proc. Natl. Acad. Sci. USA 86, 9355-9359] circumvented this problem by using ancient gene duplications that predated the last common ancestor of all living things. Their separate, reciprocally rooted gene trees for elongation factors and ATPase subunits showed Bacteria (eubacteria) as branching first from the universal tree with Archaea (archaebacteria) and Eucarya (eukaryotes) as sister groups. Given its topical importance to evolutionary biology and concerns about the appropriateness of the ATPase data set, an evaluation of the universal tree root using other ancient gene duplications is essential. In this study, we derive a rooting for the universal tree using aminoacyl-tRNA synthetase genes, an extensive multigene family whose divergence likely preceded that of prokaryotes and eukaryotes. An approximately 1600-bp conserved region was sequenced from the isoleucyl-tRNA synthetases of several species representing deep evolutionary branches of eukaryotes (Nosema locustae), Bacteria (Aquifex pyrophilus and Thermotoga maritima) and Archaea (Pyrococcus furiosus and Sulfolobus acidocaldarius). In addition, a new valyl-tRNA synthetase was characterized from the protist Trichomonas vaginalis. Different phylogenetic methods were used to generate trees of isoleucyl-tRNA synthetases rooted by valyl- and leucyl-tRNA synthetases. All isoleucyl-tRNA synthetase trees showed Archaea and Eucarya as sister groups, providing strong confirmation for the universal tree rooting reported by Iwabe et al. As well, there was strong support for the monophyly (sensu Hennig) of Archaea. The valyl-tRNA synthetase gene from Tr. vaginalis clustered with other eukaryotic ValRS genes, which may have been transferred from the mitochondrial genome to the nuclear genome, suggesting that this amitochondrial trichomonad once harbored an endosymbiotic bacterium.

Phylogenetic analysis of 70 kD heat shock protein sequences suggests a chimeric origin for the eukaryotic cell nucleus

BACKGROUND

The evolutionary relationships between archaebacteria, eubacteria and eukaryotic cells are of central importance in biology. The current view is that each of these three groups of organisms constitutes a monophyletic domain, and that eukaryotic cells have evolved fom an archaebacterial ancestor. Recent studies on a number of highly conserved protein sequences do not, however, support this view and raise important questions concerning the evolutionary relationships between all extant organisms, particularly regarding the origin of eukaryotic cells.

RESULTS

RESULTS

We have used sequences of 70 kD heat shock protein (hsp70)--the most conserved protein found to date in all species--to examine the evolutionary relationship between various species. We have obtained two new archaebacterial hsp70 sequences from the species, Thermoplasma acidophilum and Halobacterium cutirubrum. A global comparison of hsp70 sequences, including our two new sequences, shows that all known archaebacterial homologs share a number of sequence signatures with the Gram-positive group of bacteria that are not found in any other prokaryotic or eukaryotic species. In contrast, the eukaryotic homologs are shown to share a number of unique sequence features with the Gram-negative bacteria that are not present in any archaebacteria. Detailed phylogenetic analyses of hsp70 sequences strongly support a specific evolutionary relationship between archaebacteria and Gram-positive bacteria on the one hand, and Gram-negative bacteria and eukaryotes on the other. The phylogenetic analyses also indicate a polyphyletic branching of archaebacteria within the Gram-positive species. The possibility that the observed relationships are due to horizontal gene transfers can be excluded on the basis of sequence characteristics of different groups of homologs.

CONCLUSIONS

Our results do not support the view that archaebacteria constitute a monophyletic domain, but instead suggest a close evolutionary linkage between archaebacteria and Gram-positive bacteria. Furthermore, in contrast to the presently accepted view, eukaryotic hsp70s show a close and specific relationship to those from Gram-negative species. To explain the phylogenies based on different gene sequences, a chimeric model for the origin of the eukaryotic cell nucleus involving fusion between an archaebacterium and a Gram-negative eubacterium is proposed. Several predictions from the chimeric model are discussed.