Tracking Vibrio: population dynamics and ecology of Vibrio parahaemolyticus and V. vulnificus in an Alabama estuary.

Vibrio is a genus of halophilic, gram-negative bacteria found in estuaries around the globe. Integral parts of coastal cultures often involve contact with vectors of pathogenic Vibrio spp. (e.g., consuming raw shellfish). High rates of mortality from certain Vibrio spp. infections demonstrate the need for an improved understanding of Vibrio spp. dynamics in estuarine regions. Our study assessed meteorological, hydrographic, and biological correlates of Vibrio parahaemolyticus and V. vulnificus at 10 sites in the Eastern Mississippi Sound System (EMSS) from April to October 2019. During the sampling period, median abundances of V. parahaemolyticus and V. vulnificus were 2.31 log MPN/L and 2.90 log MPN/L, respectively. Vibrio spp. dynamics were largely driven by site-based variation, with sites closest to freshwater inputs having the highest abundances. The E-W wind scalar, which affects Ekman transport, was a novel Vibrio spp. correlate observed. A potential salinity effect on bacterial-particle associations was identified, where V. vulnificus was associated with larger particles in conditions outside of their optimal salinity. Additionally, V. vulnificus abundances were correlated to those of harmful algal species that did not dominate community chlorophyll. Correlates from this study may be used to inform the next iteration of regionally predictive Vibrio models and may lend additional insight to Vibrio spp. ecology in similar systems. IMPORTANCE: Vibrio spp. are bacteria found in estuaries worldwide; some species can cause illness and infections in humans. Relationships between Vibrio spp. abundance, salinity, and temperature are well documented, but correlations to other environmental parameters are less understood. This study identifies unique correlates (e.g., E-W wind scalar and harmful algal species) that could potentially inform the next iteration of predictive Vibrio models for the EMSS region. Additionally, these correlates may allow existing environmental monitoring efforts to be leveraged in providing data inputs for future Vibrio risk models. An observed correlation between salinity and V. vulnificus/particle-size associations suggests that predicted environmental changes may affect the abundance of Vibrio spp. in certain reservoirs, which may alter which vectors present the greatest vibrio risk.

Organic polymer consumption facilitates domoic acid entry into the marine food web without direct ingestion of Pseudo-nitzschia.

Domoic acid (DA) is a neurotoxin produced by diatoms from the genera Pseudo-nitzschia and Nitzschia. DA is transferred through the food web when consumed by organisms such as copepods (e.g., Acartia tonsa). DA bioaccumulates in higher trophic levels and poses a threat to human health through amnesic shellfish poisoning. Laboratory experiments using a DA reference standard demonstrated that mild turbulence facilitates formation of organic polymer aggregates >0.6 µm in-vivo that can scavenge dissolved DA (dDA). Using A. tonsa, we demonstrate that DA can be assimilated through consumption of these organic polymers which scavenged dDA -a pathway which does not require direct ingestion of the toxin-producer Pseudo-nitzschia. In filtered seawater with spiked DA, copepods accumulated 24.8 ± 4.7 pg DA copepod-1 (2.1 ppm) on average by consuming organic polymers. This was validated in one out of five experiments using ambient DA concentrations. Copepods were suspended in particle-free seawater and accumulated 14.4 ± 3.8 pg DA copepod-1 (1.20 ppm), and in particle-concentrated seawater they accumulated 40.9 ± 3.8 pg DA copepod-1 (3.42 ppm). Data from this experiment suggests that ~34% of the total assimilated DA entered via an organic polymer-bound DA pathway. This experiment had the highest Pseudo-nitzschia spp. abundance (~225,000 cells L - 1) and cellular toxin quota, up to 0.88 pg DA cell-1, relative to the other four ambient DA experiments. These results demonstrate the potential for DA to enter the marine food web through an alternate pathway and may have considerable implications to understanding the flow of DA through marine food webs, and how we monitor DA and its potential vectors into the food web.

Silicon limitation facilitates virus infection and mortality of marine diatoms.

Diatoms are among the most globally distributed and ecologically successful organisms in the modern ocean, contributing upwards of 40% of total marine primary productivity1,2. By converting dissolved silicon into biogenic silica, and photosynthetically fixing carbon dioxide into particulate organic carbon, diatoms effectively couple the silicon (Si) and carbon cycles and ballast substantial vertical flux of carbon out of the euphotic zone into the mesopelagic and deep ocean3-5. Viruses are key players in ocean biogeochemical cycles6,7, yet little is known about how viral infection specifically impacts diatom populations. Here, we show that Si limitation facilitates virus infection and mortality in diatoms in the highly productive coastal waters of the California Current Ecosystem. Using metatranscriptomic analysis of cell-associated diatom viruses and targeted quantification of extracellular viruses, we found a link between Si stress and the early, active and lytic stages of viral infection. This relationship was also observed in cultures of the bloom-forming diatom Chaetoceros tenuissimus, where Si stress accelerated virus-induced mortality. Together, these findings contextualize viruses within the ecophysiological framework of Si availability and diatom-mediated biogeochemical cycling.

Silicic acid limitation drives bloom termination and potential carbon sequestration in an Arctic bloom.

The spring diatom bloom in the Arctic Ocean accounts for significant annual primary production leading to the most rapid annual drawdown of water-column pCO2. Late-winter waters in the Atlantic Arctic & Subarctic Provinces (AASP) have lower silicic acid concentrations than nitrate, which suggests diatom blooms may deplete Si before N. Here we test a facet of the hypothesis that silicic acid limitation terminates the spring diatom bloom in the AASP and the sinking of the senescent and dead diatoms helps drive carbon sequestration. During a 6-week study, diatoms bloomed and progressively consumed silicic acid to where it limited their growth. The onset of growth limitation was concurrent with the minimum pCO2 in the surface waters and increases in both the proportion of dead diatoms and the diatom assemblage sedimentation rate. Data reanalysis within the AASP shows a highly significant and positive correlation between silicic acid and pCO2 in the surface waters, but no significant relationship with nitrate and pCO2 was observed unless data were smoothed. Therefore, understanding the future of the AASP spring diatom bloom requires models that explicitly consider changes in silicic acid supply as a driver of this process.

Diatom populations in an upwelling environment decrease silica content to avoid growth limitation.

A mix of adaptive strategies enable diatoms to sustain rapid growth in dynamic ocean regions, making diatoms one of the most productive primary producers in the world. We illustrate one such strategy off coastal California that facilitates continued, high, cell division rates despite silicic acid stress. Using a fluorescent dye to measure single-cell diatom silica production rates, silicification (silica per unit area) and growth rates we show diatoms decrease silicification and maintain growth rate when silicon concentration limits silica production rates. While this physiological response to silicon stress was similar across taxa, in situ silicic acid concentration limited silica production rates by varying degrees for taxa within the same community. Despite this variability among taxa, silicon stress did not alter the contribution of specific taxa to total community silica production or to community composition. Maintenance of division rate at the expense of frustule thickness decreases cell density which could affect regional biogeochemical cycles. The reduction in frustule silicification also creates an ecological tradeoff: thinner frustules increase susceptibility to predation but reducing Si quotas maximizes cell abundance for a given pulse of silicic acid, thereby favouring a larger eventual population size which facilitates diatom persistence in habitats with pulsed resource supplies.

Taxon-specific contributions to silica production in natural diatom assemblages.

The metabolic activity and growth of phytoplankton taxa drives their ecological function and contribution to biogeochemical processes. We present the first quantitative, taxon-resolved silica production rates, growth rates, and silica content estimates for co-occurring diatoms along two cross-shelf transects off the California coast using the fluorescent tracer PDMPO (2-(4-pyridyl)-5-((4-(2-dimethylaminoethylaminocarbamoyl)methoxy)phenyl)oxazole), and confocal microscopy. Taxon contribution to total diatom community silica production was predominantly a function of the surface area of new frustule that each taxon created as opposed to cell abundance or frustule thickness. The influential role of surface area made large diatoms disproportionately important to community silica production over short time scales (<1 d). In some cases, large taxa that comprised only ~15% of numerical cell abundance accounted for over 50% of total community silica production. Over longer time scales relevant to bloom dynamics, the importance of surface area declines and growth rate becomes the dominant influence on contribution to production. The relative importance of surface area and growth rate in relation to silica production was modeled as the time needed for a smaller, faster-growing taxon to create more surface area than a larger, slower-growing taxon. Differences in growth rate between the taxa effected the model outcome more than differences in surface area. Shifts in relative silica production among taxa are time restricted by finite resources that limit the duration of a bloom. These patterns offer clues as to how taxa respond to their environment and the consequences for both species succession and the potential diatom contribution to elemental cycling.

Patterns and regulation of silicon accumulation in Synechococcus spp.

Six clones of the marine cyanobacterium Synechococcus, representing four major clades, were all found to contain significant amounts of silicon in culture. Growth rate was unaffected by silicic acid, Si(OH)4 , concentration between 1 and 120 µM suggesting that Synechococcus lacks an obligate need for silicon (Si). Strains contained two major pools of Si: an aqueous soluble and an aqueous insoluble pool. Soluble pool sizes correspond to estimated intracellular dissolved Si concentrations of 2-24 mM, which would be thermodynamically unstable implying the binding of intracellular soluble Si to organic ligands. The Si content of all clones was inversely related to growth rate and increased with higher [Si(OH)4 ] in the growth medium. Accumulation rates showed a unique bilinear response to increasing [Si(OH)4 ] from 1 to 500 µM with the rate of Si acquisition increasing abruptly between 80 and 100 µM Si(OH)4 . Although these linear responses imply some form of diffusion-mediated transport, Si uptake rates at low Si (~1 µM Si) were inhibited by orthophosphate, suggesting a role of phosphate transporters in Si acquisition. Theoretical calculations imply that observed Si acquisition rates are too rapid to be supported by lipid-solubility diffusion of Si through the plasmalemma; however, facilitated diffusion involving membrane protein channels may suffice. The data are used to construct a working model of the mechanisms governing the Si content and rate of Si acquisition in Synechococcus.

Net biogenic silica production and the contribution of diatoms to new production and organic matter export in the Costa Rica Dome ecosystem.

We determined the net rate of biogenic silica (bSiO2) production and estimated the diatom contribution to new production and organic matter export in the Costa Rica Dome during summer 2010. The shallow thermocline significantly reduces bSiO2 dissolution rates below the mixed layer, leading to significant enhancement of bSiO2 relative to organic matter (silicate-pump condition). This may explain why deep export of bSiO2 in this region is elevated by an order of magnitude relative to comparable systems. Diatom carbon, relative to autotrophic carbon, was low (<3%); however, the contribution of diatoms to new production averaged 3 and 13% using independent approaches. The 4-old discrepancy between methods may be explained by a low average C:Si ratio (∼1.4) for the net produced diatom C relative to the net produced bSiO2. We speculate that this low production ratio is not the result of reduced C, but may arise from a significant contribution of non-diatom silicifying organisms to bSiO2 production. The contribution of diatoms to organic matter export was minor (5.7%). These results, and those of the broader project, suggest substantial food-web transformation of diatom organic matter in the euphotic zone, which creates enriched bSiO2 relative to organic matter within the exported material.

Mesozooplankton biomass and grazing in the Costa Rica Dome: amplifying variability through the plankton food web.

We investigated standing stocks and grazing rates of mesozooplankton assemblages in the Costa Rica Dome (CRD), an open-ocean upwelling ecosystem in the eastern tropical Pacific. While phytoplankton biomass in the CRD is dominated by picophytoplankton (<2-µm cells) with especially high concentrations of Synechococcus spp., we found high mesozooplankton biomass (∼5 g dry weight m-2) and grazing impact (12-50% integrated water column chlorophyll a), indicative of efficient food web transfer from primary producers to higher levels. In contrast to the relative uniformity in water-column chlorophyll a and mesozooplankton biomass, variability in herbivory was substantial, with lower rates in the central dome region and higher rates in areas offset from the dome center. While grazing rates were unrelated to total phytoplankton, correlations with cyanobacteria (negative) and biogenic SiO2 production (positive) suggest that partitioning of primary production among phytoplankton sizes contributes to the variability observed in mesozooplankton metrics. We propose that advection of upwelled waters away from the dome center is accompanied by changes in mesozooplankton composition and grazing rates, reflecting small changes within the primary producers. Small changes within the phytoplankton community resulting in large changes in the mesozooplankton suggest that the variability in lower trophic level dynamics was effectively amplified through the food web.

Quantifying diatom silicification with the fluorescent dye, PDMPO.

Diatoms require silicic acid to construct ornately detailed cell walls called frustules. The growth and geographic distribution of diatoms is often controlled by the availability of silicic acid. Analytical methods exist to assess diatom community biogenic silica (bSiO2) production, but partitioning production among taxa has been largely qualitative. We present a method for the quantitative analysis of taxa-specific silica production through labeling diatoms with the fluorescent dye PDMPO [2-(4-pyridyl)-5-((4-(2-dimethylaminoethylaminocarbamoyl)methoxy)phenyl)oxazole]. To make PDMPO a quantitative tool: diatom frustules were solubilized to assess the total diatom community incorporation by quantitation of PDMPO fluorescence using a fluorometer, and laser confocal microscopy was used to quantify the fluorescence of PDMPO in single diatom cells. We created a fluorescence standard to intercalibrate the raw fluorescence signals of the fluorometer and microscope and to determine the fluorescence per mole of PDMPO. PDMPO incorporation was converted to silica production using diatom bSiO2:PDMPO incorporation ratios which varied systematically with silicic acid concentration. Above 3 µM Si(OH)4, bSiO2:PDMPO was constant and PDMPO incorporation was converted to silica production using a mole ratio of 2,916 as determined from cultures. Below 3 µM, the ratio was a linear function of [Si(OH)4] (bSiO2:PDMPO = 912.6 × [Si(OH)4]), as determined using data from two oceanographic cruises. Field evaluation of the method showed that total community PDMPO incorporation generally agreed to within 30% of radioisotope-determined silica production. This PDMPO method has the potential to be a powerful tool for understanding physiology, silicification and resource competition among diatom taxa.

Summer diatom blooms in the North Pacific subtropical gyre: 2008-2009.

The summertime North Pacific subtropical gyre has widespread phytoplankton blooms between Hawaii and the subtropical front (∼30°N) that appear as chlorophyll (chl) increases in satellite ocean color data. Nitrogen-fixing diatom symbioses (diatom-diazotroph associations: DDAs) often increase 10(2)-10(3) fold in these blooms and contribute to elevated export flux. In 2008 and 2009, two cruises targeted satellite chlorophyll blooms to examine DDA species abundance, chlorophyll concentration, biogenic silica concentration, and hydrography. Generalized observations that DDA blooms occur when the mixed layer depth is < 70 m are supported, but there is no consistent relationship between mixed layer depth, bloom intensity, or composition; regional blooms between 22-34°N occur within a broader temperature range (21-26°C) than previously reported. In both years, the Hemiaulus-Richelia and Rhizosolenia-Richelia DDAs increased 10(2)-10(3) over background concentrations within satellite-defined bloom features. The two years share a common trend of Hemiaulus dominance of the DDAs and substantial increases in the >10 µm chl a fraction (∼40-90+% of total chl a). Integrated diatom abundance varied 10-fold over <10 km. Biogenic silica concentration tracked diatom abundance, was dominated by the >10 µm size fraction, and increased up to 5-fold in the blooms. The two years differed in the magnitude of the surface chl a increase (2009>2008), the abundance of pennate diatoms within the bloom (2009>2008), and the substantially greater mixed layer depth in 2009. Only the 2009 bloom had sufficient chl a in the >10 µm fraction to produce the observed ocean color chl increase. Blooms had high spatial variability; ocean color images likely average over numerous small events over time and space scales that exceed the individual event scale. Summertime DDA export flux noted at the Hawaii time-series Sta. ALOHA is probably a generalized feature of the eastern N. Pacific north to the subtropical front.