A comparative genomics approach to the evolution of eukaryotes and their mitochondria.

The Organelle Genome Megasequencing Program (OGMP) investigates mitochondrial genome diversity and evolution by systematically determining the complete mitochondrial DNA (mtDNA) sequences of a phylogenetically broad selection of protists. The mtDNAs of lower fungi and choanoflagellates are being analyzed by the Fungal Mitochondrial Genome Project (FMGP), a sister project to the OGMP. Some of the most interesting protists include the jakobid flagellates Reclinomonas americana, Malawimonas jakobiformis, and Jakoba libera, which share ultrastructural similarities with amitochondriate retortamonads, and harbor mitochondrial genes not seen before in mtDNAs of other organisms. In R. americana and J. libera, gene clusters are found that resemble, to an unprecedented degree, the contiguous ribosomal protein operons str, S10, spc, and alpha of eubacteria. In addition, their mtDNAs code for an RNase P RNA that displays all the elements of a bacterial minimum consensus structure. This structure has been instrumental in detecting the rnpB gene in additional protists. Gene repertoire and gene order comparisons as well as multiple-gene phylogenies support the view of a single endosymbiotic origin of mitochondria, whose closest extant relatives are Rickettsia-type alpha-Proteobacteria.     

Evolutionary relationships among "jakobid" flagellates as indicated by alpha- and beta-tubulin phylogenies

Jakobids are free-living, heterotrophic flagellates that might represent early-diverging mitochondrial protists. They share ultrastructural similarities with eukaryotes that occupy basal positions in molecular phylogenies, and their mitochondrial genome architecture is eubacterial-like, suggesting a close affinity with the ancestral alpha-proteobacterial symbiont that gave rise to mitochondria and hydrogenosomes. To elucidate relationships among jakobids and other early-diverging eukaryotic lineages, we characterized alpha- and beta-tubulin genes from four jakobids: Jakoba libera, Jakoba incarcerata, Reclinomonas americana (the "core jakobids"), and Malawimonas jakobiformis. These are the first reports of nuclear genes from these organisms. Phylogenies based on alpha-, beta-, and combined alpha- plus beta-tubulin protein data sets do not support the monophyly of the jakobids. While beta-tubulin and combined alpha- plus beta-tubulin phylogenies showed a sister group relationship between J. libera and R. americana, the two other jakobids, M. jakobiformis and J. incarcerata, had unclear affinities. In all three analyses, J. libera, R. americana, and M. jakobiformis emerged from within a well-supported large "plant-protist" clade that included plants, green algae, cryptophytes, stramenopiles, alveolates, Euglenozoa, Heterolobosea, and several other protist groups, but not animals, fungi, microsporidia, parabasalids, or diplomonads. A preferred branching order within the plant-protist clade was not identified, but there was a tendency for the J. libera-R. americana lineage to group with a clade made up of the heteroloboseid amoeboflagellates and euglenozoan protists. Jakoba incarcerata branched within the plant-protist clade in the beta- and the combined alpha- plus beta-tubulin phylogenies. In alpha- tubulin trees, J. incarcerata occupied an unresolved position, weakly grouping with the animal/fungal/microsporidian group or with amitochondriate parabasalid and diplomonad lineages, depending on the phylogenetic method employed. Tubulin gene phylogenies were in general agreement with mitochondrial gene phylogenies and ultrastructural data in indicating that the "jakobids" may be polyphyletic. Relationships with the putatively deep-branching amitochondriate diplomonads remain uncertain.     

Comprehensive multigene phylogenies of excavate protists reveal the evolutionary positions of "primitive" eukaryotes

Many of the protists thought to represent the deepest branches on the eukaryotic tree are assigned to a loose assemblage called the "excavates." This includes the mitochondrion-lacking diplomonads and parabasalids (e.g., Giardia and Trichomonas) and the jakobids (e.g., Reclinomonas). We report the first multigene phylogenetic analyses to include a comprehensive sampling of excavate groups (six nuclear-encoded protein-coding genes, nine of the 10 recognized excavate groups). Excavates coalesce into three clades with relatively strong maximum likelihood bootstrap support. Only the phylogenetic position of Malawimonas is uncertain. Diplomonads, parabasalids, and the free-living amitochondriate protist Carpediemonas are closely related to each other. Two other amitochondriate excavates, oxymonads and Trimastix, form the second monophyletic group. The third group is comprised of Euglenozoa (e.g., trypanosomes), Heterolobosea, and jakobids. Unexpectedly, jakobids appear to be specifically related to Heterolobosea. This tree topology calls into question the concept of Discicristata as a supergroup of eukaryotes united by discoidal mitochondrial cristae and makes it implausible that jakobids represent an independent early-diverging eukaryotic lineage. The close jakobids-Heterolobosea-Euglenozoa connection demands complex evolutionary scenarios to explain the transition between the presumed ancestral bacterial-type mitochondrial RNA polymerase found in jakobids and the phage-type protein in other eukaryotic lineages, including Euglenozoa and Heterolobosea.     

Mitochondrial-type hsp70 genes of the amitochondriate protists, Giardia intestinalis, Entamoeba histolytica and two microsporidians

Genes encoding putative mitochondrial-type heat shock protein 70 (mit-hsp70) were isolated and sequenced from amitochondriate protists, Giardia intestinalis, Entamoeba histolytica, and two microsporidians, Encephalitozoon hellem and Glugea plecoglossi. The deduced mit-hsp70 sequences were analyzed by sequence alignments and phylogenetic reconstructions. The mit-hsp70 sequence of these four amitochondriate protists were divergent from other mit-hsp70 sequences of mitochondriate eukaryotes. However, all of these sequences were clearly located within a eukaryotic mitochondrial clade in the tree including various type hsp70 sequences, supporting the emerging notion that none of these amitochondriate lineages are primitively amitochodrial, but lost their mitochondria secondarily in their evolutionary past.     

Oxymonads are closely related to the excavate taxon Trimastix

Despite intensive study in recent years, large-scale eukaryote phylogeny remains poorly resolved. This is particularly problematic among the groups considered to be potential early branches. In many recent systematic schemes for early eukaryotic evolution, the amitochondriate protists oxymonads and Trimastix have figured prominently, having been suggested as members of many of the putative deep-branching higher taxa. However, they have never before been proposed as close relatives of each other. We amplified, cloned, and sequenced small-subunit ribosomal RNA genes from the oxymonad Pyrsonympha and from several Trimastix isolates. Rigorous phylogenetic analyses indicate that these two protist groups are sister taxa and are not clearly related to any currently established eukaryotic lineages. This surprising result has important implications for our understanding of cellular evolution and high-level eukaryotic phylogeny. Given that Trimastix contains small, electron-dense bodies strongly suspected to be derived mitochondria, this study constitutes the best evidence to date that oxymonads are not primitively amitochondriate. Instead, Trimastix and oxymonads may be useful organisms for investigations into the evolution of the secondary amitochondriate condition. All higher taxa involving either oxymonads or Trimastix may require modification or abandonment. Affected groups include four contemporary taxa given the rank of phylum (Metamonada, Loukozoa, Trichozoa, Percolozoa), and the informal excavate taxa. A new "phylum-level" taxon may be warranted for oxymonads and Trimastix.     

Malawimonas jakobiformis n. gen., n. sp. (Malawimonadidae n. fam.): A Jakoba-like Heterotrophic Nanoflagellate with Discoidal Mitochondrial Cristae

Malawimonas jakobiformis n. gen., n. sp., is established for a bacterivorous heterotrophic nanoflagellate isolated from the Malawi shore of Lake Nyasa (eastern Africa). Trophic stages observed were anteriorly biflagellate and naked. The posterior flagellum of a trophic cell resided in a conspicuous groove on the ventral surface, and bore a prominent vane. A Golgi stack and a mitochondrion with discoidal cristae were present anterior to the nucleus. The kinetid consisted of two short, slightly separated basal bodies, four microtubular roots, and associated fibers and bands. The three microtubular roots associated with the posterior basal body were associated with the ventral groove, while the single root associated with the anterior basal body gave rise to secondary cytoskeletal microtubules. Dividing cells became rounded, with persistent flagella. Cysts were uninucleate, and had thin organic walls without clearly differentiated apertures or ornamentation but with conspicuous attachment pads. Kinetid elements were present within cysts. On the basis of microscopical features, especially those of the kinetid, the nearest relatives of M. jakobiformis are the mitochondriate “jakobid” protists (families Histionidae and Jakobidae) and the amitochondriate retortamonads. Malawimonadidae n. fam. is established to accommodate this species.     

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.     

Identification and characterization of ectosymbionts of distinct lineages in Bacteroidales attached to flagellated protists in the gut of termites and a wood-feeding cockroach

Bacterial attachments to nearly the entire surface of flagellated protists in the guts of termites and the wood-feeding cockroach Cryptocercus are often observed. Based on the polymerase chain reaction-amplified 16S rRNA gene sequences, we investigated the phylogenetic relationships of the rod-shaped, attached bacteria (ectosymbionts) of several protist species from five host taxa and confirmed their identity by fluorescence in situ hybridizations. These ectosymbionts are affiliated with the order Bacteroidales but formed three distinct lineages, each of which may represent novel bacterial genera. One lineage consisted of the closely related ectosymbionts of two species of the protist genus Devescovina (Cristamonadida). The second lineage comprised three phylotypes identified from the protist Streblomastix sp. (Oxymonadida). The third lineage included ectosymbionts of the three protist genera Hoplonympha, Barbulanympha and Urinympha in the family Hoplonymphidae (Trichonymphida). The ultrastructural observations indicated that these rod-shaped ectosymbionts share morphological similarities of their cell walls and their point of attachment with the protist but differ in shape. Elongated forms of the ectosymbionts appeared in all the three lineages. The protist cells Streblomastix sp. and Hoplonympha sp. display deep furrows and vane-like structures, but these impressive structures are probably evolutionarily convergent because both the host protists and their ectosymbionts are distantly related.     

The microsporidian Encephalitozoon

Encephalitozoon cuniculi is an attractive model system for amitochondriate intracellular eukaryotic parasites. It is characterized by a very small genome (below 3 Mbp) and a unique invasion apparatus. Furthermore, as an infectious agent, it is important in human and veterinary medicine. The compactness of its genome involves the reduction of rDNA sequences as well as of some protein-coding genes and intergenic regions. Its highly differentiated apparatus to penetrate the host cell, an extrusome-like polar tube, is composed of novel proteins and may permit various pathways of infestation. Completion of the systematic E. cuniculi sequencing project should provide an important reference system for the comparative genomics of amitochondriate and mitochondriate parasites. Further analysis of orphan genes should help to identify factors that are responsible for its intracellular parasitic way of life.     

A high density of ancient spliceosomal introns in oxymonad excavates

Background  

Certain eukaryotic genomes, such as those of the amitochondriate parasites Giardia and Trichomonas, have very low intron densities, so low that canonical spliceosomal introns have only recently been discovered through genome sequencing. These organisms were formerly thought to be ancient eukaryotes that diverged before introns originated, or at least became common. Now however, they are thought to be members of a supergroup known as excavates, whose members generally appear to have low densities of canonical introns. Here we have used environmental expressed sequence tag (EST) sequencing to identify 17 genes from the uncultivable oxymonad Streblomastix strix, to survey intron densities in this most poorly studied excavate group.     

Moramonas marocensis gen. nov., sp. nov.: a jakobid flagellate isolated from desert soil with a bacteria-like,but bloated mitochondrial genome

A new jakobid genus has been isolated from Moroccan desert soil. The cyst-forming protist Moramonas marocensis gen. nov., sp. nov. has two anteriorly inserted flagella of which one points to the posterior cell pole accompanying the ventral feeding groove and is equipped with a dorsal vane—a feature typical for the Jakobida. It further shows a flagellar root system consisting of singlet microtubular root, left root (R1), right root (R2) and typical fibres associated with R1 and R2. The affiliation of M. marocensis to the Jakobida was confirmed by molecular phylogenetic analyses of the SSU rRNA gene, five nuclear genes and 66 mitochondrial protein-coding genes. The mitochondrial genome has the high number of genes typical for jakobids, and bacterial features, such as the four-subunit RNA polymerase and Shine–Dalgarno sequences upstream of the coding regions of several genes. The M. marocensis mitochondrial genome encodes a similar number of genes as other jakobids, but is unique in its very large genome size (greater than 264 kbp), which is three to four times higher than that of any other jakobid species investigated yet. This increase seems to be due to a massive expansion in non-coding DNA, creating a bloated genome like those of plant mitochondria.     

Fructose-1,6-bisphosphate aldolases in amitochondriate protists constitute a single protein subfamily with eubacterial relationships

Sequences of putative fructose-1,6-bisphospate aldolases (FBA) in five amitochondriate unicellular eukaryotes, the diplomonads Giardia intestinalis (published earlier) and Spironucleus barkhanus, the pelobiont Mastigamoeba balamuthi,the entamoebid Entamoeba histolytica, and the parabasalid Trichomonas vaginalis all belong to Class II of FBAs and are highly similar to each other (>48% amino acid identity). The five protist sequences, however, do not form a monophyletic group. Diplomonad FBAs share a most recent common ancestor, while FBAs of the three other protist species are part of a lineage that also includes sequences from a few eubacteria (Clostridium difficile, Treponema pallidum, Chlorobium tepidum). Both clades are part of the Type B of Class II aldolases, a complex that contains at least three additional lineages (subgroups) of enzymes. Type B enzymes are distant from Type A Class II aldolases, which consists of a number of bacterial and fungal enzymes and also contains the cytosolic FBA of Euglena gracilis. Class II aldolases are not homologous to Class I enzymes, to which animal and plant enzymes belong. The results indicate that amitochondriate protists acquired their FBAs from separate and different sources, involving lateral gene transfer from eubacteria, than did all other eukaryotes studied so far and underscore the complex composition of the glycolytic machinery in unicellular eukaryotes.     

Analyses of RNA Polymerase II genes from free-living protists: phylogeny,long branch attraction,and the eukaryotic big bang

The phylogenetic relationships among major eukaryotic protist lineages are largely uncertain. Two significant obstacles in reconstructing eukaryotic phylogeny are long-branch attraction (LBA) effects and poor taxon sampling of free-living protists. We have obtained and analyzed gene sequences encoding the largest subunit of RNA Polymerase II (RPB1) from Naegleria gruberi (a heterolobosean), Cercomonas ATCC 50319 (a cercozoan), and Ochromonas danica (a heterokont); we have also analyzed the RPB1 gene from the nucleomorph (nm) genome of Guillardia theta (a cryptomonad). Using a variety of phylogenetic methods our analysis shows that RPB1s from Giardia intestinalis and Trichomonas vaginalis are probably subject to intense LBA effects. Thus, the deep branching of these taxa on RPB1 trees is questionable and should not be interpreted as evidence favoring their early divergence. Similar effects are discernable, to a lesser extent, with the Mastigamoeba invertens RPB1 sequence. Upon removal of the outgroup and these problematic sequences, analyses of the remaining RPB1s indicate some resolution among major eukaryotic groups. The most robustly supported higher-level clades are the opisthokonts (animals plus fungi) and the red algae plus the cryptomonad nm-the latter result gives added support to the red algal origin of cryptomonad chloroplasts. Clades comprising Dictyostelium discoideum plus Acanthamoeba castellanii (Amoebozoa) and Ochromonas plus Plasmodium falciparum (chromalveolates) are consistently observed and moderately supported. The clades supported by our RPB1 analyses are congruent with other data, suggesting that bona fide phylogenetic relationships are being resolved. Thus, the RPB1 gene has apparently retained some phylogenetically meaningful signal, making it worthwhile to obtain sequences from more diverse protist taxa. Additional RPB1 data, especially in combination with other genes, should provide further resolution of branching orders among protist groups within the apparently rapid early divergence of eukaryotes.     

The origin of introns and their role in eukaryogenesis: a compromise solution to the introns-early versus introns-late debate?

Background  

Ever since the discovery of 'genes in pieces' and mRNA splicing in eukaryotes, origin and evolution of spliceosomal introns have been considered within the conceptual framework of the 'introns early' versus 'introns late' debate. The 'introns early' hypothesis, which is closely linked to the so-called exon theory of gene evolution, posits that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. Under this scenario, the absence of spliceosomal introns in prokaryotes is considered to be a result of "genome streamlining". The 'introns late' hypothesis counters that spliceosomal introns emerged only in eukaryotes, and moreover, have been inserted into protein-coding genes continuously throughout the evolution of eukaryotes. Beyond the formal dilemma, the more substantial side of this debate has to do with possible roles of introns in the evolution of eukaryotes.     

On Core Jakobids and Excavate Taxa: The Ultrastructure of Jakoba incarcerata

The cellular organisation of the 'excavate' flagellate Jakoba incarcerata Bernard, Simpson and Patterson 2000 is described. Cells have one nucleus and dictyosome. The putative mitochondria lack cristae. Two flagella (anterior and posterior) insert anterior to the feeding groove. The posterior flagellum bears a dorsal vane. An 'anterior' microtubular root arises against the anterior basal body. Two main microtubular roots, left and right, and a singlet 'root' arise around the posterior basal body and support the groove. Non-microtubular fibres termed 'A', 'B', 'I', and 'composite' associate with the right root. A multilaminar 'C' fibre associates with the left root. The cytoskeleton of J. incarcerata indicates a common ancestry with other excavate taxa (i.e. diplomonads, retortamonads, heteroloboseids, 'core jakobids', Malawimonas, Carpediemonas, and Trimastix). Overall, J. incarcerata is most similar to (other) core jakobids, namely Jakoba libera, Reclinomonas, and Histiona. We regard J. incarcerata as a core jakobid and identify the group by the synapomorphy 'vanes restricted to dorsal side of the posterior flagellum'. The anterior root and position of the B fibre (and presence of dense inclusions in the cartwheels and a conscpicuous singlet root-associated fibre) in J. incarcerata are novel for core jakobids and argue for close relationships with Trimastix and/or Heterolobosea. The C fibre is similar in substructure to the costal fibre of parabasalids and it is possible that the structures are homologous.     

Molecular evolution of calmodulin-like domain protein kinases (CDPKs) in plants and protists

Many genes for calmodulin-like domain protein kinases (CDPKs) have been identified in plants and Alveolate protists. To study the molecular evolution of the CDPK gene family, we performed a phylogenetic analysis of CDPK genomic sequences. Analysis of introns supports the phylogenetic analysis; CDPK genes with similar intron/exon structure are grouped together on the phylogenetic tree. Conserved introns support a monophyletic origin for plant CDPKs, CDPK-related kinases, and phosphoenolpyruvate carboxylase kinases. Plant CDPKs divide into two major branches. Plant CDPK genes on one branch share common intron positions with protist CDPK genes. The introns shared between protist and plant CDPKs presumably originated before the divergence of plants from Alveolates. Additionally, the calmodulin-like domains of protist CDPKs have intron positions in common with animal and fungal calmodulin genes. These results, together with the presence of a highly conserved phase zero intron located precisely at the beginning of the calmodulin-like domain, suggest that the ancestral CDPK gene could have originated from the fusion of protein kinase and calmodulin genes facilitated by recombination of ancient introns. Received: 11 July 2000 / Accepted: 18 April 2001     

Origin and evolution of eukaryotic chaperonins: phylogenetic evidence for ancient duplications in CCT genes

Chaperonins are oligomeric protein-folding complexes which are divided into two distantly related structural classes. Group I chaperonins (called GroEL/cpn60/hsp60) are found in bacteria and eukaryotic organelles, while group II chaperonins are present in archaea and the cytoplasm of eukaryotes (called CCT/TriC). While archaea possess one to three chaperonin subunit-encoding genes, eight distinct CCT gene families (paralogs) have been characterized in eukaryotes. We are interested in determining when during eukaryotic evolution the multiple gene duplications producing the CCT subunits occurred. We describe the sequence and phylogenetic analysis of five CCT genes from TRICHOMONAS: vaginalis and seven from GIARDIA: lamblia, representatives of amitochondriate protist lineages thought to have diverged early from other eukaryotes. Our data show that the gene duplications producing the eight CCT paralogs took place prior to the organismal divergence of TRICHOMONAS: and GIARDIA: from other eukaryotes. Thus, these divergent protists likely possess completely hetero-oligomeric CCT complexes like those in yeast and mammalian cells. No close phylogenetic relationship between the archaeal chaperonins and specific CCT subunits was observed, suggesting that none of the CCT gene duplications predate the divergence of archaea and eukaryotes. The duplications producing the CCTdelta and CCTepsilon subunits, as well as CCTalpha, CCTbeta, and CCTeta, are the most recent in the CCT gene family. Our analyses show significant differences in the rates of evolution of archaeal chaperonins compared with the eukaryotic CCTs, as well as among the different CCT subunits themselves. We discuss these results in light of current views on the origin, evolution, and function of CCT complexes.     

Association of intron phases with conservation at splice site sequences and evolution of spliceosomal introns.

How exon-intron structures of eukaryotic genes evolved under various evolutionary forces remains unknown. The phases of spliceosomal introns (the placement of introns with respect to reading frame) provide an opportunity to approach this question. When a large number of nuclear introns in protein-coding genes were analyzed, it was found that most introns were of phase 0, which keeps codons intact. We found that the phase distribution of spliceosomal introns is strongly correlated with the sequence conservation of splice signals in exons; the relatively underrepresented phase 2 introns are associated with the lowest conservation, the relatively overrepresented phase 0 introns display the highest conservation, and phase 1 introns are intermediate. Given the detrimental effect of mutations in exon sequences near splice sites as found in molecular experiments, the underrepresentation of phase 2 introns may be the result of deleterious-mutation-driven intron loss, suggesting a possible genetic mechanism for the evolution of intron-exon structures.     

Not all are free‐living: high‐throughput DNA metabarcoding reveals a diverse community of protists parasitizing soil metazoa

Protists, the most diverse eukaryotes, are largely considered to be free‐living bacterivores, but vast numbers of taxa are known to parasitize plants or animals. High‐throughput sequencing (HTS) approaches now commonly replace cultivation‐based approaches in studying soil protists, but insights into common biases associated with this method are limited to aquatic taxa and samples. We created a mock community of common free‐living soil protists (amoebae, flagellates, ciliates), extracted DNA and amplified it in the presence of metazoan DNA using 454 HTS. We aimed at evaluating whether HTS quantitatively reveals true relative abundances of soil protists and at investigating whether the expected protist community structure is altered by the co‐amplification of metazoan‐associated protist taxa. Indeed, HTS revealed fundamentally different protist communities from those expected. Ciliate sequences were highly over‐represented, while those of most amoebae and flagellates were under‐represented or totally absent. These results underpin the biases introduced by HTS that prevent reliable quantitative estimations of free‐living protist communities. Furthermore, we detected a wide range of nonadded protist taxa probably introduced along with metazoan DNA, which altered the protist community structure. Among those, 20 taxa most closely resembled parasitic, often pathogenic taxa. Therewith, we provide the first HTS data in support of classical observational studies that showed that potential protist parasites are hosted by soil metazoa. Taken together, profound differences in amplification success between protist taxa and an inevitable co‐extraction of protist taxa parasitizing soil metazoa obscure the true diversity of free‐living soil protist communities.