skip to main content

Title: Dynamics and interactions of highly resolved marine plankton via automated high-frequency sampling

Short timescale observations are valuable for understanding microbial ecological processes. We assessed dynamics in relative abundance and potential activities by sequencing the small sub-unit ribosomal RNA gene (rRNA gene) and rRNA molecules (rRNA) of Bacteria, Archaea, and Eukaryota once to twice daily between March 2014 and May 2014 from the surface ocean off Catalina Island, California. Typically Ostreococcus, Braarudosphaera, Teleaulax, and Synechococcus dominated phytoplankton sequences (including chloroplasts) while SAR11, Sulfitobacter, and Fluviicola dominated non-phytoplankton Bacteria and Archaea. We observed short-lived increases of diatoms, mostly Pseudo-nitzschia and Chaetoceros, with quickly responding Bacteria and Archaea including Flavobacteriaceae (Polaribacter & Formosa), Roseovarius, and Euryarchaeota (MGII), notably the exact amplicon sequence variants we observed responding similarly to another diatom bloom nearby, 3 years prior. We observed correlations representing known interactions among abundant phytoplankton rRNA sequences, demonstrating the biogeochemical and ecological relevance of such interactions: (1) The kleptochloroplastidic ciliate Mesodinium 18S rRNA gene sequences and a single Teleaulax taxon (via 16S rRNA gene sequences) were correlated (Spearman r = 0.83) yet uncorrelated to a Teleaulax 18S rRNA gene OTU, or any other taxon (consistent with a kleptochloroplastidic or karyokleptic relationship) and (2) the photosynthetic prymnesiophyte Braarudosphaera bigelowii and two strains of diazotrophic cyanobacterium UCYN-A were correlated and each taxon was also correlated to other taxa, including B. bigelowii to a verrucomicrobium and a dictyochophyte phytoplankter (all r > 0.8). We also report strong correlations (r > 0.7) between various ciliates, bacteria, and phytoplankton, suggesting interactions via currently unknown mechanisms. These data reiterate the utility of high-frequency time series to show rapid microbial reactions to stimuli, and provide new information about in situ dynamics of previously recognized and hypothesized interactions.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
The ISME Journal
Medium: X Size: p. 2417-2432
["p. 2417-2432"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Tropical environments with unique abiotic and biotic factors—such as salt ponds, mangroves, and coral reefs—are often in close proximity. The heterogeneity of these environments is reflected in community shifts over short distances, resulting in high biodiversity. While phytoplankton assemblages physically associated with corals, particularly their symbionts, are well studied, less is known about phytoplankton diversity across tropical aquatic environments. We assess shifts in phytoplankton community composition along inshore to offshore gradients by sequencing and analyzing 16S rRNA gene amplicons using primers targeting the V1-V2 region that capture plastids from eukaryotic phytoplankton and cyanobacteria, as well as heterotrophic bacteria. Microbial alpha diversity computed from 16S V1-V2 amplicon sequence variant (ASV) data from 282 samples collected in and around Curaçao, in the Southern Caribbean Sea, varied more within the dynamic salt ponds, salterns, and mangroves, compared to the seemingly stable above-reef, off-reef, and open sea environments. Among eukaryotic phytoplankton, stramenopiles often exhibited the highest relative abundances in mangrove, above-reef, off-reef, and open sea environments, where cyanobacteria also showed high relative abundances. Within stramenopiles, diatom amplicons dominated in salt ponds and mangroves, while dictyochophytes and pelagophytes prevailed above reefs and offshore. Green algae and cryptophytes were also present, and the former exhibited transitions following the gradient from inland to offshore. Chlorophytes and prasinophyte Class IV dominated in salt ponds, while prasinophyte Class II, including Micromonas commoda and Ostreococcus Clade OII, had the highest relative abundances of green algae in mangroves, above-reef, off-reef, and the open sea. To improve Class II prasinophyte classification, we sequenced 18S rRNA gene amplicons from the V4 region in 41 samples which were used to interrelate plastid-based results with information on uncultured prasinophyte species from prior 18S rRNA gene-based studies. This highlighted the presence of newly described Ostreococcus bengalensis and two Micromonas candidate species. Network analyses identified co-occurrence patterns between individual phytoplankton groups, including cyanobacteria, and heterotrophic bacteria. Our study reveals multiple uncultured and novel lineages within green algae and dictyochophytes in tropical marine habitats. Collectively, the algal diversity patterns and potential co-occurrence relationships observed in connection to physicochemical and spatial influences help provide a baseline against which future change can be assessed. 
    more » « less
  2. Abstract

    Microbial associations that result in phytoplankton mortality are important for carbon transport in the ocean. This includes parasitism, which in microbial food webs is dominated by the marine alveolate group, Syndiniales. Parasites are expected to contribute to carbon recycling via host lysis; however, knowledge on host dynamics and correlation to carbon export remain unclear and limit the inclusion of parasitism in biogeochemical models. We analyzed a 4-year 18S rRNA gene metabarcoding dataset (2016–19), performing network analysis for 12 discrete depths (1–1000 m) to determine Syndiniales–host associations in the seasonally oligotrophic Sargasso Sea. Analogous water column and sediment trap data were included to define environmental drivers of Syndiniales and their correlation with particulate carbon flux (150 m). Syndiniales accounted for 48–74% of network edges, most often associated with Dinophyceae and Arthropoda (mainly copepods) at the surface and Rhizaria (Polycystinea, Acantharea, and RAD-B) in the aphotic zone. Syndiniales were the only eukaryote group to be significantly (and negatively) correlated with particulate carbon flux, indicating their contribution to flux attenuation via remineralization. Examination of Syndiniales amplicons revealed a range of depth patterns, including specific ecological niches and vertical connection among a subset (19%) of the community, the latter implying sinking of parasites (infected hosts or spores) on particles. Our findings elevate the critical role of Syndiniales in marine microbial systems and reveal their potential use as biomarkers for carbon export.

    more » « less
  3. Gilbert, Jack A. (Ed.)
    ABSTRACT Small subunit rRNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly used primers with >300 million rRNA gene sequences retrieved from globally distributed marine metagenomes. The best-performing primers compared to 16S rRNA of bacteria and archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering cyanobacterial and chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA specific; 82%)—demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R 2 = 0.98, n  = 272), indicating amplicons can produce equally accurate community composition data compared with shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly completely universal primer set that could accurately quantify biogeochemically important taxa in ecosystems ranging from the deep sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data. IMPORTANCE PCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time but will work optimally only if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally distributed short-read oceanic metagenomes. Our results demonstrate that primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is theoretically possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides, and present a software pipeline that can guide rational design of primers for any environment with available meta’omic data. 
    more » « less
  4. Jay Gan (Ed.)
    Perfluoroalkyl substances (PFAS) are of great ecological concern, however, exploration of their impact on bacteria-phytoplankton consortia is limited. This study employed a bioassay approach to investigate the effect of unary exposures of increasing concentrations of PFAS (perfluorooctane sulfonate (PFOS) and 6:2 fluorotelomer sulfonate (6:2 FTS)) on microbial communities from the northwestern Gulf of Mexico. Each community was examined for changes in growth and photophysiology, exudate production and shifts in community structure (16S and 18S rRNA genes). 6:2 FTS did not alter the growth or health of phytoplankton communities, as there were no changes relative to the controls (no PFOS added). On the other hand, PFOS elicited significant phototoxicity (p < 0.05), altering PSII antennae size, lowering PSII connectivity, and decreasing photosynthetic efficiency over the incubation (four days). PFOS induced a cellular protective response, indicated by significant increases (p < 0.001) in the release of transparent exopolymer particles (TEP) compared to the control. Eukaryotic communities (18S rRNA gene) changed substantially (p < 0.05) and to a greater extent than prokaryotic communities (16S rRNA gene) in PFOS treatments. Community shifts were concentration-dependent for eukaryotes, with the low treatment (5 mg/L PFOS) dominated by Coscinodiscophyceae (40 %), and the high treatment (30 mg/L PFOS) marked by a Trebouxiophyceae (50 %) dominance. Prokaryotic community shifts were not concentration dependent, as both treatment levels became depleted in Cyanobacteriia and were dominated by members of the Bacteroidia, Gammaproteobacteria, and Alphaproteobacteria classes. Further, PFOS significantly decreased (p < 0.05) the Shannon diversity and Pielou’s evenness across treatments for eukaryotes, and in the low treatment (5 mg/L PFOS) for prokaryotes. These findings show that photophysiology was not impacted by 6:2 FTS but PFOS elicited toxicity that impacted photosynthesis, exudate release, and community composition. This research is crucial in understanding how PFOS impacts microbial communities. 
    more » « less
  5. Eutrophication and hypoxia markedly alter trophic dynamics and nutrient cycling in estuarine water columns, but little is known about the microbial communities that drive and interact with these changes. Here we studied microbial plankton (bacteria, archaea, protists and micro-metazoans) in a large temperate estuary where bottom hypoxia occurs every summer due to warmer temperatures, stratification, and oxidation of organic matter fueled by nutrient enrichment. We used high-throughput sequencing of the 16S and 18S rRNA genes (V4 region) and quantified multiple abiotic and biotic factors in surface and bottom waters during the summer of 2019. The conditions associated with the intensification of hypoxia in bottom waters as the summer progressed were linked to significant changes in the diversity, community structure and potential functioning of microbial communities. Under maximum hypoxia (dissolved oxygen concentration: 0.9-3.1 mg l -1 ), there were increased proportions of ammonia-oxidizing archaea (AOA), bacterivorous and parasitic protists, and copepod nauplii. Sequence proportions of AOA ( Nitrosopumilus) and nitrite-oxidizing bacteria ( Nitrospinaceae) were significantly correlated with the concentration of oxidized N species (nitrite plus nitrate, which peaked at 14.4 µM) and the proportions of nauplii DNA sequences and biomass. Our data support a tight coupling of biogeochemical and food web processes, with rapid oxidation of ammonia and accumulation of oxidized N species as hypoxia intensifies during the summer. 
    more » « less