Marine gregarine genomes reveal the breadth of apicomplexan diversity with a partially conserved glideosome machinery.

Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~ 9 Mb) and coding capacity (~ 5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sisters to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions such as gliding, well known to allow access to the intracellular parasitic lifestyle in Apicomplexa.

Parasitic infections by Group II Syndiniales target selected dinoflagellate host populations within diverse protist assemblages in a model coastal pond.

Protists are integral to marine food webs and biogeochemical cycles; however, there is a paucity of data describing specific ecological niches for some of the most abundant taxa in marker gene libraries. Syndiniales are one such group, often representing the majority of sequence reads recovered from picoplankton samples across the global ocean. However, the prevalence and impacts of syndinian parasitism in marine environments remain unclear. We began to address these critical knowledge gaps by generating a high-resolution time series (March-October 2018) in a productive coastal pond. Seasonal shifts in protist populations, including parasitic Syndiniales, were documented during periods of higher primary productivity and increased summer temperature-driven stratification. Elevated concentrations of infected hosts and free-living parasite spores occurred at nearly monthly intervals in July, August, and September. We suggest intensifying stratification during this period correlated with the increased prevalence of dinoflagellates that were parasitized by Group II Syndiniales. Infections in some protist populations were comparable to previously reported large single-taxon dinoflagellate blooms. Infection dynamics in Salt Pond demonstrated the propagation of syndinian parasites through mixed protist assemblages and highlighted patterns of host/parasite interactions that better reflect many other marine environments where single taxon blooms are uncommon.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp.

BACKGROUND: Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization. RESULTS: We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions. CONCLUSION: These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

EukRef: Phylogenetic curation of ribosomal RNA to enhance understanding of eukaryotic diversity and distribution.

Environmental sequencing has greatly expanded our knowledge of micro-eukaryotic diversity and ecology by revealing previously unknown lineages and their distribution. However, the value of these data is critically dependent on the quality of the reference databases used to assign an identity to environmental sequences. Existing databases contain errors and struggle to keep pace with rapidly changing eukaryotic taxonomy, the influx of novel diversity, and computational challenges related to assembling the high-quality alignments and trees needed for accurate characterization of lineage diversity. EukRef (eukref.org) is an ongoing community-driven initiative that addresses these challenges by bringing together taxonomists with expertise spanning the eukaryotic tree of life and microbial ecologists, who use environmental sequence data to develop reliable reference databases across the diversity of microbial eukaryotes. EukRef organizes and facilitates rigorous mining and annotation of sequence data by providing protocols, guidelines, and tools. The EukRef pipeline and tools allow users interested in a particular group of microbial eukaryotes to retrieve all sequences belonging to that group from International Nucleotide Sequence Database Collaboration (INSDC) (GenBank, the European Nucleotide Archive [ENA], or the DNA DataBank of Japan [DDBJ]), to place those sequences in a phylogenetic tree, and to curate taxonomic and environmental information for the group. We provide guidelines to facilitate the process and to standardize taxonomic annotations. The final outputs of this process are (1) a reference tree and alignment, (2) a reference sequence database, including taxonomic and environmental information, and (3) a list of putative chimeras and other artifactual sequences. These products will be useful for the broad community as they become publicly available (at eukref.org) and are shared with existing reference databases.

Unveiling protist diversity associated with the Pacific oyster Crassostrea gigas using blocking and excluding primers.

BACKGROUND: Microbiome of macroorganisms might directly or indirectly influence host development and homeostasis. Many studies focused on the diversity and distribution of prokaryotes within these assemblages, but the eukaryotic microbial compartment remains underexplored so far. RESULTS: To tackle this issue, we compared blocking and excluding primers to analyze microeukaryotic communities associated with Crassostrea gigas oysters. High-throughput sequencing of 18S rRNA genes variable loops revealed that excluding primers performed better by not amplifying oyster DNA, whereas the blocking primer did not totally prevent host contaminations. However, blocking and excluding primers showed similar pattern of alpha and beta diversities when protist communities were sequenced using metabarcoding. Alveolata, Stramenopiles and Archaeplastida were the main protist phyla associated with oysters. In particular, Codonellopsis, Cyclotella, Gymnodinium, Polarella, Trichodina, and Woloszynskia were the dominant genera. The potential pathogen Alexandrium was also found in high abundances within some samples. CONCLUSIONS: Our study revealed the main protist taxa within oysters as well as the occurrence of potential oyster pathogens. These new primer sets are promising tools to better understand oyster homeostasis and disease development, such as the Pacific Oyster Mortality Syndrome (POMS) targeting juveniles.

Identification to species level of live single microalgal cells from plankton samples with matrix-free laser/desorption ionization mass spectrometry.

INTRODUCTION: Marine planktonic communities are complex microbial consortia often dominated by microscopic algae. The taxonomic identification of individual phytoplankton cells usually relies on their morphology and demands expert knowledge. Recently, a live single-cell mass spectrometry (LSC-MS) pipeline was developed to generate metabolic profiles of microalgae. OBJECTIVE: Taxonomic identification of diverse microalgal single cells from collection strains and plankton samples based on the metabolic fingerprints analyzed with matrix-free laser desorption/ionization high-resolution mass spectrometry. METHODS: Matrix-free atmospheric pressure laser-desorption ionization mass spectrometry was performed to acquire single-cell mass spectra from collection strains and prior identified environmental isolates. The computational identification of microalgal species was performed by spectral pattern matching (SPM). Three similarity scores and a bootstrap-derived confidence score were evaluated in terms of their classification performance. The effects of high and low-mass resolutions on the classification success were evaluated. RESULTS: Several hundred single-cell mass spectra from nine genera and nine species of marine microalgae were obtained. SPM enabled the identification of single cells at the genus and species level with high accuracies. The receiver operating characteristic (ROC) curves indicated a good performance of the similarity measures but were outperformed by the bootstrap-derived confidence scores. CONCLUSION: This is the first study to solve taxonomic identification of microalgae based on the metabolic fingerprints of the individual cell using an SPM approach.

Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes.

This revision of the classification of eukaryotes follows that of Adl et al., 2012 [J. Euk. Microbiol. 59(5)] and retains an emphasis on protists. Changes since have improved the resolution of many nodes in phylogenetic analyses. For some clades even families are being clearly resolved. As we had predicted, environmental sampling in the intervening years has massively increased the genetic information at hand. Consequently, we have discovered novel clades, exciting new genera and uncovered a massive species level diversity beyond the morphological species descriptions. Several clades known from environmental samples only have now found their home. Sampling soils, deeper marine waters and the deep sea will continue to fill us with surprises. The main changes in this revision are the confirmation that eukaryotes form at least two domains, the loss of monophyly in the Excavata, robust support for the Haptista and Cryptista. We provide suggested primer sets for DNA sequences from environmental samples that are effective for each clade. We have provided a guide to trophic functional guilds in an appendix, to facilitate the interpretation of environmental samples, and a standardized taxonomic guide for East Asian users.

Evidence for parasite-mediated selection during short-lasting toxic algal blooms.

Parasites play a role in the control of transient algal blooms, but it is not known whether parasite-mediated selection results in coevolution of the host and the parasites over this short time span. We investigated the presence of coevolution between the toxic dinoflagellate Alexandrium minutum and two naturally occurring endoparasites during blooms lasting a month in two river estuaries, using cross-inoculation experiments across time and space. Higher parasite abundance was associated with a large daily reduction in relative A. minutum abundances, demonstrating strong parasite-mediated selection. There was genetic variability in infectivity in both parasite species, and in resistance in the host. We found no evidence for coevolution in one estuary; however, in the other estuary, we found high genetic diversity in the two parasite species, fluctuations in infectivity and suggestion that the two parasites are well adapted to their host, as in 'Red Queen' dynamics. Thus, coevolution is possible over the short time span of a bloom, but geographically variable, and may feedback on community dynamics.

Transcriptomic profiling of Alexandrium fundyense during physical interaction with or exposure to chemical signals from the parasite Amoebophrya.

Toxic microalgae have their own pathogens, and understanding the way in which these microalgae respond to antagonistic attacks may provide information about their capacity to persist during harmful algal bloom events. Here, we compared the effects of the physical presence of the parasite Amoebophrya sp. and exposure to waterborne cues from cultures infected with this parasite, on gene expression by the toxic dinoflagellates, Alexandrium fundyense. Compared with control samples, a total of 14,882 Alexandrium genes were differentially expressed over the whole-parasite infection cycle at three different time points (0, 6 and 96 h). RNA sequencing analyses indicated that exposure to the parasite and parasitic waterborne cues produced significant changes in the expression levels of Alexandrium genes associated with specific metabolic pathways. The observed upregulation of genes associated with glycolysis, the tricarboxylic acid cycle, fatty acid ß-oxidation, oxidative phosphorylation and photosynthesis suggests that parasite infection increases the energy demand of the host. The observed upregulation of genes correlated with signal transduction indicates that Alexandrium could be sensitized by parasite attacks. This response might prime the defence of the host, as indicated by the increased expression of several genes associated with defence and stress. Our findings provide a molecular overview of the response of a dinoflagellate to parasite infection.

Evolutionary processes and cellular functions underlying divergence in Alexandrium minutum.

Understanding divergence in the highly dispersive and seemingly homogeneous pelagic environment for organisms living as free drifters in the water column remains a challenge. Here, we analysed the transcriptome-wide mRNA sequences, as well as the morphology of 18 strains of Alexandrium minutum, a dinoflagellate responsible for harmful algal blooms worldwide, to investigate the functional bases of a divergence event. Analysis of the joint site frequency spectrum (JSFS) pointed towards an ancestral divergence in complete isolations followed by a secondary contact resulting in gene flow between the two diverging groups, but heterogeneous across sites. The sites displaying fixed SNPs were associated with a highly restricted gene flow and a strong overrepresentation of nonsynonymous polymorphism, suggesting the importance of selective pressures as drivers of the divergence. The most divergent transcripts were homologs to genes involved in calcium/potassium fluxes across the membrane, calcium transduction signal and saxitoxin production. The implication of these results in terms of ecological divergence and build-up of reproductive isolation is discussed. Dinoflagellates are especially difficult to study in the field at the ecological level due to their small size and the dynamic nature of their natural environment, but also at the genomic level due to their huge and complex genome and the absence of closely related model organism. This study illustrates the possibility to identify the traits of primary importance in ecology and evolution starting from high-throughput sequencing data, even for such organisms.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing.

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Genomic insights into processes driving the infection of Alexandrium tamarense by the Parasitoid Amoebophrya sp.

The regulatory circuits during infection of dinoflagellates by their parasites are largely unknown on the molecular level. Here we provide molecular insights into these infection dynamics. Alexandrium tamarense is one of the most prominent harmful algal bloom dinoflagellates. Its pathogen, the dinoflagellate parasitoid Amoebophrya sp., has been observed to infect and control the blooms of this species. We generated a data set of transcripts from three time points (0, 6, and 96 h) during the infection of this parasite-host system. Assembly of all transcript data from the parasitoid (>900,000 reads/313 Mbp with 454/Roche next-generation sequencing [NGS]) yielded 14,455 contigs, to which we mapped the raw transcript reads of each time point of the infection cycle. We show that particular surface lectins are expressed at the beginning of the infection cycle which likely mediate the attachment to the host cell. In a later phase, signal transduction-related genes together with transmembrane transport and cytoskeleton proteins point to a high integration of processes involved in host recognition, adhesion, and invasion. At the final maturation stage, cell division- and proliferation-related genes were highly expressed, reflecting the fast cell growth and nuclear division of the parasitoid. Our molecular insights into dinoflagellate parasitoid interactions point to general mechanisms also known from other eukaryotic parasites, especially from the Alveolata. These similarities indicate the presence of fundamental processes of parasitoid infection that have remained stable throughout evolution within different phyla.

The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy.

The interrogation of genetic markers in environmental meta-barcoding studies is currently seriously hindered by the lack of taxonomically curated reference data sets for the targeted genes. The Protist Ribosomal Reference database (PR(2), http://ssu-rrna.org/) provides a unique access to eukaryotic small sub-unit (SSU) ribosomal RNA and DNA sequences, with curated taxonomy. The database mainly consists of nuclear-encoded protistan sequences. However, metazoans, land plants, macrosporic fungi and eukaryotic organelles (mitochondrion, plastid and others) are also included because they are useful for the analysis of high-troughput sequencing data sets. Introns and putative chimeric sequences have been also carefully checked. Taxonomic assignation of sequences consists of eight unique taxonomic fields. In total, 136 866 sequences are nuclear encoded, 45 708 (36 501 mitochondrial and 9657 chloroplastic) are from organelles, the remaining being putative chimeric sequences. The website allows the users to download sequences from the entire and partial databases (including representative sequences after clustering at a given level of similarity). Different web tools also allow searches by sequence similarity. The presence of both rRNA and rDNA sequences, taking into account introns (crucial for eukaryotic sequences), a normalized eight terms ranked-taxonomy and updates of new GenBank releases were made possible by a long-term collaboration between experts in taxonomy and computer scientists.

Molecular diversity of the syndinean genus Euduboscquella based on single-cell PCR analysis.

The genus Euduboscquella is one of a few described genera within the syndinean dinoflagellates, an enigmatic lineage with abundant diversity in marine environmental clone libraries based on small subunit (SSU) rRNA. The region composed of the SSU through to the partial large subunit (LSU) rRNA was determined from 40 individual tintinnid ciliate loricae infected with Euduboscquella sampled from eight surface water sites in the Northern Hemisphere, producing seven distinct SSU sequences. The corresponding host SSU rRNA region was also amplified from eight host species. The SSU tree of Euduboscquella and syndinean group I sequences from environmental clones had seven well-supported clades and one poorly supported clade across data sets from 57 to 692 total sequences. The genus Euduboscquella consistently formed a supported monophyletic clade within a single subclade of group I sequences. For most parasites with identical SSU sequences, the more variable internal transcribed spacer (ITS) to LSU rRNA regions were polymorphic at 3 to 10 sites. However, in E. cachoni there was variation between ITS to LSU copies at up to 20 sites within an individual, while in a parasite of Tintinnopsis spp., variation between different individuals ranged up to 19 polymorphic sites. However, applying the compensatory base change model to the ITS2 sequences suggested no compensatory changes within or between individuals with the same SSU sequence, while one to four compensatory changes between individuals with similar but not identical SSU sequences were found. Comparisons between host and parasite phylogenies do not suggest a simple pattern of host or parasite specificity.

Temperature Affects the Biological Control of Dinoflagellates by the Generalist Parasitoid Parvilucifera rostrata.

The increase in emerging harmful algal blooms in the last decades has led to an extensive concern in understanding the mechanisms behind these events. In this paper, we assessed the growth of two blooming dinoflagellates (Alexandrium minutum and Heterocapsa triquetra) and their susceptibility to infection by the generalist parasitoid Parvilucifera rostrata under a temperature gradient. The growth of the two dinoflagellates differed across a range of temperatures representative of the Penzé Estuary (13 to 22 °C) in early summer. A. minutum growth increased across this range and was the highest at 19 and 22 °C, whereas H. triquetra growth was maximal at intermediate temperatures (15-18 °C). Interestingly, the effect of temperature on the parasitoid infectivity changed depending on which host dinoflagellate was infected with the dinoflagellate responses to temperature following a positive trend in A. minutum (higher infections at 20-22 °C) and a unimodal trend in H. triquetra (higher infections at 18 °C). Low temperatures negatively affected parasitoid infections in both hosts (i.e., "thermal refuge"). These results demonstrate how temperature shifts may not only affect bloom development in microalgal species but also their control by parasitoids.

Intracellular development and impact of a marine eukaryotic parasite on its zombified microalgal host.

Parasites are widespread and diverse in oceanic plankton and many of them infect single-celled algae for survival. How these parasites develop and scavenge energy within the host and how the cellular organization and metabolism of the host is altered remain open questions. Combining quantitative structural and chemical imaging with time-resolved transcriptomics, we unveil dramatic morphological and metabolic changes of the marine parasite Amoebophrya (Syndiniales) during intracellular infection, particularly following engulfment and digestion of nutrient-rich host chromosomes. Changes include a sequential acristate and cristate mitochondrion with a 200-fold increase in volume, a 13-fold increase in nucleus volume, development of Golgi apparatus and a metabolic switch from glycolysis (within the host) to TCA (free-living dinospore). Similar changes are seen in apicomplexan parasites, thus underlining convergent traits driven by metabolic constraints and the infection cycle. In the algal host, energy-producing organelles (plastid, mitochondria) remain relatively intact during most of the infection. We also observed that sugar reserves diminish while lipid droplets increase. Rapid infection of the host nucleus could be a "zombifying" strategy, allowing the parasite to digest nutrient-rich chromosomes and escape cytoplasmic defense, whilst benefiting from maintained carbon-energy production of the host cell.

Dinophyceae can use exudates as weapons against the parasite Amoebophrya sp. (Syndiniales).

Parasites in the genus Amoebophrya sp. infest dinoflagellate hosts in marine ecosystems and can be determining factors in the demise of blooms, including toxic red tides. These parasitic protists, however, rarely cause the total collapse of Dinophyceae blooms. Experimental addition of parasite-resistant Dinophyceae (Alexandrium minutum or Scrippsiella donghaienis) or exudates into a well-established host-parasite coculture (Scrippsiella acuminata-Amoebophrya sp.) mitigated parasite success and increased the survival of the sensitive host. This effect was mediated by waterborne molecules without the need for a physical contact. The strength of the parasite defenses varied between dinoflagellate species, and strains of A. minutum and was enhanced with increasing resistant host cell concentrations. The addition of resistant strains or exudates never prevented the parasite transmission entirely. Survival time of Amoebophrya sp. free-living stages (dinospores) decreased in presence of A. minutum but not of S. donghaienis. Parasite progeny drastically decreased with both species. Integrity of the dinospore membrane was altered by A. minutum, providing a first indication on the mode of action of anti-parasitic molecules. These results demonstrate that extracellular defenses can be an effective strategy against parasites that protects not only the resistant cells producing them, but also the surrounding community.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

Comparative analysis between protist communities from the deep-sea pelagic ecosystem and specific deep hydrothermal habitats.

Protist communities associated with deep seawater and bivalves from six hydrothermal sites in the Pacific Ocean were characterized by microscopy and molecular rRNA gene surveys (18S rRNA) and compared with planktonic communities from Pacific deep-pelagic seawater (from 500 to 3000 m in depth). Genetic libraries from larger size fractions (>3 µm) of deep-pelagic water were mainly dominated by Dinophyceae, whereas small size fractions (<3 µm) mainly revealed radiolarians and Syndiniales. In contrast, more specific opportunistic detritivores and grazers, mostly belonging to Stramenopiles and Cercozoa, were detected from water surrounding vent chimneys. Protist communities were different in the pallial cavity of the giant hydrothermal bivalves Bathymodiolus thermophilus and Calyptogena magnifica, dominated by Ciliophora (primarily belonging to Phyllopharyngea, Oligohymenophorea and Oligotrichea) and Cercozoa. Interestingly, protist communities retrieved from the pallial cavity liquid of hydrothermal bivalves were remarkably homogeneous along the Southern East Pacific Rise, in contrast to bivalves collected on the Mid-Atlantic Ridge hydrothermal vents and cold seeps from the Gulf of Mexico. Hence, complex protist communities seem to occur inside hydrothermal bivalves, and these metazoa may constitute a stable micro-niche for micro-eukaryotes, including grazers, detritivores, symbionts and potential parasites. From these communities, new lineages within the ciliates may emerge.