Protein phylogeny gives a robust estimation for early divergences of eukaryotes: phylogenetic place of a mitochondria-lacking protozoan, Giardia lamblia

Evaluation on Elongation Factor 1 Alpha of Entamoeba histolytica Interaction with the Intermediate Subunit of the Gal/GalNAc Lectin and Actin in Phagocytosis

Entamoeba histolytica is the causative agent of amoebiasis. This disease results in 40,000 to 100,000 deaths annually. The pathogenic molecules involved in the invasion of trophozoites had been constantly being clarified. This study explored the role of elongation factor 1 alpha (EF1a) in E. histolytica pathogenicity. Biolayer interferometry binding and pull-down assays suggest that EF1a and intermediate subunit of lectin (Igl) binding are specific. Submembranous distribution of EF1a closely aligns with the localization of Igl, which appear in abundance on membranes of trophozoites. Messenger RNA (mRNA) expression of EF1a is positively correlated with trends in Igl levels after co-incubation with Chinese hamster ovary (CHO) cells in vitro, suggesting a regulatory linkage between these proteins. Erythrophagocytosis assays also imply a role for EF1a in phagocytosis. Finally, EF1a and actin are collocated in trophozoites. These results indicated elongation factor 1a is associated with E. histolytica phagocytosis, and the relationships between EF1a, Igl, and actin are worth further study to better understand the pathogenic process.

Identification and evolutionary analysis of the nucleolar proteome of Giardia lamblia

Background

The nucleoli, including their proteomes, of higher eukaryotes have been extensively studied, while few studies about the nucleoli of the lower eukaryotes – protists were reported. Giardia lamblia, a protist with the controversy of whether it is an extreme primitive eukaryote or just a highly evolved parasite, might be an interesting object for carrying out the nucleolar proteome study of protists and for further examining the controversy.

Results

Using bioinformatics methods, we reconstructed G. lamblia nucleolar proteome (GiNuP) and the common nucleolar proteome of the three representative higher eukaryotes (human, Arabidopsis, yeast) (HEBNuP). Comparisons of the two proteomes revealed that: 1) GiNuP is much smaller than HEBNuP, but 78.4% of its proteins have orthologs in the latter; 2) More than 68% of the GiNuP proteins are involved in the “Ribosome related” function, and the others participate in the other functions, and these two groups of proteins are much larger and much smaller than those in HEBNuP, respectively; 3) Both GiNuP and HEBNuP have their own specific proteins, but HEBNuP has a much higher proportion of such proteins to participate in more categories of nucleolar functions.

Conclusion

For the first time the nucleolar proteome of a protist - Giardia was reconstructed. The results of comparison of it with the common proteome of three representative higher eukaryotes -- HEBNuP indicated that the simplicity of GiNuP is most probably a reflection of primitiveness but not just parasitic reduction of Giardia, and simultaneously revealed some interesting evolutionary phenomena about the nucleolus and even the eukaryotic cell, compositionally and functionally.

The impact of history on our perception of evolutionary events: endosymbiosis and the origin of eukaryotic complexity

Evolutionary hypotheses are correctly interpreted as products of the data they set out to explain, but they are less often recognized as being heavily influenced by other factors. One of these is the history of preceding thought, and here I look back on historically important changes in our thinking about the role of endosymbiosis in the origin of eukaryotic cells. Specifically, the modern emphasis on endosymbiotic explanations for numerous eukaryotic features, including the cell itself (the so-called chimeric hypotheses), can be seen not only as resulting from the advent of molecular and genomic data, but also from the intellectual acceptance of the endosymbiotic origin of mitochondria and plastids. This transformative idea may have unduly affected how other aspects of the eukaryotic cell are explained, in effect priming us to accept endosymbiotic explanations for endogenous processes. Molecular and genomic data, which were originally harnessed to answer questions about cell evolution, now so dominate our thinking that they largely define the question, and the original questions about how eukaryotic cellular architecture evolved have been neglected. This is unfortunate because, as Roger Stanier pointed out, these cellular changes represent life's "greatest single evolutionary discontinuity," and on this basis I advocate a return to emphasizing evolutionary cell biology when thinking about the origin of eukaryotes, and suggest that endogenous explanations will prevail when we refocus on the evolution of the cell.

Putative SF2 helicases of the early-branching eukaryote Giardia lamblia are involved in antigenic variation and parasite differentiation into cysts

Background

Regulation of surface antigenic variation in Giardia lamblia is controlled post-transcriptionally by an RNA-interference (RNAi) pathway that includes a Dicer-like bidentate RNase III (gDicer). This enzyme, however, lacks the RNA helicase domain present in Dicer enzymes from higher eukaryotes. The participation of several RNA helicases in practically all organisms in which RNAi was studied suggests that RNA helicases are potentially involved in antigenic variation, as well as during Giardia differentiation into cysts.

Results

An extensive in silico analysis of the Giardia genome identified 32 putative Super Family 2 RNA helicases that contain almost all the conserved RNA helicase motifs. Phylogenetic studies and sequence analysis separated them into 22 DEAD-box, 6 DEAH-box and 4 Ski2p-box RNA helicases, some of which are homologs of well-characterized helicases from higher organisms. No Giardia putative helicase was found to have significant homology to the RNA helicase domain of Dicer enzymes. Additionally a series of up- and down-regulated putative RNA helicases were found during encystation and antigenic variation by qPCR experiments. Finally, we were able to recognize 14 additional putative helicases from three different families (RecQ family, Swi2/Snf2 and Rad3 family) that could be considered DNA helicases.

Conclusions

This is the first comprehensive analysis of the Super Family 2 helicases from the human intestinal parasite G. lamblia. The relative and variable expression of particular RNA helicases during both antigenic variation and encystation agrees with the proposed participation of these enzymes during both adaptive processes. The putatives RNA and DNA helicases identified in this early-branching eukaryote provide initial information regarding the biological role of these enzymes in cell adaptation and differentiation.

Comparative analysis of the 5S rRNA and its associated proteins reveals unique primitive rather than parasitic features in Giardia lamblia

BACKGROUND

5S rRNA is a highly conserved ribosomal component. Eukaryotic 5S rRNA and its associated proteins (5S rRNA system) have become very well understood. Giardia lamblia was thought by some researchers to be the most primitive extant eukaryote while others considered it a highly evolved parasite. Previous reports have indicated that some aspects of its 5S rRNA system are simpler than that of common eukaryotes. We here explore whether this is true to its entire system, and whether this simplicity is a primitive or parasitic feature.

METHODOLOGY/PRINCIPAL FINDINGS

By collecting and confirming pre-existing data and identifying new data, we obtained almost complete datasets of the system of three isolates of G. lamblia, two other parasitic excavates (Trichomonas vaginalis, Trypanosoma cruzi), and one free-living one (Naegleria gruberi). After comprehensively comparing each aspect of the system among these excavates and also with those of archaea and common eukaryotes, we found all the three Giardia isolates to harbor a same simplified 5S rRNA system, which is not only much simpler than that of common eukaryotes but also the simplest one among those of these excavates, and is surprisingly very similar to that of archaea; we also found among these excavates the system in parasitic species is not necessarily simpler than that in free-living species, conversely, the system of free-living species is even simpler in some respects than those of parasitic ones.

CONCLUSION/SIGNIFICANCE

The simplicity of Giardia 5S rRNA system should be considered a primitive rather than parasitically-degenerated feature. Therefore, Giardia 5S rRNA system might be a primitive system that is intermediate between that of archaea and the common eukaryotic model system, and it may reflect the evolutionary history of the eukaryotic 5S rRNA system from the archaeal form. Our results also imply G. lamblia might be a primitive eukaryote with secondary parasitically-degenerated features.

Physiological uncoupling of mitochondrial oxidative phosphorylation. Studies in different yeast species

Under non-phosphorylating conditions a high proton transmembrane gradient inhibits the rate of oxygen consumption mediated by the mitochondrial respiratory chain (state IV). Slow electron transit leads to production of reactive oxygen species (ROS) capable of participating in deleterious side reactions. In order to avoid overproducing ROS, mitochondria maintain a high rate of O(2) consumption by activating different exquisitely controlled uncoupling pathways. Different yeast species possess one or more uncoupling systems that work through one of two possible mechanisms: i) Proton sinks and ii) Non-pumping redox enzymes. Proton sinks are exemplified by mitochondrial unspecific channels (MUC) and by uncoupling proteins (UCP). Saccharomyces. cerevisiae and Debaryomyces hansenii express highly regulated MUCs. Also, a UCP was described in Yarrowia lipolytica which promotes uncoupled O(2) consumption. Non-pumping alternative oxido-reductases may substitute for a pump, as in S. cerevisiae or may coexist with a complete set of pumps as in the branched respiratory chains from Y. lipolytica or D. hansenii. In addition, pumps may suffer intrinsic uncoupling (slipping). Promising models for study are unicellular parasites which can turn off their aerobic metabolism completely. The variety of energy dissipating systems in eukaryote species is probably designed to control ROS production in the different environments where each species lives.

Evolution and diversity of the Golgi

The Golgi is an ancient and fundamental eukaryotic organelle. Evolutionary cell biological studies have begun establishing the repertoire, processes, and level of complexity of membrane-trafficking machinery present in early eukaryotic cells. This article serves as a review of the literature on the topic of Golgi evolution and diversity and reports a novel comparative genomic survey addressing Golgi machinery in the widest taxonomic diversity of eukaryotes sampled to date. Finally, the article is meant to serve as a primer on the rationale and design of evolutionary cell biological studies, hopefully encouraging readers to consider this approach as an addition to their cell biological toolbox. It is clear that the major machinery involved in vesicle trafficking to and from the Golgi was already in place by the time of the divergence of the major eukaryotic lineages, nearly 2 billion years ago. Much of this complexity was likely generated by an evolutionary process involving gene duplication and coevolution of specificity encoding membrane-trafficking proteins. There have also been clear cases of loss of Golgi machinery in some lineages as well as innovation of novel machinery. The Golgi is a wonderfully complex and diverse organelle and its continued exploration promises insight into the evolutionary history of the eukaryotic cell.

Upregulation of a novel eukaryotic translation initiation factor 5A (eIF5A) in dengue 2 virus-infected mosquito cells

BACKGROUND

Dengue virus, a mosquito-borne flavivirus, is the etiological agent of dengue fever, dengue hemorrhagic fever, and dengue shock syndrome. It generally induces apoptosis in mammalian cells, but frequently results in persistent infection in mosquito cells. That mechanism remains to be explored. In turn, a genomic survey through subtractive hybridization (PCR-select cDNA subtraction) was conducted in order to find gene(s) that may play a role in interactions between the virus and its host cells.

RESULTS

Through this technique, we identified a novel eukaryotic translation initiation factor 5A (eIF5A) which is upregulated in Aedes albopictus-derived C6/36 cells infected by the type 2 dengue (Den-2) virus. The full-length of the identified eIF5A gene consisted of 1498 bp of nucleotides with a 41.39% G+C content, and it possessed a higher similarity and shorter evolutionary distance with insects than with other organisms. Upregulation of eIF5A in response to Den-2 virus infection was validated at both the RNA and protein levels. This phenomenon was also observed by confocal microscopy. In addition, cell death obviously occurred when eIF5A activity was inhibited in C6/36 cells even when they were infected by the virus. However, viral multiplication was not obviously affected in infected C6/36 cells when eIF5A activity was reduced.

CONCLUSIONS

Taken together, we postulated that eIF5A plays a role in preventing mosquito cells from death in response to Den-2 viral infection, thus facilitating continued viral growth and potential persistent infection in mosquito cells. It would be worthwhile to further investigate how its downstream factors or cofactors contribute to this phenomenon of dengue infection.

A contiguous compartment functions as endoplasmic reticulum and endosome/lysosome in Giardia lamblia

The dynamic evolution of organelle compartmentalization in eukaryotes and how strictly compartmentalization is maintained are matters of ongoing debate. While the endoplasmic reticulum (ER) is classically envisioned as the site of protein cotranslational translocation, it has recently been proposed to have pluripotent functions. Using transfected reporter constructs, organelle-specific markers, and functional enzyme assays, we now show that in an early-diverging protozoan, Giardia lamblia, endocytosis and subsequent degradation of exogenous proteins occur in the ER or in an adjacent and communicating compartment. The Giardia endomembrane system is simple compared to those of typical eukaryotes. It lacks peroxisomes, a classical Golgi apparatus, and canonical lysosomes. Giardia orthologues of mammalian lysosomal proteases function within an ER-like tubulovesicular compartment, which itself can dynamically communicate with clathrin-containing vacuoles at the periphery of the cell to receive endocytosed proteins. These primitive characteristics support Giardia's proposed early branching and could serve as a model to study the compartmentalization of endocytic and lysosomal functions into organelles distinct from the ER. This system also may have functional similarity to the retrograde transport of toxins and major histocompatibility complex class I function in the ER of mammals.

Nucleotide substitution bias within the genus Drosophila affects the pattern of proteome evolution

The availability of complete genome sequences for 12 Drosophila species provides an unprecedented resource for large-scale studies of genome evolution. In this study, we looked for correlated shifts in the patterns of genome and proteome evolution within the genus Drosophila. Specifically, we asked if the nucleotide composition of the Drosophila willistoni genome--which is significantly less GC rich than the other 11 sequenced Drosophila genomes--is reflected in an altered pattern of amino acid substitutions in the encoded proteins. Our results show that this is indeed the case: There are large and highly significant asymmetries in the patterns of amino acid substitution between D. willistoni and Drosophila melanogaster, and they are in the direction predicted by the nucleotide biases. The implication of this result, combined with previous studies on long-term proteome evolution, is that substitutional biases at the DNA level can be a major factor in determining both the long-term and the short-term directions of proteome evolution.

Phylogenomic analyses support the monophyly of Excavata and resolve relationships among eukaryotic "supergroups"

Nearly all of eukaryotic diversity has been classified into 6 suprakingdom-level groups (supergroups) based on molecular and morphological/cell-biological evidence; these are Opisthokonta, Amoebozoa, Archaeplastida, Rhizaria, Chromalveolata, and Excavata. However, molecular phylogeny has not provided clear evidence that either Chromalveolata or Excavata is monophyletic, nor has it resolved the relationships among the supergroups. To establish the affinities of Excavata, which contains parasites of global importance and organisms regarded previously as primitive eukaryotes, we conducted a phylogenomic analysis of a dataset of 143 proteins and 48 taxa, including 19 excavates. Previous phylogenomic studies have not included all major subgroups of Excavata, and thus have not definitively addressed their interrelationships. The enigmatic flagellate Andalucia is sister to typical jakobids. Jakobids (including Andalucia), Euglenozoa and Heterolobosea form a major clade that we name Discoba. Analyses of the complete dataset group Discoba with the mitochondrion-lacking excavates or "metamonads" (diplomonads, parabasalids, and Preaxostyla), but not with the final excavate group, Malawimonas. This separation likely results from a long-branch attraction artifact. Gradual removal of rapidly-evolving taxa from the dataset leads to moderate bootstrap support (69%) for the monophyly of all Excavata, and 90% support once all metamonads are removed. Most importantly, Excavata robustly emerges between unikonts (Amoebozoa + Opisthokonta) and "megagrouping" of Archaeplastida, Rhizaria, and chromalveolates. Our analyses indicate that Excavata forms a monophyletic suprakingdom-level group that is one of the 3 primary divisions within eukaryotes, along with unikonts and a megagroup of Archaeplastida, Rhizaria, and the chromalveolate lineages.

Gain and loss of elongation factor genes in green algae

Background

Two key genes of the translational apparatus, elongation factor-1 alpha (EF-1α) and elongation factor-like (EFL) have an almost mutually exclusive distribution in eukaryotes. In the green plant lineage, the Chlorophyta encode EFL except Acetabularia where EF-1α is found, and the Streptophyta possess EF-1α except Mesostigma, which has EFL. These results raise questions about evolutionary patterns of gain and loss of EF-1α and EFL. A previous study launched the hypothesis that EF-1α was the primitive state and that EFL was gained once in the ancestor of the green plants, followed by differential loss of EF-1α or EFL in the principal clades of the Viridiplantae. In order to gain more insight in the distribution of EF-1α and EFL in green plants and test this hypothesis we screened the presence of the genes in a large sample of green algae and analyzed their gain-loss dynamics in a maximum likelihood framework using continuous-time Markov models.

Results

Within the Chlorophyta, EF-1α is shown to be present in three ulvophycean orders (i.e., Dasycladales, Bryopsidales, Siphonocladales) and the genus Ignatius. Models describing gene gain-loss dynamics revealed that the presence of EF-1α, EFL or both genes along the backbone of the green plant phylogeny is highly uncertain due to sensitivity to branch lengths and lack of prior knowledge about ancestral states or rates of gene gain and loss. Model refinements based on insights gained from the EF-1α phylogeny reduce uncertainty but still imply several equally likely possibilities: a primitive EF-1α state with multiple independent EFL gains or coexistence of both genes in the ancestor of the Viridiplantae or Chlorophyta followed by differential loss of one or the other gene in the various lineages.

Conclusion

EF-1α is much more common among green algae than previously thought. The mutually exclusive distribution of EF-1α and EFL is confirmed in a large sample of green plants. Hypotheses about the gain-loss dynamics of elongation factor genes are hard to test analytically due to a relatively flat likelihood surface, even if prior knowledge is incorporated. Phylogenetic analysis of EFL genes indicates misinterpretations in the recent literature due to uncertainty regarding the root position.

UPF1, a conserved nonsense-mediated mRNA decay factor, regulates cyst wall protein transcripts in Giardia lamblia

The Giardia lamblia cyst wall is required for survival outside the host and infection. Three cyst wall protein (cwp) genes identified to date are highly up-regulated during encystation. However, little is known of the molecular mechanisms governing their gene regulation. Messenger RNAs containing premature stop codons are rapidly degraded by a nonsense-mediated mRNA decay (NMD) system to avoid production of non-functional proteins. In addition to RNA surveillance, NMD also regulates thousands of naturally occurring transcripts through a variety of mechanisms. It is interesting to know the NMD pathway in the primitive eukaryotes. Previously, we have found that the giardial homologue of a conserved NMD factor, UPF1, may be functionally conserved and involved in NMD and in preventing nonsense suppression. In this study, we tested the hypothesis that NMD factors can regulate some naturally occurring transcripts in G. lamblia. We found that overexpression of UPF1 resulted in a significant decrease of the levels of CWP1 and cyst formation and of the endogenous cwp1-3, and myb2 mRNA levels and stability. This indicates that NMD could contribute to the regulation of the cwp1-3 and myb2 transcripts, which are key to G. lamblia differentiation into cyst. Interestingly, we also found that UPF1 may be involved in regulation of eight other endogenous genes, including up-regulation of the translation elongation factor gene, whose product increases translation which is required for NMD. Our results indicate that NMD factor could contribute to the regulation of not only nonsense containing mRNAs, but also mRNAs of the key encystation-induced genes and other endogenous genes in the early-diverging eukaryote, G. lamblia.

Identification of the major cysteine protease of Giardia and its role in encystation

Giardia lamblia is a protozoan parasite and the earliest branching clade of eukaryota. The Giardia life cycle alternates between an asexually replicating vegetative form and an infectious cyst form. Encystation and excystation are crucial processes for the survival and transmission of Giardia. Cysteine proteases in Giardia have been implicated in proteolytic processing events that enable the continuance of the life cycle throughout encystation and excystation. Using quantitative real-time PCR, the expression of twenty-seven clan CA cysteine protease genes in the Giardia genome was measured during both vegetative growth and encystation. Giardia cysteine protease 2 was the most highly expressed cysteine protease during both life cycle stages measured, with a dramatic expression increase during encystation. The mRNA transcript for Giardia cysteine protease 2 was 7-fold up-regulated during encystation and was greater than 3-fold higher than any other Giardia protease gene product. Recombinant Giardia cysteine protease 2 was expressed, purified, and biochemically characterized. The activity of the recombinant cysteine protease 2 protein was confirmed to be identical to the dominant cysteine protease activity found in G. lamblia lysates. Giardia cysteine protease 2 was co-localized with cyst wall protein in encystation-specific vesicles during encystation and processed cyst wall protein 2 to the size found in Giardia cyst walls. These data suggest that Giardia cysteine protease 2 is not only the major cysteine endoprotease expressed in Giardia, but is also central to the encystation process.

Unusually low levels of genetic variation among Giardia lamblia isolates

Giardia lamblia, an intestinal pathogen of mammals, including humans, is a significant cause of diarrheal disease around the world. Additionally, the parasite is found on a lineage which separated early from the main branch in eukaryotic evolution. The extent of genetic diversity among G. lamblia isolates is insufficiently understood, but this knowledge is a prerequisite to better understand the role of parasite variation in disease etiology and to examine the evolution of mechanisms of genetic exchange among eukaryotes. Intraisolate genetic variation in G. lamblia has never been estimated, and previous studies on interisolate genetic variation have included a limited sample of loci. Here we report a population genetics study of intra- and interisolate genetic diversity based on six coding and four noncoding regions from nine G. lamblia isolates. Our results indicate exceedingly low levels of genetic variation in two out of three G. lamblia groups that infect humans; this variation is sufficient to allow identification of isolate-specific markers. Low genetic diversity at both coding and noncoding regions, with an overall bias towards synonymous substitutions, was discovered. Surprisingly, we found a dichotomous haplotype structure in the third, more variable G. lamblia group, represented by a haplotype shared with one of the homogenous groups and an additional group-specific haplotype. We propose that the distinct patterns of genetic-variation distribution among lineages are a consequence of the presence of genetic exchange. More broadly, our findings have implications for the regulation of gene expression, as well as the mode of reproduction in the parasite.

Phylogeny and evolution of the SERA multigene family in the genus Plasmodium

The serine repeat antigen gene family of Plasmodium falciparum (Pf-SERA) consists of nine gene members. By sequence similarity search, 45 genes were identified to be homologous to the Pf-SERA genes in the ongoing seven Plasmodium genome sequencing project databases for the species: P. reichenowi, P. vivax, P. knowlesi, P. yoelii, P. berghei, P. chabaudi, and P. gallinaceum. In combination with additional PCR-based sequencing, we found that almost all SERA genes in each species were aligned in a tandem cluster and sandwiched between two conserved hypothetical protein genes, except for P. reichenowi, which could not be confirmed. The minimum and maximum numbers of clustered genes were 2 and 12 for P. gallinaceum and P. vivax, respectively. The best tree of the maximum likelihood analysis demonstrated that all Plasmodium SERA homologues, except for SERA1 of P. gallinaceum (Pg-SERA1), can be classified into four groups, represented by Pf-SERA5, Pf-SERA6, Pf-SERA7, and Pf-SERA8. Genes in the Pf-SERA8 group, although highly divergent and distantly related to the sequences of other groups, were not pseudogenes. P. berghei SERA5, the counterpart of Pf-SERA8, was expressed in the mosquito stage. P. gallinaceum lacks the orthologues to Pf-SERA5, Pf-SERA6, and Pf-SERA7, suggesting that P. gallinaceum diverged from a common ancestor of all eight Plasmodium species examined before gene duplication(s) occurred to generate these paralogous groups. Here, we reveal an evolutionary trail of SERA gene cluster in the genus Plasmodium and discuss a phylogeny of Plasmodium species from the viewpoint of the evolution of a multigene family.

Comprehensive comparative analysis of kinesins in photosynthetic eukaryotes

BACKGROUND

Kinesins, a superfamily of molecular motors, use microtubules as tracks and transport diverse cellular cargoes. All kinesins contain a highly conserved approximately 350 amino acid motor domain. Previous analysis of the completed genome sequence of one flowering plant (Arabidopsis) has resulted in identification of 61 kinesins. The recent completion of genome sequencing of several photosynthetic and non-photosynthetic eukaryotes that belong to divergent lineages offers a unique opportunity to conduct a comprehensive comparative analysis of kinesins in plant and non-plant systems and infer their evolutionary relationships.

RESULTS

We used the kinesin motor domain to identify kinesins in the completed genome sequences of 19 species, including 13 newly sequenced genomes. Among the newly analyzed genomes, six represent photosynthetic eukaryotes. A total of 529 kinesins was used to perform comprehensive analysis of kinesins and to construct gene trees using the Bayesian and parsimony approaches. The previously recognized 14 families of kinesins are resolved as distinct lineages in our inferred gene tree. At least three of the 14 kinesin families are not represented in flowering plants. Chlamydomonas, a green alga that is part of the lineage that includes land plants, has at least nine of the 14 known kinesin families. Seven of ten families present in flowering plants are represented in Chlamydomonas, indicating that these families were retained in both the flowering-plant and green algae lineages.

CONCLUSION

The increase in the number of kinesins in flowering plants is due to vast expansion of the Kinesin-14 and Kinesin-7 families. The Kinesin-14 family, which typically contains a C-terminal motor, has many plant kinesins that have the motor domain at the N terminus, in the middle, or the C terminus. Several domains in kinesins are present exclusively either in plant or animal lineages. Addition of novel domains to kinesins in lineage-specific groups contributed to the functional diversification of kinesins. Results from our gene-tree analyses indicate that there was tremendous lineage-specific duplication and diversification of kinesins in eukaryotes. Since the functions of only a few plant kinesins are reported in the literature, this comprehensive comparative analysis will be useful in designing functional studies with photosynthetic eukaryotes.

Polymorphic insertions and deletions in parabasalian enolase genes

Insertions and deletions in gene sequences have been used as characters to infer phylogenetic relationships and, like any character, the information they contain varies in utility between different levels of evolution. In one case, the absence of two otherwise highly conserved deletions in the enolase genes of parabasalian protists has been interpreted as a primitive characteristic that suggests these were among the first eukaryotes. Here, semi-environmental 3'-RACE was used to sample enolases from parabasalia in the hindgut of the termite Zootermopsis angusticolis to examine the conservation of this character within the parabasalia. Parabasalian homologues were found to be polymorphic for these deletions, and the phylogeny of parabasalian enolases shows that the deletion-possessing genes branch within deletion-lacking genes (i.e., they did not form two clearly distinct groups). Phylogenetic incongruence was detected in the carboxy-terminal third of the sequence (in the region of the deletions), but there is no unambiguous evidence for recombination. The polymorphism of this character discredits these deletions as strong evidence for the early origin of parabasalia, although the complex distribution makes it impossible to state whether parabasalian enolases were ancestrally like those of other eukaryotes. These observations stress the importance of strong corroborating evidence when considering insertion and deletion data, and raises some interesting questions about the apparent variation in degree of conservation of these deletions between different eukaryotic groups.

A class of eukaryotic GTPase with a punctate distribution suggesting multiple functional replacements of translation elongation factor 1alpha

Translation elongation factor 1alpha (EF-1alpha, or EF-Tu in bacteria) is a highly conserved core component of the translation machinery that is shared by all cellular life. It is part of a large superfamily of GTPases that are involved in translation initiation, elongation, and termination, as well as several other cellular functions. Eukaryotic EF-1alpha (eEF-1alpha) is well studied and widely sampled and has been used extensively for phylogenetic analyses. It is generally thought that such highly conserved and functionally integrated proteins are unlikely to be involved in events such as lateral gene transfer or ancient duplication and gene sorting, which would undermine phylogenetic reconstruction. Here we describe a GTPase called EF-like (EFL), which is very similar to, but also distinct from, canonical eEF-1alpha. EFL is found in a wide variety of eukaryotes (dinoflagellates, haptophytes, cercozoa, green algae, choanoflagellates, and fungi), but its distribution is punctate: organisms that possess EFL are not closely related to one another, and EFL appears to be absent from the closest relatives of organisms that do possess it. Moreover, in most genomes where EFL is present, canonical eEF-1alpha appears to be absent. Analysis of functional divergence suggests that, whereas EFL is divergent in general, putative functional binding sites involved in translation are not significantly divergent as a whole. Altogether, it appears that EFL has replaced eEF-1alpha several times independently. This finding could be an indication of an ancient paralogy or, more likely, eukaryote-to-eukaryote lateral gene transfer.

Evolution of eukaryotic transcription: insights from the genome of Giardia lamblia

The Giardia lamblia genome sequencing project affords us a unique opportunity to conduct comparative analyses of core cellular systems between early and late-diverging eukaryotes on a genome-wide scale. We report a survey to identify canonical transcription components in Giardia, focusing on RNA polymerase (RNAP) subunits and transcription-initiation factors. Our survey revealed that Giardia contains homologs to 21 of the 28 polypeptides comprising eukaryal RNAPI, RNAPII, and RNAPIII; six of the seven RNAP subunits without giardial homologs are polymerase specific. Components of only four of the 12 general transcription initiation factors have giardial homologs. Surprisingly, giardial TATA-binding protein (TBP) is highly divergent with respect to archaeal and higher eukaryotic TBPs, and a giardial homolog of transcription factor IIB was not identified. We conclude that Giardia represents a transition during the evolution of eukaryal transcription systems, exhibiting a relatively complete set of RNAP subunits and a rudimentary basal initiation apparatus for each transcription system. Most class-specific RNAP subunits and basal initiation factors appear to have evolved after the divergence of Giardia from the main eukaryotic line of descent. Consequently, Giardia is predicted to be unique in many aspects of transcription initiation with respect to paradigms derived from studies in crown eukaryotes.

The abundance of sterile transcripts in Giardia lamblia

The protozoan parasite Giardia lamblia synthesizes a diverse and surprisingly abundant array of sterile transcripts unable to code for proteins. Random sampling of cDNAs from two evolutionarily divergent Giardia strains indicates that approximately 20% of cDNAs in the libraries represent polyadenylated sterile transcripts. RNase protection analysis and northern blot hybridization of three sterile transcript loci demonstrated that both the sterile transcript and a complementary mRNA were made in each case, further categorizing these sterile transcripts as antisense transcripts. Investigation of the genomic loci for these same three sterile antisense transcripts showed typical transcription units for the sense transcripts, but still failed to reveal a usable open reading frame for the sterile antisense transcripts. 5'-RACE mapped the transcription start site for one of the sterile antisense transcripts to an AT-rich region, as is typical for GIARDIA: It is unclear whether these sterile transcripts represent errors in transcription or whether they have regulatory functions within the cell, although preliminary investigations failed to reveal evidence for a role in developmental gene regulation. In either case, the presence of such a large pool of sterile antisense transcripts is dramatic evidence of the unusual molecular machinery of the early diverging protist G.lamblia.

Biology of Giardia lamblia

Giardia lamblia is a common cause of diarrhea in humans and other mammals throughout the world. It can be distinguished from other Giardia species by light or electron microscopy. The two major genotypes of G. lamblia that infect humans are so different genetically and biologically that they may warrant separate species or subspecies designations. Trophozoites have nuclei and a well-developed cytoskeleton but lack mitochondria, peroxisomes, and the components of oxidative phosphorylation. They have an endomembrane system with at least some characteristics of the Golgi complex and encoplasmic reticulum, which becomes more extensive in encysting organisms. The primitive nature of the organelles and metabolism, as well as small-subunit rRNA phylogeny, has led to the proposal that Giardia spp. are among the most primitive eukaryotes. G. lamblia probably has a ploidy of 4 and a genome size of approximately 10 to 12 Mb divided among five chromosomes. Most genes have short 5' and 3' untranslated regions and promoter regions that are near the initiation codon. Trophozoites exhibit antigenic variation of an extensive repertoire of cysteine-rich variant-specific surface proteins. Expression is allele specific, and changes in expression from one vsp gene to another have not been associated with sequence alterations or gene rearrangements. The Giardia genome project promises to greatly increase our understanding of this interesting and enigmatic organism.

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.

An estimate of divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata by aldolase and triose phosphate isomerase clocks

Previously we suggested that four proteins including aldolase and triose phosphate isomerase (TPI) evolved with approximately constant rates over long periods covering the whole animal phyla. The constant rates of aldolase and TPI evolution were reexamined based on three different models for estimating evolutionary distances. It was shown that the evolutionary rates remain essentially unchanged in comparisons not only between different classes of vertebrates but also between vertebrates and arthropods and even between animals and plants, irrespective of the models used. Thus these enzymes might be useful molecular clocks for inferring divergence times of animal phyla. To know the divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata, the aldolase cDNAs from Ephydatia fluviatilis, a freshwater sponge, and the TPI cDNAs from Ephydatia fluviatilis and Branchiostoma belcheri, an amphioxus, have been cloned and sequenced. Comparisons of the deduced amino acid sequences of aldolase and TPI from the freshwater sponge with known sequences revealed that the Parazoa-Eumetazoa split occurred about 940 million years ago (Ma) as determined by the average of two proteins and three models. Similarly, the aldolase and TPI clocks suggest that vertebrates and amphioxus last shared a common ancestor around 700 Ma and they possibly diverged shortly after the divergence of deuterostomes and protostomes.

Presence of a mitochondrial-type 70-kDa heat shock protein in Trichomonas vaginalis suggests a very early mitochondrial endosymbiosis in eukaryotes

Molecular phylogenetic analyses, based mainly on ribosomal RNA, show that three amitochondriate protist lineages, diplomonads, microsporidia, and trichomonads, emerge consistently at the base of the eukaryotic tree before groups having mitochondria. This suggests that these groups could have diverged before the mitochondrial endosymbiosis. Nevertheless, since all these organisms live in anaerobic environments, the absence of mitochondria might be due to secondary loss, as demonstrated for the later emerging eukaryote Entamoeba histolytica. We have now isolated from Trichomonas vaginalis a gene encoding a chaperone protein (HSP70) that in other lineages is addressed to the mitochondrial compartment. The phylogenetic reconstruction unambiguously located this HSP70 within a large set of mitochondrial sequences, itself a sister-group of alpha-purple bacteria. In addition, the T. vaginalis protein exhibits the GDAWV sequence signature, so far exclusively found in mitochondrial HSP70 and in proteobacterial dnaK. Thus mitochondrial endosymbiosis could have occurred earlier than previously assumed. The trichomonad double membrane-bounded organelles, the hydrogenosomes, could have evolved from mitochondria.

Amplification, expression, and packaging of a foreign gene by giardiavirus in Giardia lamblia

Giardia lamblia is an intestinal protozoan parasite and one of the earliest eukaryotic divergents. The trophozoite multiplies via asexual binary fission and lacks all natural means of lateral gene transfer. A system is developed here for long-term expression of a foreign gene in this organism by exploiting recombinant virions derived from the giardiavirus (GLV), a double-stranded RNA virus that infects many Giardia isolates. An in vitro transcript of the cloned GLV cDNA, comprising the firefly luciferase-encoding region flanked by 5' and 3' fragments of GLV positive-strand RNA, was electroporated into GLV-infected trophozoites. Luciferase activity in electroporated cells peaked on day 2 at levels 6 orders of magnitude above background. Expression of this foreign gene remained at 80% of its peak level after 30 days in the absence of selective pressure. The chimeric RNA was replicated as double-stranded RNA and packaged into virus-like particles. The recombinant virions were partially purified from the wild-type helper virus by CsCl equilibrium density-gradient centrifugation and used to superinfect Giardia trophozoites. At multiplicities of infection of 100 or higher, these chimeric virions were able to initiate new rounds of expression of luciferase activity in the superinfected cells. Thus, the engineered virion can be successfully used to introduce and efficiently express a heterologous gene in this eukaryotic microorganism.

Unbiased estimation of symmetrical directional mutation pressure from protein-coding DNA

The most generally applicable procedure for obtaining estimates of the symmetrical, or strandnonspecific, directional mutation pressure (microD) on protein-coding DNA sequences is to determine the G+C content at synonymous codon sites (Psyn), and to divide Psyn by twice the arithmetic mean of the G+C content at synonymous codon sites of a large number of randomly generated, synonymously coding DNA sequences (Psyn). Unfortunately, the original procedure yields biased estimates of Psyn and microD and is computationally expensive. We here present a fast procedure for estimating unbiased microD values. The procedure employs direct calculation of Psyn (approximately Psyn) and two normalization procedures, one for Psyn < or = Psyn and another for Psyn > or = Psyn. The normalization removes a bias sometimes caused by codons specifying arginine, asparagine, isoleucine, and leucine. Consequently, comparison of protein-coding genes that are translated using different genetic codes is facilitated.

Protein phylogeny of translation elongation factor EF-1 alpha suggests microsporidians are extremely ancient eukaryotes

Partial regions of the mRNA encoding a major part of translation elongation factor 1 alpha (EF-1 alpha) from a mitochondrion-lacking protozoan, Glugea plecoglossi, that belongs to microsporidians, were amplified by polymerase chain reaction (PCR) and their primary structures were analyzed. The deduced amino acid sequence was highly divergent from typical EF-1 alpha's of eukaryotes, although it clearly showed a eukaryotic feature when aligned with homologs of the three primary kingdoms. Maximum likelihood (ML) analyses on the basis of six different stochastic models of amino acid substitutions and a maximum parsimony (MP) analysis consistently suggest that among eukaryotic species being analyzed, G. plecoglossi is likely to represent the earliest offshoot of eukaryotes. Microsporidians might be the extremely ancient eukaryotes which have diverged before an occurrence of mitochondrial symbiosis.

Molecular and phylogenetic analysis of PCR-amplified cyclin-dependent kinase (CDK) family sequences from representatives of the earliest available lineages of eukaryotes

Cyclin-dependent kinase (CDK) and cell division control (CDC2) sequences are strongly conserved among eukaryotes and may complement the use of other sequence families in eukaryotic phylogenetic inference. We synthesized degenerate PCR primers to amplify the catalytic region of CDK homologs in representatives of the earliest available lineages of eukaryotes. CDK family sequence-based, maximum-likelihood distance measurements with neighbor-joining, and Fitch-Margoliash least-squares analyses produced unrooted dendrograms that included protists, yeasts, and higher eukaryotes. Bootstrap confidence estimates supported CDK-based phylogenetic groupings among the protists, fungi, and vertebrates although resolution within these groups was often insignificant. However, Trichomonas vaginalis and Giardia lamblia exhibited two of the most divergent CDK-like sequences consistent with rRNA-phylogenetic inference of early divergence of these eukaryotic lineages. In the evolution from unicellular to multicellular organisms, a constellation of amino acid residues aligning with the human, CDK N-terminal beta sheet may have undergone abrupt replacement.

A nuclear gene of eubacterial origin in Euglena gracilis reflects cryptic endosymbioses during protist evolution

Genes for glycolytic and Calvin-cycle glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of higher eukaryotes derive from ancient gene duplications which occurred in eubacterial genomes; both were transferred to the nucleus during the course of endosymbiosis. We have cloned cDNAs encoding chloroplast and cytosolic GAPDH from the early-branching photosynthetic protist Euglena gracilis and have determined the structure of its nuclear gene for cytosolic GAPDH. The gene contains four introns which possess unusual secondary structures, do not obey the GT-AG rule, and are flanked by 2- to 3-bp direct repeats. A gene phylogeny for these sequences in the context of eubacterial homologues indicates that euglenozoa, like higher eukaryotes, have obtained their GAPDH genes from eubacteria via endosymbiotic (organelle-to-nucleus) gene transfer. The data further suggest that the early-branching protists Giardia lamblia and Entamoeba histolytica--which lack mitochondria--and portions of the trypanosome lineage have acquired GAPDH genes from eubacterial donors which did not ultimately give rise to contemporary membrane-bound organelles. Evidence that "cryptic" (possibly ephemeral) endosymbioses during evolution may have entailed successful gene transfer is preserved in protist nuclear gene sequences.

Phylogenetic position of Dictyostelium inferred from multiple protein data sets

The phylogenetic position of Dictyostelium inferred from 18S rRNA data contradicts that from protein data. Protein trees always show the close affinity of Dictyostelium with animals, fungi, and plants, whereas in 18S rRNA trees the branching of Dictyostelium is placed at a position before the massive radiation of protist groups including the divergence of the three kingdoms. To settle this controversial issue and to determine the correct position of Dictyostelium, we inferred the phylogenetic relationship among Dictyostelium and the three kingdoms Animalia, Fungi, and Plantae by a maximum-likelihood method using 19 different protein data sets. It was shown at the significance level of 1 SE that the branching of Dictyostelium antedates the divergence of Animalia and Fungi, and Plantae is an outgroup of the Animalia-Fungi-Dictyostelium clade.

Transient transfection and expression of firefly luciferase in Giardia lamblia

We have developed a gene transfer system for the protozoan parasite Giardia lamblia. This organism is responsible for many cases of diarrhea worldwide and is considered to be one of the most primitive eukaryotes. Expression of a heterologous gene was detected in this parasite after electroporation with appropriate DNA constructs. We constructed a series of transfection plasmids using flanking sequences of the Giardia glutamate dehydrogenase (GDH) gene to drive expression of the firefly luciferase reporter gene. The optimal construct consisted of a GDH/luciferase fusion gene in which the first 18 codons of the GDH gene immediately preceded the luciferase gene; this fusion gene was flanked by the upstream and downstream sequences of the GDH gene. Electroporation of this construct into Giardia yielded luciferase activity that was 3000- to 50,000-fold above background. Removal of either the 5' or 3' GDH flanking sequences from this construct resulted in significantly reduced luciferase activity, and removal of both flanking sequences reduced luciferase activity to background levels. Luciferase activity was proportional to the amount of DNA electroporated and was maximal at 6 hr after electroporation.

Improved dating of the human/chimpanzee separation in the mitochondrial DNA tree: heterogeneity among amino acid sites

The internal branch lengths estimated by distance methods such as neighbor-joining are shown to be biased to be short when the evolutionary rate differs among sites. The variable-invariable model for site heterogeneity fits the amino acid sequence data encoded by the mitochondrial DNA from Hominoidea remarkably well. By assuming the orangutan separation to be 13 or 16 Myr old, a maximum-likelihood analysis estimates a young date of 3.6 +/- 0.6 or 4.4 +/- 0.7 Myr (+/- 1 SE) for the human/chimpanzee separation, and these estimates turn out to be robust against differences in the assumed model for amino acid substitutions. Although some uncertainties still exist in our estimates, this analysis suggests that humans separated from chimpanzees some 4-5 Myr ago.

Phylogenetic relationships among eutherian orders estimated from inferred sequences of mitochondrial proteins: instability of a tree based on a single gene

The phylogenetic relationships among Primates (human), Artiodactyla (cow), Cetacea (whale), Carnivora (seal), and Rodentia (mouse and rat) were estimated from the inferred amino acid sequences of the mitochondrial genomes using Marsupialia (opossum), Aves (chicken), and Amphibia (Xenopus) as an outgroup. The overall evidence of the maximum likelihood analysis suggests that Rodentia is an outgroup to the other four eutherian orders and that Cetacea and Artiodactyla form a clade with Carnivora as a sister taxon irrespective of the assumed model for amino acid substitutions. Although there remains an uncertainty concerning the relation among Artiodactyla, Cetacea, and Carnivora, the existence of a clade formed by these three orders and the outgroup status of Rodentia to the other eutherian orders seems to be firmly established. However, analyses of individual genes do not necessarily conform to this conclusion, and some of the genes reject the putatively correct tree with nearly 5% significance. Although this discrepancy can be due to convergent or parallel evolution in the specific genes, it was pointed out that, even without a particular reason, such a discrepancy can occur in 5% of the cases if the branching among the orders in question occurred within a short period. Due to uncertainty about the assumed model underlying the phylogenetic inference, this can occur even more frequently. This demonstrates the importance of analyzing enough sequences to avoid the danger of concluding an erroneous tree.