*Josephine Bay-Paul Center for Comparative Molecular Biology and Evolution, Marine Biological Laboratory, Woods Hole, Massachusetts;
and
School of Biological Sciences, University of Sydney, Sydney, New South Wales, Australia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Several lines of evidence argue that jakobids are among the earliest diverging of eukaryotic lineages. First, they were implicated to be key organisms in the establishment of the mitochondrial symbiosis because of their ultrastructural similarity to the retortamonads and, indirectly, diplomonads. Diplomonads and retortamonads lack mitochondria (Brugerolle 1991
), consistent with a divergence prior to the establishment of the mitochondrial symbiosis (Cavalier-Smith 1983, 1987
; Sleigh 1989
; Sogin et al. 1989
; Sogin 1991
; Roger, Morrison, and Sogin 1999
). This hypothesis was supported by the early divergence of the diplomonads in molecular phylogenies of SSU rDNA and several protein families (Hashimoto et al. 1994
; Sogin et al. 1989
). On the basis of these data, O'Kelly suggests that the ancestral mitochondrial symbiont was first engulfed by a sister taxon to the retortamonads, producing a "jakobid-like" ancestor from which all mitochondriate eukaryotes evolved (O'Kelly 1993
). This argument has been weakened by the discovery of genes of mitochondrial origin in the diplomonad Giardia (Hashimoto et al. 1998
; Roger et al. 1998
; Roger 1999
), indicating that the symbiosis took place prior to the divergence of diplomonads and parabasalids. Nevertheless, a deeply diverging position of diplomonads among eukaryotes remains a plausible scenario and implies that mitochondriate excavate taxa such as jakobids may similarly be early-diverging eukaryotes.
Second, O'Kelly (1993)
noted that the jakobids encompass all three of the basic forms of mitochondrial cristae known in eukaryotes. Histionidae have tubular cristae, while Jakoba libera has flat cristae, and the sole studied member of Malawimonadidae, Malawimonas jakobiformis, has discoidal cristae (Patterson 1990
; O'Kelly 1993
; O'Kelly and Nerad 1999
). Mitochondrial cristae shape has been seen as a strongly conserved character and has been used to delimit the deepest evolutionary divisions within eukaryotes (Taylor 1976
; Sleigh 1989
; Patterson 1994
). Although the presence of all three forms in jakobids implies that jakobids are not monophyletic, this was also consistent with a model in which jakobid flagellates were a stem group from which all extant mitochondriate eukaryotes emerged in at least three separate radiations (O'Kelly 1993
).
Third, and most significantly, the mitochondrial genomes of Reclinomonas americana and Jakoba libera preserve ancestral bacterial features not found in any other eukaryotes studied to date. Mitochondrial DNAs (mtDNAs) from R. americana (Lang et al. 1997
) and J. libera (http://megasun.bch.umontreal.ca/ogmp/projects/individual.html) possess gene clusters that closely resemble ribosomal protein operons of eubacteria. They also code for RNase P RNAs with bacterial minimum consensus structures, as well as subunits of bacterial- type RNA polymerase genes that were ancestrally present in the proto-mitochondrial symbiont genome (Lang et al. 1997, 1999
). Indeed, all other mitochondrial genomes sequenced to date encode a subset of the proteins described in jakobid mtDNAs (Burger et al. 1999
; Gray, Burger, and Lang 1999
). Nonetheless, jakobid mtDNAs share several unique gene arrangements with mitochondrial genomes from other eukaryotic lineages to the exclusion of all known bacterial genomes. These gene arrangements, together with the monophyly of all mtDNAs (including those of jakobids) recovered in phylogenies of single and concatenated mitochondrial gene data sets, indicate that all known mitochondria evolved from a single endosymbiosis (Gray, Burger, and Lang 1999
; Lang et al. 1999
). The unique retention of many ancestral alpha-proteobacterial genes in the jakobid mtDNAs is most parsimoniously explained by "jakobids" being among the earliest diverging mitochondrion-containing taxa.
For all of their potential evolutionary significance, our phylogenetic understanding of "jakobids" is incomplete. Ultrastructural data indicate that the Jakobidae (Jakoba) and the Histionidae (Reclinomonas and Histiona) are closely related but also suggest that Malawimonas may have evolved independently from within a broader radiation of excavate taxa (O'Kelly and Nerad 1999
; Simpson and Patterson 1999
). Although phylogenies of concatenated mitochondrial genes also indicate that J. libera and R. americana form a clade, this group is not clearly allied with Malawimonas (Burger et al. 1999
). In light of the possibility of jakobid polyphyly, we will refer to Jakobidae and Histionidae collectively as the "core jakobids" (Simpson and Patterson 1999
) and treat Malawimonas as a separate lineage. Our understanding of the phylogenetic placement of both groups remains hampered by the absence of nuclear gene sequences.
Over the last decade, phylogenetic reconstructions of small-subunit ribosomal RNA (SSU rRNA) gene sequences have provided a molecular phylogenetic framework for understanding the course of early eukaryotic evolution. In these phylogenies, the amitochondriate diplomonads, microsporidia, and parabasalids consistently branch before any mitochondrion-bearing lineages (Sogin 1991
; Cavalier-Smith 1993
; Leipe et al. 1993
). However, there is doubt as to whether SSU rRNA correctly identifies the deepest diverging groups and their relative branching order (Leipe et al. 1993
; Hirt et al. 1999
). Large disparities in the rates of evolution (leading to long-branch attraction) and extreme biases in base compositions of different SSU rRNA sequences, especially among "early-branching" eukaryotic taxa, may confound accurate phylogenetic reconstruction (Loomis and Smith 1990
; Hasegawa et al. 1993
; Hasegawa and Hashimoto 1993
; Galtier and Gouy 1995
; Roger et al. 1999
; Stiller and Hall 1999
; Philippe and Germot 2000
). Highly conserved protein-coding genes might be less sensitive to some of these problems for several of the taxa identified as "deep-branches" in SSU rRNA trees (Hasegawa et al. 1993
; Hashimoto et al. 1994
; Roger et al. 1996
; Hirt et al. 1999
). Therefore, the development of these markers in combination with SSU rRNA should help in determining a more robust global eukaryotic phylogeny.
To this end, we focused on the tubulin gene family as a phylogenetic marker for early eukaryotic evolution. This gene family consists of three highly conserved subfamilies, alpha-, beta-, and gamma-tubulin, that arose from a series of gene duplications in early eukaryotic evolution (Edlind et al. 1996
; Keeling and Doolittle 1996
). Currently, alpha- and beta-tubulin genes have the widest taxonomic representation and have been identified in the putatively early diverging eukaryotic lineages (Kirk-Mason, Turner, and Chakraborty 1988
; Burns 1991
; Katiyar and Edlind 1994, 1996
; Edlind et al. 1996
; Keeling and Doolittle 1996
). In order to study the relationships among the "jakobid" flagellates and their relationship to other "excavate" protist lineages, we characterized the alpha- and beta-tubulin genes from four jakobid species: Jakoba libera, Jakoba incarcerata, Reclinomonas americana, and Malawimonas jakobiformis.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
PCR Amplification, Cloning, and Sequencing
Alpha- and beta-tubulin genes were amplified from genomic DNA. Degenerate primers were previously designed based on highly conserved N- and C-terminal regions of the alpha- and beta-tubulin gene families (Roger 1996
). BtubA was based on the amino acid motif GQCGNQ and had the sequence 5'-GCAGGNCARTGYGGNAAYCA-3'; BtubB was based on MDEMEFT and had the sequence 5'-AGTRAAYTCCATYTCRTCCAT-3'; AtubA was based on QVGNCWE with the sequence 5'-RGTNGGNAAYGCNTGYTGGGA-3'; and AtubB was based on WYVGEGM with the sequence 5'-CCATNCCYTCNCCNACRTACCA-3'. These primer combinations amplified >90% of the coding region of each gene. Amplification reactions were performed with 50100 ng of total genomic DNA under standard conditions both in the presence and in the absence of 5%6% acetamide (final concentration). Thermal cycling consisted of 35 cycles with an annealing temperature of 4855°C for 1 min, an extension temperature of 72°C for 1 min, and a denaturing temperature of 94°C for 30 s.
PCR products were purified from agarose gels using the Prep-A-Gene DNA Purification Kit (BioRad) and cloned directly into pGEM-T-easy vector (Promega) or TOPO-XL vector (Invitrogen, Inc.). Sequence data were gathered on a LI-COR 4200L automated sequencing apparatus using infrared dye-labeled T7 and M13R primers (LI-COR) and a cycle sequencing protocol (Sequitherm, Epicenter, Inc.). Single-stranded sequence of 410 clones of each PCR product was initially determined to test for the presence of multiple distinct copies of tubulin genes within each organism and to determine internal restriction sites for subcloning. Gene fragments of 100900 bp were generated by restriction enzyme digestion of plasmid templates and subsequently subcloned into pBluescript KS (-) or pUC-19 vectors. Full- length double-stranded DNA sequences of subcloned distinct gene copies were determined by combined analysis of both complete gene and subcloned template data using ALIGNIR sequence assembly software (LI-COR).
Phylogenetic Analysis
Nucleotide sequence data were submitted to NCBI BLAST (Basic Local Alignment Search Tool) (Altschul et al. 1990
) for comparison with available sequence databases to verify tubulin gene identity and locate introns. Protein sequences were compiled and manually aligned in PAUP*, version 4.0b2 (Swofford 1999
), along with tubulin sequences from representatives of diverse taxa obtained using the NCBI Entrez protein sequence search program. In total, a data set of 63 alpha tubulins and 85 beta tubulin sequences was compiled, with representatives from all available eukaryotic groups.
The alignment of intron sequences was obtained using the CLUSTAL W program, version 1.7 (Thompson, Higgins, and Gibson 1994
) under default parameter settings and manually adjusted by eye. Phylogenetic analyses of the R. americana complete beta-tubulin DNA sequences including introns were performed using the maximum-likelihood (ML) method implemented in PUZZLE, version 4.02 (Strimmer and von Haeseler 1996
), employing the Hasegawa, Kishino and Yano (1985)
model with an eight-category discrete gamma model of among-site rate variation (HKY+
). ML estimation on a neighbor-joining topology was used to derive a transition : transversion parameter of 1.21 and a gamma shape parameter (
) of 0.02. Bootstrapping (1,000 resamplings) was accomplished using the SEQBOOT and CONSENSE programs (from PHYLIP [Felsenstein 1995
]) in combination with QBOOT, a unix shell-script that allows quartet ML bootstrap analysis with PUZZLE, version 4.02 (available on request from A.J.R.).
Phylogenetic analyses of the amino acid sequences were performed using both the protein ML distance (ML-distance) and the ML quartet puzzling (ML-QP) methods implemented in PUZZLE, version 4.02 (Strimmer and von Haeseler 1996
). Data missing in more than a single sequence and regions of ambiguous alignment were removed, yielding 381 sites for the alpha-tubulin analysis, 387 sites for the beta-tubulin analysis, and 768 sites for a combined (concatenated) alpha- plus beta- tubulin analysis. The ML model employed in both distance and likelihood analyses was the Jones, Taylor, and Thornton (1992)
amino acid replacement matrix with an eight-category discrete gamma model of among-sites rate variation (JTT+
). The gamma shape parameter
was estimated by ML optimization on a neighbor-joining topology to be 0.72 for the alpha-tubulin data set, 0.52 for the beta-tubulin data set, and 0.57 for the concatenated alpha+beta-tubulin data set. Phylogenetic trees were inferred from the ML distances using the Fitch-Margoliash algorithm with global rearrangements (FITCH, version 3.57c [Felsenstein 1995
]). ML distance bootstrap values for bipartitions were calculated by analysis of 100 resampled data sets generated with the SEQBOOT program and analyzed with the PUZZLEBOOT script (http://www.tree-puzzle.de/) in conjunction with PUZZLE, version 4.02, and CONSENSE. Full protein ML analyses were performed using the quartet puzzling algorithm in PUZZLE, version 4.02, utilizing 1,000 puzzling steps.
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Another noteworthy property of the tubulin genes described here is a moderate to extreme codon usage bias toward codons ending in G or C (table 1 ). The GC- bias of alpha- and beta-tubulin genes from each organism were very similar, confirming that each set of genes derive from the same genome.
|
|
|
Placement of "Jakobid" Flagellates in Tubulin Phylogenies
ML distance/Fitch and ML-QP trees for each of the three analyses (figs. 1A, 1B, and 2
) show that the "jakobids" in the broad sense (Reclinomonas, Jakoba, and Malawimonas) (O'Kelly 1993
) and the "excavate taxa" (in this analysis, jakobids, heteroloboseids, and diplomonads) do not form clades. However, our analyses do show a close grouping between beta-tubulins from J. libera, R. americana 50283, and R. americana 50284 (99% bootstrap support in a distance analysis, and quartet puzzling support of 83%) (fig. 2B
). The extremely close relationship between R. americana and J. libera beta-tubulins (fig. 1B
) is particularly interesting because these two sequences each harbor a single nonhomologous intron that occupies different positions in the gene (data not shown). Neither intron position is found in any other beta-tubulin gene studied to date (J. M. Logsdon Jr., personal communication), thus indicating recent intron turnover in beta-tubulin in the R. americanaJ. libera lineage.
Interestingly, the J. liberaR. americana grouping is not recovered in alpha-tubulin phylogenies by either method of analysis (fig. 1A
). Lack of support for this grouping, which is well supported by other data (discussed below), is puzzling and may have several causes. For instance, placement of R. americana and J. libera alpha-tubulin sequences may be obscured by the overall poorer support for monophyletic groupings observed in the alpha-tubulin tree, especially within the plant-protist superclade (fig. 1A
). It is also possible that one of these sequences represents an aberrant paralog of alpha-tubulin which could have life-cycle-specific expression and has an altered pattern of substitution, a frequently observed phenomenon in tubulins (Kube-Granderath and Schliwa 1998
). If so, then the orthologous alpha- tubulin gene(s) from one or another of these organisms may not have been detected in our PCR approach, perhaps due to introns interrupting the primer binding sites in the genes. However, the J. liberaR. americana relationship was well supported in the combined alpha- and beta-tubulin analysis (91% bootstrap support and 90% QP support; fig. 2
), indicating that the alpha-tubulin data do not significantly conflict with the beta- tubulin data, but instead lack strong enough phylogenetic signal to resolve the relationships between these taxa.
Also unexpected was the positioning of J. incarcerata in tubulin trees. Jakoba incarcerata was assigned to Jakoba primarily on the basis of light microscope observations (Bernard, Simpson, and Patterson 2000
). Ultrastructural data support assignation of J. incarcerata to the "core jakobids" (Bernard, Simpson, and Patterson 2000
; unpublished data) but also suggest that J. incarcerata may be basal to other known "core jakobids" (i.e., Jakoba may not be monophyletic). Jakoba incarcerata would therefore be expected to group with J. libera and R. americana in a "core jakobid" clade (Simpson and Patterson 1999
), but not necessarily to group specifically with J. libera. However, in our beta-tubulin analyses, J. incarcerata did not cluster with these taxa, emerging separately within the plant-protist superclade in beta-tubulin trees. In alpha-tubulin phylogeny, it joined the animals-plus-fungiparabasalid grouping with very weak support (45% distance bootstrap and 59% QP support; fig. 1A
). Interestingly, this sequence robustly grouped with the plant-protist superclade in the combined alpha- plus beta-tubulin analysis (100% distance bootstrap and 74% QP support) and emerged with moderate support as the earliest-diverging lineage within this clade (70% distance bootstrap and QP support). Similarly, M. jakobiformis occupied a relatively basal position in the plant-protist superclade in all three analyses, varying in position only slightly (figs. 1 and 2
).
To evaluate the stability of the position of J. incarcerata in the combined analysis in light of the weakly supported affinity of J. incarcerata for the animals- plus-fungiparabasalid grouping in the alpha-tubulin analysis, we performed a separate combined alpha- plus beta-tubulin analysis that included the alpha- and beta- tubulin genes from a second diplomonad, Giardia intestinalis, and a composite trichomonad sequence composed of Trichomonas vaginalis beta-tubulin and Monocercomonas sp. alpha-tubulin (data not shown). In this analysis, there were only minor changes in support values at several nodes, and the position of the jakobids, including J. incarcerata, did not change from the combined alpha- plus beta-tubulin analysis presented here.
There is little, if any, resolution in the branching order among the major groups in the plant-protist superclade in our tubulin trees, and, in general, the jakobid lineages show inconsistent affinities to other taxa in the various analyses. One exception is the tendency in beta- tubulin trees for the J. liberaR. americana group to branch with the heteroloboseid Naegleria and the euglenozoan Euglena (fig. 1B ). This tendency is amplified, albeit weakly (44% distance bootstrap support and 56% QP support), in the combined alpha- plus beta-tubulin analysis, in which the J. liberaR. americana group became a sister group of a Heterolobosea-Euglenozoa clade. The best resolution of the major phyla and their relative branching order was obtained in our combined alpha- plus beta-tubulin analyses, yet support values for some groups (e.g., Viridiplantae) remained low. It is likely that the extreme conservation of both tubulin genes leave few phylogenetically informative sites, especially for groups in the plant-protist superclade. In this respect, ancient phylogenetic signal appears to be retained in tubulin sequences, but its overall amount is low.
Jakobid Polyphyly
Our observation of a J. liberaR. americana grouping is in agreement with concatenated mitochondrial gene phylogenies (Burger et al. 1999
; Gray, Burger, and Lang 1999
; Burger et al. 2000
), SSU rRNA phylogenies (unpublished data), and shared unique mtDNA gene composition between these taxa (e.g., their common possession of bacterial DNA-dependent RNA polymerase genes) (Lang et al. 1999
). Concatenated mitochondrial gene phylogenies also support the distinctness of Malawimonas from "core jakobids" and show them emerging from within a strikingly similar plant-protist superclade (although the monophyly of this clade is not well supported). These results are consistent with recent evaluations of ultrastructural features of excavate taxa (O'Kelly, Farmer, and Nerad 1999
; O'Kelly and Nerad 1999
; unpublished data) which indicate that "core jakobids" (including J. libera and R. americana) share several characters not found in other excavate taxa but that only one or two features of debatable importance are unique to Malawimonas and "core jakobids." Given the conflict between ultrastructural and tubulin data in the placement of J. incarcerata, further data on this taxon would be desirable. Unfortunately, mitochondrial or nuclear gene data have not yet been published.
The Cristal Hypothesis
O'Kelly's (1993)
original suggestion for jakobids as a stem group from which all three basic cristal morphologies arose requires the jakobids to be paraphyletic, giving rise to at least three clades of nonjakobid organisms. While it could be argued that this degree of paraphyly of jakobids is reconcilable with our results, it appears safe to reject O'Kelly's hypothesis from the evidence at hand. First, O'Kelly's hypothesis requires that the J. liberaR. americana grouping be interrupted by at least two nonjakobid clades (one essentially platycristate, the other tubulocristate). Second, it would predict that Malawimonas would fall adjacent to Heterolobosea and Euglenozoa to form a clade of discicristate organisms. In fact, it is the J. liberaR. americana clade which shows some affinity for Heterolobosea and Euglenozoa (figs. 1 and 2
).
Jakobids and Deep-Branching Excavate Taxa
Tubulin genes are the first molecules to be used to evaluate the affinities of jakobids to amitochondriate "excavate" protist lineages. We did not observe a close relationship between the diplomonads and any of the jakobid lineages in this study, although sequences from Malawimonas and J. incarcerata branched near the base of the plant-protist superclade and therefore emerged in a region of the tree close to the diplomonad branch. Currently, the "excavate hypothesis" can be rationalized in terms of tubulin phylogenies in the following manner. If some jakobid lineages (e.g., J. incarcerata and Malawimonas) are among the deepest branches in a plant- protist superclade, and diplomonads are the outgroup to these and all other eukaryotes, then the common ancestor of the plant-protist superclade and the common ancestor of extant eukaryotes could have had an "excavate" phenotype. Paradoxically, it is precisely R. americana and J. libera that display the most primitive mitochondrial genome organization (Lang et al. 1999
), yet they do not seem to occupy a particularly deep branch in tubulin trees. Clearly, more molecular data will be required from jakobids and other excavates, such as retortamonads, Trimastix, and Carpediemonas, to decide where these protists fit in the large-scale phylogeny of eukaryotes.
Perhaps even more problematic will be understanding the basis of conflicts between large-scale eukaryotic phylogenies inferred with the various nuclear gene markers (e.g., tubulin vs. SSU rRNA phylogenies). Substitution rate calibration analyses of "crown" eukaryotic rRNAs support the division between an animals-plus- fungi group and a large heterogeneous plant-protist clade of eukaryotes (Van de Peer et al. 1996
), as do concatenated mitochondrial genes (Burger et al. 1999
). However, the most conspicuous disagreement between SSU rRNA trees and tubulin trees is the fact that "deep" lineages in SSU rRNA trees (Leipe et al. 1993
), such as Microsporidia, Eumycetozoa, Euglenozoa, and Heterolobosea, emerge from within the these two major eukaryotic clades in tubulin trees. Several authors have suggested that the deep-branching position of these and other taxa within SSU rRNA trees could result from long-branch attraction artifacts (Embley and Hirt 1998
; Keeling, Deane, and McFadden 1998
; Philippe and Adoutte 1998
; Roger, Morrison, and Sogin 1999
) and that the alternative positioning of these taxa within the tubulin tree may be correct (Embley and Hirt 1998
; Keeling, Deane, and McFadden 1998
; Dacks and Roger 1999
). Although this seems likely to be the case for Microsporidia (Keeling and McFadden 1998
), we should also be cautious in placing too much confidence in tubulin phylogenies. Most taxa that lack flagella or centrioles, such as Microsporidia, higher fungi, Entamoeba, Dictyostelium, and red algae, tend to form extremely long branches in tubulin trees, and often group together (Edlind et al. 1996
; Keeling and Doolittle 1996
; Keeling et al. 1999
; Keeling, Luker, and Palmer 2000
; unpublished data), as do aberrant developmentally regulated tubulin isoforms (data not shown). Therefore, much of the structure of tubulin-based trees could also result from long-branch attraction effects.
![]() |
Supplementary Material |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Keywords: jakobid
alpha-tubulin
beta-tubulin
phylogeny
flagellates
1 Present address: Department of Biochemistry and Molecular Biology, Dalhousie University, Halifax, Nova Scotia, Canada
3 Address for correspondence and reprints: Mitchell L. Sogin, Josephine Bay-Paul Center of Comparative Molecular Biology and Evolution, Marine Biological Laboratory, 7 MBL Street, Woods Hole, Massachusetts 02543. sogin{at}evol5.mbl.edu
![]() |
literature cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Altschul, S. F., W. Gish, W. Miller, E. W. Myers, and D. J. Lipman. 1990. Basic local alignment search tool. J. Mol. Biol. 215:403410[ISI][Medline]
Baldauf, S. L., and J. D. Palmer. 1993. Animals and fungi are each other's closest relatives: congruent evidence from multiple proteins. Proc. Natl. Acad. Sci. USA 90:11558 11562
Bernard, C., A. G. B. Simpson, and D. J. Patterson. 2000. Some free-living flagellates from anoxic sediments. Ophelia 52:113142
Brugerolle, G. 1991. Flagellar and cytoskeletal systems in amitochondrial flagellates: Archamoeba, Metamonada and Parabasala. Protoplasma 164:7090
Burger, G., D. Saint-Louis, M. W. Gray, and B. F. Lang. 1999. Complete sequence of the mitochondrial DNA of the red alga Porphyra purpurea: cyanobacterial introns and shared ancestry of red and green algae. Plant Cell 11:1675 1694
Burger, G., Y. Zhu, T. G. Littlejohn, S. J. Greenwood, M. N. Schnare, B. F. Lang, and M. W. Gray. 2000. Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA. J. Mol. Biol. 297:365380[ISI][Medline]
Burns, R. G. 1991. Alpha-, beta-, and gamma-tubulins: sequence comparisons and structural constraints. Cell Motil. Cytoskeleton 20:181189
Cavalier-Smith, T. 1983. A 6-kingdom classification and a unified phylogeny. Pp. 10271034 in W. Schwemmler and H. E. A. Schenk, eds. Endocytobiology. Vol. 2. de Gruyter, Berlin
. 1987. Eukaryotes with no mitochondria. Nature 326: 332333
. 1993. Kingdom Protozoa and its 18 phyla. Microbiol. Rev. 57:953994[Abstract]
Dacks, J., and A. J. Roger. 1999. The first sexual lineage and the relevance of facultative sex. J. Mol. Evol. 48:779783[ISI][Medline]
Edlind, T. D., J. Li, G. S. Visvesvara, M. H. Vodkin, G. L. McLaughlin, and S. K. Katiyar. 1996. Phylogenetic analysis of beta-tubulin sequences from amitochondrial protozoa. Mol. Phylogenet. Evol. 5:359367[ISI][Medline]
Embley, T. M., and R. P. Hirt. 1998. Early branching eukaryotes? Curr. Opin. Genet. Dev. 8:624629[ISI][Medline]
Fast, N. M., J. M. Logsdon Jr., and W. F. Doolittle. 1999. Phylogenetic analysis of the TATA box binding protein (TBP) gene from Nosema locustae: evidence for a microsporidia-fungi relationship and spliceosomal intron loss. Mol. Biol. Evol. 16:14151419
Felsenstein, J. 1995. PHYLIP (phylogeny inference package). Version 3.57c. Distributed by the author, Department of Genetics, University of Washington, Seattle
Flavin, M., and T. A. Nerad. 1993. Reclinomonas americana N. G., N. Sp., a new freshwater heterotrophic flagellate. J. Eukaryot. Microbiol. 40:172179
Galtier, N., and M. Gouy. 1995. Inferring phylogenies from DNA sequences of unequal base compositions. Proc. Natl. Acad. Sci. USA 92:1131711321
Germot, A., H. Philippe, and H. Le Guyader. 1997. Evidence for loss of mitochondria in Microsporidia from a mitochondrial-type HSP70 in Nosema locustae. Mol. Biochem. Parasitol. 87:159168
Gray, M. W., G. Burger, and B. F. Lang. 1999. Mitochondrial evolution. Science 283:14761477
Hasegawa, M., and T. Hashimoto. 1993. Ribosomal RNA trees misleading? Nature 361:23
Hasegawa, M., T. Hashimoto, J. Adachi, N. Iwabe, and T. Miyata. 1993. Early branchings in the evolution of eukaryotes: ancient divergence of entamoeba that lacks mitochondria revealed by protein sequence data. J. Mol. Evol. 36:380388[ISI][Medline]
Hasegawa, M., H. Kishino, and T. Yano. 1985. Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22:160174[ISI][Medline]
Hashimoto, T., Y. Nakamura, F. Nakamura, T. Shirakura, J. Adachi, N. Goto, K. Okamoto, and M. Hasegawa. 1994. Protein phylogeny gives a robust estimation for early divergences of eukaryotes: phylogenetic place of a mitochondria-lacking protozoan, Giardia lamblia. Mol. Biol. Evol. 11:6571[Abstract]
Hashimoto, T., L. B. Sanchez, T. Shirakura, M. Muller, and M. Hasegawa. 1998. Secondary absence of mitochondria in Giardia lamblia and Trichomonas vaginalis revealed by valyl-tRNA synthetase phylogeny. Proc. Natl. Acad. Sci. USA 95:68606865
Hirt, R. P., B. Healy, C. R. Vossbrinck, E. U. Canning, and T. M. Embley. 1997. A mitochondrial Hsp70 orthologue in Vairimorpha necatrix: molecular evidence that microsporidia once contained mitochondria. Curr. Biol. 7:995998[ISI][Medline]
Hirt, R. P., J. M. Logsdon, B. Healy, M. W. Dorey, W. F. Doolittle, and T. M. Embley. 1999. Microsporidia are related to fungi: evidence from the largest subunit of RNA polymerase II and other proteins. Proc. Natl. Acad. Sci. USA 96:580585
Jones, D. T., W. R. Taylor, and J. M. Thornton. 1992. The rapid generation of mutation data matrices from protein sequences. CABIOS 8:275282
Katiyar, S. K., and T. D. Edlind. 1994. Beta-tubulin genes of Trichomonas vaginalis. Mol. Biochem. Parasitol. 64:33 42[ISI][Medline]
. 1996. Entamoeba histolytica encodes a highly divergent beta-tubulin. J. Eukaryot. Microbiol. 43:3134[ISI][Medline]
Keeling, P. J., J. A. Deane, C. Hink-Schauer, S. E. Douglas, U.-G. Maier, and G. I. McFadden. 1999. The secondary endosymbiont of the cryptomonad Guillardia theta contains alpha-, beta-, and gamma-tubulin genes. Mol. Biol. Evol. 16:13081313[Abstract]
Keeling, P. J., J. A. Deane, and G. I. McFadden. 1998. The phylogenetic position of alpha- and beta-tubulins from the Chlorarachnion host and Cercomonas (Cercozoa). J. Eukaryot. Microbiol. 45:561570[ISI][Medline]
Keeling, P. J., and W. F. Doolittle. 1996. Alpha-tubulin from early-diverging eukaryotic lineages and the evolution of the tubulin family. Mol. Biol. Evol. 13:12971305
Keeling, P. J., M. A. Luker, and J. D. Palmer. 2000. Evidence from beta-tubulin phylogeny that microsporidia evolved from within the fungi. Mol. Biol. Evol. 17:2331
Keeling, P. J., and G. I. McFadden. 1998. Origins of microsporidia. Trends Microbiol. 6:1923[ISI][Medline]
Kirk-Mason, K. E., M. J. Turner, and P. R. Chakraborty. 1988. Cloning and sequence of beta tubulin cDNA from Giardia lamblia. Nucleic Acids Res. 16:2733
Kube-Granderath, E., and M. Schliwa. 1998. Unusual tubulins for unusual cells. Protist 149:123126
Lang, B. F., G. Burger, C. J. O'Kelly, R. Cedergren, G. B. Lemieux, D. Sankoff, M. Turmel, and M. W. Gray. 1997. An ancestral mitochondrial DNA resembling a eubacterial genome in miniature. Nature 387:493497
Lang, B. F., E. Seif, M. W. Gray, C. J. O'Kelly, and G. Burger. 1999. A comparative genomics approach to the evolution of eukaryotes and their mitochondria. J. Eukaryot. Microbiol. 46:320326[ISI][Medline]
Leipe, D. D., J. H. Gunderson, T. A. Nerad, and M. L. Sogin. 1993. Small subunit ribosomal RNA+ of Hexamita inflata and the quest for the first branch in the eukaryotic tree. Mol. Biochem. Parasitol. 59:4148[ISI][Medline]
Loomis, W. F., and D. W. Smith. 1990. Molecular phylogeny of Dictyostelium discoideum by protein sequence comparison. Proc. Natl. Acad. Sci. USA 8790931097
O'Kelly, C. J. 1993. The jakobid flagellates: structural features of Jakoba, Reclinomonas and Histiona and implications for the early diversification of eukaryotes. J. Eukaryot. Microbiol. 40:627636[ISI]
O'Kelly, C. 1997. Ultrastructure of trophozoites, zoospores and cysts of Reclinomonas americana Flavin and Nerad, 1993 (Protista incertae sedis: Histionidae). Eur. J. Protistol. 33:337348[ISI]
O'Kelly, C. J., M. A. Farmer, and T. A. Nerad. 1999. Ultrastructure of Trimastix pyriformis (Klebs) Bernard et al.: similarities of Trimastix species with retortamonad and jakobid flagellates. Protist 150:14962
O'Kelly, C. J., and T. A. Nerad. 1999. Malawimonas jakobiformis n. gen., n. sp. (Malawimonadidae n. fam.): a Jakoba-like heterotrophic nanoflagellate with discoidal mitochondrial cristae. J. Eukaryot. Microbiol. 46:522531
Patterson, D. J. 1990. Jakoba libera (Ruinen 1938), a heterotrophic flagellate from deep oceanic sediments. J. Eukaryot. Microbiol. 40:172179
. 1994. Progress in protozoology. IX International Congress of Protozoology, Stuttgart, Germany
Philippe, H., and A. Adoutte. 1998. The molecular phylogeny of protozoa: solid facts and uncertainties. Pp. 2552 in G. H. Coombs, K. Vickerman, M. A. Sleigh, and A. Warren, eds. Evolutionary relationships among Protozoa. Kluwer Academic Publishers, London
Philippe, H., and A. Germot. 2000. Phylogeny of eukaryotes based on ribosomal RNA: long-branch attraction and models of sequence evolution. Mol. Biol. Evol. 17:830834
Roger, A. J. 1996. Studies on the phylogeny and gene structure of early-branching eukaryotes. Ph.D. thesis, Dalhousie University, Hallifax, Nova Scotia, Canada
. 1999. Reconstructing early events in eukaryotic evolution. Am. Nat. 154(Suppl.):S146S163
Roger, A. J., H. G. Morrison, and M. L. Sogin. 1999. Primary structure and phylogenetic relationships of a malate dehydrogenase gene from Giardia lamblia. J. Mol. Evol. 48:750755[ISI][Medline]
Roger, A. J., O. Sandblom, W. F. Doolittle, and H. Philippe. 1999. An evaluation of elongation factor 1a as a phylogenetic marker for eukaryotes. Mol. Biol. Evol. 16:218 233[Abstract]
Roger, A. J., M. W. Smith, R. F. Doolittle, and W. F. Doolittle. 1996. Evidence for the Heterolobosea from phylogenetic analysis of genes encoding glyceraldehyde-3-phosphate dehydrogenase. J. Eukaryot. Microbiol. 43:475485[ISI][Medline]
Roger, A. J., S. G. Svard, J. Tovar, C. G. Clark, M. W. Smith, F. D. Gillin, and M. L. Sogin. 1998. A mitochondrial-like chaperonin 60 gene in Giardia lamblia: evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria. Proc. Natl. Acad. Sci. USA 95:229234
Simpson, A. G. B., and D. J. Patterson. 1999. The ultrastructure of Carpediemonas membranifera: (Eukaryota), with reference to the excavate hypothesis. Eur. J. Protistol. 35: 353370
Simpson, A. G. B., C. Bernard, and D. J. Patterson. 2000. The ultrastructure of Trimastix marina Kent, 1880 (Eukaryota), an excavate flagellate. Europ. J. Protistol. 36:229 251[ISI]
Sleigh, M. 1989. Protozoa and other protists. 2nd edition. Edward Arnold, London
Sogin, M. L. 1991. Early evolution and the origin of eukaryotes. Curr. Opin. Genet. Dev. 1:457463[Medline]
Sogin, M. L., J. H. Gunderson, H. J. Elwood, R. A. Alonso, and D. A. Peattie. 1989. Phylogenetic meaning of the kingdom concept: an unusual ribosomal RNA from Giardia lamblia. Science 243:7577
Stiller, J. W., and B. J. Hall. 1999. Long-branch attraction and the rDNA model of early eukaryotic evolution. Mol. Biol. Evol. 16:12701279
Strimmer, K., and A. von Haeseler. 1996. Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies. Mol. Biol. Evol 13:964969
Swofford, D. L. 1999. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sinauer, Sunderland, Mass
Taylor, F. J. R. 1976. Flagellate phylogeny: a study in conflicts. J. Protozool. 23:2840[ISI]
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:46734680[Abstract]
Van de Peer, Y., A. Ben Ali, and A. Meyer. 2000. Microsporidia: accumulating molecular evidence that a group of amitochondriate and suspectedly primitive eukaryotes are just curious fungi. Gene 246:18
Van de Peer, Y., S. A. Rensing, U. G. Maier, and R. De Wachter. 1996. Substitution rate calibration of small subunit ribosomal RNA identifies chlorarachniophyte endosymbionts as remnants of green algae. Proc. Natl. Acad. Sci. USA 93:77327736
Voight, M. 1901. Mitteilungen aus der Biolog. Station Plön, Holstein-Uber einige bisher unbekannte Süsswasserorganismen. Zool. Az. 24:191195
Wainright, P. O., G. Hinkle, M. L. Sogin, and S. K. Stickel. 1993. Monophyletic origins of the metazoa: an evolutionary link with fungi. Science 260:340342
Weiss, L. M., T. D. Edlind, C. R. Vossbrinck, and T. Hashimoto. 1999. Microsporidian molecular phylogeny: the fungal connection. J. Eukaryot. Microbiol. 46:17S18S