Protein phylogenies robustly resolve the deep-level relationships within Euglenozoa

The deepest-level relationships amongst Euglenozoa remain poorly resolved, despite a rich history of morphological examination and numerous molecular phylogenetic studies of small subunit ribosomal RNA (SSU rRNA) data. We address this question using two nuclear-encoded proteins, the cytosolic isoforms of heat shock protein 90 (hsp90) and heat shock protein 70 (hsp70). For both proteins we examined sequences from the three primary groups within Euglenozoa (euglenids, diplonemids, and kinetoplastids), and from their close relatives, Heterolobosea. Maximum likelihood (ML) and ML distance analyses of these proteins support a close relationship between diplonemids and kinetoplastids to the exclusion of the euglenid Euglena gracilis. In hsp90 and combined protein analyses bootstrap support is very strong and alternative topologies are generally rejected by 'approximately unbiased' (AU) tests. This result is consistent with recent molecular biological and morphological data, but contradicts early structural accounts and many SSU rRNA analyses that favour a closer relationship between diplonemids and euglenids. However, a re-examination of an important SSU rRNA data set highlights the instability of the inferences from this marker. The protein analyses also suggest that bodonids are paraphyletic, with trypanosomatids grouping with 'clade 2' and 'clade 3' bodonids to the exclusion of 'clade 1' bodonids.     

PHYLOGENY OF PHAGOTROPHIC EUGLENIDS (EUGLENOZOA): A MOLECULAR APPROACH BASED ON CULTURE MATERIAL AND ENVIRONMENTAL SAMPLES1

Molecular studies based on small subunit (SSU) rDNA sequences addressing euglenid phylogeny hitherto suffered from the lack of available data about phagotrophic species. To extend the taxon sampling, SSU rRNA genes from species of seven genera of phagotrophic euglenids were investigated. Sequence analyses revealed an increasing genetic diversity among euglenid SSU rDNA sequences compared with other well‐known eukaryotic groups, reflecting an equally broad diversity of morphological characters among euglenid phagotrophs. Phylogenetic inference using standard parsimony and likelihood approaches as well as Bayesian inference and spectral analyses revealed no clear support for euglenid monophyly. Among phagotrophs, monophyly of Petalomonas cantuscygni and Notosolenus ostium, both comprising simple ingestion apparatuses, is strongly supported. A moderately supported clade comprises phototrophic euglenids and primary osmotrophic euglenids together with phagotrophs, exhibiting a primarily flexible pellicle composed of numerous helically arranged strips and a complex ingestion apparatus with two supporting rods and four curved vanes. Comparison of molecular and morphological data is used to demonstrate the difficulties to formulate a hypothesis about how the ingestion apparatus evolved in this group.     

Ribosomal RNA phylogeny of bodonid and diplonemid flagellates and the evolution of euglenozoa

Euglenozoa is a major phylum of excavate protozoa (comprising euglenoids, kinetoplastids, and diplonemids) with highly unusual nuclear, mitochondrial, and chloroplast genomes. To improve understanding of euglenozoan evolution, we sequenced nuclear small-subunit rRNA genes from 34 bodonids (Bodo, Neobodo, Parabodo, Dimastigella-like, Rhynchobodo, Rhynchomonas, and unidentified strains), nine diplonemids (Diplonema, Rhynchopus), and a euglenoid (Entosiphon). Phylogenetic analysis reveals that diplonemids and bodonids are more diverse than previously recognised, but does not clearly establish the branching order of kinetoplastids, euglenoids, and diplonemids. Rhynchopus is holophyletic; parasitic species arose from within free-living species. Kinetoplastea (bodonids and trypanosomatids) are robustly holophyletic and comprise a major clade including all trypanosomatids and most bodonids ('core bodonids') and a very divergent minor one including Ichthyobodo. The root of the major kinetoplastid clade is probably between trypanosomatids and core bodonids. Core bodonids have three distinct subclades. Clade 1 has two distinct Rhynchobodo-like lineages; a lineage comprising Dimastigella and Rhynchomonas; and another including Cruzella and Neobodo. Clade 2 comprises Cryptobia/ Trypanoplasma, Procryptobia, and Parabodo. Clade 3 is an extensive Bodo saltans species complex. Neobodo designis is a vast genetically divergent species complex with mutually exclusive marine and freshwater subclades. Our analysis supports three phagotrophic euglenoid orders: Petalomonadida (holophyletic), Ploeotiida (probably holophyletic), Peranemida (paraphyletic).     

Phylogenetic affinities of Diplonema within the Euglenozoa as inferred from the SSU rRNA gene and partial COI protein sequences

In order to shed light on the phylogenetic position of diplonemids within the phylum Euglenozoa, we have sequenced small subunit rRNA (SSU rRNA) genes from Diplonema (syn. Isonema) papillatum and Diplonema sp. We have also analyzed a partial sequence of the mitochondrial gene for cytochrome c oxidase subunit I from D. papillatum. With both markers, the maximum likelihood method favored a closer grouping of diplonemids with kinetoplastids, while the parsimony and distance suggested a closer relationship of diplonemids with euglenoids. In each case, the differences between the best tree and the alternative trees were small. The frequency of codon usage in the partial D. papillatum COI was different from both related groups; however, as is the case in kinetoplastids but not in Euglena, both the non-canonical UGA codon and the canonical UGG codon were used to encode tryptophan in Diplonema.     

Dinoflagellate,Euglenid, or Cercomonad? The Ultrastructure and Molecular Phylogenetic Position of Protaspis grandis n. sp.

Protaspis is an enigmatic genus of marine phagotrophic biflagellates that have been tentatively classified with several different groups of eukaryotes, including dinoflagellates, euglenids, and cercomonads. This uncertainty led us to investigate the phylogenetic position of Protaspis grandis n. sp. with ultrastructural and small subunit (SSU) rDNA sequence data. Our results demonstrated that the cells were dorsoventrally flattened, shaped like elongated ovals with parallel lateral sides, 32.5-55.0 mum long and 20.0-35.0 mum wide. Moreover, two heterodynamic flagella emerged through funnels that were positioned subapically, each within a depression and separated by a distinctive protrusion. A complex multilayered wall surrounded the cell. Like dinoflagellates and euglenids, the nucleus contained permanently condensed chromosomes and a large nucleolus throughout the cell cycle. Pseudopodia containing numerous mitochondria with tubular cristae emerged from a ventral furrow through a longitudinal slit that was positioned posterior to the protrusion and flagellar apparatus. Batteries of extrusomes were present within the cytoplasm and had ejection sites through pores in the cell wall. The SSU rDNA phylogeny demonstrated a very close relationship between the benthic P. grandis n. sp. and the planktonic Cryothecomonas longipes. These ultrastructural and molecular phylogenetic data for Protaspis indicated that the current taxonomy of Protaspis and Crythecomonas is in need of re-evaluation. The composition and identity of Protaspis is reviewed and suggestions for future taxonomic changes are presented. Problems within the genus Cryothecomonas are highlighted as well, and the missing data needed to resolve ambiguities between the two genera are clarified.     

Morphological convergence characterizes the evolution of Xanthophyceae (Heterokontophyta): evidence from nuclear SSU rDNA and plastidial rbcL genes

Xanthophyceae are a group of heterokontophyte algae. Few molecular studies have investigated the evolutionary history and phylogenetic relationships of this class. We sequenced the nuclear-encoded SSU rDNA and chloroplast-encoded rbcL genes of several xanthophycean species from different orders, families, and genera. Neither SSU rDNA nor rbcL genes show intraspecific sequence variation and are good diagnostic markers for characterization of problematic species. New sequences, combined with those previously available, were used to create different multiple alignments. Analyses included sequences from 26 species of Xanthophyceae plus three Phaeothamniophyceae and two Phaeophyceae taxa used as outgroups. Phylogenetic analyses were performed according to Bayesian inference, maximum likelihood, and maximum parsimony methods. We explored effects produced on the phylogenetic outcomes by both taxon sampling as well as selected genes. Congruent results were obtained from analyses performed on single gene multiple alignments as well as on a data set including both SSU rDNA and rbcL sequences. Trees obtained in this study show that several currently recognized xanthophycean taxa do not form monophyletic groups. The order Mischococcales is paraphyletic, while Tribonematales and Botrydiales are polyphyletic even if evidence for the second order is not conclusive. Botrydiales and Vaucheriales, both including siphonous taxa, do not form a clade. The families Botrydiopsidaceae, Botryochloridaceae, and Pleurochloridaceae as well as the genera Botrydiopsis and Chlorellidium are polyphyletic. The Centritractaceae and the genus Bumilleriopsis also appear to be polyphyletic but their monophyly cannot be completely rejected with current evidence. Our results support morphological convergence at any taxonomic rank in the evolution of the Xanthophyceae. Finally, our phylogenetic analyses exclude an origin of the Xanthophyceae from a Vaucheria-like ancestor and favor a single early origin of the coccoid cell form.     

Phylogenetic relationships among Syndermata inferred from nuclear and mitochondrial gene sequences

Phylogenetic relationships among Syndermata have been extensively debated, mainly because the sister-group of the Acanthocephala has not yet been clearly identified from analyses of morphological and molecular data. Here we conduct phylogenetic analyses on samples from the 4 classes of Acanthocephala (Archiacanthocephala, Eoacanthocephala, Polyacanthocephala, and Palaeacanthocephala) and the 3 Rotifera classes (Bdelloidea, Monogononta, and Seisonidea). We do so using small-subunit (SSU) and large-subunit (LSU) ribosomal DNA and cytochrome c oxidase subunit 1 (cox 1) sequences. These nuclear and mitochondrial DNA sequences were obtained for 27 acanthocephalans, 9 rotifers, and representatives of 6 phyla that were used as outgroups. Maximum parsimony (MP), maximum likelihood (ML), and Bayesian analyses were conducted on the nuclear rDNA(SSU+LSU) and the combined sequence dataset(SSU+LSU+cox 1 genes). Phylogenetic analyses of the combined rDNA and cox 1 data uniformly provided strong support for a clade including rotifers plus acanthocephalans (Syndermata). Strong support was also found for monophyly of Acanthocephala in analyses of the combined dataset or rDNA sequences alone. Within the Acanthocephala the monophyletic grouping of the representatives of each class was strongly supported. Our results depicted Archiacanthocephala as the sister-group to the remaining acanthocephalans. Analyses of the combined dataset recovered a sister-group relationship between Acanthocephala and Bdelloidea by parsimony, likelihood, and Bayesian methods. Support for this clade was generally strong. Alternative topologies that depicted a different rotifer sister-group of Acanthocephala (or monophyly of Rotifera) were significantly worse. In this paraphyletic assemblage of rotifers, the relative positions of Seisonidea and Monogononta to the clade Bdelloidea+Acanthocephala were inconsistent among trees based on different inference methods. These results indicate that Bdelloidea is the free-living sister-group to acanthocephalans, which should prove key for comparative investigations of the morphological, molecular, and ecological changes accompanying the evolution of parasitism.     

New taxa refresh the phylogeny and classification of pleurostomatid ciliates (Ciliophora,Litostomatea)

A high diversity of pleurostomatid ciliates has been discovered in the last decade, and their systematics needs to be improved in the light of new findings concerning their morphology and molecular phylogeny. In this work, a new genus, Protolitonotus gen. n., and two new species, Protolitonotus magnus sp. n. and Protolitonotus longus sp. n., were studied. Furthermore, 19 novel nucleotide sequences of SSU rDNA, LSU rDNA and ITS1‐5.8S‐ITS2 were collected to determine the phylogenetic relationships and systematic positions of the pleurostomatid ciliates in this study. Based on both molecular and morphological data, the results demonstrated that: (i) as disclosed by the sequence analysis of SSU rDNA, LSU rDNA and ITS1‐5.8S‐ITS2, Protolitonotus gen. n. is sister to all other pleurostomatids and thus represents an independent lineage and a separate family, Protolitonotidae fam. n., which is defined by the presence of a semi‐suture formed by the right somatic kineties near the dorsal margin of the body; (ii) the families Litonotidae and Kentrophyllidae are both monophyletic based on both SSU rDNA and LSU rDNA sequences, whereas Amphileptidae are non‐monophyletic in trees inferred from SSU rDNA sequences; and (iii) the genera Loxophyllum and Kentrophyllum are both monophyletic, whereas Litonotus is non‐monophyletic based on SSU rDNA analyses. ITS1‐5.8S‐ITS2 sequence data were used for the phylogenetic analyses of pleurostomatids for the first time; however, species relationships were less well resolved than in the SSU rDNA and LSU rDNA trees. In addition, a major revision to the classification of the order Pleurostomatida is suggested and a key to its families and genera is provided.     

Phylogeny of phagotrophic euglenids (Euglenozoa) as inferred from hsp90 gene sequences

Molecular phylogenies of euglenids are usually based on ribosomal RNA genes that do not resolve the branching order among the deeper lineages. We addressed deep euglenid phylogeny using the cytosolic form of the heat-shock protein 90 gene (hsp90), which has already been employed with some success in other groups of euglenozoans and eukaryotes in general. Hsp90 sequences were generated from three taxa of euglenids representing different degrees of ultrastructural complexity, namely Petalomonas cantuscygni and wild isolates of Entosiphon sulcatum, and Peranema trichophorum. The hsp90 gene sequence of P. trichophorum contained three short introns (ranging from 27 to 31 bp), two of which had non-canonical borders GG-GG and GG-TG and two 10-bp inverted repeats, suggesting a structure similar to that of the non-canonical introns described in Euglena gracilis. Phylogenetic analyses confirmed a closer relationship between kinetoplastids and diplonemids than to euglenids, and supported previous views regarding the branching order among primarily bacteriovorous, primarily eukaryovorous, and photosynthetic euglenids. The position of P. cantuscygni within Euglenozoa, as well as the relative support for the nodes including it were strongly dependent on outgroup selection. The results were most consistent when the jakobid Reclinomonas americana was used as the outgroup. The most robust phylogenies place P. cantuscygni as the most basal branch within the euglenid clade. However, the presence of a kinetoplast-like mitochondrial inclusion in P. cantuscygni deviates from the currently accepted apomorphy-based definition of the kinetoplastid clade and highlights the necessity of detailed studies addressing the molecular nature of the euglenid and diplonemid mitochondrial genome.     

Unique mitochondrial genome structure in diplonemids, the sister group of kinetoplastids

Kinetoplastid flagellates are characterized by uniquely massed mitochondrial DNAs (mtDNAs), the kinetoplasts. Kinetoplastids of the trypanosomatid group possess two types of mtDNA molecules: maxicircles bearing protein and mitoribosomal genes and minicircles specifying guide RNAs, which mediate uridine insertion/deletion RNA editing. These circles are interlocked with one another to form dense networks. Whether these peculiar mtDNA features are restricted to kinetoplastids or prevail throughout Euglenozoa (euglenids, diplonemids, and kinetoplastids) is unknown. Here, we describe the mitochondrial genome and the mitochondrial ultrastructure of Diplonema papillatum, a member of the diplonemid flagellates, the sister group of kinetoplastids. Fluorescence and electron microscopy show a single mitochondrion per cell with an ultrastructure atypical for Euglenozoa. In addition, DNA is evenly distributed throughout the organelle rather than compacted. Molecular and electron microscopy studies distinguish numerous 6- and 7-kbp-sized mitochondrial chromosomes of monomeric circular topology and relaxed conformation in vivo. Remarkably, the cox1 gene (and probably other mitochondrial genes) is fragmented, with separate gene pieces encoded on different chromosomes. Generation of the contiguous cox1 mRNA requires trans-splicing, the precise mechanism of which remains to be determined. Taken together, the mitochondrial gene/genome structure of Diplonema is not only different from that of kinetoplastids but unique among eukaryotes as a whole.     

Species boundaries in gregarine apicomplexan parasites: a case study-comparison of morphometric and molecular variability in Lecudina cf. tuzetae (Eugregarinorida, Lecudinidae)

Trophozoites of gregarine apicomplexans are large feeding cells with diverse morphologies that have played a prominent role in gregarine systematics. The range of variability in trophozoite shapes and sizes can be very high even within a single species depending on developmental stages and host environmental conditions; this makes the delimitation of different species of gregarines based on morphological criteria alone very difficult. Accordingly, comparisons of morphological variability and molecular variability in gregarines are necessary to provide a pragmatic framework for establishing species boundaries within this diverse and poorly understood group of parasites. We investigated the morphological and molecular variability present in the gregarine Lecudina cf. tuzetae from the intestines of Nereis vexillosa (Polychaeta) collected in two different locations in Canada. Three distinct morphotypes of trophozoites were identified and the small subunit (SSU) rDNA was sequenced either from multicell isolates of the same morphotype or from single cells. The aim of this investigation was to determine whether the different morphotypes and localities reflected phylogenetic relatedness as inferred from the SSU rDNA sequence data. Phylogenetic analyses of the SSU rDNA demonstrated that the new sequences did not cluster according to morphotype or locality and instead were intermingled within a strongly supported clade. A comparison of 1,657 bp from 45 new sequences demonstrated divergences between 0% and 3.9%. These data suggest that it is necessary to acquire both morphological and molecular data in order to effectively delimit the "clouds" of variation associated with each gregarine species and to unambiguously reidentify these species in the future.     

DNA barcoding in autotrophic euglenids: evaluation of COI and 18s rDNA

Autotrophic euglenids (Euglenophyceae) are a common and abundant group of microbial eukaryotes in freshwater habitats. They have a limited number of features, which can be observed using light microscopy, thus species identification is often problematic. Establishing a barcode for this group is therefore an important step toward the molecular identification of autotrophic euglenids. Based on the literature, we selected verified species and used a plethora of available methods to validate two molecular markers: COI and 18S rDNA (the whole sequence and three fragments separately) as potential DNA barcodes. Analyses of the COI gene were performed based on the data set of 43 sequences (42 obtained in this study) representing 24 species and the COI gene was discarded as a DNA barcode mainly due to a lack of universal primer sites. For 18S rDNA analyses we used a data set containing 263 sequences belonging to 86 taxonomically verified species. We demonstrated that the whole 18S rDNA is too long to be a useful marker, but from the three shorter analyzed variable regions we recommend variable regions V2V3 and V4 of 18S rDNA as autotrophic euglenid barcodes due to their high efficiency (above 95% and 90%, respectively).     

Identity and systematic position of Paradinium poucheti and other Paradinium-like parasites of marine copepods based on morphology and nuclear-encoded SSU rDNA

Paradinium and Paradinium-like parasites were detected in various copepod hosts collected in the NW Mediterranean Sea, the North Atlantic Ocean, and the Godth?bsfjord (Greenland). The identity and systematic position of the parasitic, plasmodial protist Paradinium was investigated on the basis of SSU rDNA and morphology. SSU rDNA sequences were obtained from 3 specimens of Paradinium poucheti isolated from their cyclopoid copepod host, Oithona similis. In addition, a comparable sequence was obtained from a hitherto undescribed species of Paradinium from the harpactacoid copepod Euterpina acutifrons. Finally, SSU rDNA sequences were acquired from 2 specimens of a red plasmodial parasite (RP parasite) isolated from Clausocalanus sp. Both morphological and SSU rDNA sequence data supported that P. poucheti and Paradinium sp. are closely related organisms. In phylogenetic analyses based on SSU rDNA sequences, Paradinium spp. clustered with sequences from an uncultured eukaryote clone from the Pacific Ocean and two sequences from haplosporidian-like parasites of shrimps, Pandalus spp. This Paradinium clade branched as a sister group to a clade comprising the Haplosporidia and the Foraminifera. The RP parasite had a superficial morphological resemblance to Paradinium and has previously been interpreted as a member of this genus. However, several morphological characters contradict this and SSU rDNA sequence data disagree with the RP parasite and Paradinium being related. The phylogenetic analyses suggested that the RP parasite is a fast-evolved alveolate and a member of the so-called marine alveolate Group I (MAGI) and emerging data now suggest that this enigmatic group may, like the syndinian dinoflagellates, consist of heterotrophic parasites.     

Combining Molecular Data with Classical Morphology for Uncultured Phagotrophic Euglenids (Excavata): A Single‐Cell Approach

Phagotrophic euglenids are one of the most diverse and important forms of heterotrophic flagellates in sediment systems, and are key to understanding the evolution of photosynthetic euglenids and ‘primary osmotrophs’, yet relatively little is known about their biodiversity and phylogenetic relationships. A wealth of light microscopy‐based information is available, but little progress has been made in associating this with molecular sequence data. We established a protocol to obtain light microscopy data and molecular data from single euglenid cells isolated from environmental samples. Individual cells from freshwater and marine benthic samples were isolated and rinsed by micropipetting, documented using high‐resolution photomicroscopy, then subjected to single‐cell nested PCR using taxon‐specific primers in combination with universal eukaryotic primers, generating > 75% or full‐length SSU rDNA sequences. As a proof‐of‐principle eight individuals were characterised and subjected to phylogenetic analyses. Many of these cells were identified as Anisonema or Dinema, and grouped with existing sequences assigned to these taxa, and with a ‘Peranema sp.’ sequence that we could now clearly demonstrate was misidentified or misannotated. Another cell is Heteronema c.f. exaratum, the first ‘skidding heteronemid’ for which sequence data are available. This is not closely related to Heteronema scaphurum, and intriguingly, branches as the sister group to primary osmotrophs. A cell similar to Ploeotia vitrea (the type of this genus), shows no particular phylogenetic affinity to Ploeotia costata, the best studied Ploeotia species. Our experimental protocol provides a useful starting point for future analyses on euglenid biodiversity (including environmental sequence surveys), and their evolution and systematics.     

EVOLUTION OF PHACUS (EUGLENOPHYCEAE) AS INFERRED FROM PELLICLE MORPHOLOGY AND SSU rDNA

This research integrates a large morphological data set into a molecular context. Nineteen pellicle characters and 62 states from 13 euglenid taxa were analyzed cladistically. The pellicle morphology of Euglena tripteris (Klebs), Lepocinclis ovata (Conrad), Phacus brachykentron (Pochmann), P. oscillans (Klebs), P. pyrum (Stein), and P. triqueter (Dujardin) is described comprehensively. These data are compared with new information on the pellicle morphology of Euglena acus (Ehrenberg), E. stellata (Mainx), and Peranema trichophorum (Stein) in addition to published data on Entosiphon sulcatum (Dujardin), Euglena gracilis (Klebs), Distigma proteus (Pringsheim), and Petalomonas cantuscygni (Cann and Pennick). Nuclear small subunit (SSU) rDNA sequences provided an independent test for establishing a robust organismal pedigree of the same taxa. A synthetic tree derived from the combined phylogenetic analyses of pellicle morphology and SSU rDNA enabled us to parsimoniously map morphological character states. This approach demonstrated the utility of pellicle morphology for inferring phylogenetic relationships of euglenids and establishing apomorphy-based clade definitions. Three robust clades with unambiguous pellicle-based apomorphies can be recognized within taxa traditionally classified as Phacus : (1) L. ovata and P. pyrum , (2) E. tripteris and P. triqueter , and (3) P. brachykentron and P. oscillans. Taxonomic concerns that emerged from these results are discussed.     

Compartmentalization of a glycolytic enzyme in Diplonema, a non-kinetoplastid euglenozoan

Glycosomes are peroxisome-related organelles containing glycolytic enzymes that have been found only in kinetoplastids. We show here that a glycolytic enzyme is compartmentalized in diplonemids, the sister group of kinetoplastids. We found that, similar to kinetoplastid aldolases, the fructose 1,6-bisphosphate aldolase of Diplonema papillatum possesses a type 2-peroxisomal targeting signal. Western blotting showed that this aldolase was present predominantly in the membrane/organellar fraction. Immunofluorescence analysis showed that this aldolase had a scattered distribution in the cytosol, suggesting its compartmentalization. In contrast, orotidine-5'-monophosphate decarboxylase, a non-glycolytic glycosomal enzyme in kinetoplastids, was shown to be a cytosolic enzyme in D. papillatum. Since euglenoids, the earliest diverging branch of Euglenozoa, do not possess glycolytic compartments, these findings suggest that the routing of glycolytic enzymes into peroxisomes may have occurred in a common ancestor of diplonemids and kinetoplastids, followed by diversification of these newly established organelles in each of these euglenozoan lineages.     

Higher-level phylogeny of Foraminifera inferred from the RNA polymerase II (RPB1) gene

Macroevolutionary relations among main lineages of Foraminifera have traditionally been inferred from the small subunit ribosomal genes (SSU rDNA). However, important discrepancies in the rates of SSU rDNA evolution between major lineages led to difficulties in accurate interpretation of SSU-based phylogenetic reconstructions. Recently, actin and beta-tubulin sequences have been used as alternative markers of foraminiferal phylogeny and their analyses globally confirm results obtained with SSU rDNA. In order to test new protein markers, we sequenced a fragment of the largest subunit of the RNA polymerase II (RPB1), a nuclear encoded single copy gene, for 8 foraminiferal species representing major orders of Foraminifera. Analyses of our data robustly confirm previous SSU rDNA and actin phylogenies and show (i) the paraphyly and ancestral position of monothalamid Foraminifera; (ii) the independent origin of miliolids; (iii) the monophyly of rotaliids, including buliminids and globigerinids; and (iv) the polyphyly of planktonic families Globigerinidae and Candeinidae. Additionally, the RPB1 phylogeny suggests Allogromiidae as the most ancestral foraminiferal lineage. In the light of our study, RPB1 appears as a valuable phylogenetic marker, particularly useful for groups of protists showing extreme variations of evolutionary rates in ribosomal genes.