Protistan community composition and metabolism in the North Pacific Subtropical Gyre: Influences of mesoscale eddies and depth.

Marine protists and their metabolic activities are intricately tied to the cycling of nutrients and the flow of energy through microbial food webs. Physiochemical changes in the environment, such as those that result from mesoscale eddies, may impact protistan communities, but the effects that such changes have on protists are poorly known. A metatranscriptomic study was conducted to investigate how eddies affected protists at adjacent cyclonic and anticyclonic eddy sites in the oligotrophic ocean at four depths from 25 to 250 m. Eddy polarity impacted protists at all depths sampled, although the effects of eddy polarity were secondary to the impact of depth across the depth range. Eddy-induced vertical shifts in the water column yielded differences in the cyclonic and anticyclonic eddy protistan communities, and these differences were the most pronounced at and just below the deep chlorophyll maximum. An analysis of transcripts associated with protistan nutritional physiology at 150 m revealed that cyclonic eddies may support a more heterotrophic community, while anticyclonic eddies promote a more phototrophic community. The results of this study indicate that eddies alter the metabolism of protists particularly in the lower euphotic zone and may therefore impact carbon export from the euphotic zone.

Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities.

The 24-h cycle of light and darkness governs daily rhythms of complex behaviors across all domains of life. Intracellular photoreceptors sense specific wavelengths of light that can reset the internal circadian clock and/or elicit distinct phenotypic responses. In the surface ocean, microbial communities additionally modulate nonrhythmic changes in light quality and quantity as they are mixed to different depths. Here, we show that eukaryotic plankton in the North Pacific Subtropical Gyre transcribe genes encoding light-sensitive proteins that may serve as light-activated transcription factors, elicit light-driven electrical/chemical cascades, or initiate secondary messenger-signaling cascades. Overall, the protistan community relies on blue light-sensitive photoreceptors of the cryptochrome/photolyase family, and proteins containing the Light-Oxygen-Voltage (LOV) domain. The greatest diversification occurred within Haptophyta and photosynthetic stramenopiles where the LOV domain was combined with different DNA-binding domains and secondary signal-transduction motifs. Flagellated protists utilize green-light sensory rhodopsins and blue-light helmchromes, potentially underlying phototactic/photophobic and other behaviors toward specific wavelengths of light. Photoreceptors such as phytochromes appear to play minor roles in the North Pacific Subtropical Gyre. Transcript abundance of environmental light-sensitive protein-encoding genes that display diel patterns are found to primarily peak at dawn. The exceptions are the LOV-domain transcription factors with peaks in transcript abundances at different times and putative phototaxis photoreceptors transcribed throughout the day. Together, these data illustrate the diversity of light-sensitive proteins that may allow disparate groups of protists to respond to light and potentially synchronize patterns of growth, division, and mortality within the dynamic ocean environment.

The Mixoplankton Database (MDB): Diversity of photo-phago-trophic plankton in form, function, and distribution across the global ocean.

Protist plankton are major members of open-water marine food webs. Traditionally divided between phototrophic phytoplankton and phagotrophic zooplankton, recent research shows many actually combine phototrophy and phagotrophy in the one cell; these protists are the "mixoplankton." Under the mixoplankton paradigm, "phytoplankton" are incapable of phagotrophy (diatoms being exemplars), while "zooplankton" are incapable of phototrophy. This revision restructures marine food webs, from regional to global levels. Here, we present the first comprehensive database of marine mixoplankton, bringing together extant knowledge of the identity, allometry, physiology, and trophic interactivity of these organisms. This mixoplankton database (MDB) will aid researchers that confront difficulties in characterizing life traits of protist plankton, and it will benefit modelers needing to better appreciate ecology of these organisms with their complex functional and allometric predator-prey interactions. The MDB also identifies knowledge gaps, including the need to better understand, for different mixoplankton functional types, sources of nutrition (use of nitrate, prey types, and nutritional states), and to obtain vital rates (e.g. growth, photosynthesis, ingestion, factors affecting photo' vs. phago' -trophy). It is now possible to revisit and re-classify protistan "phytoplankton" and "zooplankton" in extant databases of plankton life forms so as to clarify their roles in marine ecosystems.

Daily dynamics of contrasting spring algal blooms in Santa Monica Bay (central Southern California Bight).

Protistan algae (phytoplankton) dominate coastal upwelling ecosystems where they form massive blooms that support the world's most important fisheries and constitute an important sink for atmospheric CO2 . Bloom initiation is well understood, but the biotic and abiotic forces that shape short-term dynamics in community composition are still poorly characterized. Here, high-frequency (daily) changes in relative abundance dynamics of the metabolically active protistan community were followed via expressed 18S V4 rRNA genes (RNA) throughout two algal blooms during the spring of 2018 and 2019 in Santa Monica Bay (central Southern California Bight). A diatom bloom formed after wind-driven, nutrient upwelling events in both years, but different taxa dominated each year. Whereas diatoms bloomed following elevated nutrients and declined after depletion each year, a massive dinoflagellate bloom manifested under relatively low inorganic nitrogen conditions following diatom bloom senescence in 2019 but not 2018. Network analysis revealed associations between diatoms and cercozoan putative parasitic taxa and syndinean parasites during 2019 that may have influenced the demise of the diatoms, and the transition to a dinoflagellate-dominated bloom.

Biological composition and microbial dynamics of sinking particulate organic matter at abyssal depths in the oligotrophic open ocean.

Sinking particles are a critical conduit for the export of organic material from surface waters to the deep ocean. Despite their importance in oceanic carbon cycling and export, little is known about the biotic composition, origins, and variability of sinking particles reaching abyssal depths. Here, we analyzed particle-associated nucleic acids captured and preserved in sediment traps at 4,000-m depth in the North Pacific Subtropical Gyre. Over the 9-month time-series, Bacteria dominated both the rRNA-gene and rRNA pools, followed by eukaryotes (protists and animals) and trace amounts of Archaea. Deep-sea piezophile-like Gammaproteobacteria, along with Epsilonproteobacteria, comprised >80% of the bacterial inventory. Protists (mostly Rhizaria, Syndinales, and ciliates) and metazoa (predominantly pelagic mollusks and cnidarians) were the most common sinking particle-associated eukaryotes. Some near-surface water-derived eukaryotes, especially Foraminifera, Radiolaria, and pteropods, varied greatly in their abundance patterns, presumably due to sporadic export events. The dominance of piezophile-like Gammaproteobacteria and Epsilonproteobacteria, along with the prevalence of their nitrogen cycling-associated gene transcripts, suggested a central role for these bacteria in the mineralization and biogeochemical transformation of sinking particulate organic matter in the deep ocean. Our data also reflected several different modes of particle export dynamics, including summer export, more stochastic inputs from the upper water column by protists and pteropods, and contributions from sinking mid- and deep-water organisms. In total, our observations revealed the variable and heterogeneous biological origins and microbial activities of sinking particles that connect their downward transport, transformation, and degradation to deep-sea biogeochemical processes.

SSU-rRNA Gene Sequencing Survey of Benthic Microbial Eukaryotes from Guaymas Basin Hydrothermal Vent.

Microbial eukaryotes have important roles in marine food webs, but their diversity and activities in hydrothermal vent ecosystems are poorly characterized. In this study, we analyzed microbial eukaryotic communities associated with bacterial (Beggiatoa) mats in the 2,000 m deep-sea Guaymas Basin hydrothermal vent system using 18S rRNA gene high-throughput sequencing of the V4 region. We detected 6,954 distinct Operational Taxonomic Units (OTUs) across various mat systems. Of the sequences that aligned with known protistan phylotypes, most were affiliated with alveolates (especially dinoflagellates and ciliates) and cercozoans. OTU richness and community structure differed among sediment habitats (e.g. different mat types and cold sediments away from mats). Additionally, full-length 18S rRNA genes amplified and cloned from single cells revealed the identities of some of the most commonly encountered, active ciliates in this hydrothermal vent ecosystem. Observations and experiments were also conducted to demonstrate that ciliates were trophically active and ingesting fluorescent bacteria or Beggiatoa trichomes. Our work suggests that the active and diverse protistan community at the Guaymas Basin hydrothermal vent ecosystem likely consumes substantial amounts of bacterial biomass, and that the different habitats, often defined by distances of just a few 10s of cm, select for particular assemblages and levels of diversity.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone.

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Ecological Patterns Among Bacteria and Microbial Eukaryotes Derived from Network Analyses in a Low-Salinity Lake.

Microbial communities are comprised of complex assemblages of highly interactive taxa. We employed network analyses to identify and describe microbial interactions and co-occurrence patterns between microbial eukaryotes and bacteria at two locations within a low salinity (0.5-3.5 ppt) lake over an annual cycle. We previously documented that the microbial diversity and community composition within Lake Texoma, southwest USA, were significantly affected by both seasonal forces and a site-specific bloom of the harmful alga, Prymnesium parvum. We used network analyses to answer ecological questions involving both the bacterial and microbial eukaryotic datasets and to infer ecological relationships within the microbial communities. Patterns of connectivity at both locations reflected the seasonality of the lake including a large rain disturbance in May, while a comparison of the communities between locations revealed a localized response to the algal bloom. A network built from shared nodes (microbial operational taxonomic units and environmental variables) and correlations identified conserved associations at both locations within the lake. Using network analyses, we were able to detect disturbance events, characterize the ecological extent of a harmful algal bloom, and infer ecological relationships not apparent from diversity statistics alone.

Effect of light and prey availability on gene expression of the mixotrophic chrysophyte, Ochromonas sp.

BACKGROUND: Ochromonas is a genus of mixotrophic chrysophytes that is found ubiquitously in many aquatic environments. Species in this genus can be important consumers of bacteria but vary in their ability to perform photosynthesis. We studied the effect of light and bacteria on growth and gene expression of a predominantly phagotrophic Ochromonas species. Axenic cultures of Ochromonas sp. were fed with heat-killed bacteria (HKB) and grown in constant light or darkness. RNA was extracted from cultures in the light or in the dark with HKB present (Light + HKB; Dark + HKB), and in the light after HKB were depleted (Light + depleted HKB). RESULTS: There were no significant differences in the growth or bacterial ingestion rates between algae grown in light or dark conditions. The availability of light led to a differential expression of only 8% of genes in the transcriptome. A number of genes associated with photosynthesis, phagotrophy, and tetrapyrrole synthesis was upregulated in the Light + HKB treatment compared to Dark + HKB. Conversely, the comparison between the Light + HKB and Light + depleted HKB treatments revealed that the presence of HKB led to differential expression of 59% of genes, including the majority of genes involved in major carbon and nitrogen metabolic pathways. Genes coding for unidirectional enzymes for the utilization of glucose were upregulated in the presence of HKB, implying increased glycolytic activities during phagotrophy. Algae without HKB upregulated their expression of genes coding for ammonium transporters, implying uptake of inorganic nitrogen from the culture medium when prey were unavailable. CONCLUSIONS: Transcriptomic results agreed with previous observations that light had minimal effect on the population growth of Ochromonas sp. However, light led to the upregulation of a number of phototrophy- and phagotrophy-related genes, while the availability of bacterial prey led to prominent changes in major carbon and nitrogen metabolic pathways. Our study demonstrated the potential of transcriptomic approaches to improve our understanding of the trophic physiologies of complex mixotrophs, and revealed responses in Ochromonas sp. not apparent from traditional culture studies.

UniEuk: Time to Speak a Common Language in Protistology!

Universal taxonomic frameworks have been critical tools to structure the fields of botany, zoology, mycology, and bacteriology as well as their large research communities. Animals, plants, and fungi have relatively solid, stable morpho-taxonomies built over the last three centuries, while bacteria have been classified for the last three decades under a coherent molecular taxonomic framework. By contrast, no such common language exists for microbial eukaryotes, even though environmental '-omics' surveys suggest that protists make up most of the organismal and genetic complexity of our planet's ecosystems! With the current deluge of eukaryotic meta-omics data, we urgently need to build up a universal eukaryotic taxonomy bridging the protist -omics age to the fragile, centuries-old body of classical knowledge that has effectively linked protist taxa to morphological, physiological, and ecological information. UniEuk is an open, inclusive, community-based and expert-driven international initiative to build a flexible, adaptive universal taxonomic framework for eukaryotes. It unites three complementary modules, EukRef, EukBank, and EukMap, which use phylogenetic markers, environmental metabarcoding surveys, and expert knowledge to inform the taxonomic framework. The UniEuk taxonomy is directly implemented in the European Nucleotide Archive at EMBL-EBI, ensuring its broad use and long-term preservation as a reference taxonomy for eukaryotes.

Estimating Protistan Diversity Using High-Throughput Sequencing.

Sequencing hypervariable regions from the 18S rRNA gene is commonly employed to characterize protistan biodiversity, yet there are concerns that short reads do not provide the same taxonomic resolution as full-length sequences. A total of 7,432 full-length sequences were used to perform an in silico analysis of how sequences of various lengths and target regions impact downstream ecological interpretations. Sequences that were longer than 400 nucleotides and included the V4 hypervariable region generated results similar to those derived from full-length 18S rRNA gene sequences. Present high-throughput sequencing capabilities are approaching protistan diversity estimation comparable to whole gene sequences.

Investigating microbial eukaryotic diversity from a global census: insights from a comparison of pyrotag and full-length sequences of 18S rRNA genes.

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.

Development, comparison, and validation using ELISAs for the determination of domoic acid in California sea lion body fluids.

Mortalities of California sea lions (Zalophus californianus) attributed to the neurotoxin domoic acid (DA) produced by the diatom Pseudo-nitzschia have occurred repeatedly along the U.S. west coast since the late 1990s. Quantifying the amount of DA in these animals and correlating this information with the presence of DA in phytoplankton and the local food web has become a research focus for many scientists. However, differences in materials, equipment, technical capability, budgets, and objectives of the various groups and/or agencies involved in this work have influenced the DA quantification platforms used. The goal of the present study was to compare the performance of two commercially available ELISAs for the determination of DA in a spectrum of California sea lion body fluids and to compare the results with LC/MS of the same samples. The results indicated differences among these approaches, presumably owing to matrix effects (particularly urine) and antibody reactivities. This information implies that care should be taken in attempting to compare datasets generated using different analytical platforms and interpreting the results of published studies.

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Seasonality and disturbance: annual pattern and response of the bacterial and microbial eukaryotic assemblages in a freshwater ecosystem.

High-throughput pyrosequencing of SSU rDNA genes was used to obtain monthly snapshots of eukaryotic and bacterial diversity and community structure at two locations in Lake Texoma, a low salinity lake in the south central United States, over 1 year. The lake experienced two disturbance events (i) a localized bloom of Prymnesium parvum restricted to one of the locations that lasted from January to April, and (ii) a large (17 cm), global rain event in the beginning of May, overlaid onto seasonal environmental change. Eukaryotic species richness as well as both eukaryotic and bacterial community similarity exhibited seasonal patterns, including distinct responses to the rain event. The P. parvum bloom created a natural experiment in which to directly explore the effects of an Ecosystem Disruptive Algal Bloom (EDAB) on the microbial community separated from seasonal changes. Microbial species richness was unaffected by the bloom, however, the eukaryotic community structure (evenness) and the patterns of both eukaryotic and bacterial community similarity at bloom and non-bloom sites were statistically distinct during the 4 months of the bloom. These results indicate that physical and biological disturbances as well as seasonal environmental forces contribute to the structure of both the eukaryotic and bacterial communities.

Towards a molecular taxonomy for protists: benefits, risks, and applications in plankton ecology.

The increasing use of genetic information for the development of methods to study the diversity, distributions, and activities of protists in nature has spawned a new generation of powerful tools. For ecologists, one lure of these approaches lies in the potential for DNA sequences to provide the only immediately obvious means of normalizing the diverse criteria that presently exist for identifying and counting protists. A single, molecular taxonomy would allow studies of diversity across a broad range of species, as well as the detection and quantification of particular species of interest within complex, natural assemblages; goals that are not feasible using traditional methods. However, these advantages are not without their potential pitfalls and problems. Conflicts involving the species concept, disagreements over the true (physiological/ecological) meaning of genetic diversity, and a perceived threat by some that sequence information will displace knowledge regarding the morphologies, functions and physiologies of protistan taxa, have created debate and doubt regarding the efficacy and appropriateness of some genetic approaches. These concerns need continued discussion and eventual resolution as we move toward the irresistible attraction, and potentially enormous benefits, of the application of genetic approaches to protistan ecology.

Diversity and distributional patterns of ciliates in Guaymas Basin hydrothermal vent sediments.

Little is known about protists at deep-sea hydrothermal vents. The vent sites at Guaymas Basin in the Gulf of California are characterized by dense mats of filamentous pigmented or nonpigmented Beggiatoa that serve as markers of subsurface thermochemical gradients. We constructed 18S rRNA libraries to investigate ciliate assemblages in Beggiatoa mats and from bare sediments at the Guaymas vent site. Results indicated a high diversity of ciliates, with 156 operational taxonomic units identified in 548 sequences. Comparison between mat environments demonstrated that ciliate and bacterial assemblages from pigmented mats, nonpigmented mats, and bare sediments were significantly different and highly correlated with bacterial assemblages. Neither bacterial nor ciliate assemblages were correlated with environmental factors. The most abundant ciliates at Guaymas were more likely to be represented in clone libraries from other hydrothermal, deep-sea, and/or anoxic or microaerophilic environments, supporting the hypothesis that these ciliate species are broadly distributed. The orange mat environment included a higher proportion of ciliate sequences that were more similar to those from other environmental studies than to cultured ciliate species, whereas clone libraries from bare sediments included sequences that were the most highly divergent from all other sequences and may represent species that are endemic to Guaymas.

The Effect of pH and Salinity on the Toxicity and Growth of the Golden Alga, Prymnesium parvum.

Bioassays using cultures of the toxic haptophyte Prymnesium parvum and the ciliate Cyclidium sp. as prey were conducted to test the effect of pH (range = 6.5 - 8.5), salinity (range = 1.50 - 7.50‰), and a combination of pH and salinity on the toxicity of P. parvum. pH had a significant effect on P. parvum toxicity. Toxicity was rapidly (within 24 hr) induced by increasing pH of the medium, or reduced by lowering pH. Conversely, lowering salinity reduced toxicity, albeit less effectively compared to pH, and P. parvum cells remained toxic at the lowest values tested (1.50‰ at pH 7.5). An additional effect between pH and salinity was also observed: low salinity combined with low pH led to not only decreased toxicity, but also resulted in lower P. parvum growth rates. Such effects of pH and salinity on P. parvum growth and toxicity provide insight into the environmental factors supporting community dominance and toxic blooms of the alga.