Presence of elongation factor 1 in nuclei and nucleoli of rat liver

Purified rat liver nuclei contain 11% of the total cellular translation elongation factor 1 activity. Ninety percent of the nuclear EF-1 activity was in the nucleoplasm and 10% was nucleolar. The specific activities of the nuclear and nucleolar EF-1 were 2 to 3 times higher, respectively, than EF-1 activity of the liver homogenate. The presence of EF-1 in the purified nuclei did not result from cytoplasmic contamination since only 0.14% of the cellular lactate dehydrogenase was present in the nuclei. These results provide the first evidence for the presence in the cell nucleus of translational factors of protein synthesis.     

Identification of the major Mr 100,000 substrate for calmodulin-dependent protein kinase III in mammalian cells as elongation factor-2

The major substrate for Ca2+/calmodulin-dependent protein kinase III in mammalian cells is a species of Mr 100,000 that has a primarily cytoplasmic localization. This substrate has now been identified as elongation factor-2 (EF-2), a protein that catalyzes the translocation of peptidyl-tRNA on the ribosome. The amino acid sequence of 18 residues from the N-terminal of the Mr 100,000 CaM-dependent protein kinase III substrate purified from rat pancreas was found to be identical to the N-terminal sequence of authentic rat EF-2 as previously deduced from nucleic acid sequencing of a cDNA (Kohno, K., Uchida, T., Ohkubo, H., Nakanishi, S., Nakanishi, T., Fukui, T., Ohtsuka, E., Ikehara, M., and Okada, Y. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 4978-4982). CaM-dependent protein kinase III phosphorylated EF-2 in vitro with a stoichiometry of approximately 1 mol/mol on a threonine residue. Amino acid sequencing of the purified tryptic phosphopeptide revealed that this threonine residue lies within the sequence: Ala-Gly-Glu-Thr-Arg-Phe-Thr-Asp-Thr-Arg (residues 51-60 of EF-2). The Mr 100,000 protein was stoichiometrically ADP-ribosylated in vitro by the addition of diphtheria toxin and NAD. The Mr 100,000 protein was photoaffinity labeled with a GTP analog and the protein had an endogenous GTPase activity that could be stimulated by the addition of salt-washed ribosomes. These properties are all characteristic of EF-2. Dephospho-EF-2 could support poly(U)-directed polyphenylalanine synthesis in a reconstituted elongation system when combined with EF-1. In the same system, phospho-EF-2 was virtually inactive in supporting polypeptide synthesis; this effect could be reversed by dephosphorylation of phospho-EF-2. These results suggest that intracellular Ca2+ inhibits protein synthesis in mammalian cells via CaM-dependent protein kinase III-catalyzed phosphorylation of EF-2.     

Identification of the gene encoding the mitochondrial elongation factor G in mammals.

Protein synthesis in cytosolic and rough endoplasmic reticulum associated ribosomes is directed by factors, many of which have been well characterized. Although these factors have been the subject of intense study, most of the corresponding factors regulating protein synthesis in the mitochondrial ribosomes remain unknown. In this report we present the cloning and initial characterization of the gene encoding the rat mitochondrial elongation factor-G (rEF-Gmt). The rat gene encoding EF-Gmt (rMef-g) maps to rat chromosome 2 and it is expressed in all tissues with highest levels in liver, thymus and brain. Its DNA sequence predicts a 752 amino acid protein exhibiting 72% homology to the yeast Saccharomyces cerevisiae mitochondrial elongation factor-G (YMEF-G), 62% and 61% homology to the Thermus thermophilus and E. coli elongation factor-G (EF-G) respectively and 52% homology to the rat elongation factor-2 (EF-2). The deduced amino acid sequence of EF-G contains characteristic motifs shared by all GTP binding proteins. Therefore, similarly to other elongation factors, the enzymatic function of EF-Gmt is predicted to depend on GTP binding and hydrolysis. EF-Gmt differs from its cytoplasmic homolog, EF-2, in that it contains an aspartic acid residue at amino acid position 621 which corresponds to the EF-2 histidine residue at position 715. Since this histidine residue, following posttranslational modification into diphthamide, appears to be the sole cellular target of diphtheria toxin and Pseudomonas aeruginosa endotoxin A, we conclude that EF-Gmt will not be inactivated by these toxins. The severe effects of these toxins on protein elongation in tissues expressing EF-Gmt suggest that EF-Gmt and EF-2 exhibit nonoverlapping functions. The cloning and characterization of the mammalian mitochondrial elongation factor G will permit us to address its role in the regulation of normal mitochondrial function and in disease states attributed to mitochondrial dysfunction.     

Characterization of the statin-like S1 and rat elongation factor 1 alpha as two distinctly expressed messages in rat.

Previously, we reported a rat S1 protein that is antigenically related to statin, a nonproliferating cell-specific marker; however, it shares high homology with the known human elongation factor-1 alpha (EF-1 alpha). To differentiate S1 from rat EF-1 alpha and to study their respective regulation for expression, a rat EF-1 alpha cDNA clone was isolated and characterized. The nucleotide and deduced amino acid sequences of this partial rat EF-1 alpha cDNA are compared with that of human and mouse as well as with rat S1. Both their messages were detected in rat brain by EF-1 alpha- or S1-specific probes. However, the mRNA encoding EF-1 alpha is more abundant than that encoding S1. S1 and EF-1 alpha expression were investigated in the parotid and submandibular glands of untreated rats and those treated with isoproterenol, a proliferation-inducing catecholamine. Quantitative solution hybridization demonstrated a dramatic reduction (approximately 68%) in the S1 mRNA following isoproterenol injection in proliferation-responsive parotid glands and a mild reduction (approximately 20%) of S1 steady-state messages in the proliferation-refractile submandibular glands. A slight increase or no changes of EF-1 alpha levels in both parotid and submandibular glands following isoproterenol treatment are also observed. Therefore, the EF-1 alpha and S1 genes are different genes, both expressed and regulated in vivo, but in differential quantitative and qualitative patterns.     

Cloning and expression of Bombyx mori silk gland elongation factor 1gamma in Escherichia coli

Elongation factor 1 (EF-1) from the silk gland of Bombyx mori consists of alpha-, beta-, gamma-, and delta-subunits. EF-1alpha GTP catalyzes the binding of aminoacyl-tRNA to ribosomes concomitant with the hydrolysis of GTP. EF-1betagammadelta catalyzes the exchange of EF-1alpha-bound GDP for exogenous GTP and stimulates the EF-1alpha-dependent binding of aminoacyl-tRNA to ribosomes. EF-1gamma cDNA, which contains an open reading frame (ORF) encoding a polypeptide of 423 amino acid residues, was amplified and cloned by PCR from a silk gland cDNA library. The calculated molecular mass and predicted pI of the product were 48,388 Da and 5.84, respectively. The silk gland EF-1gamma shares 67.3% amino acid identity with Artemia salina EF-lgamma. The N-terminal domain (amino acid residues 1-211) of silk gland EF-lgamma is 29.3% identical to maize glutathione S-transferase. We demonstrated that silk gland EF-lgamma bound to glutathione Sepharose, suggesting that the N-terminal domain of EF-1gamma may have the capacity to bind to glutathione.     

Protein synthesis in yeast. Purification of elongation factor 3 from temperature-sensitive mutant 13-06 of the yeast Saccharomyces cerevisiae

An altered form of the elongation factor 3 (EF-3) has been purified to near homogeneity from a thermolabile yeast mutant ts 13-06. The isolation procedure involved chromatography on DEAE-Sephadex, CM-Sepharose, and hydroxylapatite columns. The final purification of this protein was obtained by affinity chromatography on an ATP-Sepharose column. Because of the extreme lability of the mutant protein, the yield was very poor. Silver stain analysis of the sodium dodecyl sulfate electrophoretograms indicated that the affinity-purified protein was better than 90% pure. From the studies of the physical and biochemical properties, the following characteristics of the purified wild type and the mutant protein have been established. The two proteins were indistinguishable by their molecular weight, amino acid composition, and isoelectric point. Purified mutant EF-3 was rapidly inactivated between 37 and 39 degrees C. Under this condition, wild type EF-3 was completely stable. Ribosome-dependent GTPase and ATPase activities of the mutant EF-3 were heat sensitive; GTPase activity was more labile than the ATPase activity. Mutant EF-3, after exposure to a nonpermissive temperature, failed to stimulate binding of the ternary complex of EF-1 X GTP X aminoacyl-tRNA to ribosome. The wild type protein was fully active under this condition. Other biochemical and physical properties of these two proteins are under current investigation.     

Molecular cloning,biochemical and structural analysis of elongation factor-1 alpha from Leishmania donovani: comparison with the mammalian homologue

The Src-homology 2 domain containing protein tyrosine phosphatase-1 (SHP-1) is involved in the pathogenesis of infection with Leishmania. Recently, we identified elongation factor-1 alpha (EF-1 alpha) from Leishmania donovani as a SHP-1 binding and activating protein [J. Biol. Chem. 277 (2002) 50190]. To characterize this apparent Leishmania virulence factor further, the cDNA encoding L. donovani EF-1 alpha was cloned and sequenced. Whereas nearly complete sequence conservation was observed amongst EF-1 alpha proteins from trypanosomatids, the deduced amino acid sequence of EF-1 alpha of L. donovani when compared to mammalian EF-1 alpha sequences showed a number of significant changes. Protein structure modeling-based upon the known crystal structure of EF-1 alpha for Saccharomyces cerevisiae-identified a hairpin loop present in mammalian EF-1 alpha and absent from the Leishmania protein which corresponded to a 12 amino acid deletion. Consistent with these structural differences, the sub-cellular distributions of L. donovani EF-1 alpha and host EF-1 alpha were strikingly different. Interestingly, infection of macrophages with L. donovani caused redistribution of host as well as pathogen EF-1 alpha. Since EF-1 alpha is essential for survival, the distinct biochemical and structural properties of Leishmania EF-1 alpha may provide a novel target for drug development.     

Human elongation factor 1 beta: cDNA and derived amino acid sequence.

From a cDNA library in lambda gt11 derived from poly (A+)RNA of human ovarian granulosa cells a cDNA clone lambda HGP34, containing an EcoRI insert of 829 bp, was identified. After subcloning of the insert into pUC18, the clone pHGP34 was obtained and sequenced. The derived amino acid sequence, corresponding to a protein of 225 amino acids, shows a high degree of homology to elongation factor 1 beta (EF-1 beta) of Artemia salina (57%) and known peptide sequences of Xenopus laevis EF-1 beta (86%). We therefore assume that the protein coded for by pHGP34 represents human EF-1 beta. Northern analysis reveals an EF-1 beta specific mRNA of 900 bp. Southern analysis indicates that EF-1 beta in the human genome, like EF-1 alpha, appears to be specified by more than one gene. A high degree of sequence homology for EF-1 beta specific sequences is observed for bovine, rat and mouse species.     

Molecular cloning of the black tiger shrimp (Penaeus monodon) elongation factor 2 (EF-2): sequence analysis and its expression on the ovarian maturation stage

The techniques of homology cloning and anchored PCR were used to clone the elongation factor 2 (EF-2) gene from black tiger shrimp (Penaeus monodon). The full length cDNA of black tiger shrimp EF-2 (btsEF-2) contained a 5' untranslated region (UTR) of 73 bp, an ORF of 2541 bp encoding a polypeptide of 846 amino acids with an estimated molecular mass of 95 kDa, and a 3( UTR of 112 bp. The searches for protein sequence similarities with BLAST analysis indicated that the deduced amino acid sequence of btsEF-2 was homological to the EF-2 of other species and even the mammalians. The conserved signature sequence of EF-2 gene family, GTPase effector domain and ADP-ribosylation domain were found in the btsEF-2 deduced amino acid sequence. The temporal expressions of gene in the different ovarian stages were measured by real time PCR. The mRNA expressions of the gene were constitutively expressed in ovary and different during the maturation stages. The result indicated that EF-2 gene was constitutively expressed and could play a critical role in the ovarian maturation stage.     

The primary structure and the functional domains of an elongation factor-1 alpha from Mucor racemosus

We have determined the complete nucleotide sequence for TEF-1, one of three genes coding for elongation factor (EF)-1 alpha in Mucor racemosus. The deduced EF-1 alpha protein contains 458 amino acids encoded by two exons. The presence of an intervening sequence located near the 3' end of the gene was predicted by the nucleotide sequence data and confirmed by alkaline S1 nuclease mapping. The amino acid sequence of EF-1 alpha was compared to the published amino acid sequences of EF-1 alpha proteins from Saccharomyces cerevisiae and Artemia salina. These proteins shared nearly 85% homology. A similar comparison to the functionally analogous EF-Tu from Escherichia coli revealed several regions of amino acid homology suggesting that the functional domains are conserved in elongation factors from these diverse organisms. Secondary structure predictions indicated that alpha helix and beta sheet conformations associated with the functional domains in EF-Tu are present in the same relative location in EF-1 alpha from M. racemosus. Through this comparative structural analysis we have predicted the general location of functional domains in EF-1 alpha which interact with GTP and tRNA.     

Isolation and characterization of eft-1, an elongation factor 2-like gene on chromosome III of Caenorhabditis elegans.

A gene (eft-1) encoding an elongation factor 2-like protein was isolated from a region adjacent to the polyubiquitin gene, ubq-1, of Caenorhabditis elegans. Sequence analysis of genomic and cDNA clones revealed that the deduced amino acid sequence of the protein (EFT-1) is 38% identical to that of mammalian and Drosophila elongation factor 2 (EF-2). The entire eft-1 gene is approximately 3.8 kb in length and contains 5 exons separated by short introns of 46-75 bp. The 2,547-bp open reading frame predicts a protein of 849 amino acid residues (calculated Mr, 96,151). Conserved sequences shared among a variety of GTP-binding proteins including EF-2 are found in the amino-terminal third of EFT-1. The carboxy-terminal half contains regions with 40-57% similarity (including conservative changes) with segments characteristic of EF-2 and its prokaryotic homolog, EF-G. However, the histidyl residue target for ADP-ribosylation of EF-2 by diphtheria toxin is replaced by tyrosine in EFT-1. Southern and Northern blot analyses indicate that eft-1 is a single-copy gene that is expressed at all stages of nematode development. Amplification of fragments encoding highly conserved regions of EF-2 using the polymerase chain reaction led to the isolation of a fragment encoding the modifiable histidyl residue and which likely represents part of the C. elegans EF-2 gene (eft-2). This suggests that EFT-1 is not the C. elegans homolog of EF-2, but a closely related protein.     

Structure and expression of elongation factor 2 gene during development of Dictyostelium discoideum

A cDNA library constructed from poly(A)+ RNA isolated from Dictyostelium discoideum cells at 12 h of development was screened with the hamster elongation factor 2 (EF-2) cDNA. Several different cDNA clones which hybridized were isolated after a second screening. A cDNA clone representing the 5'-end of the mRNA was obtained by primer extension. By comparing the amino acid sequence deduced from the nucleotide sequences of these clones with that of hamster EF-2, we found enough homology between them to conclude that the isolated clones were complementary to the mRNA of D. discoideum EF-2. The N terminus which is the GTP-binding domain and the C-terminal half where it interacts with a ribosome showed a high degree of homology. The amino acid sequence of the carboxyl half includes that it contain a site of ADP-ribosylation by diphtheria toxin. From the Northern blotting analysis, the size of the mRNA was estimated to be 2.6 kilobases. The expression of the mRNA was high in vegetative cells, became maximal at the aggregation stage, and decreased thereafter through development. Upon differentiation of prespore and prestalk cells, the mRNA was highly enriched in the former over the latter. ADP-ribosylation assay of EF-2 protein by diphtheria toxin showed nearly the same developmental changes for the protein as the mRNA. However, prestalk cells were found to contain the same amount of the protein as prespore cells. The Southern blot analyses indicated that the gene encoding EF-2 is unique.     

Reduced ribosomal binding of eukaryotic elongation factor 2 following ADP-ribosylation. Difference in binding selectivity between polyribosomes and reconstituted monoribosomes

The biological activity of elongation factor 2 (EF-2) following NAD+ - and diphtheria-toxin-dependent ADP-ribosylation was studied (i) in translation experiments using the reticulocyte lysate system and (ii) in ribosomal binding experiments using either reconstituted empty rat liver ribosomes or programmed reticulocyte polysomes. Treatment of the lysates with toxin and NAD+ at a NAD+/ribosome ratio of 4 resulted in a 90% inhibition of the amino acid incorporation rate. The inhibition was overcome by the addition of native EF-2. At this level of inhibition more than 90% of the EF-2 present in the lysates was ADP-ribosylated and the total ribosome association of EF-2 was reduced by approx. 50%. All of the remaining unmodified factor molecules were associated with the ribosomes, whereas only about 3% of the ribosylated factor was ribosome-associated. The nucleotide requirement for the binding of EF-2 to empty reconstituted rat liver ribosomes and programmed reticulocyte polysomes was studied together with the stability of the resulting EF-2 X ribosome complexes using purified 125I-labelled rat liver EF-2. With both types of ribosomes, the complex formation was strictly nucleotide-dependent. Stable, high-affinity complexes were formed in the presence of the non-hydrolysable GTP analogue guanosine 5'-(beta, gamma-methylene)triphosphate (GuoPP[CH2]P). In contrast to the reconstituted ribosomes, GTP stimulated the formation of high-affinity complexes in the presence of polysomes, albeit at a lower efficiency than GuoPP[CH2]P. The formation of high-affinity complexes was restricted to polysomes in the pretranslocation phase of the elongation cycle. Low-affinity post-translocation complexes, demonstrable after fixation, were formed in the presence of GTP, GuoPP[CH2]P and GDP. In polysomes, these complexes involved a different population of particles than did the high-affinity complexes. In the binding experiments using reconstituted or programmed ribosomes, the pretranslocation binding of EF-2 observed in the presence of GuoPP[CH2]P was reduced by approx. 50% after ADP-ribosylation, whereas the post-translocation binding in the presence of GDP was unaltered. The data indicate that the inhibition of translocation caused by diphtheria toxin and NAD+ is mediated through a reduced affinity of the ADP-ribosylated EF-2 for binding to ribosomes in the pretranslocation state.     

Molecular cloning and characterization of the Caenorhabditis elegans elongation factor 2 gene (eft-2)

A Caenorhabditis elegans lambda ZAP cDNA library was screened using a fragment amplified from highly conserved regions of the mammalian and Drosophila elongation factor 2 (EF-2). Two types of cDNA clones were obtained, corresponding to two mRNA species with 3'-untranslated regions of 60 and 115 nucleotides, both encoding identical polypeptides. Sequence analysis of these clones and comparisons with hamster and Drosophila EF-2 sequences suggests that they encode C. elegans EF-2. Clone pCef6A, encoding the entire C. elegans EF-2 mRNA sequence including 45 nucleotides of 5'-untranslated region, contains a 2,556-bp open reading frame which predicts a polypeptide of 852 amino acid residues (Mr 94,564). The deduced amino acid sequence is greater than 80% identical to that of mammalian and Drosophila EF-2. Conserved sequence segments shared among a variety of GTP-binding proteins are found in the amino-terminal region. The carboxy-terminal half contains segments unique to EF-2 and its prokaryotic homolog, EF-G, as well as the histidyl residue which is ADP-ribosylated by diphtheria toxin. The C. elegans protein contains a 12-amino-acid insertion between positions 90 and 100, and a 13-amino-acid deletion between positions 237 and 260, relative to hamster EF-2. Partial sequencing of a genomic clone encoding the entire C. elegans EF-2 gene (named eft-2) has so far revealed two introns of 48 and 44 bp following codons Gln-191 and Gln-250, respectively. Southern and Northern blot analyses indicate that eft-2 is a single-copy gene and encodes a 3-kb mRNA species which is present throughout nematode development.     

Modification of the reactivity of three amino-acid residues in elongation factor 2 during its binding to ribosomes and translocation

The accessibility of three amino acids of EF-2, located within highly conserved regions near the N- and C-terminal extremities of the molecule (the E region and the ADPR region, respectively) to modifying enzymes has been compared within nucleotide-complexed EF-2 and ribosomal complexes that mimic the pre- and posttranslocational ones: the high-affinity complex (EF-2)-nonhydrolysable GTP analog GuoPP[CH2]P ribosome and the low-affinity (EF-2)-GDP-ribosome complex, EF-2 and ribosomes being from rat liver. We studied the reactivity of two highly conserved residues diphthamide-715 and Arg-66, to diphtheria-toxin-dependent ADP-ribosylation and trypsin attack, and of a threonine that probably lies between residues 51 and 60, to phosphorylation by a Ca2+/calmodulin-dependent protein kinase. Diphthamide 715 and this threonine residue were unreactive within the high-affinity complex but seemed fully reactive in the low-affinity complex. Arg-66 was resistant to trypsin in both complexes. The possible involvement of the E and ADPR regions of EF-2 in the interaction with ribosome in the two complexes is discussed.     

Characterization of the endogenous ADP-ribosylation of wild-type and mutant elongation factor 2 in eukaryotic cells.

Anti-[ADP-ribosylated elongation factor 2 (EF-2)] antiserum has been used to immunoprecipitate the modified form of EF-2 from polyoma-virus-transformed baby hamster kidney (pyBHK) cells [Fendrick, J. L. & Iglewski, W. J. (1989) Proc. Natl Acad. Sci. USA 86, 554-557]. This antiserum also immunoprecipitates a 32P-labelled protein of similar size to EF-2 from a variety of primary and continuous cell lines derived from many species of animals. One of these cell lines, chinese hamster ovary CHO-K1 cells was further characterized. The time course of labelling of ADP-ribosylated EF-2 with [32P]orthophosphate was similar in pyBHK cells and in CHO-K1 cells. The kinetics of labelling were more rapid for cells cultured in 2% serum than 10% serum, with incorporation of 32P reaching a maximum at 6 h and 10 h, respectively. EF-2 mutants of pyBHK and CHO-K1 cells resistant to diphtheria-toxin-catalyzed ADP-ribosylation of EF-2 remain sensitive to cellular ADP-ribosylation of EF-2. The 32P-labelled moiety of ADP-ribosylated EF-2 was digested by snake venom phosphodiesterase and the product was identified as AMP. The same 32P-labelled tryptic peptide was modified by toxin in wild-type EF-2 and by the cellular transferase in mutant EF-2. When purified EF-2 from pyBHK cells was incubated with [carbonyl-14C]nicotinamide and diphtheria toxin fragment A, under conditions for reversal of the ADP-ribosylation reaction, [14C]NAD was generated. The results suggest that cellular ADP-ribosylated EF-2 exists in a variety of cell types, and the ribosylated product is identical to that produced by toxin ADP-ribosylation of EF-2, except in diphthamide mutant cells. Studies with the mutant cell lines indicate that the toxin and the cellular transferase, however, recognize different determinants at the ADP-ribose acceptor site in EF-2. The cellular transferase does not require the diphthamide modification of the histidine ring in the amino acid sequence of EF-2 for the transfer of ADP-ribose to the ring. Therefore, we would expect the cellular transferase active site to be similar to, but not identical to, the critical amino acids demonstrated in the active site of diphtheria toxin and Pseudomonas exotoxin A.     

Phosphorylation of the elongation factor 2: the fifth Ca2+/calmodulin-dependent system of protein phosphorylation

Elongation factor 2 (EF-2) has been recently shown to be extensively phosphorylated in a Ca2+/calmodulin-dependent manner in extracts of mammalian cells (A. G. Ryazanov (1987) FEBS Lett. 214, 331-334). In the present study, we partially purified the protein kinase which phosphorylates EF-2 from rabbit reticulocytes. The molecular weight of the enzyme determined by gel filtration was about 140,000. Unlike the substrate, the EF-2 kinase had no affinity for RNA and therefore could be separated from EF-2 by chromatography on RNA-Sepharose. After chromatography on hydroxyapatite, the kinase activity became calmodulin-dependent. Two-dimensional separation of the phosphorylated EF-2 according to O'Farrell's technique revealed that there were two phosphorylation sites within the EF-2 molecule; in both cases, the phosphorylated amino acid was threonine. The EF-2 kinase differed from the four known types of Ca2+/calmodulin-dependent protein kinases. Thus, the system of EF-2 phosphorylation represents the novel (fifth) Ca2+/calmodulin-dependent system of protein phosphorylation. This system is supposed to be responsible for the regulation of the elongation rate of protein biosynthesis in eukaryotic cells.     

Phylogenetic place of mitochondrion-lacking protozoan, Giardia lamblia, inferred from amino acid sequences of elongation factor 2

Partial regions of the mRNA encoding a major part of translation elongation factor 2 (EF-2) from a mitochondrion-lacking protozoan, Giardia lamblia, were amplified by polymerase chain reaction, and their primary structures were analyzed. The deduced amino acid sequence was aligned with other eukaryotic and archaebacterial EF-2's, and the phylogenetic relationships among eukaryotes were inferred by the maximum likelihood (ML) and the maximum parsimony (MP) methods. The ML analyses using six different models of amino acid substitutions and the MP analysis consistently suggest that among eukaryotic species being analyzed, G. lamblia is likely to have diverged from other higher eukaryotes on the early phase of eukaryotic evolution.